Бесплатный автореферат и диссертация по наукам о земле на тему
Термодинамическое моделирование формирования морского ледяного покрова в Арктике
ВАК РФ 25.00.28, Океанология

Автореферат диссертации по теме "Термодинамическое моделирование формирования морского ледяного покрова в Арктике"

АРКТИЧЕСКИЙ И АНТАРКТИЧЕСКИЙ НАУЧНО-ИССЛЕДОВАТЕЛЬСКИЙ ИНСТИТУТ

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

Андреев Олег Михайлович

ТЕРМОДИНАМИЧЕСКОЕ МОДЕЛИРОВАНИЕ ФОРМИРОВАНИЯ МОРСКОГО ЛЕДЯНОГО ПОКРОВА В АРКТИКЕ

Специальность 25.00.28 - океанология

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Санкт-Петербург - 2005

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

Научный руководитель:

кандидат географических наук, с.н.с Иванов Борис Вячеславович

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

доктор физико-математических наук, с.н.с. Макштас Александр Петрович

доктор физико-математических наук, профессор Царев Валерий Анатольевич

Ведущая организация:

Санкт-Петербургский государственный университет, факультет географии и геоэкологии.

Защита состоится 20 октября 2005 г. в 15 час на заседании диссертационного совета Д 327.002.01 при Арктическом и антарктическом научно-исследовательском институте по адресу: 199397, Санкт-Петербург, ул. Беринга, 38.

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

Автореферат разослан «20 » 2005 г.

Ученый секретарь диссертационного совета

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

Актуальность проблемы

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

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

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

Первые попытки математического моделирования промерзания среды относятся к середине XIX века, и наиболее серьезным исследованием в этом направлении следует признать работу И.О.Стефана (1891), в которой были заложены основы современных методов термодинамического моделирования взаимодействующих сред с движущейся границей раздела фаз. Однако наибольшая исследовательская активность началась в 40-60 годах XX века, когда появились работы А.А.Шепелевского, Л.И.Рубинштейна и В.Г.Меламеда, продемонстрировавшие возможности термодинамического подхода применительно к моделированию морского льда. Отдельно следует выделить классическую работу Ю.П. Доронина (1963). Эта термодинамическая модель морского льда и в настоящее время является одной из наиболее часто используемых. Появившиеся несколько позже модели зарубежных авторов Мэйкута-Унтерштейнера (1971) и Семтнера (1976) стали основой для современных термодинамических моделей морского льда, используемых для различных теоретических исследований.

Однако многие аспекты в термодинамическом моделировании формирования и эволюции морского льда до сих пор являются не решенными. Эволюция вертикального распределения солености морского льда, осолонение подледного слоя

изменения его толщины, причем равновесие этой системы поддерживается и

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

Цель диссертационной работы:

Исследовать основные механизмы формирования морского льда и воспроизвести сезонный ход основных физических параметров морского льда в Арктике путем математического моделирования, в рамках термодинамического подхода.

(\.нпишля чяляии исследования

1. Разработать численную термодинамическую модель морского льда.

2. Определить вклад различных факторов в формирование морского льда, оценить скорость промерзания и таяния тороса.

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

4. Исследовать применимость модели к различным физико-географическим условиям Арктики.

Научная и практическая значимость

Диссертационная работа выполнялась в соответствии с плановой тематикой ААНИИ. Полученные результаты могут быть использованы:

- при исследованиях изменений климата Арктики;

- при разработке прогностических моделей эволюции морского ледяного покрова в Арктическом бассейне;

- в виде отдельного самостоятельного термодинамического блока в совместных динамико-термодинамических моделях морского ледяного покрова;

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

- для оценки прочностных характеристик морского льда.

«М1<1>.

Научная новизна диссерташга

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

2. Впервые предложена, основанная на гипотезе об автомодельности вертикального профиля, параметризация для годового цикла солености однолетнего морского льда.

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

4. В работе, на одной расчетной базе с моделью ровного льда, создана оценочная термодинамическая модель торосистого образования.

Дпрпбяци» работы

Основные результаты работы докладывались на итоговых сессиях ученого совета ААНИИ (1999-2004) и РГГМУ (1999), на семинарах отдела взаимодействия океана и атмосферы ААНИИ (1999-2005), на конференции молодых ученых, посвященной 140 летаю Н.М. Книповича (ПИНРО, 2002), на совместных российско-норвежских семинарах по проекту «Транспорт и судьба загрязнений в Северных морях» (ААНИИ-НПИ, 1999-2001), и в рамках совместной российско-норвежской лаборатории по исследованию климата Арктики им. «Фрама». Некоторые, полученные в ходе исследования результаты опубликованы в трудах ААНИИ и журнале «Метеорология и гидрология»

Публикации

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

Основные положения, выносимые и? тяп^п-у

1. Усовершенствованная термодинамическая модель морского льда с учетом реального профиля солености.

2. Метод усвоения притока тепла солнечной радиации в численной термодинамической модели морского льда с учетам снежного покрова.

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

4. Модельные оценки скорости промерзания и таяния торосистого образования в зависимости от коэффициента заполнения.

5. Результаты тестирования предложенной термодинамической модели к различным физико-географическим условиям.

Структура ди^-гергятрга

Диссертационная работа состоит из введения, четырех глав и заключения. Объем диссертации составляет 166 стр. и включает 33 рис., 11 табл. Список цитируемой литературы содержит 110 наименований.

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

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

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

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

Параграф 1.2 посвящен описанию общих закономерностей таяния морского льда. Указываются основные природные механизмы, определяющие как таяние льда, так и его интенсивность.

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

В параграфе 1.4 рассмотрен вопрос о фазовом составе морского льда и методах его описания. Указываются источники экспериментальных данных, по описанию соотношения фаз в морских льдах и приводится, по-видимому, наиболее удобная для использования в расчетах, связанных с учетом фазовых переходов во льду аппроксимация соотношения фаз предложенная Ю.Л. Назинцевым.

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

В параграфе 2.1 приводится определение структуры и текстуры морского льда, рассматривается структурно-генетическая классификация морских льдов Н.В.Черепанова.

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

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

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

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

В настоящее время, широкое применение нашел следующий способ расчета коротковолнового радиационного баланса поверхности (Г) и проникающей в толщу среды солнечной радиации (/о):

р = р0а-«)а-«, а)

/0 =

(2)

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

Затухание солнечной радиации с глубиной описывается в соответствии законом ослабления Буггера-Ламберта:

'г = 10е'кг, (3)

где 1г - солнечная радиация, проникающая на глубину г; к - коэффициент ослабления солнечной радиации в слое от поверхности до глубины равной г.

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

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

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

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

где г в этом случае равно шагу расчетной сетки по вертикали.

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

где а - альбедо поверхности льда, выраженное в процентах, у - шаг сетки по вертикали, м.; С - константа, имеющая размерность - м (зависящая, по-видимому, от структуры и/или текстуры льда). В первом приближении численное значение константы С получилось равным 0,003722 м.

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

где анЬ- эмпирические функции плотности снега.

При подстановке плотности снега в кг/м3 (в интервале изменения плотности снега (А) от 250 до 550 кг/м3 и альбедо от 85 до 40 %), эти функции имеют следующий вид:

(l-io) = 1-е ь или ig =e~b ,

(4)

(5)

а = -(45,0 - 0,06 р,), Ъ = -(2,35 - 0,003 -Pj),

(7)

(8)

Проверка полученного соотношения осуществлялась по результатам расчетов ослабления в верхнем слое снега. Для этого использовались экспериментальные

значения, полученные в ходе экспедиционных исследований автора на НЭС "Михаил Сомов" (май 2001 г., дрейфуюший лед в районе Земли Франца-Иосифа), а также на припайных льдах и в тундре о. Западный Шпицберген (апрель - май 2002, 2003, 2004 гг.). Показано, что рассчитанные по выражению (6) и измеренные значения величины коэффициента пропускания для верхнего слоя согласуются вполне удовлетворительно. Предложенная функция описывает 59 и 74 % дисперсии указанных экспериментальных выборок для слоев 5 и 10 см соответственно.

Третья глава посвящена формулировке и описанию одномерной термодинамической модели морского льда.

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

ЗГХ ЭГ2ЭГ, ^

; 0<г1<А(<)

(9)

<?Г, Э <л ог2

{ Эг2

; 0 <г2<Н(0

Со стандартными для задач, связанных с морским ледяным покровом граничными условиями:

21=0 Л,

= Ф;

(10)

г, =0

¿Г,

(П)

г2=0

г2=0 Г1(Г,А) = Г2(Г,0);

(12)

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

г2=Я Т2(1,Н) =0;

(13)

дН _ 3 <?Г2

ел <71г

"Фа

(14)

г2=И

где: индекс 1 - обозначает характеристики снега; 2 - характеристики морского льда; г - время; Т - температура; г - вертикальная координата; Ь - эффективная теплота плавления; Н - толщина льда; к - толщина снега; с - коэффициент эффективной теплоемкости; р - плотность; X - теплопроводность; Ф*, - поток тепла от воды; Ф -суммарный поток тепла на границе снег (лед)- атмосфера; 1о - поток коротковолновой солнечной радиации, проникающей в среду; © - температура замерзания морской воды.

При таянии снега, верхнее граничное условие приобретает вид:

г, =0 Г,=0,

(15)

и добавляется уравнение, описывающее таяние снежного покрова:

(

Эй

1

Э* 1лР\

А

¿1=0

(16)

С началом таяния льда (снег полностью растаял) система уравнений принимает

вид:

ЭГ2 Э

(17)

г2=0 Г2 = 0;

(18)

г2=Я Г2(Г,Я) = 0;

(19)

дН = 1 Э/

Л

дт2 Эг2

-Ф«

г2=Н

ЬгРг

дгг

(20)

гг=0

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

ЭГ , д2Т, ЭХ дТ: Э/„

(21)

Уравнение теплопроводности в виде (21) отличается от ранее используемых при моделировании уравнений двумя дополнительными нелинейными слагаемыми дХ дГ Э/0

——— и ——. Таким образом уравнении теплопроводности морского льда, 01 Л Эг

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

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

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

В параграфе 3.4 формулируется разработанная параметризация эволюции профиля солености однолетнего морского льда, основанная на гипотезе об автомодельности профиля солености в морском льду.

Так, во время активного роста морского льда, профиль солености можно описать зависимостью вида:

5 = (0,981 -1,472• Т] + 3,751-т]2 -5,662• т]3 + 3,461 V). (22)

гДе $тах - соленость на нижней границе морского льда; Т] = — - безразмерная

Н

вертикальная координата.

Для получения величины солености морского льда на его нижней границе, в зависимости от поставленной задачи, можно воспользоваться известной формулой В.Л. Цурикова или иными выражениями.

Во время начала таяния снега (лед продолжает расти под действием запаса холода):

5 = • (0.68 + 5,74 • ц - 34,97 -т? + 81,66 Г}3 - 84,3 • Т]* + 32,19 У), (23)

В период таяния ледяного покрова (толщина льда более 0,7 м):

5 = 5^- (0,02 - 0,12 • ?7 - 4,44 • 772 + 50,78 • 7]ъ -90,99 У + 45,44 У), (24)

здесь Бпех - соленость среднего слоя льда (слоя максимальной солености).

Соленость срединного слоя, видимо, необходимо выбирать, опираясь на характерные значения максимальной солености льда для данного периода и данного географического объекта. Однако следует отметить, что для арктических морей, как показывают измерения, типичная максимальная соленость, характерная для описанного профиля, находится в пределах величины 3,0-4,0 %с.

Значение коэффициентов корреляции значений солености вычисленных по формулам (22, 23 и 24) с приведенными в работе В.Л.Цурикова более 0,99. Это говорит о качественном воспроизведении степенными многочленами исходных данных.

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

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

В параграфе 3.5. как частный случай термодинамической модели - приводится формулировка оценочной модели промерзания и таяния торосистого образования.

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

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

Система уравнений, для вычисления промерзания торосистого образования (при отсутствии слоя снега на поверхности), в одномерном представлении, имеет следующий вид:

ЭГ Эг

= (а 2 ■у1+аа

1 э/„

Эг2 с2р2 Эг

О2г<й;

(25)

ЭГ Эг

Э 2Г 1 (дЛгдТ Э/0

Эг2 сгРг \ дг Эг Эг

(26)

Граничные условия: ЭГ

2=0 Лз

Эг

= Ф;

г=0

(27)

г = Л (Л-у.+Я^а-^))^

г = А Г(г, гАЧ)) = Г(г,гЛ+о); г = Я Гг(г,Я) = 0;

г=А-0

г=А+0

(28)

(29)

(30)

дН =_1__

1 3 ЭГ

-Фи

(31)

где у - коэффициент заполнения надводной (1) или подводной (2) части тороса; Л -толщина надводной части тороса; Н - толщина консолидированного слоя; а -коэффициент температуропроводности льда или воздуха; Л - коэффициент теплопроводности льда или воздуха.

При наличии слоя снега, в систему добавляется уравнение теплопроводности снега как в (9), граничное условие (27) записывается в виде (10) и в систему добавляется условие склейки потоков на границе снег - надводная часть тороса в виде, подобном (11-12) с учетом коэффициента заполнения. Таяние торосистого образования сверху вычисляется аналогично таянию ровного морского льда с учетом коэффициента заполнения паруса тороса.

В параграфе 3.6 описывается специально поставленный лабораторный опыт по намораживанию морского льда из раствора №С1 в бассейне, по результатам которого, на начальном этапе, тестировалась разработанная термодинамическая модель морского льда. Приведено сравнение толщин лабораторного морского льда, измеренных в ходе эксперимента, с рассчитанными по модели (коэффициент корреляции превысил 0,99).

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

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

В параграфе 4.2 приводятся результаты сравнения модельных расчетов и натурных измерений вертикального распределения температуры морского льда по данным СП-13. Показано хорошее соответствие вычисленных и измеренных значений для холодного периода года (коэффициент корреляции 0,97 - 0,99), и неплохое для периода таяния (коэффициент корреляции 0,75 - 0,80).

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

Параграф 4.4 посвящен оценке влияния некоторых природных характеристик, учитываемых моделью, на изменение толщины морского льда. На примере припайных льдов районов о. Визе и Кандалакшского залива показано влияние гидрометеорологических условий ледообразования на скорость роста морского льда. Получено, что лед вертикально-волокнистой структуры, формирующийся при относительно спокойных погодных условиях растет на 13-15% быстрее льда зернистой структуры, формирующегося в динамичных погодных условиях. Отмечено, что это влияние сказывается лишь в период роста морского льда. На примере расчетов для СП-13 и Кандалакшского залива установлено, что при расчетах морского льда в период таяния важнейшую роль играет величина коэффициента пропускания. Показано, что при некорректном учете его величины

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

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

В параграфе 4.5 приводится сравнение скоростей роста ровного льда и промерзания торосистого образования. Показано, что торос промерзает быстрее, чем ровный лед. Причем с ростом коэффициента заполнения надводной и подводной частей скорость промерзания увеличивается. Сравнение скоростей таяния паруса тороса и ровного льда показало, что торос тает быстрее и скорость таяния увеличивается при уменьшении коэффициента заполнения паруса.

В заключении сформулированы основные результаты диссертационной работы:

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

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

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

коэффициента теплопроводности и наличие внутреннего источника тепла). Выполнена модельная оценка чувствительности параметров морского льда к изменениям внешних условий, выделен ряд наиболее влияющих на изменение толщины морского льда внешних факторов. Осуществлена проверка работоспособности и адекватности воспроизведения моделью термических условий и динамики изменения толщины морского льда и снега для различных географических объектов Арктики (Кандалакшский залив Белого моря, районы о. Визе и о. Диксон в Карском море, залив Гренфиорд (арх. Шпицберген), дрейфующая станция СП-13).

4. Показано влияние гидрометеорологических условий ледообразования на скорость роста морского льда. Получено, что лед вертикально-волокнистой структуры, формирующийся при относительно спокойных погодных условиях растет на 13-15% быстрее льда зернистой структуры, формирующегося в динамичных погодных условиях. Отмечено, что это влияние сказывается лишь в период роста морского льда.

5. Установлено, что при расчетах морского льда в период таяния важнейшую роль играет величина коэффициента пропускания. Показано, что при некорректном учете его величины в численной термодинамической модели, разница в величинах толщин тающего толстого морского льда может достигать 20-25 %. При относительно тонких, сезонных льдах, серьезно смещаются даты схода снежного покрова и взлома припая.

6. Оценена величина потока соли от растущих льдов различной толщины. Получено, что поток соли ото льдов толщиной 50-70 см более чем вдвое превосходит поток соли ото льдов 150-170 см и более чем вчетверо, чем ото льдов 280-300 см толщины.

7. Разработана простая термодинамическая модель промерзания и таяния торосистого образования. На ее основании проведены оценки скоростей промерзания тороса в зависимости от коэффициента заполнения его надводной (парус) и подводной (киль) частей. Показано, что скорость промерзания тороса с коэффициентом заполнения подводной части 0,6 в 1,4 раза, а с коэффициентом 0,9 в 3,4 раза превышают скорость роста ровного льда при одинаковых условиях. Оценена скорость таяния надводной части тороса в зависимости от коэффициента заполнения и альбедо поверхности. Получено, что надводная часть тороса тает быстрее даже при одинаковых значениях альбедо поверхности льда и тороса, чем ровный лед, причем эта скорость растет по мере уменьшения коэффициента заполнения.

По теме диссертации опубликованы следующие работы:

1. Параметризация радиационных процессов в модели ледяного покрова.-Метеорология и гидрология. -2001. -№ 2. -С. 81-88 (вместе с Б.В.Ивановым).

2. О решении задачи усвоения коротковолновой радиации в толще снега.-Метеорология и гидрология. -2001. -№ 12. -С. 65-69 (вместе с Б.В.Ивановым).

3. Параметризация переноса коротковолновой солнечной радиации в снежно-ледяном покрове.- Метеорология и гидрология. -2003. -№ 2. -С. 54-58 (вместе с Б.В.Ивановым).

4. Предварительные результаты исследования свойств снежного покрова архипелага Шпицберген в мае 2002 г. Комплексные исследования природы Шпицбергена. - Изд. КНЦ РАН, г. Аппатиты. -2003. -вып.З. -С. 124-127 (вместе с Б.В.Ивановым, А.М.Безгрешновым).

5. Рост толщины льда в Кандалакшском заливе Белого моря // Материалы Ш международного семинара "Рациональное использование прибрежных районов полярных морей", 17 июля 1998. Кандалакша. -С-Петербург. -1999. -С. 81-95 (вместе с Н.В.Кубышкиным).

6. Термодинамическое моделирование морского ледяного покрова в Арктическом бассейне / Тезисы докладов "Всероссийской конференции молодых ученых, посвященной 140-летию со дня рождения Н.М.Книповича", 23-25 апреля 2002 Мурманск. -Мурманск. -2002. -С. 10-11.

7. Особенности термической структуры арктических снежниц. - Труды ААНИИ т.446. -СПб. -Гидрометеоиздат. -2003. -С. 176-184. (вместе с Б.В.Ивановым)

8. Исследования радиационных и теплофизических свойств снежно-ледяного покрова в районе ледника Альдегонда. В сб.: «IV Международная конференция -Комплексные исследования природы Шпицбергена». -Изд. КНЦ РАН г. Апатиты. -2004. -вып.4. -С. 199-206 (совместно с Б.ВИвановым, П.Н.Священниковым, Р.Е.Власенковым)

9. Radiation and thermophisical characteristics of snow and sea ice Covers in the Arctic areas / Arctic-alpine ecosystems and people in a changing environment. Polar Environmental Centre, Tromso, Norway. 24.02-01.03.2003. -Book of abstracts. -P. 79. (with B.V. Ivanov)

Ротапринт ААНИИ. Заказ № 17. Тираж 100 экз.

№166 i 8

РНБ Русский фонд

2006-4 13559

Содержание диссертации, кандидата физико-математических наук, Андреев, Олег Михайлович

Введение

1. Закономерности роста и таяния морского льда

1.1 Нарастание морского льда

1.2 Таяние морского льда

1.3 Соленость морского льда

1.4 Характеристики фазовых переходов в морском льду

2. Свойства морского льда и снега с учетом особенностей энергомассообмена в высоких широтах

2.1 Особенности строения морского льда

2.1.1 Структура и текстура льда

2.1.2 Структурно-генетическая классификация Н. В. Черепанова

2.2 Теплофизические характеристики морского ледяного покрова

2.3 Структурные и текстурные особенности, а также теплофизические характеристики снежного покрова в Арктическом бассейне

2.4 Энергообмен между атмосферой и океаном в Арктическом бассейне при наличии снежно-ледяного покрова

3. Одномерная нестационарная термодинамическая модель эволюции морского снежно-ледяного покрова

3.1 Постановка задачи

3.2 Разностная аппроксимация

3.3 Расчет теилофизических характеристик и потоков тепла

3.3.1 Расчет теилофизических характеристик

3.3.2 Расчет потоков тепла

3.4 Параметризация профиля солености однолетнего морского льда

3.5 Оценочная термодинамическая модель эволюции торосистого образования

3.6 Тестирование одномерной нестационарной термодинамической модели морского льда по данным лабораторного эксперимента

4. Результаты моделирования

4.1 Проверка работоспособности блоков модели и оценка чувствительности модели к изменению внешних параметров

4.1.1 Проверка работоспособности термической части снегового блока модели

4.1.2 Проверка работоспособности термического блока морского льда

4.1.3 Оценка чувствительности модели снежно-ледяного покрова к вариациям внешних параметров

4.2 Тестирование одномерной термодинамической модели морского льда по данным дрейфующей станции СП

4.3 Проверка адекватности воспроизведения толщины морского льда по натурным данным

4.3.1 Расчеты для растущего припайного льда залива Гренфиорд о. Западный Шпицберген)

4.3.2 Расчеты по модели для припайных льдов бухты о. Диксон

4.4 Оценки влияния некоторых природных характеристик, учитываемых моделью, на изменение толщины морского льда

4.4.1 Оценки влияния различных гидрометеорологических условий образования морского льда на скорость его роста

4.4.2 Модельные оценки влияния выбора коэффициента пропускания на скорость таяния морского льда

4.4.3 Модельные оценки потока тепла от подледного слоя морской воды и оценки потока соли ото льдов различной толщины

4.4.4 Оценка скорости роста соленого и пресного льдов по результатам модельных расчетов

4.5 Модельные оценки скоростей промерзания и таяния торосистого образования

Введение Диссертация по наукам о земле, на тему "Термодинамическое моделирование формирования морского ледяного покрова в Арктике"

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

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

Первые попытки математического моделирования промерзания среды относятся к середине XIX века, и наиболее серьезным исследованием в этом направлении следует признать работу И.О. Стефана (1891) /1/, в которой были заложены основы современных методов термодинамического моделирования взаимодействующих сред с движущейся границей раздела фаз. Однако наибольшая исследовательская активность началась в 40-60 годах XX века, когда появились работы Л.Л. Шепелевского /2/, Л.И. Рубинштейна /3/ и В.Г. Меламеда /4/, продемонстрировавшие возможности термодинамического подхода применительно к моделированию морского льда. Отдельно следует выделить классическую работу Ю.П. Доронина (1963) /5/. Эта термодинамическая модель морского льда и в настоящее время является одной из наиболее часто используемых. Появившиеся несколько позже модели зарубежных авторов Мэйкута-Унтерштейнера (1971) /6/ и Семтнера (1976) III стали основой для современных термодинамических моделей морского льда, используемых для различных теоретических исследований, например /8/.

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

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

Основные задачи исследования

1. Разработать численную термодинамическую модель морского льда.

2. Определить вклад различшлх факторов в формирование морского льда, оценить скорость промерзания и таяния тороса.

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

4. Исследовать применимость модели к различным физико-географическим условиям Арктики.

Научная новизна диссертации

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

2. Впервые предложена, основанная на гипотезе об автомоделыюсти вертикального профиля, параметризация для годового цикла солености однолетнего морского льда.

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

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

Научная и практическая значимость

Диссертационная работа выполнялась в соответствии с плановой тематикой ААНИИ. Полученные результаты могут быть использованы:

- при исследованиях изменений климата Арктики;

- при разработке прогностических моделей эволюции морского ледяного покрова в Арктическом бассейне;

- в виде отдельного самостоятельного термодинамического блока в совместных динамико-термодинамических моделях морского ледяного покрова;

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

- для оценки прочностных характеристик морского льда.

Апробация работы

Основные результаты работы докладывались на итоговых сессиях ученого совета ААПИИ (1999-2004) и РГГМУ (1999), на семинарах отдела взаимодействия океана и атмосферы ААНИИ (1999-2005), на конференции молодых ученых, посвященной 140 летию Н.М. Книповича (ПИНРО, 2002), на совместных российско-норвежских семинарах по проекту «Транспорт и судьба загрязнений в Северных морях» (ААНИИ-НПИ, 1999-2001), и в рамках совместной российско-норвежской лаборатории по исследованию климата Арктики им. «Фрама». Некоторые, полученные в ходе исследования результаты опубликованы в трудах ААПИИ и журнале «Метеорология и гидрология»

Публикации

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

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

1. Усовершенствованная термодинамическая модель морского льда с учетом реального профиля солености.

2. Метод усвоения притока тепла солнечной радиации в численной термодинамической модели морского льда с учетом снежного покрова.

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

4. Модельные оценки скорости промерзания и таяния торосистого образования в зависимости от коэффициента заполнения.

5. Результаты тестирования предложенной термодинамической модели к различным физико-географическим условиям.

Структура диссертации

Диссертационная работа состоит из введения, четырех глав и заключения.

Заключение Диссертация по теме "Океанология", Андреев, Олег Михайлович

ЗАКЛЮЧЕНИЕ.

Завершая описание проведенного исследования по термодинамическому моделированию морского льда отметим некоторые основные моменты, которые удалось решить в рамках данной работы:

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

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

3. Рассмотрены преимущества и недостатки существующих ранее термодинамических моделей морского ледяного покрова. Разработана усовершенствованная термодинамическая модель морского льда, учитывающая три источника нелинейности (зависимость теплофизичсских характеристик ог температуры и солености, учет вертикального градиента коэффициента теплопроводности и наличие внутреннего источника тепла). Выполнена модельная оценка чувствительности параметров морского льда к изменениям внешних условий, выделен ряд наиболее влияющих на изменение толщины морского льда внешних факторов. Осуществлена проверка работоспособности и адекватности воспроизведения моделью термических условий и динамики изменения толщины морского льда и снега для различных географических объектов Арктики (Кандалакшский залив Белого моря, районы о. Визе и о. Диксон в Карском море, залив Гренфиорд (арх. Шпицберген), дрейфующая станция СП-13).

4. Показано влияние гидрометеорологических условий ледообразования на скорость роста морского льда. Получено, что лед вертикально-волокнистой структуры, формирующийся при относительно спокойных погодных условиях растет на 13-15% быстрее льда зернистой структуры, формирующегося в динамичных погодных условиях. Отмечено, что это влияние сказывается лишь в период роста морского льда.

5. Установлено, что при расчетах морского льда в период таяния важнейшую роль играет величина коэффициента пропускания. Показано, что при некорректном учете его величины в численной термодинамической модели, разница в величинах толщин тающего толстого морского льда может достигать 20-25 %. При относительно тонких, сезонных льдах, серьезно смещаются даты схода снежного покрова п взлома припая. Произведен опосредованный учет влияния снежниц на таяние морского льда в одномерной термодинамической модели. Получено, что разница в величинах толщин тающего морского льда с учетом и без учета влияния снежниц может превышать 15 %.

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

7. Оценена величина потока соли от растущих льдов различной толщины. Получено, что поток соли ото льдов толщиной 50-70 см более чем вдвое превосходит поток соли ото льдов 150-170 см и более чем вчетверо, чем ото льдов 280-300 см толщины.

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

9. Разработана простая термодинамическая модель промерзания и таяния торосистого образования. На ее основании проведены оценки скоростей промерзания тороса в зависимости от коэффициента заполнения его надводной (парус) и подводной (киль) частей. Показано, что скорость промерзания тороса с коэффициентом заполнения подводной части 0.6 в 1.4 раза, а с коэффициентом 0.9 в 3.4 раза превышают скорость роста ровного льда при одинаковых природных условиях. Оценена скорость таяния надводной части тороса в зависимости от коэффициента заполнения и альбедо поверхности. Получено, что надводная часть тороса тает быстрее даже при одинаковых значениях альбедо поверхности льда и тороса, чем ровный лед, причем эта скорость растет по мере уменьшения коэффициента заполнения.

Библиография Диссертация по наукам о земле, кандидата физико-математических наук, Андреев, Олег Михайлович, Санкт-Петербург

1. Stefan J. Uber die Theorie der Eisbildung, insbesondere iiber Eisbildung im Polarmeere // Annln Phys.- band 42.- N2.- P. 269-286.

2. Шепелевекий A.A. Скорость нарастания пресного льда при заданной температуре верхней его поверхности // Проблемы Арктики.- 1940.- № 10.

3. Рубинштейн Л.И. К вопросу о единственности решения одномерной задачи Стефана в случае однофазного начального состояния теплопроводящей среды //Докл. АН СССР.- 1957.- т. 117.- №3.

4. Меламед В.Г. Сведение задачи Стефана к системе обыкновенных дифференциальных уравненпй.-Изд. АН СССР, сер. геофизика.- 1958.- №7.

5. Доронин Ю.П. Численный метод расчета толщины и температуры льда // Проблемы Арктики и Антарктики.- 1963.- №14.- с. 17-25.

6. Maykut G.A., Untersteiner N. Some results from a time dependent thermodynamic model of sea ice // J.Geophys. Res. 1971. -Vol.76, N 6. - P. 1550-1575.

7. Semtner A.J. A model for the thermodynamic growth of sea ice in numerical investigations of climate // J. Phisic. Oceanogr.-1976.- N6.- P 379-389.

8. Ebert E.E., Curry J.A. An intermediate one-dimensional thermodynamic sea ice model for investigating ice-atmosphere interactions // J.Geophys. Res. 1993. -Vol.98, NC6.-P. 10085-10109.

9. Доронин Ю.П. Тепловое взаимодействие атмосферы и гидросферы в Арктике. Л.: Гидрометеоиздат, 1969. - 299 с.

10. Доронин Ю.П. Рост и таяние морского льда // Морской лед. Сбор и анализ данных наблюдений, физические свойства и прогнозирование ледовых условий: справочное пособие. СПб.: Гидрометеоиздат, 1997. - С. 107-125.

11. Лыков А.В. Теория теплопроводности.-М.: Гостехиздат, 1952.-357 с.

12. Савельев Б.А. Строение, состав и свойства ледяного покрова морских и пресных водоемов. М.: Изд. МГУ, 1963. - 541 с.

13. Доронин Ю.П., Хейсин Д.Е. Морской лсд.- JI.: Гидромстеоиздат, 1975. -320 с.

14. Зубов H.H. Льды Арктики.- М.: Изд. Главсевмориути, 1945.-424 с.

15. Гаврило В.П., Сухоруков К.К. Исследование динамики образования морского льда гидроакустическим методом // Труды ААНИИ.- 1980, т.374.-с. 85-96.

16. Кузнецов И.М., Гарманов А.Л. О практическом использовании формул для расчета толщины морского льда // Метеорология и гидрология.- 1987.-Ко 11.-е. 116-120.

17. Доронин Ю.П. К вопросу о нарастании морского льда // Проблемы Арктики и Антарктики.- 1959.- вып. 1.- с.78-83.

18. Эфроимсон В.О. К вопросу нарастания льда // Метеорология и гидрология.- 1977.- №8.- с. 43-54.

19. Андреев О. М., Иванов Б. В. Параметризация радиационных процессов в модели ледяного покрова // Метеорология и гидрология.- 2001.- № 2.- с. 81-88.

20. Доронин Ю.П., Кубышкин И.В. Рост и таяние морского льда,- СПб.: Гидромстеоиздат.- 2001.- 44 с.

21. Морской лед. Справочное пособие // под ред. И. Е. Фролова, В. П. Гаврило,- СПб., Гидромстеоиздат.- 1997.- 402 с.

22. Цуриков В.Л. Жидкая фаза в морских льдах. М.: Наука, 1976.- 210 с.

23. Мальмгрен Ф. О свойствах морского льда. М.: Изд. Гидрогр. упр. и Гидрометеокомитета, 1930. - 90 с. - (Пер. с англ.).

24. Назинцев ЮЛ., Панов В.В. Фазовый состав и теплофизические характеристики морского льда.- СПб.: Гидрометеоиздат.-2000.- 84 с.

25. Untersteiner N. Natural desalination and equilibrium salinity profile of perennial sea ice//J. Geophys. Res.- 1968.- Vol. 73(4).- P. 1251-1257.

26. Weeks W.F., Ackley S.F. The growth, structure and properties of sea ice // U.S.Army CRREL Monograh. 1982. - 82-1. - 130 p.

27. Рывлии Ф.Я. Метод прогнозирования предела прочности ледяного покрова на изгиб // Проблемы Арктики и Антарктики.-1974.- Вып. 45.- С. 7986.

28. Cox G.F.N., Weeks W.F. Salinity variations in sea ice // J.Claciol. -1974. -Vol.23, N67.-P. 109-120.

29. Шестсперов И.А. Проверка формулы Цурнкова но результатам наблюдений // Океанология. 1969. -Т.9.- Вып. 4. - С. 616-618.

30. Цыцарип А.Г. Гидрохимические аспекты ледообразования. Опыт экспериментального моделирования // Исследование океанов и морей: Юбилейный вып.2. СПб.: Гидрометеоиздат, 1995. - С.224-249.

31. Assur A., Weeks W.F. Growth, structure and strength of sea ice // U.S. Army CRREL Res.Rep. 1964. - N 135. - P. 1-19.

32. Weeks W.F., Lofgren G. The effective solute distribution coefficient during the freezing of NaCl solutions // Phys. snow and ice. 1967. - 1, Pt. 1. - P.579 -597.

33. Burton I.A., Prim R.C., Slichter W.P. The distribution of solute in crystals growth from the melt. Part I. Theoretical // J.Chemical Physics. 1953. - Vol.21, N 11.-P. 1987-1991.

34. Даричева JI.B., Молчанова Л.А., Чупрынин В.И. Модель расчета солености морского льда / РАН. Тихоокеан. институт географии. -Владивосток, 1996. 16 с. - Деп. в ВИНИТИ 20.09.96, № 3346-В96.

35. Лялягин Н.А., Тышко К.П. Влияние кристаллической структуры льда на формирование и изменчивость вертикального распределения солености в однолетнем ледяном покрове // Метеорология и гидрология.- 2002.- № 2 с. 64-71.

36. Черепанов Н.В. Классификация льдов природных водоемов // Тр.ААНИИ. 1978. - Т. 331. - С. 77-99.

37. Назннцев Ю.Л., Дмитраж Ж.А., Моисеев В.И. Теплофизические свойства морского льда,- Л.: Изд.ЛГУ, 1988.- 260 с.

38. Сох G.F.N., Weeks W.F. Equations for determining the gas and brine volumes in sea-ice samples //J. Glaciol. 1983. - Vol. 29, N 102. - P. 306-316.

39. Ringer W.E. Changes in the composition of sea water upon freezing (in Dutch) // Chem. Weekblad. 1906. - Vol. 3. - P. 223-249.

40. Assur A. Composition of sea ice and its tensile strength // Arctic sea ice. -1958.-N598.-P. 106-138.

41. Schwerdtferger P. The thermal properties of sea ice // J.Glaciol. 1963.-Vol. 4, N36. -P. 789-807.

42. Nelson K.II. A study of the freezing of sea water // Ph. D. Thesis. University of Washington, 1953.- 129 p.

43. Nelson K.H., Thompson T.G. Deposition of salts from sea water by frigid concentration // J.Mar. Res. 1954. - Vol 13, N2. - P. 166-182.

44. Тышко К.П., Федотов В.И., Черепанов Н.В. Кристаллическое строение морского ледяного покрова.- СПб.: Гидрометеоиздат.- 2000.- 67 с.

45. Lepparanta М, Manninen Т. The brine and gas content of sea ice with attention to low salinities and high temperatures // Internal report Finn. Inst of Mar. Res., Helsinki.- 1988.

46. Назпнцев IO.JI. Модели теплопроводности природных льдов // Тр.ЛЛНИИ.- 1976.-Т. 331.-С. 104-115.

47. Назпнцев IO.JI. Экспериментальные определения теплоемкости и теплопроводности морского льда // Пр. Арктики и Антарктики. 1959. -Вып. 1. - С.65-72.

48. Радионов В. Ф., Брязгин И. II., Александров Е. И., Снежный покров в арктическом бассейне.- СПб., Гидрометеоиздат, 1996.- 124 с.

49. Снег. Справочник // под ред. Г. М. Грея, Д. X. Мейла, пер. с англ.- Л., Гидрометеоиздат, 1986.- 752 с.

50. Anderson Е.А. A point energy and mass balance model of snow cover // NOAA Tech. Rep.- NWS 19, Washington D.C.- 1976.

51. Yen Y.-C. Review of thermal properties of snow, ice and sea ice.- Cold Regions Research Engineering Laboratory (CRREL) Report 81-10.- Hanover, NII.-26 pp.

52. Осокин Н.И., Самойлов P.C. и др. К оценке влияния изменчивости характеристик снежного покрова на промерзание грунтов // Криосфера земли.-1999.- т 3.- №1.- с 3-10.

53. Федосеева В.И., Макаров B.II., Федосеев Н.Ф. О минерализации снежного покрова // Метеорология и гидрология.- 1986.- №4.- с. 75-79.

54. Иванов Б.В., Макштас Л.П. Квазистационарная нульмерная модель морского льда//Тр. ЛЛНИИ.- 1990.- т.420.- с. 18-31.

55. Макштас Л.П. Тепловой бюджет арктических льдов в зимний период.- Д.: Гидрометеоиздат, 1984.- 147 с.

56. Brunt D. Notes on radiation in the atmosphere // Q.J.R. Meteorol. Soc. 58.- P. 389-420.

57. Маршупова M.C. Основные закономерности радиационного баланса подстилающей поверхности и атмосферы в Арктике // Труды ААНИИ.-1961.- т.229.- с.56-67.

58. Maykut G.A., Church Р.Е. Radiation climate of Barrow, Alaska, 1962-66 // J. Appl. Meteorol.- 12.- 1973.- P. 620-628.

59. Satterlund D.R. An improved equation for estimating long-wave radiation from the atmosphere // Water Ressour. Res. -15.- 1979,- P. 1649-1650.

60. Кузьмин П.П. Таяние снежного покрова.- JL: Гидрометеоиздат, 1961.

61. Ефимова Н.А. Метод вычисления месячных величин длинноволнового радиационного баланса // Метеорология и гидрология.- 1961.- №10.- с. 2833.

62. Idso S.B., Jackson R.D. Thermal radiation from the atmosphere // J. Geoph. Res.-74(23).- 1969,- P. 5397-5403.

63. Idso S.B. A set of equations for full spectrum and 8-11 microns and 10.5-12.5 microns thermal radiation from cloudless skies // Water Resource Res.- 17.- 1981.-P. 295-304.

64. Guest P.S. Surface radiation conditions in the Eastern Weddell Sea during winter//J. Geoph. Res.- 103(C13).- 1998,- P. 30761-30771.

65. Makshtas A. P., Andreas E. L., Svyashchennikov P. N., Timachev V. F. Accounting for Clouds in Sea Ice Models.- Cold Regions Research and Engineering Laboratory.- vol. 98-9.- 30 p.

66. Zillman J.W. Л study of some aspects of the radiation and heat budgets of the southern hemisphere oceans // Meteorol. Stud. Rep 26, Bur. of Mcteorol., Dep of the Inter., Canberra, Л.С.Т.- 1972.

67. Shine K.P. Parameterization of short-wave flux over high albedo surfeccs as a function of cloud thickness and surface albedo // Q.J.R. Meteorol. Soc.- 110.1984.- P. 747-764.

68. Андреев О. M., Иванов Б. В. О решении задачи усвоения коротковолновой радиации в толще снега // Метеорология и гидрология.- 2001.- № 12.- с. 65-69.

69. Brandt R. Е., Warren S. G. Solar heating rates and temperature profile in Antarctic snow and ice//J. Glaciology.- 1993.-vol. 39.-№ 131.-P. 99-110.

70. Zeebe R. E., Eicken H, et al. Modeling the heating and melting of sea ice through light absorption by microalgae // J. Geophus. Res.- 1996.- vol. 101.- P. 1163-1181.

71. Maykut G. A. The surface heat and mass balance. / In: The Geophysics of Sea Ice (N. Untersteiner, Ed.).- Plenum Press, NY.- 1986.- P. 385-465.

72. Макштас А. П., Тимачев В. Ф. Чувствительность термодинамической модели морского льда к параметризации коротковолновой и длинноволновой радиации //Труды ААНИИ.- 1992.- т. 430.- с. 116-137.

73. Богословский В. И., Кузнецов М. А. Структура радиационного потока на поверхности и в толще снега и льда // Труды советской антарктической экспедиции.- I960.-т. 10.- с. 101-106.

74. Бузуев А. Я., Шестериков II. П. Некоторые радиационные свойства однолетнего льда в период таяния // Проблемы Арктики и Антарктики.- 1964.-вып. 18.-с. 87-91.

75. Дубовцев В. Ф., Тпмерев А. А. Проникновение тепла солнечной радиации через снежно-ледяной покров Северного Ледовитого океана // Тр. ААНИИ.-1976.- т.319.-с. 81-95.

76. Клюева II.В. Проникновение солнечной радиации в снежно-ледяной покров в низовьях и устьях арктических рек // Проблемы Арктики и Антарктики.- 1978.- вып. 53.- с. 66-70.

77. Шестериков П.П. Поглощение солнечной радиации морским льдом // Труды ДАНИИ.- 1962.- т. 257.- с. 78-83.

78. Коптев Л.П., Пятпенков Б.Л. О поглощении и проникновении солнечной радиации в снег и лед в Арктике // Проблемы Арктики и Антарктики.- 1962.-вын.10.- с.71-76.

79. Богословский В.Н. Кузнецов М.Л. Структура радиационного потока на поверхности и в толще снега и льда // Труды Советской Антарктической эксп.- I960,-т. 10,-с. 101-106.

80. Gerland S., Winter J-G., Ivanov В., et al. Physical and optical properties of snow covering Arctic tundra on Svalbard // J. Hydrol. Process.- 1999.- vol. 13.- P. 2331-2343.

81. Winter J-G., Gerland S., Ivanov В., et al. Spectral reflectance of melting snow in a high Arctic watershed on Svalbard: some implications for optical satellite remote sensing studies // J. Hydrol. Process.- 1999.- vol. 13,- P. 2033-2049. .

82. Шестериков Н.П. О поглощении радиации льдом под снежницей // Проблемы Арктики и Антарктики.- 1967.- вып. 25.- с.66-70.

83. Спичкин В.А. Учет альбедо при таянии льда // Труды ААНИИ.- 1966.-Т.269.- с.71-78.

84. Шестериков Н.П. О некоторых закономерностях изменения альбедо льда в летний период в Арктике // Проблемы Арктики и Антарктики,- 1963.- вып. 14.-с. 27-31.

85. Налимов Ю.В., Тимерев А.А. Величины альбедо снежно-ледяного покрова в низовьях и устьях арктических рек // Метеорология и гидрология.- 1974.-№5.- с. 64-68.

86. Barry R. The parameterization of surface albedo for sea ice and its snow cover // Progress in Phis. Geogr.- 1996.- 20(1).- P. 63-79.

87. Нагурный А.П. Моделирование континентального и морского льда в моделях климата.- М.: Изд. АН СССР, сер. Метеорология и климатология, т.13.- 1986.- 79 с.

88. Light В., Eickcn II., Maykut G.A., Grenfell T.C. The effect of included particulates on the spectral albedo of sea ice // J. Geoph. Res.- 1998.- vol. 103.- N C12.- P. 27739-27752.

89. Perovich D.,Roesler С., Pegau W-S. Variability in arctic sea ice optical properties//J. Geoph. Res.- 1998.-vol 103.-NCI.-P. 1193-1208.

90. Красс M.C., Мсрзликин В.Г. Радиационная теплофизика снега и льда. -JL: Гидрометеоиздат, 1990.- 261 с.

91. Вабшдевич П.Н. Численные методы решения задач со свободной границей. -М.: Изд. МГУ, 1987. 164 с.

92. Самарский A.A. и др. Численное моделирование конвективно-диффузионных процессов прн фазовых переходах / A.A. Самарский, П.Н. Вабшцевич, O.JI. Илиев, А.Г. Чурбанов. М., 1992. - 34 с.

93. Самарский A.A. Теория разностных схем. М.: Наука, 1983. - 616 с.

94. Осокин Н.И., Самойлов P.C., Сосновский A.B., Жидков В.А. Влияние параметров снега на температурный режим снежной толщи при внутрисуточпых колебаниях температуры воздуха // материалы гляциологических исследований.- 2001.- № 91.-е. 71-74.

95. Сосновский A.B. Математическое моделирование промерзания грунта с учетом пространственно-временной изменчивости параметров снежного покрова // Материалы гляциологических исследований.- 2000.- №89.- с. SO-SS.

96. Sturm М., Holmgren J., Konig M., Morris К. The thermal conductivity of seasonal snow // J. of Glaciol.- 1997.- vol. 43.- N 143.- P. 26-42.

97. Gerland S., Haas C., Nicolaus M., Winter J-G. Seasonal development of structure and optical properties of fast ice in Kongsfjorden, Svalbard // reports on Polar and Marine Res.- 2004.- N 492.- P. 26-34.

98. Svendsen II. et. al. The physical environment of Kongsfjorden-Krossfjorden, an arctic fjord system in Svalbard // Polar Res.- 2002.- N 21(1).- P. 133-166.

99. Ilamre В et. al. Modeled and measured optical transmittance of snow-covered first-year sea icc in Kongsijorden, Svalbard // J. Geoph. Res.- 2004.- vol. 109.-C10006, doi:10.1029/2003JC001926.

100. Бузуев Л.Я., Шестериков П.П. Зависимость средней толщины припайного льда от торосистости // Проблемы Арктики и Антарктики.-1969,- Bi.in.32.- с. 30-36.

101. Schramm J., Flato G., Curry J. Toward the modeling of enhanced basal melting in ridge keels//J. Geoph. Res.- 2000.- vol. 105.- N Сб.-P. 14081-14092.

102. Lepparanta M., Lensu M., Kosloff P., Veitch B. The life story of a first-year sea ice ridge // Cold Reg. and Technology.- 1995.- N 23.- P. 279-290.

103. Марченко A.B. Влияние консолидации торосов на тепловые потоки из океана в атмосферу // Труды ААНИИ.- 2003.- т. 446.- с. 150-164.

104. Астафьев В.П., Сурков B.II., Трусков П.А. Торосы и стамухи Охотского моря.- СПб.: Прогресс-Погода, 1997.-184 с.

105. Hoyland K.V. Simulations of the consolidation process in first-year ice ridges // Cold Reg. Sci. Technol.-2002.- N34.- P. 143-158.

106. Петров И.Г., Федотов В.И., Черепанов II.В. О методике приготовления лабораторного льда с заданными физическими свойствами // Тр. ААНИИ. -1988.-Т.401.-С. 77-94.

107. Коптев А.П. О роли снежного покрова в процессах теплообмена // Проблемы Арктики и Антарктики,- 1966.- вып.22.- с. 82-89

108. Андреев О. М., Иванов Б. В. Параметризация радиационных процессов в модели ледяного покрова // Метеорология и гидрология.- 2001.- № 2.- с. 81-88.

109. Launiainen J., Cheng В. Modeling of ice thermodynamics in natural waterbodies // Cold Region. Sci. Tech.- 1998.- vol. 27.- P. 153-178.

110. Бузуев А.Я., Спичкин В. А. Роль снежниц в формировании пространственной неравномерности толщины многолетних льдов в зимний период // Проблемы Арктики и Антарктики.- 1977.- вып. 49.- с. 53-58.

111. Аппель И.Л., Гудкович З.М. Численное моделирование и прогноз эволюции ледяного покрова арктических морей в пероид таяния.- С-Пб: Гидрометеоиздат.-1992,- 144 с.