Бесплатный автореферат и диссертация по географии на тему

Моделирование собственных колебаний и ветровой циркуляции западной части Черного моря

ВАК РФ 11.00.08, Океанология

Автореферат диссертации по теме "Моделирование собственных колебаний и ветровой циркуляции западной части Черного моря"

ГОСУДАРСТВЕННЫЙ КОМИТЕТ РСФСР ПО ДЕЛАМ НАУКИ И ВЫСШЕЙ ШКОЛЫ ЛЕНИНГРАДСКИЙ ГИДРОМЕТЕОРОЛОГИЧЕСКИЙ ИНСТИТУТ

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

УДК 551.465 : 551.468

ДЕМИРОВ Енчо Койчев

МОДЕЛИРОВАНИЕ СОБСТВЕННЫХ КОЛЕБАНИЙ И ВЕТРОВОЙ ЦИРКУЛЯЦИИ ЗАПАДНОЙ ЧАСТИ ЧЕРНОГО МОРЯ

Специальность 11.00.08 — ОКЕАНОЛОГИЯ

АВТОРЕФЕРАТ

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

ЛЕНИНГРАД — 1991

/ У

/

(

Г

/

Работа выполнена в Ленинградском гидрометеорологическом институте.

Научный руководитель — доктор физико-математических наук, профессор Б. А. КАГАН.

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

доктор физико-математических наук Л. А. ТИМОХОВ, кандидат физико-математических наук И. А. НЕЕЛОВ.

Ведущая организация: Ленинградское отделение Государственного океанографического института.

Защита состоится « <¿¿¿1/»¿¿¿¿^ 1991 г.

в /#* час. мин. на заседании Специализированного

совета Д. 063.19.01 в Ленинградском гидрометеорологическом институте по адресу: Ленинград, Малоохтинский пр., 98.

С диссертацией можно ознакомиться в библиотеке Ленинградского гидрометеорологического института.

Отзывы; & двух экземплярах, заверенных печатью, просьба направлять по адресу: 198196 Ленинград, Малоохтинский пр., 98, Ленинградский гидрометеорологический институт.

Автореферат разослан

2*? ¿рь^^ъгл^ 1991

Ученый секретарь Специализированного совета, доктор географических наук,

профессор 10. И. ЛЯХИН

г

г; " " 3 "

4 ОИЦАЯ ХАРШЕРИСИШ РАБОТУ

1. Л..: Отдел сврт;.".

Актуальность проблем».Р последнее время в связи со все бо-

'лёе'возраствюпцш загрязнением Черного моря п нем сложилась сложная экологическая обстановка. Особенно тревожна ситуация в при-бреглой лоно,где из-за относительно слабой циркуляции и пороче-щи?аняя,с одной стороны,л непрекращчгшэгося сброса нрсмьглленных отходов,с другой .кногие мелководные районы окапались на грани экологической катастроф. ?д'соъ все чаща наблюдается так называемый "красжЛ прилив",уменьшается видовое разнообразие флоры и фаун« и резко ухудшается качество поди. Дальнейшее ухудпение качестве. годы мотне предотвратить лишь в том случае,если при рт-гаенин хозяйственных вопросов будут приписаться во внимание особенное::' происходя®»:', в прибрежной части моря процессов.

ОбЪ'Л1 сучестпу^тих сейчас экспериментальных дянньге в гель-форсЧ зоне западно!! част Черного мор» не достаточен для реле-пня практические проблем. Поэтому для получения количественных сценой остается одмчотпенная возможность - привлечь результаты численного моделирования.

Цели и задали работы. Цель работы состояла в воспроисчл«-нпи спектра собственник кт'сбаний и ветровой циркуляции в зг'/а-дчой части Черного моря и шяснеиии механизмов их формирования.

Достижение этой цели требует решения спектральной задачи и задачи о ветровой циркулями в ограниченной облаеп*, гклпчающр?' мелководный »гли*!,резки!', оп->д глубин н окрестности внес«:"!! кринки иелгфл и глубок о во дну р о^ятть моря. Использование в такой области тр'а^ип'снных численных методов интегрирования уравнений гид| елинсн'и!'.м гг'.5г"Н!ч';< ьп »й1?1'ии рт-лярньге сеток,еопря

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

Научная новизна. На основе решения спектральной задачи впервые определены периоды и пространственная структура собственных колебаний Черного моря в интервале от 2 ч. до 20 сут.; истолкована природа этих колебаний; показано,что максимальный за последние 100 лет штормовой нагон в западной части Черного моря,имевший место с 17 по 21 февраля 1979 г..является следствием резонансного усиления захваченных береговых волн. Воспроизведены особенности ветровой циркуляции в районе банки Кокет-ряйс,в районе прибрежного апвеллинга и в окрестности мыса Кали-акра.где происходит формирование фронтальных зон.

Практическая ценность работы. Разработанная численная модель ветровой циркуляции реализована в одном из наиболее вяжнык в прикладном отношении районов - районе банки Кокетрайс. Доказана возможность её использования для оценки потенциальных изменений циркуляции, выз ванных планируемой ¡эксплуатацией этой области .

Установленный в результате решения спектральной задачи вы-

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

Аппробация работы. Основные результаты докладывались на семинаре секции физики моря Института океанологии ЕАН (октябрь 1990 г.). Полностью диссертационная работа обсуждалась на семинаре кафедры динамики океана ЛГШ (ноябрь 1990 г.) и совместном семинаре Отдела взаимодействия океана и атмосферы и Отдела океанологи.: МШШ (декабрь 1990 г.).

Публикации. По теме диссертации опубликованы 3 статьи.

Автор выносит на защиту:

1. численную модель расчета собственных колебаний и результат её реализации применительно к условиям Черного моря;

2. установленный в реэутьтате решения спектральной задачи и аналг. зэ спектров скорости летра и колебаний уровня моря вывод о резонансном усилении захваченных береговых волн и преобладающем влиянии резонансных эффектов в формировании экстремальных штормовых нагонов на северо-западном мелководном шельфе Черн-'^о моря;

3. модель ветровой циркуляции в западной части Черного моря в Ларотропном и двухслойном приближениях,использующую криволинейную систему координат;

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

Структура и oSi.€M_Fç,tfpru. Диссертация состоит из введения, четырех глав и заключения. ОбишЙ об'ьеь! работа состаапяет 152 стр. количество иллюстраций - 60,количество таблиц - 7, Список литературы содержит 65 наименования.

СОДЕШАШЕ РЛБОШ

Do jbe^fibmi обосновывается актуальность темы. Здесь форму лирук-тся также :.члачи работа .кратко излагаются её структура и содержание, укал нва-тся степень научной новизны полученных результатов и перечисляется осио raine положения,выносимые на рацкту.

Первая глава. В раздело I.I рассматриваются существуйте методы расчета собственных колебаний морских бассейнов. Из честно,что в простейием случае неиршчамцегося бассейна постоянной глубины в спектре собственных колебаний присутствуют только гравитационные моды,имеющие форму чисто стоя'мх вилн. Учет вращения I' топографии дна приводит к модификации гравитационных мод и к возникновению нового класса волн - волн Россби. Наклон дна в пределах шельфа и континентального склона обусловливает существование береговых захваченных волн,распространяющихся параллельно берегу, Физический механипы воглшкнопении етих волн,как и топографических волн Россби, - растяжение вихревых нитей в поле планетарного вихря. Экспериментальные исследования в разных районах Мирового океана свидетельствует о том,что энергетика меао-масгатабннх движений не шельфе э гиачитсяишЯ норе определяется захваченными волнами. Их песпроиоведение в рамках численных мо-дспей возможно лишь при задании достаточно тонкого пространственного рр зргшенпя.

Ииве.с.тнц дче группы численных методов решения споктраwott

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

В разделе 1.2 содержится обзор суцествртцих методов расчета вынужденны.!; движений п ограниченной области. Эти методы делятся на три группы: к первой принадлежат методы,использующие те-лескопиэации модельний области,ко второй - методы,использующие на открытом контуре исследуемой области условие кмпедпнснсго типа,н к третьей - методы,используквде нерегулярные сетки. Недостатком методов эторсй группы является отказ от учета обратной связи между решениями в ограниченной области (скажем,в прибрежной зоне) и во всем бассейне,тогда как недостатками методов третьей группы являются искажение профиля волны,а также фиктивное отражение и преломление волн на границе областей с различным пространственным раэреиением. Указанные негативные последствия можно ослабить введением криволинейной системы координат. В отом случае решение исходной задачи сводится к реализации двух процедур - построению криволинейной сетки и численному интегрированию уравнений модели в новых координатах.

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

дани,не превосходящими 24 ч. Между тем экспериментальные данные свидетельствуют о существовании в Черном море россбиевских мод с более продолжительными периодами.

Расчеты штормовых нагонов в западной части Черного моря производились на неравномерных сетках с резким изменением пространственного шага. Лри этом рассматривался только верхний 200-метровый слой,а нижний слой считался неподвижным. Ясно.что такое приближение исключает движения,обусл^зчнваемые нерегулярностями топографии дна,и фильтрует или искажает захваченные береговые волны.

В разделе 1.4 рассматриваются основные особенности синоптических процессов и астрового режима Черного моря. Отмечается, что сильные штормовые ветры над западной частью Черного моря имеют преимущественно северо-восточное и восточное направления. Их скорость в отдельных случаях достигает 40 м/с. Чаще всего штормовые ветры над Черным морем связаны с процессами формирования и перемещения средиземноморских циклонов и глубоких антициклонов, сформировавшихся в районах Азорского и Сибирского максимумов.

Во второй главе обсоздаются постановка спектральной задачи, разностный метод её решения и результаты,полученные применительно к условиям Черного моря.

Если аппроксимировать плотностную стратификацию Черного моря двухслойным приближвнием,то уравнения гидродинамики в криво-

линейной системе координат вид

^ (-С > У ) принимают

(I)

ЭЖ . ШША ЭУсК п {3)

ЪЬ - и >

где LC¿ , - составляющие контравариантной скорости; ¿. = 1,2 - индексы,указывающие на принадлежность к нижнему и верхнему слоям соответственно; - плотность с-го слоя; Н- его то-

ящина; =./>г (р( "Рг)£ ^ " £г "ЯЦ^г ~ Прения;

^ ^ ^ ~ отклонения поверхности раздела и свободной поверхности моря от их не возмущенных; положений; = ^ ; - ^;

«-Зи/У 'РЗчР' - ком-

риантные компоненты метрического тензора; У = Р(*>У). - якобиан преобразования координет -^(£>1)

Криволинейные координат ^ = ^ (*,%)> ЧСХ>-*) находятся из решения краевой зада1®

+ 0.(1,1),

где X. , у - декартовые координаты, В и - управляющие функции; граничное условие (6) определяется соответствием значений Х/ Ч & Э& и £ ) 2, ^ считается также,что одна

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

Уравнение спектральной задачи следует из (П-(3) ггри ^ = - •

= О, = О, с. = 1,2. На границе область з^л-эстся ус-

ловие )..»протекания. Собственные колебания В'г.'буз-.'дп-'иси печаль-

ним возмущением уровня свободной поверхности. Временные производные в (I)—(3) аппроксимируются методом Кранка-Николсона,пространственные - центральными разностями. Полученная в результате решения последовательность значений функционала

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

Описанная вше процедура была реализована в трех случаях: I) для глубоководной части Черного моря в баротропном приближении, 2) для всего Черного моря (включая мелководный северо-западный шельф) в баротропном приближении и 3) для всего Черного ¡.¡оря в двухслойном приближении. Сравнение этих трех решений ггс зволяет оценить вклад северо-западного шельфа и плотнос.тной стратификации г формирование спектра собственных колебаний Черного моря.

Как оказалось,низшая мода собственных колебаний глубоково-

дной части Черного моря имеет период 4,62 ч. Такой же период низшей гравитационной моды дает формула И"риана для прямоугольного невращающегося бассейна. Совпадение периодов указывает на то,что вращение Земли не оказывает сколько-нибудь значительного влияния на формирование низшей гравитационной моды. Соизмеримое!., длины большой полуоси Черного моря (~1Ю0 км) и баротрогшого радиуса деформации Россби 1200 км) подтверждает этот вывод. Характерная особенность пространственной структуры низшей гравитационной моды - узловая зона,ориентированная в меридиональном направлении и проходящая через середину Черного моря.

Вторая мода глубоководной части Черного моря имеет период 2,3 час. Её пространственная структура определяется существованием двух положительных (циклонических) амфидромкй,одна из которых располагается в восточной,другая - в западной > э.стях моря.

Учет северо-западного мелководного шельфа приводит к обогащению спектра собственных колебаний Черного моря и появлению дополнительных мод с периодами 5.53; 6,66; 9,82 ч. Поскольку ба-ротропнпй родиус деформации Россби 250 км) на северо-западном мелководном шельфе заметно меньше его пространственной протяженности, вращение Земли оказывает существенное влияние на формирование моды с периоде..! 9,82 ч. Эта мода локализована преимущественно в шельфовой зоне и имеет пространственную с тру муру запоминающую волну Кельвина. При учете северо-эяшдного мелководного шельфа несколько изменяются также-и'собственные колебания глубоководной части Черного моря. 3.частности.период низшей гравитационной моды увеличивается с 4,62' до 4,74 ч.

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

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

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

Влияние плотностной стратификации способствует дальнейшему обогащению спектра собственных колебаний Черного моря и появлению в нем бароклинных мод с периодами 16,65 ч., 17,07 ч.,17,73 ч., 19,23 ч., 23,95 ч.,25,28 ч.,31,75 ч.,40,16 ч.,44,04 ч., 40,76 ч. Существование бароклинной модн с периодом 31,75 ч. подтверждается экспериментальными данными,в соответствии с которыми эта мода имеет вид внутренней волны Кельвина.

Вторая глава заканчивается обсуждением механизма формирования катастрофического штормового нагона,имевшего место в западной части Черного моря с 17 по 31 феврапя 1979 г. .когда подъем уровня в пункте Ахтопол достиг отметки 2,0 м.

Спектральный анализ рядов измерении скорости ветра в этот период обнаруживает существование '•ре.«. низших иод с периода!.!;! ~64; 36,6 и 27 ч. Эти периедч близки к периодам = 62,53, 77= 36,57 и ^ = 27,31 ч трех собственных мод Черного моря.Отсюда следует,что катастрофический штормовой нагон может быть Еызван резонансны.! усилением собственных колебаний Черного моря

Кроме мод с периодами,близкими к и (в силу ограниченности выборки,выявить гармонику с периодом не удьется) спектр колебаний уровня содержит также моды с периодами 10,3; 12,8; 0,5; 6,1 и 4,3 ч. Их появление связано с нелинейным взаимодействием собственных колебаний Черного моря. Подтверждением тому могут служить оценки периодов 19,00; 12,50; 8,58; 5,91 и 4,01 ч. комбинационных мод,обуславливаемых взаимодействием мод с периодами 62,53 ; 36,57 и 27,31 ч.

3 третьей главе рассматривается задача о ветровой циркуляции в западной частч Чгрного моря. Задача решается с использованием криволинейной системы координат и с учетом эффектов адвекции импульса.горизонтальной турбулентной вязкости и придонного трения.

Уравнения модели в двухслойном приближении получаются подстановкой в (1)-(3) выражений

...........(8)

о); =

(9) (10)

г ; jrн¿

^ Т-ГГ * ^ >

(Я*.- Г!|, Ц К*(Щ >%)**?, ^

** > (13)

(&гг\ АТ?) - - т-/} гг* <15)

С, (16)

(10)

где , - составляющие скорости с. -го слоя в декартовых координатах; , , - коэЭДшдаенты трения на поверхности моря,на поверхности раздела двух слоев жидкости и на днелсоответственно; ^ , Ц/у - составляющие скорости ветра;- плотность воздуха; й - коэффициент горизонтальной турбулентной вязкости. На береговом контуре задается условие прилипания

' Щ = = О на Г-и } (19)

в начальный момент гремени - состояние покоя

£¿ = 0 ' при ¿=0. (20)

Временные производные в (1)-(3) ,(8)-(20) аппроксимируются по схеме Кранка-Николсона,пространственные производные - центральными разностями. Как и в задаче о собственных колебания*,полученная система разностных уравнений решается итерационным методом. Решение получается е баротропном и двухслойном приближениях

В баротропной задаче скорость ветра считается неизменной во времени и в пространстве и равной 20 м/с. Такое же предположение делается и в отношении направления ветра: оно задается либо северны.!, либо южным.

Показано,что в установившемся состоянии колебания уровня в прибрежной зоне не превышают 15-20 см. При этом изолинии уровня практически совпадают с изобатами. Такое распределение нарушает-

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

Используемая в расчете нерегулярная сетка,хотя и имеет относительно хорошее пространственное разрешение в мелководной зоне зпппг.ней части Черного моря,тем не ненее не позволяет произ-втги детальный учет мелкомасштабных особенностей рельефа дна и очертаний береговой линии (например,белки Кокетрайс в окрестности болгарского побережья). Поэтому задача о гэтровой циркуляция в районе банки Кокетрайс решалась па более меляей нерегулярно;! сетке с граничшии условиями на открытом контуре »заимствуем юм из решения соответствующей задачи для всего Черного меря.

Полученное решение воспроизвело основные особенности циркуляции, известные по данным измерений. Именно,в восточной части рассматриваемой области скорости течения принимают наибольшие значения (до 33-4Э с«/е), воспроизводится также тенденция появления антициклонической завихренности в восточной и циклонической в западней частях области. 1С северу от Несербского полуострова скорости течения составл,гат 20-25 си/с.

Четвертая глава посвящена воспроизведении особенностей циркуляции в окрестности мнеа Калиакра,вызванных существованием здесь апвеллинга. С этой цельп била использована двухслойная модель (1)-(3),(8)-(20). Скорость ветра принималась ряпчсЛ 10 м/с вбливи бг-пгчрсяого побережья и убирощей по экспонате к

центру моря.направление ветра - южным,юго-западным (такое направление соответствует сгонным ветра),1). Интегрирование уравнений модели начиналось с состояния покоя и производилось на срок,равный 3 сут.

- Полученное нестационарное решение содержит четко выраженные колебания с периодом 31-32 ч.,близкие к периоду внутренней волны Кельвина. Решение воспроизводит также две области с максимальным подъемом поверхности раздела. Одна из них располагается в окрестности пролива Босфор,другая - у мыса Калиакра. Существование области подъема поверхности раздела у мыса Калиакра подтверждается данными измерений температуры и солености.

В заключении подводятся основные итоги работы. Они сводят-"п к следующему:

1. Определены периоды и пространственная структура собственных колебаний Черного моря. Периоды собственных колебаний сосредоточены в двух интервалах спектра - ниже 10 ч. и выше 27 ч. В интервале от 10 ч. до 27 ч. баротропных собственных колебаний Черного моря нет. Этим объясняется отсутствие значимее приливов в Черном море. Гравитационная мода с периодом 9,82 ч. имеет пространственную структуру,напоминающую волну Кельвина. Моды с периодами 1С,7 сут и 10 сут. - это низшие захваченные волны,а мода с периодом 4,3 сут. - низшая топографическая велна Россби. Найдены 'периоды внутренних гравитационных волн. Они равны 16,65, 17,07, 17,73, 19,23, 23,95, 25,28, 31,75, 40,16, 44,04 и 48,76 ч. Согласно экспериментальным данным ь:ода с периодом 31,75 ч. является внутренней волной Кельвина.

2. Показано,что экстремальный штормовой нагон,имевший место на болгарском побережье Че рного моря в фспрплс 1979 г.^ яшш-

ется следствием резонансного усиления захваченных береговых мод с периодши 62,33, 36,57 и 27,31 ч. Анализ спектра колебаний уровня,измеренных в пункте Ахтопол в указанный период,свидетельствует о том,что обнаруживаемые в спектре колебаний уровня коря максимумы на периодах 10,3, 12,6, 8,5, 6,1 и 4,3 ч. отвечает комбинационным модам,сбуслоплинаемым взаимодействием мод с периодами 62,53, 36,57 и 27,31 ч.

3. 3 рамках баротропноП модели ветровой циркуляции,записанной в криволинейной системе координат,воспроизведено поле течений в западной части Черного моря. Показано,что резкое по;'чх;ение уровня в районе мысов является следствием конвергенции течений

и искривлений линий тока.обусловленных резким изменением очертаний береговой линии.

4. Баротропная модель ветровой циркуляции с использованием криволинейной системы координат и телескошзацпи области реализована применительно к условия).! банки Кокетрайс. Воспроизведены основны» особенности циркуляции в этой области .известные пЬ-йнс-периментальным данный. В частности,установлено,что в восточной части рассматриваемой области поло скорости обладай? антпцпкло-нической.а з западной - циклонической тглхрснностьв;' ¡МЧсямаль-ные скорости течения отмечаются в восточной части облас® и севернее Несербского полуострова.

5. 3 рамках двухслойной модели ветровой циркуляции воспроизведено поле течений в западной части Черного моря. Показам;, что при задании постоянной скорости ветра процесс установления циркуляции сопровождается появлением внутренней волны с периодом 31-32 ч..близким к периоду внутренней волны Кельвина. Максимальные отклонения поверхности раздела от её невозмущенного

- ю -■

положения обнаруживаются в окрестности мыса Калиакра и в районе пролива Босфор.

По теме диссертации опубликованы следукхцие статьи

1. Дамиров Е.К. Численное решение задачи о собственных колебаниях. Черного моря // Океанология. - 1907. - Т.27, вып.5. -С.721-726.

2. Dimitrov p.s., Stcyanov A.S., Eemirov E.K., SoJakov d.P. On the phenomenon of upu«lling in the western part of Bluet sea and its connection with the blocun cf plankton // Докл.БАН. -1937. - 40, Jf 10. - C. 127-130.

3. Димитров П.С., Демиров Е.К., Чернева К.И. Геоэкологические аспекты эксплуатации прибреашых аккумулятивных тел (на примере песчаной банки КокетраИз - Черное море) // Геологический мониторинг и проблемы геоэкологии Балтийского и Черного морей: Тез.докл. - Ленинград, 1990. - С.67-68.

'!одвиоз>!й к печати 13.02.у1г. иоьем 1 п. п. Тиран 100. Зак.7П. '«'.платно

Типография ВАС