Авторефераты по всем темам  >>  Авторефераты по техническим специальностям

На правах рукописи

УДК 639.2.055: 639.2.053.7

БУЛГАКОВА Татьяна Ивановна РЕГУЛИРОВАНИЕ МНОГОВИДОВОГО РЫБОЛОВСТВА НА ОСНОВЕ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ

Специальность - 05.18.17 Цпромышленное рыболовство

АВТОРЕФЕРАТ

диссертации на соискание ученой степени доктора технических наук

Москва - 2009 г.

Работа выполнена в лаборатории Системного анализа промысловых биоресурсов Всероссийского научно-исследовательского института рыбного хозяйства и океанографии (ФГУП "ВНИРО")

Официальные оппоненты:

доктор технических наук, профессор Мельников Виктор Николаевич доктор технических наук, профессор Юдович Юрий Борисович доктор биологических наук, профессор, член-корреспондент РАН Криксунов Евгений Аркадьевич

Ведущая организация: Федеральное государственное унитарное предприятие Полярный научно-исследовательский институт морского рыбного хозяйства и океанографии им. Н.М. Книповича (ПИНРО)

Защита состоится " 24 " декабря 2009 г. в 11 часов на заседании диссертационного совета Д 307.004.03 во Всероссийском научноисследовательском институте рыбного хозяйства и океанографии (ФГУП "ВНИРО") по адресу: 107140, Москва, В. Красносельская, д.17.

Факс: (499) 264-91-87, e-mail: fishing@vniro.ru

С диссертацией можно ознакомиться в библиотеке Всероссийского научно-исследовательского института рыбного хозяйства и океанографии

Автореферат разослан "___"_____________2009 г.

Отзывы просим направлять в двух экземплярах, заверенных печатью, по адресу: 107140, Москва, В. Красносельская, 17а.

Ученый секретарь диссертационного совета, кандидат технических наук В.А. Татарников

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Настоящая работа относится к проблеме рационального использования водных экологических систем, и как показывает отечественный и международный опыт, для построения теоретических основ решения этой проблемы необходимо применение методов математического моделирования. Ф.И. Баранов (1918, 1925) заложил основы теории регулирования рыболовства, в дальнейшем в ее развитии участвовали многие исследователи: А.Н. Державин, Г.В. Никольский, Г.Н. Монастырский, Н.Н.

Андреев, А.В. Засосов, У. Рикер, М. Шефер, Д. Бивертон, С. Холт, Д. Галланд, В.Н. Мельников, И.В. Никоноров, Дж. Шепард, Дж. Поуп, Р.Г. Бородин, В.К.

Бабаян, Д.А. Васильев, П.С. Гасюков, Ю.Н. Ефимов, В.Л. Третьяк и др. В большинстве случаев эти авторы исследовали одновидовые промысловые запасы. Математические модели (ММ) изолированного запаса строятся на основе предположения, что промысел является важным фактором, влияющим на динамику популяции, хотя для многих промысловых видов влияние хищничества вполне соизмеримо, а иногда и значительно превышает влияние промысла.

Теория управления многовидовыми запасами водных животных до сих пор достаточно не разработана, в ее основе должны лежать многовидовые ММ.

Начало использования ММ для анализа многовидового сообщества связывают с работами В. Вольтерра (Volterra, 1931; Volterra and d'Ancona, 1935), в которых модели динамики популяций в виде систем дифференциальных уравнений построены с учетом разных типов взаимодействий видов в сообществе.

Непосредственным развитием такого подхода для анализа динамики промыслового сообщества послужили многовидовые продукционные модели промысловых сообществ, причем одной из первых была работа автора (Булгакова, 1970; Pope and Harris,1975; Horwood,1976). Отличительные черты этой группы простых моделей сложных экосистем - относительная простота их исследования, возможность оценки параметров по промысловым данным и возможность анализа различных правил регулирования (ПР) промысла. В данной работе показано, как следует учитывать межвидовые отношения при выборе ПР промысла и при формулировке цели такого регулирования.

Многовидовая когортная модель MSVPA (Pope, 1979; Helgason and Gislason, 1979) разработана для сообществ рыб Северного моря и Балтики, и встала задача модификации этой модели для сообщества Баренцева моря, которая решена нами в соавторстве с Д.А.Васильевым и специалистами ПИНРО.

Всемирная Конференция по устойчивому развитию (Йоганнесбург, 2002) приняла решение перейти к применению экосистемного подхода для оценки морских ресурсов и регулирования их промысла уже к 2010 г. Настоящая работа напрямую посвящена анализу и развитию методологии экосистемного подхода с учетом многовидовых отношений и при включении некоторых факторов среды в процедуры оценки запасов и регулирования промысла.

Объектом исследований является промысловое многовидовое сообщество, популяции которого взаимодействуют трофически и находятся под влиянием промысла и абиотических факторов.

Цель работы. Разработать теоретические основы регулирования многовидовых промысловых сообществ c учетом межвидовых отношений посредством математического моделирования.

Основные задачи работы:

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

Х найти область допустимых режимов эксплуатации для сообщества простой структуры;

Х разработать продукционную модель двухвидового сообщества промысловых популяций: методы оценки параметров только по промысловой информации, настройка модели, анализ разных схем регулирования такой системы запас-промысел;

Х разработать многовидовую когортную модель MSVPA промыслового сообщества четырех видов Баренцева моря:

- провести модификацию модели MSVPA для сообщества Баренцева моря, позволяющую дополнительно вводить переменные параметры, - разработать алгоритмы автоматической обработки информации базы данных по питанию рыб, в том числе для определение возраста хищников и их жертв, - применить геостатистические процедуры для учета неравномерности распределения станций по сбору желудков и для заполнения пробелов в данных по питанию, применить эти процедуры к данным не только пространственной, но и временной структуры;

Х проанализировать современное состояние экосистемного моделирования промысловых сообществ;

Х разработать математическую модель формирования пополнения с учетом влияния экосистемных факторов на примере трески Баренцева моря;

Х разработать метод тестирования различных правил регулирования промысла на примере двух промысловых запасов.

Основные положения, выносимые на защиту:

- устойчивость модели к внешним воздействиям является одним из критериев ее адекватности, поскольку как правило моделируется сообщество, существующее в течение долгого времени;

- при регулировании многовидового промысла точечные граничные ориентиры управления следует заменять на границы допустимой области;

- разработанная автором динамическая продукционная модель промыслового сообщества (на примере ставриды и мерлузы ЮВА) позволяет оценить тип и интенсивность взаимодействия видов и является эффективным инструментом для сравнения разных стратегий промысла;

- когортная многовидовая модель, адаптированная к сообществу 4-х видов Баренцева моря, оценивает динамику промыслового сообщества и дает количественные оценки взаимоотношений видов;

- предлагаемый принцип сценарного моделирования, основанный на ретроспективных оценках запаса, эффективен для тестирования ПР запаса.

Научная новизна. Впервые автором проведено исследование на устойчивость по методу А.М. Ляпунова (1950) математических моделей (ММ) сообществ простой структуры (Булгакова, 1966,1968), получены условия устойчивости аналитических ММ сообщества с учетом трофических взаимоотношений разного типа. Показаны стабилизирующее влияние косного компонента на динамику экосистемы, влияние на динамику модельных популяций ограничения рациона хищника (по В.С. Ивлеву, 1955) и др. Доказана устойчивость элементарной ячейки модели с лимитирующими факторами.

Построена область допустимого управления и найдены режимы оптимальной эксплуатации двухвидового промыслового сообщества конкурирующих за общий ресурс популяций и популяций, взаимодействующих по типу "хищникжертва", причем как по критерию максимального вылова, так и по критерию максимальной экономической прибыли.

Впервые разработана динамическая продукционная ММ многовидового промысла, реализованная на примере двухвидового сообщества ЮгоВосточной Атлантики. По данным промысловой статистики оценены параметры модели и получены рекомендации по наилучшей стратегии промысла.

Модифицирована многовидовая когортная модель MSVPA для применения к промысловому сообществу Баренцева моря, исследована динамика этого сообщества, ориентированного на треску (северо-восточная арктическая треска, мойва, атлантическо-скандинавская весенне-нерестующая сельдь и северная креветка). На основе совместной базы данных (БД) по питанию хищников, созданной учеными России (ПИНРО) и Норвегии, разработаны новые компьютерные программы обработки данных по питанию трески. Среди них программы по определению возраста жертв в каждом из желудков, что удается делать и в ряде случаев, когда жертвы частично переварены (предварительно оценены коэффициенты реконструкции массы). Эта процедура позволила увеличить размер выборок желудков хищника.

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

С помощью модели MSVPA получены количественные оценки влияния многовидовых отношений на динамику популяций сообщества, которые необходимо учитывать при регулировании промысла.

Разработана ММ пополнения трески в зависимости от индекса структуры половозрелого запаса трески Баренцева моря, ее каннибализма и устойчивого притока Атлантических вод в Баренцево море, использующаяся как для прогнозирования численности пополнения, так и как часть имитационной модели популяции трески.

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

Предложен метод построения сценарных моделей популяции для тестирования различных правил регулирования промысла, основанный на анализе динамики параметров популяции в ретроспективе, на примере важных промысловых запасов - трески северо-восточной Атлантики и атлантическоскандинавской весенне-нерестующей сельди. Эти ММ позволяют выбрать наилучшее ПР и оценить вероятность риска нарушения биологических ориентиров (БО) управления по биомассе и промысловой смертности.

Достоверность результатов обусловлена четкой формулировкой гипотез, на основе которых строятся ММ, проверкой алгоритмов и программ на конкретных примерах, использованием на практике полученных автором результатов, использованием предложенных автором методов в международных организациях (ИКСЕАФ, рабочие группы ИКЕС).

Теоретическое и практическое значение работы. Теоретическое значение имеют результаты качественного исследования моделей экосистем на устойчивость, динамическая продукционная модель двухвидового сообщества.

Теоретическое значение имеют утверждение, что в многовидовом случае точечные БО управления следует заменять на границы области допустимого управления, и предложенный автором принцип построения сценарных моделей.

Практическое значение работы имеет модификация модели MSVPA, позволяющая использовать большее количество переменных параметров, предложенные автором алгоритмы и программы обработки информации базы данных (БД) по питанию хищников, а также результаты тестирования некоторых ПР для важных промысловых запасов, полученные с помощью сценарной модели, и модель формирования пополнения трески - для прогнозирования запаса.

Реализация результатов работы. Результаты, изложенные в диссертации, использовались на Рабочих группах Международного Совета по исследованию моря (ИКЕС): по многовидовым оценкам в 1983 (г. Гдыня, Польша), в 1995 г.

(г. Берген, Норвегия) и 1997 г. (г. Копенгаген, Дания), по арктическому рыболовству в 2000-2009 гг. (гг. Берген, Мурманск, Копенгаген, Виго, СанСебастьян); на Рабочей группе Международной Комиссии по рыболовству в ЮВА (г. Мадрид, 1987), в отраслевых НИИ (КамчатНИРО, МагаданНИРО, КаспНИРХ) при разработке годовых прогнозов, при работе с моделью MSVPA (1990-1996, ПИНРО), а также в учебном процессе - включены в Учебное пособие для студентов биологических специальностей ВУЗов, использованы на занятиях Отраслевого семинара по изучению методических основ рационального использования промысловых биоресурсов (п. Рыбное Московской обл., гг. Волгоград, Магадан, в 2000-2003 гг. и в 2009 г.) для повышения квалификации научных сотрудников отрасли. Автор была руководителем международного гранта ИНТАС-РФФИ 95-1204 "Многовидовой анализ сообщества Баренцева моря, ориентированного на треску" (1996-1998).

Имеются справки - акты о внедрении результатов исследований автора от КамчатНИРО, ПИНРО, КаспНИРХ, МагаданНИРО. Опубликованы методические рекомендации (в соавторстве).

Апробация работы. Результаты диссертации были представлены и обсуждены:

Х на заседаниях Ученого совета ВНИРО (многократно) и Ученого совета ПИНРО;

Х на заседаниях Ихтиологической Комиссии АН РАН (г. Москва, 2000, 2003);

Х на отраслевых, всероссийских и международных конференциях и совещаниях:

на V Всесоюзном совещании по морским млекопитающим (г. Махачкала, 1972); на XV Конвенции по котикам Северной части Тихого океана (г.

Москва, 1972); на отраслевой школе-семинаре по применению математических методов для прогнозирования уловов (п. Рыбное, 1983-1989, 2000-2003, 2009); на Всесоюзной школе-семинаре по методам анализа эксплуатируемых гидробионтов ВГБО АН СССР (г. Петрозаводск, 1987); на III Съезде советских океанологов (г. Ленинград, 1987); на Международной Конференции "Проблемы изучения и рационального использования биологических ресурсов окраинных и внутренних морей СНГ" (г. Ростов-на Дону, 1992); на Международной конференции "Современные проблемы гидроэкологии" (г. Санкт-Петербург, 1995); на Ежегодной научной Конференции ИКЕС (1990, 1995, 2003, 2004 - гг. Копенгаген, Таллинн, Виго);

на XVI Симпозиуме "Экосистемный подход к регулированию рыболовства" (г.Анкоридж, США, 1998); на VII Всероссийской конференции по промысловому прогнозированию (Мурманск, 1998); на Всероссийских конференциях по промысловой океанологии (гг. Калининград, 1999 и Светлогорск, 2002); на VIII съезде гидробиологического общества РАН (г.

Светлогорск, 2001); на Международной конференции Рациональное использование биологических ресурсов Мирового океана (г. ПетропавловскКамчатский, 2001г.); на встречах ученых России и Норвегии по теме "Многовидовые модели и управление промыслом в Баренцевом море" (г.

Берген, Норвегия, 1998) и по вопросам рыболовства в Норвежском и Баренцевом морях (г. Сванховд, Норвегия, 2004); на Международном Форуме по проблемам науки, техники и образования (г. Москва, 2002 и 2006); на Российско-Норвежском Симпозиуме (г. Тромсё, Норвегия, 2007); на двух Международных научных конференциях в г. Астрахани (2006 г. и 2007г.), а также на расширенном коллоквиуме подразделений ВНИРО: отделов методов исследований, промысловой ихтиологии, международного рыбопромыслового сотрудничества (2009 г.).

Материалы и методика В работе использованы методы качественного анализа решений систем дифференциальных уравнений А.М. Ляпунова; методы и модели теории рыболовства; многовидовое моделирование; новые разработанные автором алгоритмы и программы обработки базы данных по питанию рыб;

статистические и регрессионные методы оценки параметров моделей;

геостатистические методы обработки больших массивов данных; сценарное моделирование. Поскольку диссертация в большой своей части касается методических вопросов исследования биологических сообществ, подробное описание применяемых методов приводится в соответствующих главах.

Исходными материалами для главы 1 являются разработанные В.

Вольтерра (Volterra, 1931) модели взаимодействующих популяций по типу хищник-жертва и конкуренция, а также модели конкуренции видов, предложенные Г. Гаузе (1933,1935). Основой построения многовидовой динамической продукционной модели (гл.3 и 4) стала продукционная модель одной промысловой популяции в состоянии равновесия (Sсhaefer, 1954). При определении параметров многовидовой модели использованы данные промысловой статистики с 1974 г. по 1988 г. международного промысла в районе Намибии (Статистические сведенияЕ, 1975-1984; ICSEAF, 1976-1985;

ICSEAF, 1983), а именно: временные ряды общих международных уловов южноафриканской ставриды (Trachurus tr. capensis) и капской мерлузы (Merluccius capensis) ЮВА, ряды уловов и промыслового усилия на промысле ставриды, выбранного в качестве стандартного промыслового комплекса - советского судна типа БМРТ и разноглубинного трала, аналогичные данные выбранного в качестве стандартного промыслового комплекса на промысле мерлузы в том же районе Цсоветского судна типа БМРТи донного трала.

Исходной моделью для многовидового анализа сообщества Баренцева моря (гл.5 и 6) была когортная многовидовая модель MSVPA (Sparre, 1984; Gislason and Sparre, 1987; Sparre, 1991) в виде программного пакета в среде Фортран, которая модифицирована нами путем введения дополнительной возможности использовать переменные во времени значения рационов хищника и лостаточной смертности. Уловы каждого из четырех видов сообщества и их возрастной состав для каждого квартала каждого года из интервала 1980-20гг. рассчитаны на основе суммарных международных уловов, долей национальных уловов в общем вылове и возрастной структуры национальных уловов. При расчете возрастного состава квартальных уловов трески до 1994 г.

величина международного вылова распределялась по возрастам согласно процентному составу российских квартальных уловов, а начиная с 1994 г. - по данным российской и норвежской промысловой статистики. Распределение уловов сельди по кварталам получено из ПИНРО, возрастной состав уловов - из отчетов Рабочей группы (ICES 1996b; 1997b), его квартальная изменчивость не принималась в расчет. Данные по уловам мойвы и креветки и их возрастному составу получены из ПИНРО. Возрастной состав уловов креветки рассчитан по их размерному составу, с помощью условно принятого соответствия размерных и возрастных групп.

Входной информацией для MSVPA являются также данные по средней массе особи по возрастам, годам и кварталам, рассчитанные по размерно-возрастным ключам, полученным из ПИНРО (для трески и мойвы). Для сельди масса особи по возрастам принята одинаковой для всех кварталов года, причем масса особи в улове и в запасе тоже принята одинаковой. Для креветки приняты неизменные значения по годам и кварталам.

Среднемноголетние значения доли половозрелых особей по возрастам для каждого вида взяты из материалов Рабочих групп ИКЕС (ICES 1994; 1996a;

1997b; 1998). Для креветки половозрелая часть запаса не имеет большого значения и носит формальный характер из-за особенностей биологического развития этого вида.

Для оценки параметров питания трески использованы материалы РоссийскоНорвежской компьютерной базы данных по питанию рыб Баренцева моря (БД), начиная с 1984 г. Нами разработаны алгоритмы, оценивающие для каждого возраста хищников возраст жертв разных видов в желудках, средний за квартал вес пищи в желудке и средний состав пищи. Рационы оценивали по методу Б.

Богстада и С. Меля (1992), учитывая разную скорость переваривания разных видов жертв в желудке, полученную экспериментально для трески (dos Santos, 1990). Для улучшения качества данных по питанию применена методология кригинга (пакет WINSURF), параметры кригинга получены при анализе вариограмм (пакет GEO-EAS).

Настройка модели MSVPA (гл.6) для трески и сельди проводилась на результаты одновидовых оценок запасов (ICES, 1997; 1997b; 1998), при этом искали максимум корреляции между рядами численности возрастной группы популяции, полученными по одновидовой и по многовидовой моделям.

Настройку мойвы и креветки проводили на результаты учетных съемок.

Для модели формирования пополнения трески (гл.7) популяционная плодовитость трески рассчитана по методу М.Бондаренко и др.(2003), а данные по среднемесячной температуре воды в слое воды 0-200 м на станциях 3-7 на Кольском разрезе получены с сайта ПИНРО. Для построения сценарных моделей трески и сельди (гл.7) в качестве исходных данных взяты результаты прогонов одновидовых моделей на соответствующих рабочих группах ИКЕС (ICES, 2003; 2005; 2007; 2008).

Публикации. По теме диссертации опубликованы 69 работ в отечественных и международных изданиях, из них 33 статьи, в том числе 8 статей в журналах, рекомендованных ВАК.

Структура и объем диссертации. Диссертация состоит из введения, 7 глав, заключения и списка литературы из 295 наименований. Общий объем работы 340 cтр., в том числе 79 рисунков и 25 таблиц. Обзор литературы приведен в соответствующих главах.

СОДЕРЖАНИЕ РАБОТЫ

Во Введении приводится общая характеристика работы, использованные материалы и методы, даются определения основных понятий и перечень сокращений и обозначений.

Глава 1. Математические модели динамики экосистемы и исследование их устойчивости. С 20-х годов ХХ столетия начали развиваться теоретические основы современных моделей многовидового биологического сообщества (Volterra, 1928; Volterra & DТAncona, 1935), принимающие в расчет трофические взаимодействия между популяциями. Это модели однородных популяций, не дифференцированных по полу, возрасту, размерам особей.

Модель представляла собой систему обыкновенных дифференциальных уравнений, при этом задача определения параметров модели не ставилась.

Иногда удавалось найти точное решение такой системы, но чаще исследование модели ограничивалось качественным исследованием решения. Эти модели относятся к классу аналитических моделей.

Возможны два подхода к изучению экосистемы (ЭС) с помощью моделирования. При первом подходе, не углубляясь в строение каждой популяции, изучают влияние межпопуляционных связей на динамику ЭС. При втором подходе, условно оборвав эти связи, подробно изучают структуру отдельной популяции: с разной степенью детализации учитывают половую и возрастную структуру популяции или даже качественный состав биомассы одной особи (содержание в ней белка и жира), влияние кормовой базы на популяцию (Булгакова, 1978), процессы естественной и промысловой смертности и размножения, влияние абиотических факторов на динамику различных процессов (Булгакова, 2006) и т.д. Подробное описание всех этих факторов одновременно может привести к настолько сложной модели, что ее невозможно будет исследовать. Такую модель сложно обеспечить необходимой информацией, а поскольку каждый параметр определяется с некоторой ошибкой, то при большом количестве параметров ошибки накапливаются, и надежность полученных выводов может быть низкой.

В данной главе рассматриваем ММ, допускающие аналитическое исследование, популяции считаются однородными, при этом изучаем влияние на динамику их численности (биомассы) межпопуляционных связей разного типа. Одним из способов анализа адекватности аналитической модели является исследование на устойчивость ее решения. Как правило, реальная ЭС существует достаточно долго (в масштабах времени наблюдений). Если одна из компонент модели через некоторое время исчезает (ее численность стремится к нулю), т.е. нарушается структура экосистемы, то можно считать, что модель не соответствует оригиналу. Аналогично, если численность какого-нибудь вида сообщества в модели со временем неограниченно возрастает, это тоже говорит о несоответствии модели экосистеме, поскольку в природе всегда существуют факторы, ограничивающие рост численности. Согласно концепциям термодинамики, любая замкнутая система с протекающим сквозь нее потоком энергии (будь то биосфера или маленькое озеро) с большой вероятностью развивается в сторону некоторого устойчивого состояния, и в ней с необходимостью должны вырабатываться саморегулирующие механизмы (Свирежев и Логофет, 1978).

Устойчивость здесь понимается в смысле А.М.Ляпунова (1950) и означает, что при малых отклонениях системы от состояния равновесия система стремится к нему вернуться. Таким образом, существование устойчивого равновесия обеспечивает и структурную устойчивость ЭС (сохранение всех видов сообщества).

Рассматриваются математические модели (ММ) трофической цепи (трофической пирамиды), когда особи каждого последующего звена (численностью xi) питаются только особями предыдущего звена (xi-1), а сами являются пищей только для особей последующего (xi+1), здесь i =1,2,Еn; n - число звеньев (рис.1).

Рис.1 Блок схемы трех трофических цепей, соответствующих: модели Вольтерра (Volterra, 1931) (а), модели n видов с обратной связью на косную компоненту типа запас (b), с косной компонентой типа поток (c) (b и c - по Булгаковой,1966) Анализ на устойчивость нескольких моделей трофических цепей показал следующее (Булгакова, 1966).

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

В предложенной нами модели сообщества, где запас пищи для организмов первого трофического уровня ограничен (количество биогенов в морской системе или количество гумуса в почве для наземного сообщества), количество биогенов x0 составляет нижнее звено, которое пополняется за счет потоков трупов и выделений из всех остальных звеньев, и эти потоки пропорциональны n n биомассе соответствующих звеньев: + si )xi = pixi. Модель имеет вид:

( i i=1 i=dx0 n = pi xi - 1x0x dt i=dxi = (ixi-1 - si - i+1xi+1)xi, где i = 2,3,...,n; n+1 = 0.

dt Показано, что такая обратная связь потоков вещества на нижнее звено стабилизирует динамику сообщества, а модель с учетом косных компонент является более реалистичным отражением природной ЭС. При ослаблении прямой связи первой компоненты с косной (уменьшении параметров межвидовой связи) усиливаются колебания численностей, тот же эффект получается при ослаблении обратной связи, особенно при ослабления потока вещества с первого вида (трофического уровня) на косный компонент.

Рассмотрена структура модели ЭС (рис.1с), когда для первой компоненты ЭС (например, для фитопланктона) лимитирующим ее развитие фактором является поток некоторого необходимого вещества (например, биогенов) или поток солнечной энергии постоянной интенсивности P (Булгакова, 1966). Это означает, что при увеличении численности или биомассы первой компоненты (потребителя этого вещества или энергии) на единицу ее плотности будет приходиться меньшая доля вещества или энергии. Абиотический фактор типа поток тоже стабилизирует динамику ЭС.

Введение временной задержки - длительности времени разложения вещества, поступающего со всех трофических уровней на косный компонент, до биогенов, приводит к задержке по фазе затухающих колебаний относительно варианта с =0. (Численные расчеты проведены при n=3 и =3).

В модель двух видов хищник-жертва Вольтера вместо пропорциональной зависимости рациона хищника от плотности корма нами введена обобщенная эмпирическая зависимость по В.С. Ивлеву (1955). Условие устойчивости точки равновесия такой модели выполняется только при наличие внутривидовой конкуренции первого вида. При ее отсутствии ограничение рациона хищника не позволяет ему регулировать численность жертв так, как при неограниченном рационе в модели Вольтерра.

ММ конкуренции видов исследовались неоднократно (Volterra, 1931; Гаузе 1933, 1934; Hutchinson, 1947, 1948; Canningham, 1955; Rescigno and Richardson, 1965; Логофет и Свирежев,1975). В этих MM фактор, лимитирующий развитие популяций, не рассматривается в явном виде, а в нашей модели (Булгакова, 1968а) конкуренция между двумя видами-потребителями ведется за корм, пополняемый за счет постоянного притока. Условие устойчивости равновесия этой модели - существование внутривидовой конкуренции.

Глава 2. Управление промысловым сообществом в условиях равновесия.

Цель этой главы - на простых ММ показать, как именно учет трофических взаимодействий между промысловыми популяциям влияет на выбор режима промысла двух видов. Разработанный в рамках ИКЕС в последнее десятилетие предосторожный подход к управлению природными запасами (Jakobsen, 1992, 1993; ICES 1996c; ICES 1997d; Бабаян, 2000) относится к одновидовым запасам.

Отметим, что большая часть наших исследований, освещенных в этой главе, проведена раньше, чем появилось понятие предосторожного подхода (Булгакова, 1970, 1977, 1998).

2.1. Анализ допустимой области управления в системе хищник-жертва Естественно предположить, что эксплуатируется ЭС, реально существующая в природе в течение долгого времени, потому выбрана простейшая модель хищник-жертва, имеющая устойчивую точку равновесия, а именно dx= x1(a1 - bx1 - cx2 ) - F1xdt, (1) dx= x2(a2 x1 - M ) - F2xdt где x1 - численность вида-жертвы, x2 -численность вида-хищника, a1 и b - параметры роста популяции жертв в условиях лимитирующей среды, коэффициент смертности от хищников cx2 представлен УвольтерровскимФ членом; коэффициент воспроизводства хищников a2x1, пропорционален численности вида жертв, M - постоянный коэффициент естественной смертности хищников, F1 и F2 - коэффициенты промысловой смертности.

Согласно этой модели второй вид при отсутствии жертв вымирает.

Если сообщество не эксплуатируется, точка равновесия системы (1) лежит в положительной области плоскости {x1,x2} при условии D = a1a2 - bM > 0.

Считаем, что это условие выполняется. Если промысел ведется с неизменными во времени параметрами F1 и F2, координаты точки равновесия системы (1) принимают вид M + F2 D - (a2 F1 + bF2 ) X10 = ; X =. (2) a2 20 a2c Точка (2) лежит в положительной области, если коэффициенты F1 и F(параметры управления сообществом) удовлетворяют условию (3):

a2F1 + bF2 < D. (3) Рис.2. Область допустимого управления промыслом на плоскости {F1,F2} для модели системы хищник-жертва (1).

При выполнении условия (3) и при b>0 (существование внутривидовой конкуренции) эта точка асимптотически устойчива (Bulgakova, 1999). При b=система (1) превращается в простейшую систему уравнений Вольтерра, решение которой имеет вид незатухающих колебаний.

Область допустимого управления сообществом (рис.2) построена на плоскости {F1,F2} Цнезаштрихованный треугольник MON вместе с границами ON и OM, но без границы MN. Уравнение прямой MN получается из (3) при замене знака неравенства на знак равенства: F1=D/a2-F2b/a2.

Рассмотрены разные варианты регулирования промысла. Пусть жертвы - ценный промысловый вид, и цель промысла - максимизировать ее равновесный вылов (критерий MSY). Если менеджер пренебрегает влиянием хищника, то он оценивает граничный БО по промысловой смертности по принятой терминологии (ICES,1996) как F1lim=a1 и рекомендует выбрать F*1MSY=a1/2, соответствующий точке S на рис.2, ожидая получить вылов MSY*=a12/4a2 (знак * показывает, что формула относится к модели изолированной популяции). Но поскольку хищники существуют, величину граничного ориентира следует снизить: при отсутствии промысла хищников до уровня F1lim = D/a2 < a1, а при увеличении F2 величина F1lim линейно уменьшается. Линия MN соответствует соотношениям между параметрами F1 и F2, при которых популяция хищников будет уничтожена, т.е. линия MN - пример граничной линии в области допустимого управления сообществом и соответствует управлению, при котором структура ЭС разрушается.

Функция уравновешенного вылова первого вида за год Y10 (F1,F2) = F1X10 = F1(M+F2)/aне имеет максимума и возрастает с ростом F1 и F2 до пересечения с плоскостью, перпендикулярной плоскости {F1,F2} и проходящей через линию MN. При этом возможны 2 варианта:

1)Если a1/2

2)Если a1/2>D/a2, т.S лежит за границей допустимой области выше т.М, т.е. попытки регулировать первую популяцию по критерию MSY для изолированной популяции приведут к исчезновению хищника даже при отсутствии изъятия хищников! Этот случай - самый опасный для многовидовой системы запас-промысел. Наибольший вылов первого вида (D*M/a22) соответствует т.М при F1=D/a2, F2=0. Но для сохранения ЭС следует выбрать точку внутри допустимой области.

Пусть менеджер заинтересован в получении вылова хищников.

Уравновешенный годовой вылов хищников является функцией F1 и F2:

F2 ( D - a2 F1 - bF2 ) Y20 = F2 X = a2c.

При фиксированном F1>0 наибольший вылов хищников достигается на прямой (D - a2F1)ML (медиане угла OMN) при F2=(D-a2F1)/2b и равен. Максимальный 4a2bc же вылов хищников можно получить при отсутствии промысла жертв (F1=0) в т. L, которую можно выбрать в качестве целевого БО при промысле хищников.

При предосторожном подходе выбираем точку, лежащую на ON левее т. L.

Проблема оптимального управления сообществом в целом приводит к необходимости выбора вида критериальной функции, например, в виде взвешенной суммы уловов двух видов. Более общий вид критериальной функции получается при выборе ее в виде зависимости годовой экономической прибыли промысла двух видов в равновесном режиме от значений F1 и Fвнутри области OMN. В простейшем случае, если два специализированных промысла ведутся независимо, экономическая прибыль выражается формулой P(F1,F2) = A1F1(M + F2)/a2 + A2F2(D-a2F1-bF2)/ca2 ЦB1F1-B2F2 (4) Параметры A1 и A2 - стоимость единицы вылова I и II популяции, а B'1 и B'2 - затраты на единицу промыслового усилия E1 и E2 соответственно; тогда при Fi=qiEi (i=1,2) при постоянных qi параметры B = Bi' / qi BiЦпостоянные величины.

Показано (Булгакова, 1977; 1999), что функция прибыли (4) не имеет максимума в допустимой области, и в зависимости от соотношений параметров модели ее наибольшее значение достигается на одной из границ области.

Приводится алгоритм его нахождения. В частном случае при A1=A2 и B1=B2=функция (4) представляет собой сумму уловов двух видов.

2.2 Анализ допустимой области управления в системе двух конкурирующих видов. Конкуренция возможна лишь в среде с лимитирующими ресурсами, потому в качестве модели популяции берется модель логистического роста:

dxi/dt = xi ai/Ki ( Ki -xi - bi xj ) - Fi xi, i =1,2 (5) где Ki - максимальная численность i-го вида в отсутствии конкурирующего вида, Fi Цмгновенные коэффициенты промысловой смертности, ai - коэффициенты прироста численности популяции. Параметры конкуренции bi имеют следующий смысл: xj oсобей j-го вида потребляют столько же ограниченного ресурса, сколько bixj особей i-го вида. Координаты точки равновесия (5) зависят от промысловой смертности обоих видов, причем равновесная численность каждого из видов увеличивается с ростом промысловой смертности конкурента.

aFR Rpr U A F2 aРис.3. Область допустимого управления промыслом на плоскости {F2,F1} для системы двух конкурирующих за общий ресурс популяций.

Область допустимого управления (рис.3) ограничена осями координат и двумя прямыми граничными линиями, на которых система разрушается. Граничная линия для первого вида соответствует x10=0, т.е. при F1 и F2, соответствующих точкам на этой линии, первый вид исчезает, а с ростом F2 увеличивается допустимое значение F1 и наоборот. Если бы виды развивались независимо, допустимая область управления была бы шире, ограничиваясь для первого вида величиной F1lim=a1, а для второго - F2lim=a2. Равновесный вылов каждого вида зависит от промысловой смертности обоих видов. При фиксированном Fмаксимум вылова первого вида достигается на MSY-линии (рис.3) и растет с ростом F2, вплоть до т. R на граничной линии второго вида. Это означает, что управление, направленное на получение максимального вылова одного из видов, приводит к вытеснению второго, т.е. к разрушению ЭС. Чтобы сохранить ЭС, нужно выбирать точку управления внутри допустимой области, например, на пунктирной линии AR (т.Rpr на рис.3), аналогично тому, как это делается для одновидовых моделей промыслового запаса (Бабаян, 2000).

Функция экономической прибыли при промысле двух конкурирующих видов имеет тот же вид, что (4). В данном случае, поскольку равновесные выловы Y10 и Y20 - полиномы второй степени от F, поиск максимума функции прибыли внутри области допустимого управления представляет собой задачу нелинейного программирования. Если ценным является только один вид, наибольшее значение прибыли следует искать на границе области, если же оба вида ценные, возможно существование максимума прибыли внутри области управления. Сформулированы условия существования этого максимума.

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

В данной главе рассматриваются упрощенные модели, структура реальных ЭС гораздо сложнее. Например, в реальной ЭС при недостатке одного вида корма хищник переключается на другие, т.е. опасность его исчезновения не так велика. Включение в модель альтернативной жертвы (Spencer, Collie, 1995) делает модель более реалистичной, но значительно усложняет анализ.

Глава 3. Исследование многовидового промысла с помощью продукционных моделей. При решении задач оценки запасов и регулирования промысла возникает необходимость в таких ММ, параметры которых могут быть реально оценены на основе доступной информации. Продукционные модели одной изолированной популяции предназначались для исследования равновесных состояний запаса (Schaefer, I954, Fox, 1970), для оценки параметров модели в качестве исходной используют только промысловую информацию: данные об уловах Y и промысловом усилии E за год. По этим данным рассчитывается величина улова на единицу промыслового усилия Y/E (часто обозначается СPUE и считается индексом промыслового запаса). Позже появились одновидовые динамические продукционные модели (Butterworth and Andrew, 1984; Babayan and Kizner, 1988). На базе модели Шефера построены и модели многовидового промысла: модель промысла, эксплуатирующего одновременно запасы двух биологически не взаимодействующих видов (Horwood, I976) или двух взаимодействующих видов (Pope and Harris, 1975). В последнем случае виды эксплуатируются независимыми промыслами, и в качестве входной информации использовались временные ряды, полученные по одновидовому ВПА при условии, что система находится в состоянии равновесия. Такой метод оценки параметров не совсем корректен, поскольку метод ВПА определяет динамику запасов, что противоречит принятому условию равновесия системы, а с другой стороны, ВПА не учитывает межвидовых отношений и дает смещенные оценки запаса.

Нами разработана динамическая продукционная модель двух взаимодействующих эксплуатируемых запасов P1и P2 (Булгакова, Кизнер, 1986, 1987a; 1987б; Bulgakova, Kizner, I986), которая записана в виде dP1/dt =P1(t)(A1- B1P1 - C1P2)-Y1 (6) dP2/dt =P2(t)(A2- B2P2 + C2P1)-Y2, где константы Ai, Bi 0 (i = 1,2) характеризуют внутрипопуляционные процессы, а параметры C1 и C2 - межвидовые отношения. При Сi> 0 второй вид является хищником по отношению к первому, а при условии (С1> 0, С2< О) виды являются конкурентами. Могут существовать отношения между видами сразу двух типов, тогда знаки C1 и C2 говорят о том, какое из них превалирует. С помощью соотношений ui=Yi/Ei=qiPi от (6) переходим к уравнениям du= a1 - b1u1 - c1u2 - q1Edtuотносительно индексов запасов: (7) du= a2 - b2u2 + c2u1 - q2Edtuгде ai=Ai, bi=Bi/qi, с1=C1/q2, c2=C2/q1. Оценив параметры ai, bi, ci, qi, можно восстановить и величины параметров Аi, Вi; Сi исходной системы (6). В (7) E1 и E2 Цстандартизированное усилие промыслов первого и второго вида, u1 и u2 - уловы на единицу такого усилия, соответственно. В качестве входной информации для параметризации модели используется информация по трем временным рядам для каждого вида: ui(t), Ei(t) и относительное изменение dui ui (tk ) - ui (tk -1) = улова на единицу усилия за год. (8) uidt ui Именно использование временных рядов типа (8) делает наш подход динамическим. Все три указанных временных ряда рассчитываются до параметризации модели. Поскольку исходные данные содержат случайные ошибки, фильтрация данных является необходимым этапом первичной обработки данных.

При параметризации модели следует минимизировать функционалы n k = i (9) [ui (t ) - u (t )]j j i=1 j= n k 1 = i или (10) [vi (t ) - v (t )]2, j j i=1 j=где n - количество видов в сообществе, k - длина рядов (количество лет), ui(tj) и u (tj ) i ряды заданных и расчетных значений уловов на единицу усилия, vi=dui/dtui и v (t ) - ряды значений относительных скоростей изменения уловов i j на единицу усилия. Расчетные значения получаются при решении системы уравнений (7) с уже оцененными параметрами в процессе поиска минимумов.

Задача настройки модели, состоящая в минимизации функционалов (9,10), разрешима при условии, что kmn (mЦколичество неизвестных параметров модели), т.е. нужны достаточно длинные ряды исходных данных. Поэтому предложено (Булгакова, Кизнер,1987) заменить функционалы (9) и (10) на n функционалов вида k k i = [ui (t ) - u (t )]2 или i ' = [vi (t ) - v (t )]2 при i=1,2,Е,n, (11) i i j j j j j=1 j=и решать задачу параметризации для каждого уравнения отдельно. Тогда условие ее разрешимости имеет вид km. Если же оценить параметры qi независимо, число неизвестных параметров уменьшится, тогда в качестве третьего ряда удобно рассматривать величины dui/dtui+qiEi вместо (8).

Описанные приемы настройки можно перенести и на сообщество более сложной структуры.

Глава 4. Рекомендации по регулированию промысла ставриды и хека в ЮВА, основанные на двухвидовой динамической продукционной модели.

Динамическая модель (10) использована для описания промыслового сообщества южноафриканской ставриды (Trachurus tr. capensis) и капской мерлузы (Merluccius capensis) (хека) ЮВА (Булгакова и Кизнер, 1985; 1987а).

Промысловая статистика использована с 1974 по 1983г., когда ставриду ловили главным образом пелагическими тралами (с приловом хека), а хека - донными (с приловом ставриды). Рассмотрены вопросы подготовки входных данных, в том числе стандартизация промыслового усилия по Галланду (Gulland, 1969), оценка параметров, анализ чувствительности к изменениям ее параметров и вопросы прогнозирования при разных вариантах регулирования промысла. В качестве стандартных выбраны исходные ряды советских данных о промысле ставриды комплексами БМРТ - разноглубинный трал и хека - комплексами БМРТ - донный трал за 1974-1983 гг.

Входные ряды данных сглаживали, затем аппроксимировали полиномами от времени. По этим полиномам рассчитаны значения каждого ряда с полугодовым шагом, затем для каждой точки оценены величины уловов на усилие. Настройка модели проводилась в несколько этапов. Оценки запасов ставриды и хека в районах 1.3+1.4 ИКСЕАФ (Babayan and Bulgakova, 1983a;

1983b; Bulgakova et al., I984; Babayan et al, 1985) сопоставлены с уловами на усилие, в результате получены приближенные оценки: q1= 0,9*10-6 и q2 = 2,5*10-6, что позволило уменьшить число неизвестных параметров модели.

Результаты последнего этапа подготовки данных служат для оценки начальных приближений параметров модели (10). Процедуру оценки параметров можно выполнить строго только в случае независимости троек "гладких" переменных {u1, u2, E1) и {u1, u2, E2}. Однако существует достаточно сильная линейная связь между усилием и уловом на усилие. Начальные приближения параметров при заданных q1 и q2 определены из (10) с помощью программы линейной множественной регрессии:

a1= I,76; b1=0,00; c1=1,00; a2= 0,56; b2=0,34; c2= 0,18.

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

Коэффициент внутривидовой плотностной регуляции запаса отчетливо выражен для хека, что подтверждает биологический смысл параметров - ведь каннибализм хека в этом районе широко известен (Lleonart et al., 1985;

Konchina, 1986; Crawford et al., 1987). Параметры c1 и c2 получены положительными, что означает, что хищничество хека по отношению к ставриде является более сильным фактором, чем конкуренция, хотя некоторые ученые (Krzeptowski, 1982) указывали на совпадение состава пищи капского хека и ставриды, основой которой являются эвфаузииды и миктофиды, и на этом основании делали вывод о конкуренции между ними. По данным же Ю.В.

Кончиной (Konchina, 1986), в 1985 г. ставрида составляла до трети пищи хека длиной 24-44 см и до 100% пищи хека длиной 32-40 см, что говорит в пользу более сильного влияния хищничества.

Тестирование модели. В (10) подставлены зависимости E1 и E2 от времени в виде полиномов, приняты начальные значения ui, соответствующие середине 1976г., затем уравнения (10) решались методом Рунге-Кутта. Варьировали значения оценок коэффициентов уравнений в некоторых диапазонах и для каждой их комбинации численно определяли величины u1 и u2. Варианты сравнивали по среднеквадратическому отклонению расчетных величин du1/duв исходных точках от сглаженных фактических значений. Наилучшей оказалась комбинация: a1=I,70, b1=0,05, c1= 0,9, a2=0,5, b2=0,35, c2=0,2.

Результаты тестирования показали (рис.4), что модель способна описывать главные закономерности в динамике запасов. В окончательном варианте параметр b1 больше нуля, хотя и значительно меньше b2.

Анализ чувствительности показал, что самая высокая чувствительность результатов к изменениям a1 и a2, меньше - к b2.

Исследование различных стратегий управления двухвидовым промыслом показало (Булгакова, Кизнер 1987), что данная динамическая система при выборе стратегии постоянных годовых уловов является неустойчивой. Это означает, что если в некоторый момент времени система не находится в равновесном состоянии, а промыслом берется один и тот же вылов из года в год, возникают колебания численности популяций с возрастающей амплитудой, которые через некоторый промежуток времени могут разрушить ЭС. Такая стратегия промысла опасна для популяций - под действием внешних факторов система запас-промысел выходит из состояния равновесия, а в реальной ситуации только устойчивая система может существовать длительное время. Если же промысел ведется с постоянным промысловым усилием (которое соответствует области допустимого управления), то всякое равновесное состояние системы будет устойчивым.

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

Рис. 4. Фазовая диаграмма сообщества ставрида-хек: 1 - сглаженные фактические данные; 2 - результат тестирования модели (7); тонкие пунктирные линии - изолинии равновесных уловов (в I06 т в год), крестики - исходные данные.

Глава 5. Многовидовая модель промысла, построенная на основе возрастного состава уловов. Основное положение многовидового когортного анализа (MSVPA) заключается в том, что для каждого вида сообщества коэффициент естественной смертности подразделяется на две компоненты Mia = M 2ia + M1ia, (i Циндекс вида, а - индекс возраста). Первое слагаемое - коэффициент смертности от хищников, второе - смертность от всех других причин или остаточная естественная смертность. Для каждого вида сообщества параметры M1 должны быть оценены вне модели, а параметры M2 Цв рамках модели. Модель MSVPA (Daan, 1973; Pope, 1979; Helgason&Gislason, 1979;

Sparre, 1980, 1984, 1991) применялась сначала к сообществу Северного моря. В модель включены промысловые виды, находящиеся в трофических отношениях: виды-хищники и виды-жертвы. Входная информация для первых двух категорий дифференцирована по возрастам, годам и кварталам и включает возрастной состав уловов, массу особи, темп полового созревания, параметры М1. Для хищников дополнительно требуется подробная информация по количеству и составу пищи, переменная по возрастам и временным интервалам.

Один и тот же вид в модели может выступать и как жертва в младших возрастах и как хищник.

В главе описаны состав входной и выходной информации модели, уравнения и алгоритмы расчетов. Модель была нами несколько модифицирована для применения к ЭС Баренцева моря, в этом в рамках международного гранта принял участие Х. Гисласон (Дания). В традиционной модели не рассматривались сезонные (межквартальные) различия параметра M1, а рационы хищников брались разными для разных кварталов, но неизменными по годам. Значительная изменчивость рационов трески по годам (Булгакова и др., 1995; Bulgakova et al., 1995a) заставила ввести в программу возможность задавать их переменными по годам и кварталам. Для мойвы, с ее высокой смертностью после нереста, M1 берем разными для разных кварталов.

Глава 6. Исследование функционирования промыслового сообщества Баренцева моря. В состав моделируемого сообщества включены виды, составляющие "ближайшее окружение" трески (Gadus morhua L.): помимо трески, вошли мойва (Mallotus villosus), северная креветка (Pandalus borealis) и атлантическо-скандинавская сельдь (Clupea harengus L.) - основные компоненты питания трески. Кроме того, моделируется каннибализм трески (молодь трески в возрасте от 0+ до 4 лет потребляется особями старших возрастов). В модели возраст трески меняется от 0+ до 15 лет, сельди от 0+ до 12 лет (треской потребляется сельдь в возрасте от 0+ до 4 лет), мойвы и креветки Цот 0+ до 6 лет. Для креветки приняты условные возраста, соответствующие разным размерным группам. Ретроспективные расчеты проводились для интервала лет с 1980 г. по 1996 г. с квартальным шагом.

В модели используется информация двух типов по питанию особей каждого возраста трески: квартальные рационы и средняя за квартал масса содержимого желудка и его состав (вид и возраст жертв), все эти данные изменяются по годам и кварталам. Информация по составу желудков регулярно, начиная с 1984 г., заносится в совместную Российско-Норвежскую базу данных по питанию (БД). Эта информация была любезно передана нам сотрудниками ПИНРО. Для расчета входных параметров модели по питанию нами разработаны специальные алгоритмы и программы (Bulgakova et al., 1995b;

Bulgakova and Vasilyev, 1995; Булгакова и др., 1998). В БД возраст хищника был определен не для всех особей, потому разработана программа в среде PARADOX, которая определяет возраст хищника с помощью размерновозрастных ключей трески. Поскольку для всех записей БД первоначально неизвестен возраст жертв в желудках, для его определения заранее рассчитаны коэффициенты реконструкции массы разных видов жертв в зависимости от степени переваривания, а затем специальная программа позволяет оценить возраст жертв для всех записей БД, за исключением тех, для которых код длины жертвы не определен и одновременно либо неизвестно количество жертв данной категории в желудке, либо неизвестна степень их переваривания. На следующем этапе для каждого года и квартала оценивали средний за квартал вес содержимого желудка трески и его состав, под составом понимается вес каждой категории видов-жертв (категория - это сочетание вида и возраста).

Для увеличения надежности оценок состава пищи хищников применена методология кригинга, позволяющая учитывать пространственные аспекты информации, сглаживать информацию и заполнять пробелы в данных.

(Булгакова и Васильев, 1998; Bulgakova, Vasilyev, Daan, 2001; Булгакова и др., 2001). Впервые эта методология применена нами к массивам данных не только пространственной, но и временной природы (год, возраст).

Расчет суточных рационов трески проводился по методу (Bogstad and Mehl, 1992), в основу которого положены предположения об экспоненциальной модели переваривания, когда скорости процесса переваривания каждого вида жертв постоянны и оценены для трески Баренцева моря экспериментально (dos Santos, 1990). Детально результаты расчетов рационов обсуждаются в работе Булгаковой и др.(1995). Изменчивость рациона по кварталам и по годам значительна и сильно зависит от состояния кормовой базы. Из четырех кварталов хуже всего в БД представлен четвертый. Доля пустых желудков самая большая в I и II кварталах. Полученные нами оценки относительного суточного рациона лежат в интервале от 0,4% до 2,0% и близки к оценкам А.В.

Долгова с соавторами (Dolgov et al., 1992) - 0,54-2,2%. В ряде работ, где применялись другие методы оценки рационов или расчет велся для отдельного района или определенного периода откорма трески, оценки рациона иногда превышали 4% (Тарвердиева и Ярагина, 1989; Orlova, 1992).

По аналогии с одновидовым ВПА, в модели MSVPA для каждого вида требуется задать терминальные значения параметра F по возрастам для IV квартала последнего расчетного (терминального) года, а для старшей возрастной группы - для всех расчетных лет. Оценка этих параметров и является целью настройки многовидовой модели. Настройка MSVPA проводилась на оценки запасов, полученные по ВПА или с помощью съемок.

Предложен и новый метод настройки, не требующий дополнительной информации (Vasilyev, Bulgakova et al., 1999; Булгакова и Васильев, 1999б;

Васильев, 2000). Терминальные F и коэффициенты М1 оцениваются при поочередном применении MSVPA и квартальной версии одновидовой сепарабельной когортной модели QISVPA. Этот метод применим только для сельди и трески, у которых только младшие возрастные группы подвержены хищничеству. Мойву и креветку настраивали на результаты учетных съемок, варьируя M1 и терминальные F так, чтобы минимизировать медиану квадратов отклонений логарифмов оценок биомассы, полученных из съемок и из MSVPA.

Результаты основного прогона модели (рис.5) показали, что потребляемая треской биомасса мойвы и креветки в некоторые периоды много выше промыслового изъятия этих видов. Велик пресс каннибализма - параметр Mдля молоди трески усиливается в годы низкого запаса мойвы и при росте запаса трески начиная с 1992г., для возрастной группы 1 достигает 0,6 в 1996 г..

Видна главная роль мойвы в питании трески: при низком ее запасе усиливается потребление других видов и другой пищи.

треска креветка треска 800 30800 45662030441015220 0 1980 1984 1988 1992 1996 1980 1984 1988 1992 19мойва сельдь мойва 80140161206012100804086020404201980 1984 1988 1992 191980 1984 1988 1992 19Рис. 5. Результаты основного прогона MSVPA для каждого вида: динамика биомассы общего запаса (TSB, линия с ромбиками), биомассы, потребленной за год треской (пустые квадратики) и вылова (сплошная линия), 1000 т.

Смертность мойвы от хищника особенно велика для возрастных групп 3, 4 и 5 в годы низкого запаса мойвы и при росте запаса трески, M2 мойвы в возрасте 4 г. достигает 1,7 [1/год]).

Для проверки чувствительности модели к выбранному набору данных по составу желудков проводились прогоны с разными такими наборами: для всего указанного периода (1984-1995) (основной прогон), для периода 1986-1987гг., когда мойва была в депрессии, и для периода 1990-1992гг., когда запас мойвы восстановился. В основном прогоне рационы трески варьировали по кварталам и годам, а в остальных прогонах в расчетах использованы переменные по кварталам, но усредненные за выбранный период рационы. Оказалось, что TSB TSB вылов вылов Потребление хищником и Потребление хищником и TSB и вылов и вылов Потребление хищником TSB, потребление хищником усреднение квартальных рационов по годам сильно сглаживает динамику потребленной хищником биомассы каждого вида жертв. Например, в годы с низким запасом мойвы отклонения потребленной биомассы мойвы от основного расчета достигают 80-90%. Для креветки наблюдается та же картина.

Если взять за основу данные по составу желудков в годы, когда было много мойвы, то для всех лет расчетного периода возрастают оценки суммарной биомассы потребленной пищи и доля мойвы в ней. Оценки запаса мойвы получаются сильно завышенными по сравнению с основным расчетом, например, для 1991 г. - в 5 раз. Если же взять за основу данные по составу желудков за годы, когда было мало мойвы, получаем еще более сложную картину. Период 1986-1987гг. характеризуется не только низким запасом мойвы, но и выделяется составом пищи трески, потому распространение данных по составу желудков этих лет на весь период исследований дает сильные и неожиданные смещения оценок.

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

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

Для совершенствования эксплуатации многовидовых систем необходимы следующие этапы работы, представляющие собой экосистемный подход к оценке запасов и регулированию промысла:

1. определить, какие факторы среды оказывают наиболее значительное влияние на ЭС или на отдельные ее элементы, разработать ММ, включающие влияние факторов среды, прогнозировать их значения на несколько лет вперед, 2. сформулировать цель регулирования промысла основных компонентов сообщества, определить границы области допустимого управления в координатах коэффициентов промысловой смертности, разработать схему регулирования промысла для сообщества, которая должна включать в себя план восстановления некоторых единиц запасов, находящихся в угнетенном состоянии, и учитывать неопределенность параметров управления и запаса.

Одним из важнейших процессов, влияющих на динамику промысловых запасов рыб, является процесс формирования пополнения. В качестве примера включения факторов среды в ММ популяции разработана модель пополнения баренцевоморской трески (Булгакова и Борисов, 2001; Булгакова, 2006). За основу взята модель запас-пополнение У. Рикера (Ricker, 1954), включающая каннибализм, но зависимость численности пополнении от биомассы нерестового запаса заменяется на его зависимость от популяционной плодовитости, что позволяет учесть влияние возрастной структуры нерестового запаса. Делается предположение, что на выживание молоди влияет интенсивность и продолжительность притока атлантических вод, так как образующиеся при этом вихри вдоль основных потоков создают благоприятные условия для развития кормового планктона и личинок трески (по рекомендации А.А.Елизарова, ВНИРО). Индекс притока атлантических вод N(Y) рассчитан по данным Кольского разреза в слое 0-200м как количество месяцев в год рождения поколения Y с положительными температурными аномалиями.

Уравнение модели имеет вид:

Rmod(Y+3) =IndPF(Y)*EXP[ - a - Can*SSB(Y)+ny*N(Y)], (12) где Rmod (Y+3)Цчисленность пополнения в возрасте 3г., IndPF(Y) - индекс популяционной плодовитости; SSB(Y) - биомасса нерестового запаса.

Параметры a, Can и ny - оцениваемые константы. Выражение под знаком экспоненты - коэффициент естественной смертности поколения до возраста 3 г.

со знаком минус. Параметры модели оценены с помощью временных рядов численности пополнения и нерестового запаса трески для поколений 1984-20гг., полученных при оценке запаса трески на рабочей группе ИКЕС по арктическому рыболовству, рядов популяционной плодовитости (Бондаренко и др., 2003) и N(Y). Модель позволяет объяснить более 70% изменчивости пополнения, она используется при сценарном моделировании популяции трески и прогнозировании пополнения.

Разработаны имитационные модели для тестирования разных ПР для двух запасов СВА - трески и сельди (Bulgakova, 2003, 2004; 2008). Расчеты проводим по ретроспективному периоду, для которого известны многие параметры, в том числе факторы среды, и моделируя введение тестируемого ПР, получаем разные "сценарии" динамики характеристик запаса. Модель включает блоки динамики численности популяции, формирования пополнения, схемы регулирования, абиотических факторов, датчика случайных чисел.

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

По модели для трески тестировали два ПР. Первым было предосторожное правило ИКЕС (ICES, 1997c) со следующими БО: Blim=220 тыс. т, Bpa=460 тыс.

т, Flim=0,74 и Fpa=0,40. Согласно этому ПР рекомендуемое значение промысловой смертности Frec(t) для года t рассчитывается по формуле:

SSB(t) - Blim Frec (t) = Fpa. (13) Bpa - Blim Другое ПР, называемое трехлетним ПР, предложено Совместной РоссийскоНорвежской Комиссией в 2002 г. с целью стабилизации ОДУ и вылова трески.

Порядок применения этого ПР следующий: в год t оценить ОДУ для трех последующих лет, основываясь на (13). Оценка ОДУ для первого года приравнивается к средней из полученных. В следующем году повторяется процедура со сдвигом на один год, но ОДУ не должно отличаться от предыдущего более чем на 10% по модулю. Если SSB(t+1)

Критериями качества ПР являются: средний вылов за период, вероятность риска нарушения БО по биомассе и относительное изменение вылова от года к году. В модели случайный шум добавляется к численности запаса трески в стартовый год и к численности пополнения для каждого года (остатки после описания пополнения моделью описываются логнормальным распределением).

Поскольку фактические уловы не равны ОДУ, после расчета ОДУ по трехлетнему правилу рассчитывается величина улова - к оценкам ОДУ добавляется логнормальный шум (Булгакова, 2008, 2009).

Показано, что при использовании ПР (13) можно на 30 000 т в год увеличить среднемноголетний вылов, если увеличить Fpa до 0,5 (риск при этом равен нулю). При использовании трехлетнего ПР введение случайного шума в величину вылова резко увеличивает доверительные интервалы SSB, и вероятность риска пересечь граничный ориентир Blim достигает 14%.

Увеличение допустимого лимита изменения ОДУ трески до 15 % приведет к уменьшению риска до 8%.

Для сельди протестированы ПР в виде простой и двойной ступеньки (Bulgakova, 2008). Использование ступенчатой функции ПР вместо постоянного коэффициента промысловой смертности уменьшает величину риска падения запаса ниже уровня Blim и позволяет в прогнозные годы поддерживать вылов на высоком уровне.

Далее рассмотрены общие проблемы регулирования многовидового промысла. В первую очередь экосистемное регулирование должно обеспечивать нормальное существование ЭС - сохранение ее структуры, и по мере возможности - увеличение ее продуктивности. При разработке ПР многовидового промысла нужна четкая формулировка цели и определение БО управления аналогично тому, как определялось ПР предосторожного подхода в одновидовом случае. Если при выборе цели управления сообществом пренебрегать межвидовыми отношениями, легко разрушить структуру ЭС. Для эксплуатируемого сообщества одновидовой критерий MSY вообще теряет смысл, поскольку максимальные уловы разных видов достигаются в различных точках области допустимого управления на плоскости {F1,F2}(см. гл.2).

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

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

С помощью простых моделей показано (ICES, 1997b; Rice, 1997), что хищничество и каннибализм влияют на оценки БО, принятые в одновидовом анализе: пренебрежение многовидовыми отношениями приводит к занижению зависимости равновесного улова на пополнение Y/R(F) и завышению зависимости нерестовой биомассы на пополнение SSB/R(F) от коэффициентов промысловой смертности, и эти искажения зависят от биомассы хищников. При нестабильном промысле одновидовой подход может привести к переоценке величины ОДУ.

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

При оценивании БО для вида-жертвы следует учитывать значительное влияние на этот вид хищников высших трофических уровней (Furness, 2003).

Чтобы добытая пища скомпенсировала затраты на поиск и поимку жертв и на воспроизводство, как показали исследования по морским птицам, запас корма должен на 2-3 порядка превышать потребляемое количество пищи.

При выработке ПР промыслового сообщества необходимо учитывать как многовидовые биологические отношения, так и смешанный характер промысла, когда в одном районе одновременно работают разные промысловые комплексы (ПК) и каждый ПК вылавливает несколько видов. Изымаемые промыслом виды могут взаимодействовать биологически и/или через промысел. Приведен краткий анализ существующих моделей смешанного промысла. В главах 3 и приводится пример смешанного промысла двух видов рыб ЮВА, где промысел велся многими ПК, а виды взаимодействовали также как хищник и жертва.

Задача решалась путем выбора для каждого вида стандартного ПК и проведения стандартизации промыслового усилия.

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

ЗАКЛЮЧЕНИЕ В диссертационной работе описаны различные подходы к анализу многовидовых промысловых сообществ, разработаны и исследованы математические многовидовые модели разной сложности, предложены новые методы подготовки входной информации для многовидового моделирования, рассмотрены многовидовые аспекты регулирования промысла. Основные результаты, полученные диссертантом, заключаются в следующем:

1. Исследование на устойчивость к внешним возмущениям простых моделей сложных ЭС с разными трофическими отношениями между видами позволило получить ряд теоретических выводов о влиянии разных биологических процессов внутри популяции и взаимодействий между популяциями на стабильность сообщества и динамику его компонент:

Х модель биогеоценоза является более устойчивой к внешним возмущениям, чем модель биоценоза, а именно - обратная связь в системе в виде потоков вещества на косную компоненту (биогены в морской ЭС или количество гумуса в почве для наземного сообщества) стабилизирует динамику сообщества, и его устойчивое равновесие может существовать при любом количестве видов; стабилизирует динамику сообщества и ограниченный приток вещества (или солнечной энергии), необходимого для продуцентов;

Х временные задержки при круговороте вещества (время разложения вещества, поступающего со всех трофических уровней на нижний уровень, до биогенов, потребляемых первым трофическим уровнем) приводят к сдвигу по фазе колебаний численности разных трофических уровней.

2. С помощью простых моделей промыслового сообщества двух видов получены следующие теоретические выводы и результаты:

Х трофические отношения между видами резко сужают область допустимого управления в координатах коэффициентов промысловой смертности видов;

Х в модели хищник-жертва максимумы вылова популяций жертв и хищников достигаются на разных границах области допустимого управления сообществом;

Х пренебрежение конкурентами данного вида приводит к завышению рекомендаций по величине его ОДУ; для увеличения возможного вылова одной из популяций следует вылавливать и ее конкурента;

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

3. Разработанная автором динамическая продукционная модель двух промысловых видов, исходными данными для которой служит только промысловая информация, способна описывать промысловые сообщества двух и более видов с различными типами внутри- и межвидовых отношений и при смешанном характере промысла. На примере реального промысла ставриды и хека ЮВА в I974-I983 гг. проиллюстрирована специально разработанная процедура настройки модели, проведены расчеты по фактическим данным, получены рекомендации по регулированию промысла. Промысел, при котором годовой улов каждого вида сообщества постоянен и соответствует выбранной оптимальной точке, приводит к неустойчивости системы. Если же для перехода системы в ее оптимальное состояние изымать постоянную долю биомассы каждого вида, переходный период получается длинный, но система устойчива к такому воздействию. Приведена и более эффективная схема перехода к оптимальному состоянию сообщества.

4. Многовидовая когортная модель MSVPA, позволяющая исследовать довольно сложную структуру сообщества взаимодействующих видов, адаптирована к сообществу 4-х видов Баренцева моря. Результаты работы позволяют сделать следующие выводы:

Х Разработанные алгоритмы компьютерной обработки данных позволяют определить возраст хищников и возраст жертв в каждом желудке БД, что удалось сделать и в ряде случаев, когда жертвы были частично переварены.

Х Методология кригинга, апробированная нами при работе с БД по питанию, оказалась эффективной для учета неоднородного распределения хищников по ареалу. Эта методология может быть применена и к двумерным массивам с координатами времени{возраст хищника, год} или {возраст жертвы, возраст хищника} для экстраполяции и заполнения пробелов в данных.

Х Новый метод настройки модели MSVPA - последовательное применение одновидовых моделей класса ISVPA и MSVPA - является эффективным для видов, у которых только особи младших возрастов подвержены хищничеству.

Х Исследование динамики промыслового сообщества Баренцева моря по MSVPA показало, как велико влияние многовидовых отношений на динамику популяций. Например, в среднем популяция трески съедает в год около 1 млн т мойвы, в некоторые годы - до 3 млн т., от 200 до 400 тыс.т креветки, при этом смертность мойвы и креветки от хищничества часто намного превышает ее промысловую смертность. Каннибализм является мощным фактором, регулирующим численность поколений трески (только в 1996г. съедено 2тыс.т. ее молоди).

Х В Баренцевом море (в отличие от Северного и Балтийского морей, для которых тоже применялась модель MSVPA) велика изменчивость состава пищи трески, потому для получения реальной картины динамики популяций в ЭС следует собирать данные по питанию в течение длинного ряда лет, включающего годы подъема и падения численности как хищников, так и основных объектов их питания.

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

5. Для выбора ПР многовидового сообщества и его эффективного использования необходимы следующие этапы работы:

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

Х Выбрать критерий управления. Критерий MSY в многовидовом случае теряет смысл, поскольку максимальный вылов разных компонент сообщества достигается в разных точках области допустимого управления.

Х Внутри допустимой области найти точку, в которой критериальная функция принимает оптимальное значение.

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

Х Для каждого вида сообщества рассчитать величину ОДУ в оптимальной точке области допустимого управления.

Х При смешанном характере промысла найти распределение промыслового усилия между промысловыми комплексами, позволяющее как можно ближе подойти к оптимальной точке управления.

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

7. Предложенные автором стохастические сценарные модели, основанные на ретроспективных оценках запаса и на модели пополнения, учитывающей экосистемные факторы, позволяют оценивать эффективность разных ПР.

Основные публикации по теме диссертации Статьи в журналах, рекомендованных ВАК 1. Владимиров В.А., Булгакова Т.И. О смертности детенышей северных морских котиков// Рыбное хозяйство.Ц1973. №6. - С.13-15.

2. Андреев В.А., Булгакова Т.И. Челноков Ф.Г. Оценка вероятности потери меток и влияние мечения на выживание морских котиков// Рыбное хозяйство. - 1976. №8. ЦС.31-33.

3. Булгакова Т.И., Ефимов Ю.Н. Метод расчета величины возможного улова с учетом зависимости естественной смертности от возраста// Вопросы ихтиологии. Ц1982. Том 22 №2. - С.200-206.

4. Булгакова Т.И., Васильев Д.А., Долгов А.В., Фролова З.Н.. Сезонная и межгодовая динамика рационов трески Gadus Morhua Morhua Баренцева моря// Вопросы ихтиологии, - 1995. Tом 35 №2. - С.206-218.

5. Булгакова Т.И. Использование многовидовой модели для оценки естественной смертности видов-жертв// Вопросы рыболовства. - 2000. Т.1, №2-3. - С.71-73.

6. Васильев Д.А., Булгакова Т.И. Альтернативная оценка запаса баренцевоморской трески с использованием модели TISVPA// Рыбное хозяйство. - 2007. № 5. - С.54-60.

7. Булгакова Т.И. Сценарное моделирование, направленное на тестирование правила регулирования промысла северо-восточной арктической трески// Рыбное хозяйство. - 2009. N 4. - С.77-80.

8. Bulgakova T., Vasilyev D., Daan N. Weighting and smoothing of stomach content data as input for MSVPA with particular reference to the Barents Sea// ICES Journal of Marine Science. - 2001. Vol. 58, N 6 - p.1208-1218.

Статьи, опубликованные в других научных изданиях 9. Булгакова Т.И. (Эман) О некоторых математических моделях биогеоценозов// Проблемы кибернетики. Ц1966. Вып.16, Наука, М.: - С.191-202.

10. Булгакова Т.И. О моделях конкуренции видов// Проблемы кибернетики. - 1968. Вып.20, Наука, М.: - С. 263-270.

11. Булгакова Т.И. Об устойчивости простейшей модели биогеоценоза// Проблемы кибернетики. Ц1968. Вып. 20, Наука, М.: - С. 271-276.

12. Булгакова Т.И. К проблеме оптимального вылова для случая двух конкурирующих популяций (математическая модель)// Труды молодых ученых ВНИРО, - 1970. Вып.4. ОНТИ ВНИРО, М.: - С.143-156.

13. Булгакова Т.И. Применение математических методов для расчета оптимального выбоя морских котиков о.Тюленьего// Сб. трудов по промышл.

рыболовству, Ц1973. Т.2. ЦНИИТЭИРХ, М.: - С. 30-52.

14. Булгакова Т.И., Засосов А.В., Шатуновский М.И. О моделировании некоторых хозяйственных систем рыболовства и рыбоводства с учетом эколого-физиологических параметров// Труды ВНИРО - 1973. Т.94. - С.9-23.

15. Булгакова Т.И. Модель популяции типа "запас-пополнение", учитывающая влияние кормовой базы производителей на выживание молоди// Труды ВНИРО Ц1978. Т.128. М.: ВНИРО. - С.16-22.

16. Булгакова Т.И., Андреев В.А., Челноков Ф.Г. Метод оценки некоторых параметров популяции морских котиков по материалам мечения// Труды ВНИРО Ц1978. Т.128, М.: ВНИРО. - С.23-32.

17. Булгакова Т.И., Кизнер З.И. Математическое моделирование двухвидового промысла// Информац. и математ. обеспечение исследования сырьевой базы.

Труды ВНИРО. - 1985. М.: - С.97-104.

18. Булгакова Т.И. Программа "Расчет возможного улова с учетом зависимости естественной смертности от возраста"// ВНТИЦЕНТР гос. регистр.

50850000888 Ц1985.

19. Булгакова Т.И., Бабаян В.К. Программа реализации анализа размерных когорт// ВНТИЦЕНТР гос. регистр. 50850000889. - 1985.

20. Булгакова Т.И., Грауман Г.Б. Модель эмбрионального периода развития трески Балтийского моря// Fischerei-Fortschung,Ц1991. B.29(2). Rostock - s.61-67.

21. Булгакова Т.И., Васильев Д.А., Долгов А.В., Фролова З.Н. Подготовка входной информации для многовидовой модели промыслового сообщества рыб Баренцева моря на основе базы данных по питанию хищника// Вопросы теории и практики промышленного рыболовства. Поведение гидробионтов в зоне действия орудий лова. Ц1998. М.: Изд-во ВНИРО. - С. 88-96.

22. Булгакова Т.И., Васильев Д.А. Разработка методологии оценки переменных параметров многовидовой модели с учетом их пространственной и временной изменчивости// Деп. ВНИЭРХ 1325 98; ВИНИТИ №6, Ц1998. Ц52 с.

23. Бабаян В.К., Булгакова Т.И., Васильев Д.А., Фролова З.Н.

Информационное обеспечение рационального использования запасов осетровых //Рыбоводство и рыболовство. - 2001. № 1. - С.42-44.

24. Булгакова Т.И. Экосистемный подход к оценке запасов и регулированию промысла // Вопросы промысловой океанологии. - 2004. В.1. М.: ВНИРО. - С.275-295.

25. Булгакова Т.И., Экосистемные аспекты рационального рыболовства// Морехозяйственный комплекс России: эколого-географические проблемы, Сб. научных трудов, - 2005. СПб. Российское геогр. общество. - С. 19-24.

26. Булгакова Т.И., Дьяков Ю.П. Проблемы оценки запаса желтоперой камбалы у западного побережья Камчатки в условиях недостатка информационного обеспечения // Методические аспекты исследования рыб Дальнего Востока:

Труды ВНИРО. Ц2006. Т. 146, М.: ВНИРО.ЦС.191-204.

27. Babayan V., Bulgakova T. Assessment of natural mortality and growth parameter of Cape horse mackerel (Tr. trachurus capensis) in Div.1.3+1.4// ICSEAF Colln.

scient. papers.Ц1983. V. 10. Madrid: - pp.49-53.

28. Babayan V., Bulgakova T. Estimation de lТetat du stock at de la capture admissible de chinchard du Cap (Tr. trachurus capensis) dans les Div.1.3 et 1.4 a partir de lТanalise de la composition par longuer des captures// ICSEAF Colln. scient.

papers. Ц1983. V.10. Madrid: - pp.37-47.

29. Babayan V., Bulgakova T. Comparison of some biostatistical methods applied for stock assessment and catch forecasting of Cape horse mackerel in Div.1.3+1.4// ICSEAF Colln. scient. papers.Ц1984. V.11(1). Madrid. - pp.5-13.

30. Bulgakova T., Kolarov P. et al. Stock assessment and catch projection for Cape horse mackerel in ICSEAF// ICSEAF Colln. scient. pap.Ц1984. V.11(1). Madrid. - p.15-27.

31. Bulgakova T., Daan N., Gislason H., Vasilyev D. Final Report INTAS-RFBR project 95-1204 УMultispecies Analysis for the cod-related community of the Barents SeaФ 1996-1998. RIVO-DLO report N C009/99, Ц1999. Ijmuiden: Ц24 pp.

32. Babayan V.K., Bulgakova T. et al. Caspian sturgeon TAC foundation in modern conditions// VNIRO Publishing. - 2006. - M.: Ц27 p.

33. Vasilyev D., Bulgakova T. North-East Arctic cod stock assessment by means of the TISVPA model// Towards development of optimal strategy for sustainable fishery in the Barents and Norwegian seas (for the 36-th session of the joint RussianNorwegian Fisheries committee), - 2007. VNIRO. Цp.7-22.

Работы, опубликованные в материалах конференций и симпозиумов 34. Булгакова Т.И., Модели многовидового промысла в рыбном хозяйстве// Имитационное моделирование в решение задач прогнозирования состояния запасов и оптимизации промысла рыб// Препринт докладов 4-8 мая 1987г. - 1987. Петрозаводск: Цс.30-38.

35. Булгакова Т.И., Бородин Р.Г., Васильев В.А., Коржев В.А. Многовидовой анализ промыслового сообщества// Труды междунар. форума по проблемам науки, техники и образования. 2- 6 декабря 2002. М.: Том 3. - С. 21-22.

36. Бабаян В.К., Булгакова Т.И., Васильев Д.А. и др. Методический подход к оценке запасов и ОДУ каспийских осетровых// Материалы междунар. конф.

Современное состояние и пути совершенствования научных исследований в Каспийском бассейне. Ц2006. Изд-во КаспНИРХ, Астрахань: - С. 89-96.

37. Бабаян В.К., Булгакова Т.И., Васильев Д.А. и др. Реализация модельного подхода к обоснованию допустимого промыслового изъятия русского осетра Acipenser gueldenstaedtii в росс. зоне Каспийского моря. Там же. - С.97-103.

38. Бабаян В.К., Булгакова Т.И., Васильев Д.А. и др. Альтернативный подход к обоснованию ОДУ каспийских осетровых в современных условиях// Труды междунар. форума по проблемам науки, техники и образования под ред. В.П.

Савиных и В.В. Вишневского. Ц2006. Том 2. М.: - с. 138-140.

39. Булгакова Т.И. Имитационное моделирование для тестирования разных вариантов правила регулирования промысла северо-восточной арктической трески // Материалы Междунар. научно-практич. конф. Комплексный подход к проблеме сохранения и восстановления биоресурсов Каспийского бассейна (1316 окт. 2008 г., Астрахань) - 2008. Изд-во КаспНИРХ: ЦС.61-40. Бабаян В.К., Булгакова Т.И., Васильев Д.А. Модельный подход к оценке неучтенного вылова каспийских осетровых// Там же: 2008. - С.36-40.

41. Bulgakova T. Grauman G.B. A model of embryonic period of development in Baltic Sea cod.// ICES. C.M. 1990/J:10. Ц21 p.

42. Bulgakova T. Vasilyev D. The relationship between the initial meal weight and weight in stomach for various regimes of nutrition// ICES CM 1995/D:15 Ref.

G,H,J,M - 8p.

43. Bulgakova T., Vasilyev D. et al. The results of multispecies analysis for the Barents Sea fishery community (cod, capelin, shrimp and herring)// ICES CM 1995/D:14 Ref. G,H,K. - 24p.

44. Bulgakova T., Vasilyev D. et al. Some special algorithms for the Barents Sea fish stomach content data base processing// ICES CM 1995/D:13 Ref. G,H,K. Ц12 p.

45. Bulgakova T., Vasilyev D. Seasonal and yearly dynamics of the Barents Sea cod rations// ICES CM 1995/G:32 Ref. D,H,K. Ц28 p.

46. Bulgakova T., Vasilyev D. Stomach data base analysis with special attention to spatial aspects of feeding parameter estimation and to the recovery of gaps in the data// ICES. Multispecies Assessment WG 1997. WP 2. Ц48p.

47. Bulgakova T., Vasilyev D. Estimation of the MSVPA input parameters characterizing by temporal and spatial variability// Joint Russian-Norw. Meeting on Multisp. Models of Relevance to Management of the Barents Sea, Norway, 9-June 1998. Ц16 p.

48. Bulgakova T., Vasilyev D. The results of application of 4-species MSVPA to the Barents Sea fish community// Joint Russ.-Norw. Meeting on Multispecies Models of Relevance to Management of the Barents Sea, Norway, 9-12 June 1998. Ц21p.

49. Vasilyev D., Bulgakova T. Tuning of MSVPA for cod-oriented fish community of the Barents Sea by means of Quarter-based Instantaneous Separable VPA (QISVPA)// Joint Russian-Norwegian Meeting on Multispecies Models of Relevance to Management of the Barents Sea, Norway, 9-12 June 1998. Ц12p.

50. Bulgakova T. Optimum control in a predator-prey model based on a two-species exploited ecosystem// Ecosystem Approaches for Fisheries Management. Alaska Sea Grant College Program. Ц1999. AK-SG-99-01, - p. 149-162.

51. Bulgakova T., Borodin R., Kotenev B. Ecological role of marine mammals in the north-eastern Atlantic//The 5th Atlantic Fisheries Ministers Confer.(NAFMC). - 2000.

52. Bulgakova T. NEA cod dynamics simulation to test various management scenarios// ICES ASC 24-27 Sept. Tallinn, Estonia. CM 2003/X:08. Ц2003.Ц17p.

53. Bulgakova T. A simulation testing of various management regimes for the NEA cod stock// ICES ASC 20-25 Sept. Vigo, Spain. CM 2004/FF:08. Ц2004. Ц14 p.

54. Bulgakova T. A model of Norwegian spring spawning herring dynamics as tool for testing for harvest control rule variants// Proceed. of 12th Norw.-Russian Symp.

Troms, 21-22 Aug. 2007 // Edited by T. Hauge et al.Ц2008. Bergen: - P.103-113.

Методические рекомендации и пособия 55. Методические рекомендации. Применение математических методов и моделей для оценки запасов рыб// Авторы составители: Бабаян В.К., Бородин Р.Г., Булгакова Т.И., Ефимов Ю.Н. Ц1984. ОНТИ ВНИРО М.: - с.155.

56. Методические рекомендации по математическому моделированию двухвидового промысла// Авторы составители: Булгакова Т.И., Кизнер З.И. - 1987. ВНИРО. М.: Ц39 с.

57. Многовидовой анализ промыслового сообщества: Методическое пособие// Под редакцией Т.И. Булгаковой. Авторы-составители: Булгакова Т.И., Бородин Р.Г., Васильев Д.А., Коржев В.А. 2001. ВНИРО, М.: - 113с.

58. Методические рекомендации по обоснованию общих допустимых уловов (ОДУ) каспийских осетровых// Изучение экосистем рыбохоз. водоемов, обработка данных о водных биол. ресурсах, техника и технология их добычи и переработки. Вып. 4. Авторы составители: Бабаян В.К., Булгакова Т.И., и др. - 2006. ВНИРО, М.: - 59 с.

   Авторефераты по всем темам  >>  Авторефераты по техническим специальностям