CC BY
31
ИССЛЕДОВАНИЯ ВОДНЫХ БИОЛОГИЧЕСКИХ РЕСУРСОВ КАМЧАТКИ И СЕВЕРО-ЗАПАДНОЙ ЧАСТИ ТИХОГО ОКЕАНА, 2017, вып. 45
УДК 597.553.2 DOI: 10.15853/2072-8212.2017.45.5-23
К ВОПРОСУ ОБ ИСПОЛЬЗОВАНИИ СТРУКТУРЫ ЧЕШУИ АЗИАТСКОЙ ГОРБУШИ (ОМСОЯНУЖНШ вОЯБШСНА, WALBAUM), С ЦЕЛЬЮ ИДЕНТИФИКАЦИИ ВНУТРИВИДОВОГО СОСТАВА МОРСКИХ УЛОВОВ
В.А. Савин, Р.А. Шапорев*, В.С. И, А.В. Бугаев
Н. с.; инж.-исслед.; зав. отд., док. биол. наук.; Камчатский научно-исследовательский институт рыбного хозяйства и океанографии
683000 Петропавловск-Камчатский, Набережная, 18 Тел., факс: (4152) 41-27-01, 20-16-08. E-mail: savin.v.a@kamniro.ru *Вед. инж.-программист, Правительство Камчатского края 683000 Петропавловск-Камчатский, пр. Ленина, 1
ГОРБУША, ЧЕШУЯ, СКЛЕРИТЫ, РЕГИОНАЛЬНЫЙ КОМПЛЕКС ЛОКАЛЬНЫХ СТАД, ДИФФЕРЕНЦИАЦИЯ СТАД, ИДЕНТИФИКАЦИЯ СТАД, СМЕШАННЫЕ МОРСКИЕ УЛОВЫ Представлены результаты исследования структуры чешуи азиатской горбуши и оценка возможности ее использования в качестве дифференцирующего маркера при идентификации смешанных морских уловов. На основе анализа реперной базы данных за 2010 г. выделены пять региональных комплексов стад, характеризующихся своеобразием формирования склеритов в первой годовой зоне на чешуе: 1) Япония (о. Хоккайдо); 2) о. Сахалин; 3) северокурильские популяции (о. Парамушир); 4) Западная Камчатка; 5) Восточная Камчатка. Проведен ряд симуляционных тестов, в результате которых выявлена достаточно высокая региональная разрешающая способность созданной модели (в среднем 93,6%). Выполнена идентификация выборок производителей из дрифтерных уловов РМС «Кадет-701» и СРТМ-К «Леонид Ходыко», отобранных в июне-августе 2010 г., и показано, что ее результаты вполне соответствуют общей промысловой обстановке в исследуемом году. Проанализированы зависимости коэффициентов вариации региональных бутстреп-оценок от количества используемых склеритов ПГЗ. Выявлено, что устойчивость модели достигается при использовании набора из 11-15 склеритов, что дает основание говорить о возможности проводить идентификацию регионального происхождения молоди азиатской горбуши в морских уловах.
TO THE ISSUE OF THE USE OF ASIAN PINK SALMON (ONCORHYNCHUS GORBUSCHA, WALBAUM) SCALE STRUCTURE FOR IDENTIFICATION IN MIXED MARINE CATCHES
Vladimir А. Savin, Roman А. Shaporev*, Viktor S. I, Аlexander V. Bugaev
Researcher; Researcher Engineer; Head of Dep., Dr. Sc. (Biology); Kamchatka Research Institute of Fisheries and Oceanography
683000 Petropavlovsk-Kamchatsky, Naberezhnaya, 18 Tel., fax: (4152) 41-27-01, 20-16-08. E-mail: savin.v.a@kamniro.ru *Leading software engineer, The Government of the Kamchatka Territory 683000 Petropavlovsk-Kamchatsky, Lenin Square, 1
PINK SALMON, SCALE, SCLERITES, REGIONAL COMPLEX OF LOCAL STOCKS, STOCK DIFFERENTIATION, STOCK IDENTIFICATION, MIXED MARINE CATCHES
Results of analysis of Asian pink salmon scale structure are demonstrated and estimation of the possible use of the structure as a differentiation marker in identification of mixed marine catches is made. Five regional complexes of stocks are figured out on the base of the analysis of baselines collected in 2010, each having specific formation of sclerites in the first annual zone on the scale: 1) of Japan (Hokkaido); 2) of Sakhalin; 3) Northern Kuril populations (Paramushir); 4) of West Kamchatka; 5) of East Kamchatka. A complex of simulation tests is applied and results demonstrate quite high regional resolution ability of the model suggested (93,6% averaged). Identification is made for the samples from the drift net catches collected by RMV "Kadet-701" and SRTM-K "Leonid Khodyko" in June-August of 2010, and it is demonstrated that the results illustrate general commercial conjuncture in the year of the research. Correlations between coefficients of variation of the regional boot-stripe values and the number of sclerites used (in the zone of the 1st year) were analyzed. It is demonstrated that the model is stable when the number of the sclerites used is 11-15, and that allows regional identification of Asian juvenile pink salmon in marine catches possible.
Необходимость организации и регулирования промысла ставит перед современной рыбохозяйствен-ной наукой ряд задач, где одной из основных является дифференциация локальных стад тихоокеанских лососей и идентификация их на путях нагульных и нерестовых миграций. Особенно
остро этот вопрос стоит для горбуши (Oncorhynchus gorbuscha, Walbaum) — наиболее массового вида тихоокеанских лососей, которому свойственны не только многолетние тренды изменения численности, но и значительные колебания запаса в смежные годы, что связано с обычным для многих рай-
онов доминированием одного из поколений четных или нечетных лет.
Современные популяционные исследования лососей, связанные с выявлением внутривидовой подразделенности, в основном опираются на два методологических подхода: выявление фенотипи-ческих и генетических дифференцирующих маркеров. Каждый из методов обладает своими достоинствами и недостатками, которые необходимо учитывать для достижения оптимального результата при выполнении поставленных рыбохозяй-ственных задач. Настоящее исследование относится к первому направлению, т. к. в качестве дифференцирующего критерия в работе использована чешуя азиатской горбуши.
Использование чешуи горбуши в качестве дискриминирующего маркера, в отличие от других видов лососей, связано с определенными специфическими особенностями. К ним, в первую очередь, относится отсутствие пресноводного кольца, анализ структуры которого достаточно успешно применяется, в частности, при дифференциации нерки (Бугаев В.Ф., 1978, 1995). Это связано с коротким периодом пребывания горбуши в пресной речной воде: закладка чешуи происходит при достижении рыбой длины 4-6 см (Pearson, 1966), когда она уже скатилась и находится в прибрежных морских водах. Тем не менее к настоящему времени опубликовано уже достаточное количество работ, которые позволяют говорить о наличии региональных особенностей при формировании склеритов у данного вида, которые можно использовать в качестве дифференцирующих маркеров. В качестве примера ниже рассмотрим основные методы анализа изменчивости чешуи горбуши и некоторые результаты, полученные отечественными учеными.
В начале восьмидесятых годов прошлого века анализ различных сочетаний типов роста чешуи, ширины четырех зон из пяти склеритов и количества склеритов в зоне первого года жизни позволил выделить три крупные азиатские группировки стад: восточнокамчатскую, североохотоморскую и южноохотоморскую (Грачев, 1983). Позднее при исследовании чешуи горбуши из основных районов воспроизводства и промысла на Сахалине был впервые применен метод склеритограмм (Иванова, Октябрьский, 1990). Авторами показано, что значения величины склеритов в зимнем кольце и первых склеритов в плюсовой зоне позволяют диф-
ференцировать отдельные локальные группировки (япономорского и тихоокеанского происхождения) горбуши в смешанных скоплениях.
Комплексное исследование изменчивости сахалинской горбуши, включающее анализ структуры чешуи и морфометрических данных, позволило выявить отличия локальных стад Юго-Восточного Сахалина от популяций юго-западной части острова и залива Анива (Иванков и др., 1996). И.М. Иванова, изучив пять критериев первой морской зоны чешуи (число склеритов в первой и второй половинах этой зоны, расстояние от центра чешуи до внешнего края зоны, расстояние от центра до внешнего края пятого склерита и расстояние между внешними краями пятого и десятого склеритов), показала достоверные фено-типические различия между локальными стадами р. Амур, Восточной Камчатки, о-вов Сахалин, Кунашир и Хоккайдо (Иванова, 1997). На основе полученных результатов была выполнена идентификация смешанных морских выборок, выловленных в Японском и Охотском морях, а также в при-курильских водах Тихого океана.
Цикл работ, посвященный изучению изменчивости структуры чешуи горбуши и идентификации на ее основе траловых уловов, опубликован О.С. Темных и соавторами (Темных, 1997, 1998, 2001а, б, 2004; Темных и др., 1997, 2002). В данных исследованиях использован метод склеритограмм и показано, что популяции, относящиеся к крупным региональным комплексам, имеют выраженные различия в строении первой годовой зоны на чешуе (ПГЗ). Так, склеритограмма ПГЗ горбуши Восточносахалинского и Южнокурильского регионов характеризуется наличием зоны локального минимума, соответствующей зоне эстуарного кольца, в то время как для горбуши, нерестящейся в реках Камчатки и северного побережья Охотского моря, распределение межсклеритных расстояний на склеритограмме ПГЗ более равномерное.
Указанная особенность позволяет разделять эти два крупных («южный» и «северный») миграционных потока в смешанных скоплениях в морской период жизни и, по мнению автора, имеет наибольшую целесообразность применения для Охотоморского региона (Темных, 1998).
Н.П. Антонов с соавторами, изучив изменчивость чешуи горбуши, нерестящейся в реках Камчатки, северного побережья Охотского моря, а
также островов Хоккайдо, Сахалин и Итуруп в 1997-2000 гг., пришли к выводу, что вероятность идентификации отдельных локальных стад достаточно низка и колеблется в пределах 0-85% (Антонов, Балуева, 2000; Антонов и др., 2005). Отмечено, что популяции четных лет дифференцируются, в среднем, с большей точностью, чем нечетных. Для увеличения точности идентификации авторы объединяют близко расположенные локальности по степени сходства структуры чешуи, на основании чего выделяют следующие группы популяций: Японская, Сахалинская, Иту-рупская, Магаданская, Западнокамчатская и Вос-точнокамчатская. Также, на основе созданной ре-перной базы данных, была выполнена идентификация смешанных морских выборок из дрифтер-ных уловов в Беринговом море и северо-западной части Тихого океана (Антонов, Балуева, 2000).
Метод сопоставления склеритограмм позволил выявить высокодостоверные различия горбуши внутри Сахалино-Курильского региона (Каев, Ро-масенко, 2001). Анализируя первую и вторую годовые зоны роста чешуи, авторы дифференцируют региональные комплексы, жизненный цикл которых связан с обитанием в Японском и Охотском морях. Позднее А.М. Каев, используя распределения склеритограмм ПГЗ на чешуе, пришел к выводу, что превышение в 4,6 раза вылова горбуши над прогнозными ожиданиями у восточного побережья Сахалина в 2006 году связано с необычно высокой выживаемостью данного вида за годы исследований в этом регионе (Каев, 2007).
Подводя итог вышесказанному, можно сделать вывод, что чешуя горбуши, при относительной доступности сбора материала и невысокой себестоимости работ, является достаточно информативным маркером. Несмотря на то, что используемые в описанных работах методы исследования дифференцирующих особенностей формирования склеритов не позволяли получать стабильные надежные оценки на уровне отдельных популяций, при анализе крупных региональных группировок вероятность их идентификации, как правило, была на уровне, удовлетворяющем требованиям, предъявляемым к рыбохозяйственной науке в текущий момент.
В настоящее время развитие компьютерных технологий позволяет существенно расширить возможности исследования в данном направлении. Автоматизация процесса выбора и вычисления
дифференцирующих критериев дает возможность рассчитывать на значительное сокращение времени обработки при максимально полном и детальном анализе строения чешуи. Такой подход, с большой долей вероятности, должен способствовать улучшению результатов дифференциации как региональных комплексов популяций, так и отдельных локальных стад, что, в свою очередь, должно отразиться и на более точных оценках идентификации смешанных морских скоплений на путях нагульных и преднерестовых миграций.
Последние несколько лет в КамчатНИРО проводились исследования, целью которых являлись разработка и усовершенствование методики автоматизированного получения дискриминирующих критериев на основе подготовленного заранее пакета сканированных имиджей чешуи. Первоначально эта методика была применена для выявления фенотипической изменчивости кеты (Шапорев и др., 2007, 2014; Бугаев А.В. и др., 2012а, б, 2014; Bugaev et al., 2012; Савин и др., 2015), а позднее — адаптирована для анализа чешуи горбуши. Полученные предварительные результаты, представленные в ряде публикаций (Savin et al., 2012; Bugaev et al., 2012; Бугаев А.В. и др., 2013; Бугаев А.В., 2015), показали перспективность данного подхода и необходимость продолжения научных изысканий в выбранном направлении.
Таким образом, настоящая работа представляет собой как обобщение ранее опубликованных нами результатов, так и новые данные об изменчивости структуры чешуи азиатской горбуши, а ее цель — оценка возможности применения особенностей строения чешуи горбуши в качестве дифференцирующего маркера при региональной идентификации смешанных морских уловов.
МАТЕРИАЛ И МЕТОДИКА
Для создания реперной базы данных были использованы чешуйные пробы, собранные из речных уловов горбуши западного и восточного побережий Камчатки, Сахалина, о-вов Парамушир (Курильские о-ва) и Хоккайдо (Япония) сотрудниками КамчатНИРО, СахНИРО, Севвострыбвода, а также национального лососевого центра (г. Саппоро, Япония) в 2010 г. (рис. 1). Всего в качестве реперов исследовано 3486 особей горбуши из 33 азиатских локальных стад, а общее количество обработанных образцов чешуи составило 11 825 экземпляров (табл. 1).
В качестве смешанных выборок для данного исследования был использован материал, собранный в ходе проведения дрифтерного промысла в водах ИЭЗ РФ в июне-августе 2010 г. на промысловых судах РМС «Кадет-701» и СРТМ-К «Леонид Ходыко» (рис. 2). Лов осуществлялся в следующих промысловых районах: Карагинская подзона — 61.02.1 (период взятия проб: 12.0623.06; 4.07), Петропавловск-Командорская подзона — 61.02.2 (24.06-19.07), Тихоокеанская подзона — 61.03.1 (10.06-23.06; 13.07; 20.07-23.07; 8.08-16.08), Камчатско-Курильская подзона —
61.05.4 (25.07-6.08). Производителей горбуши отбирали из стандартных промысловых дрифтер-ных сетей (длина одной сети ~ 50 м, высота стенки ~ 9 м, ячея — 65 мм), после чего анализировали согласно стандартной методике, принятой для данных экспедиционных исследований (Карпенко и др., 1997; Ерохин, 2007).
За все время исследований был выполнен полный биологический анализ (ПБА) для 726 особей горбуши: 376 экз. на борту РМС «Кадет-701» и 350 экз. — на СРТМ-К «Леонид Ходыко». После камеральной обработки и отбраковки чешуи, не-
Рис. 1. Локализация сбора проб горбуши в нерестовых водоемах западного и восточного побережий Камчатки, Сахалина, о-вов Парамушир (Курильские о-ва) и Хоккайдо (Япония), использованных для формирования че-шуйной реперной базы 2010 г. Номера соответствуют таковым в табл. 1
Fig. 1. The map of sampling pink salmon in spawning rivers on the west and east coasts of Kamchatka, Sakhalin, Para-mushir and Hokkaido islands use for making scale baseline 2010. The numbers like in the table 1
пригодной для дальнейшей работы (с разрушенной структурой склеритов), общее количество исследуемых особей немного уменьшилось (табл. 2).
Создание цифровых имиджей чешуи из репер-ных и смешанных выборок проводили с помощью визуально-аналитической системы LEICA DMLS DC 100. Склеритограммы по первым 20 склеритам получали по всей чешуйной пластинке с помощью авторского программного пакета Р.А. Шапорева и др. (2007). При этом, по возможности, использовали по несколько образцов чешуи (максимум
восемь) от каждой особи в целях получения более объективных и устойчивых оценок.
Дифференциацию комплексов стад горбуши осуществляли методом многомерного шкалирования на основе квадрата расстояний Махалано-биса в статистическом пакете Statistica 6.0. Выбор метода многомерного шкалирования, как способа представления структурной организации изучаемого массива данных, определен нечетким или вероятностным характером кластеризации реперов в соответствии с их региональной принад-
Таблица 1. Локализация, период сбора и объем обработанного материала, использованного при формировании реперной базы данных за 2010 г.
Table 1. The area and the time of sampling, the sample size analyzed for making the baseline 2010
Регион № Локальность Период сбора Time of sampling ^личество I Number
Region Sampling area особей individuals чешуи scales
Япония (о. Хоккайдо) 1 р. Ичани / Ichani 15.09 187 187
Japan (Hokkaido) 2 р. Токор / Tokoro 24.08-27.09 239 239
3 пр. оз. Айнское / Aynskoye 16.08 31 107
4 р. Новоселка / Novoselka 1.07 93 270
5 р. Урюм / Urum 27.08 98 450
6 р. Очепуха / Ochepukha август 50 188
О. Сахалин 7 р. Бахура / Bakhura 6.08 50 187
Sakhalin 8 р. Найба / Naiba 6.07 20 71
9 р. Фирсовка / Firsovka 4.08 50 198
10 р. Дудинка / Dudinka 27.07-20.08 207 834
11 р. Лазовая / Lazovaya 19.07 59 165
12 р. Поронай / Poronay 21.08 25 74
13 р. Лангери / Langeri 13.08 98 364
Курильские о-ва (о. Парамушир) 14 р. Шелеховка / Shelekhovka 3.08 33 50
Kuril Islands (Paramushir) 15 р. Океанская / Okeanskaya 11.08 37 60
Западная Камчатка West Kamchatka
16 р. Палана / Palana 30.07-18.08 68 270
17 р. Хайрюзова / Khairyuzova 4.08-13.08 120 478
18 р. Ича / Icha 5.07-11.08 171 663
19 р. Крутогорова Krutogorova 30.07-2.09 165 640
20 р. Колпакова / Kolpakova 24.07-4.09 165 642
21 р. Б. Воровская / B. Vorovskaya 25.06-21.08 174 658
22 р. Кихчик / Kikhchik 18.08-3.09 152 521
23 р. Утка / Utka 23.07-13.08 167 593
24 р. Большая / Bolshaya 24.07-4.08 150 578
2 5_р. Опала / Opala_15.07-2.08_155_612
26 р. Паратунка I Paratunka - 75 285
27 р. Авача I Avacha 22.06 48 183
28 р. Налычева I Nalycheva 2.07 93 364
Восточная Камчатка 29 р. Жупанова I Zhupanova 11.07-28.07 149 575
East Kamchatka 30 р. Хайлюля I Khailulya 17.07-30.07 144 527
31 р. Вывенка I Vyvenka 19.06-27.06 67 227
32 р. Апука I Apuka 5.07 69 266
33 р. Аничкланваям Anichklanvayam 20.07-22.07 77 299
Всего / In total 3486 11 825
Таблица 2. Локализация, период сбора и объем обработанного материала из дрифтерных уловов РМС «Кадет-701» и СРТМ-К «Леонид Ходыко» в 2010 г.
Table 2. The area and the time of sampling, the sample size analyzed in the drift net catches of the RMV "Kadet-701" and the SRTM-К "Leonid Khodyko" in 2010_
Подзоны Subzone Период взятия проб Time of sampling Число особей Fish sample size Кол-во образцов чешуи Scale sample size
61.02.1 61.02.2 61.03.1 61.05.4 12.06-23.06; 04.07 24.06-19.07 10.06-23.06; 13.07; 20.07-23.07; 8.08-16.08 25.07-06.08 37 302 167 211 87 886 526 605
ИЭЗ РФ RF EEZ 12.06-12.08 717 2104
Рис. 2. Схема промрайонов и локализации мест сбора проб горбуши из дрифтерных уловов РМС «Кадет-701» и СРТМ-К «Леонид Ходыко» в июне-августе 2010 г.
Fig. 2. The scheme of the commercial fishery districts and of the sampling points from pink salmon drift net catches by the RMV "Kadet-701" and the SRTM-К "Leonid Khodyko" in June-August of 2010
лежностью. Классификация (идентификация) смешанных выборок произведена методом максимального правдоподобия. Реализация EM-алгоритма по MLE-оценкам долей гауссовских смесей с процедурой бутстрепа, а также расчет линейных дискриминантных функций выполнены в среде статистической обработки R 2.14.1.
РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ
Реперная база данных
График распределения центроидов реперов в пространстве многомерного шкалирования представлен на рисунке 3, анализ которого позволил определить следующие основные региональные комплексы горбуши четного поколения: 1) японские популяции о. Хоккайдо; 2) о. Сахалин; 3) се-
верокурильские популяции о. Па-рамушир; 4) Западная Камчатка и 5) Восточная Камчатка (представлена двумя группами стад, нерестящихся в юго-восточной и северовосточной частях полуострова). Характер агломерации выборок хорошо соответствует их географической подразделенности и в целом подобен региональным комплексам, выделенным ранее на основе анализа строения чешуи (Антонов и др., 2005), а также генетическими методами (Варнавская, 2006). Необходимо отметить, что отсутствие в реперной базе стад материкового побережья Охотского и Японского морей, а также довольно ограниченная информация о своеобразии популяций Курильских о-вов (в частности, отсутствие реперов крупных и важных, в промысловом значении, южно-курильских стад) и Японии не позволяют нам в полной мере представить схему региональной организации азиатской горбуши четной линии.
Для оценки качества построенной модели, в частности для определения ее разрешающей способности, был проведен ряд симуляци-онных тестов, результатом которых являются распределения оценок и их соответствующие параметры различных, т. н. «нулевых» выборочных сценариев. Значения региональных долей получены как результат суммирования MLE-долей реперов определяющих данный регион. Выборки из рек Шелеховка и Океанская объединены в единый репер Курильских островов в силу их недостаточной репрезентативности. Результаты симуляции показали, что региональная разрешающая способность созданной модели составляет от 78,6% для Курильских островов до 99,9% для Восточной Камчатки (в среднем 93,6%) (табл. 3). На уровне отдельных реперов региональная разрешающая способность составляла в среднем 91,0%. В серии тестов Хотеллинга прецедентов идентичности оценок векторов средних по выборкам из исследуемых локальностей выявлено не было. Таким
образом, оцененный уровень точности является наилучшим при идентификации горбуши четного поколения на основе изменчивости структуры чешуи и сопоставим или превышает аналогичные значения, полученные в ходе предыдущих исследований изменчивости изозимов (Варнавская, 2006) и митохондриальной ДНК (Шпигальская и др., 2012). Применение определенных подходов при формировании реперной базы данных, а именно расчет склеритограммы по всей площади че-шуйной пластинки и использование дублирующих образцов чешуи с каждой особи, позволило частично элиминировать источники вариабельности классов, и, как следствие, в значительной мере снизить величину ошибки модели.
Необходимо отметить, что используемый при построении классификатора принцип максимального правдоподобия является одним из самых распространенных подходов при идентификации
смешанных выборок лососей, как по структуре чешуи (Бугаев А.В., 2003а, 2005, 2011; Бугаев А.В. и др., 2006, 2009, 2012а, б), так и по различным системам генетических маркеров (Варнавская, 2006; Шпигальская и др., 2012). Идентификация по указанному принципу, являясь оптимальной с точки зрения минимизации ее ошибок, тем не менее приводит к несколько оптимистичным результатам в оценке ее качества, поскольку налагает достаточно жесткие требования на соответствие выборок представлению плотностных распределений генеральных совокупностей идентифицируемых локальностей. При идентификации реальных смешанных морских выборок, в условиях отсутствия полной информации по характеристикам всех без исключения стад, определяющих исследуемые регионы, а также существования явления расовой подразделенности изучаемого вида, следует ожидать несколько большей доли
• Япония / Japan
Ф О. Сахалин / Sakhalin
• Курильские о-ва / Kuril Islands
• Восточная Камчатка / East Kamchatka
• Западная Камчатка / West Kamchatka
1—Р.Утка/UtkaR.
2 — Р. Ича / Icha R.
3 — P. Большая / Bolshaya R.
4 — P. Опала / Opala R.
5 — P. Хайрюзова / Khaiiyuzova R.
Рис. 3. Распределение чешуйных центроидов локальных стад горбуши Западной и Восточной Камчатки, Сахалина, о-вов Парамушир (Курильские о-ва) и Хоккайдо (Япония) 2010 г. возврата, построенное методом многомерного шкалирования на основе расстояний Махаланобиса
Fig. 3. The distribution of the sclae centroids of the local stocks of pink salmon of West and East Kamchatka, Sakhalin, Paramushir (Kuril Islands) and Hokkaido (Japan) in returns of 2010, boult by the method of the multi-dimensional scaling based on Mahalanobis distances
неверно классифицированных особей, чем указано в таблице 3.
В противоположность вышесказанному, наиболее консервативными являются оценки, полученные с помощью метода классических линейных дискриминантных функций, основная идея которого заключается в простом маркировании обла-
стей признакового пространства. Результаты идентификации региональных массивов данных с использованием этого метода показали, что оценки качества, несмотря на общее их снижение (в среднем с 93,6 до 80,5%), остаются достаточно высокими (табл. 4). Вероятность региональной идентификации отдельных стад также немного снизи-
Таблица 3. Значения оценок разрешающей способности региональных комплексов стад горбуши 2010 г. возврата, где: MLE — средняя оценка по максимальному правдоподобию; SD — стандартное отклонение; CI — 95% -й симметричный доверительный интервал; N — число особей
Table 3. The values of the resolution ability values for the complexes of pink salmon stocks in the return in 2010, where: MLE - the average value from the maximum credibility; SD - the standard deviatiob; CI - 95% symmetric confedence interval; N - fish sample size
N, экз. Западная Камчатка West Kamchatka Восточная Камчатка East Kamchatka Курильские острова* Kuril Islands* Сахалин Sakhalin Япония Japan
1487 722 70 781 426
MLE SD CI low high 0,994 0,007 0,980 1,000 0,005 0,007 0,000 0,018 о о о о ООО о о о о о о о о о о5оо1 о^оз о^оо о.оо8 о^оо о;оо1 о^оо о.оо2
MLE SD CI low high 0,001 0,002 0,000 0,005 0,999 0,003 0,994 1,000 о о о о 'о о о о о о о о о о о о о5оо1 оооо2 о^оо о.ооз о^оо о;оо1 о^оо о.оо2
MLE SD CI low high 0,015 0,019 0,000 0,052 0,080 0,056 0,000 0,191 о;786 о^бо о;668 о^оз оД17 о^55 о;ою оД24 о^оз о;оо8 о^оо очо18
MLE SD CI low 0,005 0,006 0,000 0,003 0,004 0,000 о^оо о^оо о^оо о5991 оооо7 о;977 оооо2 о^оз о^оо
high 0,016 0,009 о^оо 15ооо о;оо8
MLE SD CI low high 0,005 0,004 0,000 0,013 0,007 0,005 0,000 0,018 о^оо о^оо о^оо о.оо! о^77 о^22 оД34 о.Ш о;912 о^22 о;868 о.956
В среднем, % Average, % 93о6%
Примечание. Средние MLE оценки получены методом бутстрепа с числом повторов R = 100; бутстреп произведен по реперам и по смешанной выборке; * — оценки по данному региону приведены по объединенным выборкам в единый репер. Note. The average MLE values obtained by the bootstripe method with the number of repetitions R = 100; the bootstripe was made on the baselines and the mixed sample; * - the values for this region were transformed from the unitied samples into one baseline.
Таблица 4. Значения оценок разрешающей способности региональных комплексов стад горбуши 2010 г. возврата на основе линейных дискриминантных функций
Table 4. The values of the resolution ability for the regional stock complexes of pink salmon in return 2010 based on the linear discriminant functions
Регион Region Западная Камчатка West Kamchatka Восточная Камчатка East Kamchatka Курильские острова Kuril Islands Сахалин Sakhalin Япония Japan
N, экз. 1487 722 70 781 426
Западная Камчатка West Kamchatka Восточная Камчатка East Kamchatka Курильские острова Kuril Islands Сахалин / Sakhalin Япония/Japan 0,855 0,163 0,029 0,078 0,042 0,094 0,759 0,043 0,046 0,026 0,000 0,001 0,871 0,003 0,002 о^46 о^72 о^57 о;846 о.237 о;оо5 оооо4 о^оо о^26 о.693
В среднем по регионам, % Average by regions 8о;5%
Северные регионы Northern regions Южные регионы Southern regions
Северные регионы Northern regions 0,936 0,064
Южные регионы Southern regions 0,088 0,912
лась, составив в среднем 84,1%. Очевидно, что минимально возможные значения ошибок, которые могут быть достигнуты при этом подходе, вне зависимости от рассматриваемой для идентификации системы маркеров, определяются в конечном итоге уровнем существующего межрегионального стреинга. Таким образом, полученные значения ошибочной классификации особей могут служить некоторым ориентиром для предосторож-ной верхней оценки его величины. Например, в нашем случае, при рассмотрении двух крупных региональных комплексов стад, «северного» и
«южного», эти оценки составляют 6,4 и 8,8% соответственно (табл. 4), что вполне сопоставимо с результатами, полученными О.С. Темных (1998).
Идентификация смешанных выборок
Оценка регионального состава дрифтерных уловов в 2010 году показала, что в большинстве рассматриваемых промысловых районов преобладали стада горбуши, нерестящиеся в реках западного побережья Камчатки (табл. 5, рис. 4). Максимальная доля западнокамчатской рыбы выявлена в Камчатско-Курильской подзоне (87,5%), а
Таблица 5. Значения оценок регионального состава стад горбуши из дрифтерных уловов вИЭЗ РФ в 2010 г., где MLE — точечная оценка по максимальному правдоподобию состава наблюдаемой смеси; MLE — средняя оценка по максимальному правдоподобию; SD — стандартное отклонение; CI — 95%-й симметричный доверительный интервал
Table 5. The values of the regional composition of pink salmon stocks in the drift net catches in the EEZ of Russia in 2010, where MLE - a point value from the maximum likelihood of the composition of the observed mixture; MLE - the average value from the maximum likelihood; SD - the standard deviation; CI - 95% symmetrical confidence interval
Подзоны Subzones Число особей Number CI
Регион / Region MLE MLE SD low high
of fish
Западная Камчатка West Kamchatka Восточная Камчатка 61.02.1 37 East Kamchatka
Курильские острова / Kuril Islands Сахалин / Sakhalin _Япония / Japan
0,031
0,915
0,000 0,000 0,054
0,332
0,558
0,000 0,080 0,030
0,206
0,238
0,003 0,076 0,038
0,000
0,092
0,000 0,000 0,000
0,736
1,000
0,006 0,229 0,104
61.02.2
302
Западная Камчатка West Kamchatka Восточная Камчатка East Kamchatka Курильские острова Kuril Islands Сахалин / Sakhalin Япония/Japan
0,748 0,118
0,000
0,128 0,006
0,695 0,195
0,000
0,107 0,003
0,073 0,100
0,552 0,837
0,000 0,390
0,000 0,000 0,000
0,054 0,005
0,001 0,000
0,213 0,014
61.03.1
Западная Камчатка West Kamchatka Восточная Камчатка 167 East Kamchatka
Курильские острова / Kuril Islands Сахалин Sakhalin _Япония Japan
0,728
0,179
0,000 0,093 0,000
0,663
0,224
0,001 0,112 0,000
0,113
0,109
0,003 0,053 0,001
0,442 0,884
0,009 0,438
0,000 0,008 0,000
0,007 0,217 0,002
Западная Камчатка West Kamchatka Восточная Камчатка 61.05.4 211 East Kamchatka
Курильские острова / Kuril Islands Сахалин / Sakhalin _Япония / Japan
0,875
0,024
0,000 0,073 0,028
0,845
0,055
0,000 0,079 0,021
0,054
0,051
0,002 0,041 0,014
0,739 0,950
0,000
0,000 0,000 0,000
0,155
0,004 0,160 0,049
Суммарно по всем подзонам In all subzones in total
Западная Камчатка West Kamchatka Восточная Камчатка 717 East Kamchatka
Курильские острова / Kuril Islands Сахалин / Sakhalin Япония/Japan
0,788
0,100
0,000 0,102 0,010
0,740
0,147
0,000 0,105 0,009
0,047
0,066
0,000 0,041 0,005
0,648 0,832
0,017
0,000 0,024 0,000
0,276
0,001 0,186 0,019
Примечание. Средние MLE оценки получены методом бутстрепа с числом повторов R = 100; веден по реперам и по смешанной выборке.
Note. The average MLE values were obtained by the bootstripe method with the number of repetitions stripe was carried out on the baselines and the mixed sample.
бутстреп произ-R = 100; the boot-
^ / ч \
О. Сахалин / Sakhalin у
Япония / Japan
OQ
Восточная Камчатка / East Kamchatka t
Западная Камчатка / West Kamchatka ^ / \ \ О. Сахалин / Sakhalin /
Япония / Japan
•О
о<а
L
Рис. 4. Результаты региональной идентификации стад горбуши из дрифтерньк_уловов РМС «Кадет-701» и СРТМ-К «Леонид Ходыко» в июне-августе 2010 г: А — по промысловым подзонам; Б — в общем в ИЭЗ РФ Fig. 4. The results of the regional identification of pink salmon socks from the droft net catches by the RMV "Kadet-701" and the SRTM-К "Leonid Khodyko" in June-August с" """"
of Russia
; of 2010: А - by commercial subzones; Б - in the total in the EEZ
в среднем по всем промрайонам она составила 78,8%. Восточнокамчатская горбуша доминировала в уловах в Карагинской подзоне (91,5%). В других районах количество особей, идентифицируемое по принадлежности к Восточной Камчатке, колебалось от 2,4 до 17,9%. Суммарно же, по всей исследуемой акватории, доля рыбы, нерестящейся на восточном побережье Камчатки, составила около 10,0%. Практически таким же оказалось и общее количество горбуши, идентифицируемой по происхождению как сахалинская (10,2%) (рис. 4, табл. 5).
В трех (из четырех) промысловых районах было также отмечено небольшое количество особей, принадлежащих популяциям о. Хоккайдо, однако общая их численность не превысила 1,0%, что находится в пределах статистической ошибки модели. Проведенный анализ не выявил в дриф-терных уловах северокурильских популяций.
Анализируя результаты идентификации, возникает естественный вопрос: насколько полученные оценки соответствуют реальной промысловой обстановке в текущем году? Как известно, основными региональными центрами, определяющими численность дальневосточной горбуши четного поколения, являются стада рек Западной Камчатки, о. Сахалин и Южных Курильских островов (суммарно 88,4% от фактического прибрежного
вылова по Дальневосточному бассейну (Шунтов, Темных, 2010) и 83,9% по Азиатско-Тихоокеанскому региону в 2010 г.). Нельзя не отметить численность японских популяций, которые также вносят достаточно весомый вклад в общую численность вида в азиатской части ареала, составляя около 5% в общем прибрежном вылове 2010 г. (рис. 5, табл. 6).
Необходимо понимать, что данные прибрежных уловов, величина которых в немалой степени зависит от материально-технической базы рыбо-добывающих предприятий и других сопутствующих путине нюансов, не могут служить абсолютным ориентиром, на который можно безоглядно ссылаться, как на показатель региональной численности вида. Тем не менее они вполне могут использоваться для получения общих представлений о мощности нерестовых подходов. Таким образом, сопоставив данные идентификации дриф-терных уловов с относительными величинами выловов в соответствующих регионах в 2010 г., была получена следующая картина (рис. 6). Очевидно, что выявленные в смешанных морских выборках суммарные доли «северной» (Западная и Восточная Камчатка) рыбы существенно превышают соответствующие относительные величины прибрежного вылова, в то время как доли «южной» (Сахалин и Япония) — значительно ниже.
50 45 40 35 30 ^ 25 20 15 10 5 0
I
л
о
+-> И
&
Я
3
ев Л и
и cd
W
а
J-
I £
О
о PQ
.д ТЗ
«3
Я Н
О О
о %
I
Ш
3 х а я
<L>
О
1 О (Я
I
О
3
&
« о-
3 з
1'Е
Рч U еГ +* VO ¿4
<Я оЗ
и а я л
ц Я О Л
« -t?
йё ой
li
о £
О
«
к us о
к *
§
0
1
О
4>
г?
0
1
Рн
О
Л &
s
в &
a о
'БЬ
3
Рн
■з gp
в
ё о
а
и
St, 9
ГЗ
& i—>
« a
§ к
1 о в
5
3
Рис. 5. Относительные величины прибрежных уловов горбуши в промысловых районах Азиатско-Тихоокеанского бассейна в 2010 г. (Шунтов, Темных, 2010; данные NPAFC (http://www.npafc.org/new/science statistics.html)) Fig. 5. The relative values of the coastal catches of pink salmon in the commercial fishery zones in the Asian Pacific basin in 2010 (Shuntov, Temnykh, 2010; the data by NPAFC (http://www.npafc.org/new/science_statistics.html))
В полученных результатах идентификации, по нашему мнению, нет существенных противоречий, если учитывать ограниченность реперной базы и, в большей степени, специфику дрифтер-ного промысла. Так, в очередной раз отметим отсутствие реперов из целого ряда регионов, популяции которых хоть и являются в своем большинстве малочисленными, в сумме составляют более 20% общего прибрежного вылова Азиатско-Тихоокеанского бассейна в 2010 г. (в том числе крупные южнокурильские стада (14,2%)) (рис. 5, табл. 6). Таким образом, при проведении MLE-анализа, когда в смешанной выборке попадается особь из региона, не представленного в базе данных, модель относит ее к реперной выборке наиболее (максимально) схожей с ней, и, соответственно, ошибочно увеличивает (уменьшает) какую-либо оценку. Однако мы считаем, что непосредственно этот фактор не очень сильно влияет на итоговый результат, т. к. речь идет об идентификации полово-
зрелой части стад, идущей на нерест, в отличие от посткатадромных скоплений горбуши в Охотском море, когда происходит смешение рыб из разных регионов и отсутствие каких-либо реперов может существенно повлиять на получаемые оценки.
По нашему мнению, на результаты идентификации более значительное влияние оказывает непосредственно специфика отечественного дриф-терного промысла, который, как известно, нацелен, в основном, на вылов наиболее ценного вида — нерки. Вследствие этого дрифтерные порядки стараются выставлять таким образом, чтобы как можно более полно обловить преднерестовые скопления двух основных промысловых стад Камчатки — рек Озерная и Камчатка, среднегодовая численность нерестовой части которых относительно всей камчатской нерки составляет около 50 и 25% соответственно (Бугаев В.Ф., 1995; Варнавская и др., 2005). Таким образом, в середине июня - начале июля работы в Карагинской под-
Таблица 6. Фактический прибрежный вылов горбуши в Азиатско-Тихоокеанском регионе в 2010 г. (по: Шунтов, Темных, 2010; данные NPAFC (http://www.npafc.org/new/science_statistics.html))
Table 6. The infact coastal catch of pink salmon in the Asia-Pacific region in 2010 (from: Shuntov, Temnykh, 2010; the data by NPAFC)
Район / District I Вылов (тыс. т) / Catch (th. t) %
Западная Камчатка / West Kamchatka 90,1 42,9
Восточная Камчатка / East Kamchatka 6,9 3,3
О. Сахалин / Sakhalin 56,3 26,8
Северные Курильские о-ва / North Kuril Islands 0,3 0,1
Южные Курильские о-ва / South Kuril Islands 29,8 14,2
Амурский бассейн / The Amur basin 6,4 3,0
Татарский пролив (Хабаровский край) 6,0 2,9
Tatar Strait (Khabarovsk territory)
Охотско-Аянский / Okhotsk-Ayan 0,7 0,3
Приморье / Primorye 2,3 1,1
Магаданская область / Magadan region 0,4 0,2
Анадырский / Anadyr 0,1 0,0
Япония/Japan 10,8 5,1
Всего / In the total
210,1
100,0
ox
80
70
60
50
40
30
20
10
0
Рис. 6. Региональное соотношение производителей горбуши из дрифтер-ных уловов в ИЭЗ РФ (I) и общий прибрежный вылов вида в соответствующих районах промысла (относительно всего Азиатско-Тихоокеанского бассейна) (II) в 2010 г.
Fig. 6. The regional ratio for pink salmon in the drift net catches in the EEZ of Russia (I) and the total coastal catch of pink salmon in the relevant commercial fishery districts (in sense of the Asian Pacific basin in the whole) (II) in 2010
зоне проводились фактически на ее южной границе (56-я параллель), а в Петропавловск-Командорской подзоне — южнее командорской островной группы (рис. 2, табл. 2). В этот период на рассматриваемой акватории происходит, главным образом, облов основной и арьергардной частей нерестовых скоплений нерки р. Камчатки, а также второстепенных восточнокамчатских стад этого вида (Бугаев А.В., 20036). В связи с этим, вероятность попадания в дрифтерные сети производителей восточнокамчатской (в основном западнобе-ринговоморской) горбуши достаточно велика, т. к. именно здесь проходит один из путей ее анадром-ных миграций (Шунтов, Темных, 2011).
С началом июля дрифтерные суда начинают постепенно смещаться на юг, с целью облова мигрирующего озерновского стада, и к концу июля достигают оконечности полуострова (мыс Лопатка) (Бугаев А.В., 2003б). В этот период проходит вылов горбуши из многочисленных нерестовых
скоплений западнокамчатских и, в меньшей степени, сахалино-курильских популяций, идущих на нерест через северокурильские проливы (Шунтов, Темных, 2010, 2011). В Охотское море (Кам-чатско-Курильская подзона) суда заходят в конце июля - начале августа и выставляют промысловые порядки в непосредственной близости от 12-мильной зоны на траверзе р. Озерной, в связи с чем доминирование в уловах западнокамчатской горбуши выглядит более чем вероятно.
Еще одним фактором, влияющим на итоговые оценки в нашей работе, является характеристика орудий лова. По результатам летней траловой съемки НИС «Профессор Кагановский», проведенной с 1 июня по 8 июля в 2010 г. в прикуриль-ских водах Тихого океана, было выделено два района повышенных концентраций горбуши («южный» и «северный»), разделенных полосой пониженных уловов от Средних Курил на юго-восток (рис. 7). Авторы показывают, что к юго-западу от
Рис. 7. Пространственное распределение траловых уловов горбуши на фоне поверхностных изотерм в верхней эпипелагиали прикурильских вод Тихого океана (НИС «Профессор Кагановский», 1.06-8.07.2010 г.): 1 — 0; 2 — 1-10; 3 — 11-50; 4 — 51-100; 5 — 101-370 экз./час. трал. (Шунтов, Темных, 2010)
Fig. 7. The spatial distribution of the trawl catches of pink salmon against the background of the surface isotherms in the upper epipelagic zone of the Pacific Ocean waters near the Kurils (the R/V "Professor Kaganovsky", 1.06-8.07.2010): 1 - 0; 2 - 1-10; 3 - 11-50; 4 - 51-100; 5 - 101-370 fishes per an hour-long trawling (Shuntov, Temnykh, 2010)
этой полосы в большинстве тралений средний вес особей был несколько ниже и составил в среднем 0,97 кг, в то время как в северо-восточном направлении, на подходах к району дрифтерного промысла, средний вес был равен 1,10 кг (Шунтов, Темных, 2010). Как упоминалось ранее, работа дрифтерных судов осуществлялась промысловыми сетями с размером ячеи 65 мм, которые практически не облавливают мелкоразмерных рыб. В результате этого исследуемые в настоящей работе смешанные выборки горбуши представляли собой особей (более 85% самцы) со средним размером 50,0 см (SD±2,9 см) и весом 1,49 кг (SD±0,26 кг) (рис. 8). Таким образом, используемые дрифтерные сети производили облов менее половины мигрирующего размерного состава рыб, с акцентом в сторону североохотоморских группировок (в основном камчатских), в силу их больших
размерно-весовых характеристик (Шунтов, Темных, 2010, 2011).
Подводя итог вышесказанному, очевидно, что полученные в данной работе результаты, учитывая специфику промысла, вполне укладываются в современные представления о нерестовых миграциях азиатской горбуши. Используемые в настоящем исследовании методы и модели вполне заслуживают внимания в качестве инструмента для решения аналогичных рыбохозяйственных задач. Другой вопрос — насколько представленная методика актуальна и востребована в каждом конкретном случае. Если возвратиться непосредственно к нашей работе, то промышленный дрифтерный селективный вылов крупноразмерных особей горбуши в любом случае не дает полноценного представления о составе смешанных нерестовых скоплений, и следовательно, такая информация име-
сг-
20-| 181614 12 H 10 864 2-1 0
N=726
,пп
см / sm
20-| 181614
12-^ ю-864 2j
XL
Я
ГШ,
ооооооооооооооо ооооооооооооооо
ts m
о о о о
«") ЧО
г/g
M (S NN M (S M
Рис. 8. Гистограммы размерно-весовых показателей уловов горбуши, собранные в ходе проведения дрифтерного промысла 2010 г. (размер ячеи 65 мм)
Fig. 8. The histograms of the body length and weight of pink salmon, collectes during drift net fishery campaign in 2010 (the mesh size 65 mm)
ет довольно ограниченную область применения. Кроме того, с 1 января 2016 г. данный вид промысла запрещен Федеральным законом от 29.06.2015 г. № 208-ФЗ, и это на обозримую перспективу закрывает возможность проведения подобных работ. Что же касается траловых учетных съемок, где селективность орудий лова сведена к минимуму, то получение данных о региональном соотношении стад в смешанных морских уловах на путях миграций имеет большую актуальность, т. к. способствует заблаговременному и оперативному прогнозированию численности горбуши, возвращающейся на нерест к районам размножения.
Согласно более ранним исследованиям, как уже упоминалось выше, анализ изменчивости структуры чешуи (как отдельно, так и совместно с экологическими и биологическими признаками) позволяет с достаточной долей вероятности выделять в траловых уловах «северные» (камчатско-североохотоморские) и «южные» (сахалино-ку-рильские) группировки стад. При этом авторами показано, что региональные комплексы популяций, составляющие данные группировки, не обладают достаточным своеобразием по исследуемым чешуйным критериям, и, следовательно, их идентификация в смешанных уловах затруднена. Нет сомнения в том, что выявление в морских скоплениях горбуши (особенно в Охотском море и СЗТО) даже этих двух составляющих, относящихся к разным частям Азиатско-Тихоокеанского бассейна, позволяет получить важнейшую информацию о путях миграции и численности нагульных и преднерестовых стад. Тем не менее, возможность получения достоверных аналогичных данных по каждому из регионов воспроизводства могла бы, в перспективе, существенно повысить уровень разного рода прогнозов и, как следствие, благотворно отразиться на рыбопромышленной отрасли в целом. В настоящее время эту задачу пытаются решить с применением генетических методов, анализируя, в частности, изменчивость структуры митохондриальной ДНК (Шпигальская и др., 2011; Шевляков и др., 2012, 2014).
Не претендуя на оценку эффективности такого подхода, мы бы, в свою очередь, хотели представить наше предварительное заключение о перспективах использования особенностей структуры чешуи для идентификации траловых уловов с применением вышеописанной методики.
Как известно, региональные различия в строении первой годовой зоны на чешуе горбуши формируются, главным образом, в летний период. Следовательно, к моменту проведения осенних учетных траловых съемок в Охотском море основные дифференцирующие критерии должны быть уже сформированы. По опубликованным данным, количество склеритов на чешуе горбуши, выловленной в ноябре-декабре 1995 г. в тихоокеанских водах Курильских островов и западной части Охотского моря, может колебаться от 15 до 19 (Темных, 1998). Анализ чешуи нагульной молоди из траловых уловов НИС «Профессор Каганов-ский» в южной части Охотского моря в период 7-8 ноября 2009 г. показал, что количество склеритов варьировало в пределах 11-17, составляя в среднем около 14 (трал № 184 (48°51Т' с. ш. 148°46'1" в. д.), трал № 185 (49°32'0" с. ш. 150°00'0" в. д.)). Учитывая, что оценка разрешающей способности применяемой модели проводилась на основе изменчивости первых двадцати склеритов ПГЗ, важно понять, набор из какого количества склеритов можно использовать для идентификации горбуши из траловых уловов, чтобы не допустить качественных потерь. Чтобы получить ответ на этот вопрос, была исследована зависимость коэффициентов вариации региональных бутстреп-оценок от числа используемых первых склеритов ПГЗ. Выполненный анализ показал, что устойчивость бут-стреп-оценок наблюдалась на склеритораммах из 11-15 склеритов, что, в свою очередь, дает все основания говорить о возможности без потери качества проводить идентификацию постката-дромной молоди горбуши на путях ее нагульных миграций (рис. 9).
ЗАКЛЮЧЕНИЕ
Получение объективных данных по региональному составу морских скоплений лососей является важной составной частью информации, необходимой для составления долгосрочных и оперативных промысловых прогнозов. Представленная работа является своеобразной демонстрацией возможностей и перспектив использования изменчивости структуры чешуи горбуши в качестве дифференцирующего маркера при идентификации смешанных уловов на путях нагульных и преднерестовых миграций.
В результате проведенного анализа показано, что своеобразие формирования склеритов в ПГЗ
на чешуе позволяет выделять следующие региональные комплексы стад горбуши четного поколения: 1) Япония (о. Хоккайдо); 2) о. Сахалин;
3) северокурильские популяции (о. Парамушир);
4) Западная Камчатка; 5) Восточная Камчатка. При этом разрешаюшая способность реперной базы на уровне отдельных реперов составляла в среднем 91,0%, а на уровне регионов — 93,6%. Несомненно, что отсутствие реперного материала из крупных популяций южнокурильских островов, а также стад материкового побережья Охотского и Японского морей в определенной степени ограничивает полноту представления географической подраз-деленности вида. Тем не менее в использованной для настоящего исследования реперной базе данных представлены почти все основные региональные комплексы популяций, формирующие общую численность азиатской горбуши в четные годы. Кроме того, благодаря используемым методикам, полученный уровень точности разрешающей способности созданной модели является наилучшим при идентификации данного вида по изменчивости структуры чешуи и, как минимум, сопоставим
с аналогичными значениями, полученными при анализе полиморфизма генетических маркеров.
Оценка регионального состава дрифтерных уловов в ИЭЗ РФ в 2010 году показала, что в большинстве рассматриваемых промысловых районах доминировала западнокамчатская горбуша. Доли сахалинской и восточнокамчатской рыбы были существенно ниже и составили в сумме около 10,0% каждая. Полученные результаты, учитывая специфику отечественного «дрифтера», вполне соответствуют общей промысловой обстановке в исследуемом году.
Анализ зависимости коэффициентов вариации региональных бутстреп-оценок от числа используемых склеритов ПГЗ выявил, что применяемая модель проявляет устойчивость при использовании набора из 11-15 склеритов, что дает основание говорить о возможности без потери качества проводить идентификацию горбуши на путях ее нагульных миграций.
Подводя итог, хотелось бы акцентировать внимание на том, что несмотря на набирающие в последнее время популярность популяционно-гене-
10 ' и
Восточная Камчатка / East Kamchatka Западная Камчатка / West Kamchatka — Курильские о-ва / Kuril Islands
О. Сахалин / Sakhalin Япония/ Japan
12 ' is
Число склеритов/Numb,
er of sclerites
19 ' 20
Рис. 9. График зависимости коэффициентов вариации (CV) региональных бутстреп-оценок от количества используемых склеритов ПГЗ
Fig. 9. The graphic correlation between the coefficient of variations (CV) of the regional bootstripe results and the number of sclerites in the first year zone used
тические методы, потенциал изменчивости структуры чешуи горбуши, как маркера внутривидовой дифференциации, еще существенно недооценен. Представленные результаты исследования наглядно демонстрируют, что анализ строения чешуи, с использованием возможностей, предоставляемых научно-техническим прогрессом, позволяют получать достаточно высокие оценки разрешающей способности, достаточные, на наш взгляд, для идентификации преднерестовых скоплений. Также выявлены и показаны предпосылки для решения вопроса о возможности региональной идентификации посткатадромной молоди горбуши. Таким образом, с большой долей вероятности можно говорить о том, что уровень информативности структуры чешуи горбуши достаточно высок для успешного решения подобных задач, что, в комплексе с относительной доступностью сбора материала и невысокой себестоимостью работ, делает представленное направление исследований довольно перспективным.
СПИСОК ЛИТЕРАТУРЫ
АнтоновН.П., БалуеваЕ.С. 2000. Идентификация горбуши Oncorhynchus gorbuscha ^а1Ьаиш) из смешанных морских уловов по структуре чешуи // Исслед. водн. биол. ресурсов Камчатки и сев.-зап. части Тихого океана: Сб. науч. тр. Камчат. НИИ рыб. хоз-ва и океанографии. Вып. 5. С. 51-55. Антонов Н.П., Зорбиди Ж.Х., Балуева Е.С. 2005. Сравнительный анализ внутривидового разнообразия горбуши Oncorhynchus gorbuscha ^а1Ь.) и кижуча O. kisutch ^а1Ь.) северо-западной части Тихого океана по структуре чешуи // Популяцион-ная биология, генетика и систематика гидробион-тов. Сб. науч. тр. Т. 1. Петропавловск-Камчатский: КамчатНИРО. С. 160-167.
Бугаев А.В. 2003а. Идентификация локальных стад нерки Oncorhynchus nerka по чешуйным критериям в юго-западной части Берингова моря и сопредельных водах Тихого океана в период преднерестовых миграций. Сообщение 1 (Формирование реперных баз данных) // Изв. ТИНРО. Т. 132. С. 154-177. Бугаев А.В. 2003б. Идентификация локальных стад нерки Oncorhynchus nerka по чешуйным критериям в юго-западной части Берингова моря и сопредельных водах Тихого океана в период преднерестовых миграций. Сообщение 2 (Пространственное и темпоральное распределение) // Изв. ТИНРО. Т. 132. С.178-203.
Бугаев А.В. 2005. Идентификация локальных стад нерки Oncorhynchus nerka в дрифтерных уловах в 2001-2002 гг. // Вопр. ихтиологии. Т. 45. № 1. С. 41-54.
Бугаев А.В. 2011. Идентификация и промысловое изъятие азиатских стад нерки Oncorhynchus nerka в экономической зоне России // Изв. ТИНРО. Т. 167. С. 3-31.
Бугаев А.В. 2015. Преднерестовые миграции тихоокеанских лососей в экономической зоне России. Петропавловск-Камчатский: КамчатНИРО. 416 с. Бугаев А.В., Заволокина Е.А., Заварина Л.О., Шубин А.О., Золотухин С.Ф., Капланова Н.Ф., Во-лобуев М.В., Киреев И.Н. 2006. Идентификация локальных стад кеты Oncorhynchus keta в западной части Берингова моря по данным траловых съемок НИС «ТИНРО» в сентябре-октябре 2002-2003 гг. // Изв. ТИНРО. Т. 146. С. 3-34. Бугаев А.В., Заволокина Е.А., Заволокин А.В., Заварина Л.О., Киреев И.Н., Шубин А.О., Игнатьев Ю.И., Золотухин С.Ф., Капланова Н.Ф., Во-лобуевМ.В. 2009. Происхождение и распределение локальных стад кеты Oncorhynchus keta в западной части Берингова моря по данным траловых съемок НИС «ТИНРО» в 2004 и 2006 гг. // Изв. ТИНРО. Т. 157. С. 3-33.
Бугаев А.В., Шапорев Р.А., Заволокин А.В. 2014. Происхождение и распределение стад кеты Oncorhynchus keta в западной части Берингова моря и северо-западной части Тихого Океана в 2009 и 2010 гг. // Изв. ТИНРО. Т. 179. С. 177-203. Бугаев А.В., ШапоревР.А., КовальМ.В., Огородников В.С., Шубин А.О., КойновА.А., ХохловЮ.Н., Золотухин С.Ф., Кульбачный С.Е., ПодорожнюкЕ.В. 2012а. Распределение региональных комплексов азиатских стад кеты Oncorhynchus keta в период преднерестовых миграций в экономической зоне России // Изв. ТИНРО. Т. 171. С. 3-25. Бугаев А.В., Шапорев Р.А., Савин В.А., Чистякова А.И. 2012б. Результаты камчатских феноти-пических исследований тихоокеанских лососей Oncorhynchus spp. в 1990-2010 гг. // Исслед. водн. биол. ресурсов Камчатки и сев.-зап. части Тихого океана: Сб. науч. тр. Камчат. НИИ рыб. хоз-ва и океанографии. Вып. 25. С. 37-57. Бугаев А.В., Шапорев Р.А., Чистякова А.И., Савин В.А., Савенкова Е.В., ИВ.С. 2013. Результаты исследований, проведенных лабораторией морфологических исследований лососей // Мат-лы отчетной сессии ФГУП «КамчатНИРО» по итогам
научно-исследовательских работ в 2012 г. Петропавловск-Камчатский: КамчатНИРО. С. 106-121. Бугаев В.Ф. 1978. К экологии пресноводного периода жизни и дифференциации популяций нерки бассейна р. Камчатка / В кн.: Биология лососевых. Владивосток. С. 35-36.
Бугаев В.Ф. 1995. Азиатская нерка (пресноводный период жизни, структура локальных стад, динамика численности). М.: Колос. 464 с. Варнавская Н.В. 2006. Генетическая дифференциация популяций тихоокеанских лососей. Петропавловск-Камчатский: КамчатНИРО. 488 с. Варнавская Н.В., Маслов А.В., Савин В.А. 2005. Региональные комплексы популяций дальневосточных лососей // Популяционная биология, генетика и систематика гидробионтов. Сб. науч. тр. Т. 1. Петропавловск-Камчатский: КамчатНИРО. C. 8-27. Грачев Л.Е. 1983. Дифференциация азиатских стад горбуши // Биологические основы развития лососевого хозяйства в водоемах СССР. М.: Наука. С. 84-97.
Ерохин В.Г. 2007. Дрифтерные исследования морского периода жизни тихоокеанских лососей в при-камчатских водах ИЭЗ России в 1993-2006 гг. // Вопр. рыболовства. Т. 8. № 3 (31). С. 484-524. Иванков В.Н., Добрицкий О.Ю., Скуба Н.С., Карпенко А.И. 1996. Дифференциация популяций горбуши Oncorhynchus gorbuscha Южного Сахалина // Биология моря. Т. 22. № 3. С. 167-173. Иванова И.М. 1997. Идентификация стад горбуши в северо-западной части Тихого океана, Охотского и Японского морей по структуре чешуи // Тез. докл. I Конгресса ихтиологов России (Астрахань, сентябрь 1997). М.: ВНИРО. С. 39.
Иванова И.М., Октябрьский Г.А. 1990. Применение метода склеритограмм для дифференциации локальных стад сахалинской горбуши // Тез. докл. Межд. симпозиума по тихоокеанским лососям (Южно-Сахалинск, 9-17 сентября 1989). Владивосток: ТИНРО. С. 82-84.
Каев A.M. 2007. Чья горбуша на Сахалине? // Рыбное хозяйство. № 2. С. 52-54. Каев А.М., Ромасенко Л.В. 2001. Дифференциация пространственно-временных группировок горбуши Сахалино-Курильского региона на основании изучения склеритограмм // Вопр. рыболовства. Т. 2. № 4 (8). С. 638-652.
Карпенко В.И., Коваленко М.Н., Василец П.М., Багин Б.Н., Кондрашенков Е.Л., Ерохин В.Г., Адамов А.А., Смородин В.П., Максименков В.В., Яков-
лев В.М. 1997. Методика морских исследований тихоокеанских лососей. Метод. пособие. Петропавловск-Камчатский: КамчатНИРО. 64 с. Савин В.А., Шапорев Р.А., Бугаев А.В., Заволо-кин А.В., ИВ.С. 2015. Происхождение и распределение стад кеты Oncorhynchus keta в западной части Берингова моря и северо-западной части Тихого океана в летне-осенние периоды 2013 и 2014 гг. // Бюл. № 10 изучения тихоокеанских лососей на Дальнем Востоке. Владивосток: ТИНРО-Центр. С.80-95.
Темных О.С. 1997. Пространственная дифференциация азиатской горбуши в период морских миграций // Тез. докл. I Конгресса ихтиологов России (Астрахань, сентябрь 1997). М.: ВНИРО. С. 52. Темных О.С. 1998. Региональная изменчивость склеритограмм чешуи азиатской горбуши // Изв. ТИНРО. Т. 124. С. 375-390. Темных О.С. 2001а. Опыт прогнозирования подходов горбуши в Охотское море по данным морских траловых съемок // Вопр. рыболовства. Т. 2. № 1 (5). С. 140-153.
Темных О.С. 2001б. Горбуша 2001 года: где, сколько? // Рыбное хозяйство. № 2. С. 27-28. Темных О.С. 2004. Азиатская горбуша в морской период жизни: биология, пространственная дифференциация, место и роль в пелагических сообществах. Автореф. дис. ... докт. биол. наук. Владивосток: ТИНРО-Центр. 47 с.
Темных О.С., МалининаМ.Е., ПодлесныхА.В. 1997. Дифференциация анадромных миграционных потоков горбуши четных поколений в Охотском море в 90-е гг. // Изв. ТИНРО. Т. 122. С. 131-151. Темных О.С., Мельников И.В., Вайтукевич Д.В. 2002. Учет сеголетков горбуши в Охотском море и перспективы путины 2002 г. // Рыбное хозяйство. № 2. С. 31-32.
ШапоревР.А., Варнавская Н.В., Савин В.А., Степанов В.Г., Бажин А.Г. 2007. Методика комплексного анализа изменчивости строения слоистых структур чешуи тихоокеанских лососей на примере кеты Камчатки // Исслед. водн. биол. ресурсов Камчатки и сев.-зап. части Тихого океана: Сб. науч. тр. Камчат. НИИ рыб. хоз-ва и океанографии. Вып. 9. С. 131-142.
Шапорев Р.А., Савин В.А., И В.С., Бугаев А.В. 2014. Предварительные данные по идентификации стад кеты Oncorhynchus keta в западной части Берингова моря и северо-западной части Тихого океана в период нагульных и преднерестовых миграций
2013 г. // Бюл. № 9 изучения тихоокеанских лососей на Дальнем Востоке. Владивосток: ТИНРО-Центр. С.95-102.
Шевляков Е.А., Дубынин В.А., Ерохин В.Г., Бугаев В.Ф., Зорбиди Ж.Х., Заварина Л.О., Шубкин С.В., Зикунова О.В., Попова Т.А., Шпигальская Н.Ю. 2014. Характеристика прибрежного промысла тихоокеанских лососей в Камчатском крае в 2014 г. // Бюл. № 9 изучения тихоокеанских лососей на Дальнем Востоке. Владивосток: ТИНРО-Центр. С.11-36.
Шевляков Е.А., Ерохин В.Г., Дубынин В.А., Шпигальская Н.Ю., ЗорбидиЖ.Х. 2012. Характеристика прибрежного промысла лососей в Камчатском крае в 2012 г. // Бюллетень № 7 реализации «Концепции Дальневосточной бассейновой программы изучения тихоокеанских лососей». Владивосток. ТИНРО-Центр. С. 10-28.
Шпигальская Н.Ю., Брыков В.А., Кухлевский А.Д., Сараванский О.Н., Климов А.В., Четвертак А.А., Шевляков Е.А. 2011. Региональная идентификация смешанных морских скоплений молоди горбуши ^п^^п^ш gorbuscha Wa1baum) на основе изменчивости фрагмента СуШ/0-1оор митохондри-альной ДНК // Изв. ТИНРО. Т. 165. С. 89-103. ШпигальскаяН.Ю., Пильганчук О.А., СавенковВ.В., Кустова А.С., Муравская У.О., Сараванский О.Н. 2012. Результаты исследований внутривидовой
структуры промысловых видов рыб методами по-пуляционной генетики // Исслед. водн. биол. ресурсов Камчатки и сев.-зап. части Тихого океана: Сб. науч. тр. Камчат. НИИ рыб. хоз-ва и океанографии. Вып. 25. С. 69-97.
Шунтов В.П., Темных О.С. 2010. Лососевая путина — 2010: снова рекордный вылов, но только для четных лет // Бюл. № 5 реализации «Концепции дальневосточной бассейновой программы изучения тихоокеанских лососей». Владивосток: ТИНРО-Центр. С. 3-11.
Шунтов В.П., Темных О.С. 2011. Тихоокеанские лососи в морских и океанических экосистемах. Владивосток: ТИНРО-Центр. Т. 2. 473 с. Bugaev А.У., Shpigalskaya N.Yu., Pilganchuk O.A., Shaporev R.A., Savin V.A., Chistyakova A.I., Muravs-kaya U.O., Saravansky O.N. 2012. Review of some results of Russian population studies of Pacific salmon in 2010 and 2011. // NPAFC Doc. No. 1397. 18 p. Pearson R.E. 1966. Report on the investigations by the United States for the International North Pacific Fisheries Commission - 1964: Pink salmon scale studies // INPFC Annu. Rep. P. 118-124. Savin V.A., Shaporev R.A., Savenkova E.V., Varnavs-kaya N.V. 2012. Differentiation of local pink salmon stocks on the basis of variations in their scale structure // NPAFC Tech. Rep. No 8. P. 67-70.
