Содержание к диссертации
Введение
Глава 1. Физико-географические и климатические особенности Каспийского моря 18
Глава 2. Модель гидродинамики внутреннего моря 22
2.1. Уравнения и граничные условия 23
2.2. Разностная аппроксимация уравнений и граничных условий 25
2.3. Подмодель взаимодействия моря и приводного слоя атмосферы 36
2.4. Параметризация вертикального турбулентного обмена 42
Глава 3. Постановка численных экспериментов 51
3.1. Начальные условия 52
3.2. Пространственная сетка. Топография дна 52
3.3. Граничные условия 53
3.4. Пространственно-временные масштабы описываемых процессов 60
Глава 4 . Крупномасштабные характеристики решения 63
4.1. Термохалинные поля 63
4.2. Годовой ход уровня моря 65
4.3. Циркуляция бассейнового и суббассейнового масштаба 66
Глава 5. Мезомасштабная изменчивость 84
5.1. Виды мезомасштабной изменчивости 84
5.2. Высокочастотная изменчивость по данным наблюдений 85
5.3. Спектральные характеристики модельного решения 88
5.4. Инерционно-гравитационные волны 90
5.5. Сейшевые колебания 106
5.6. Суточные колебания 109
Глава 6. Размерностный анализ модельного решения 111
Заключение 120
Литература 124
- Разностная аппроксимация уравнений и граничных условий
- Граничные условия
- Термохалинные поля
- Высокочастотная изменчивость по данным наблюдений
Введение к работе
Изучение океанов и морей позволяет выявить закономерности изменчивости гидродинамических характеристик в широком пространственно-временном диапазоне. Определенная часть спектра пространственно-временной изменчивости состояния моря может быть непосредственно измерена. Например, достаточно полную информацию о синоптической динамике в настоящее время могут дать спутниковые измерения, хотя эти данные ограничены поверхностными характеристиками. Тем не менее, количественное описание циркуляции вод моря, понимание его динамики и взаимозависимости происходящих в море процессов не могут быть получены только с помощью измерительных средств. Численные модели необходимы как средство для интерпретации данных наблюдений и как средство для понимания динамики моря.
Сказанное относится и к изучению циркуляции вод Каспийского моря. Накопленные за десятилетия измерений океанографические данные позволяют составить представление о динамике вод моря и его термохалинной структуре в масштабах всего бассейна. Описание процессов синоптического и меньших масштабов как правило основывается на результатах локальных измерений. Из имеющихся данных о состоянии моря выделим следующие. Существуют статистически обеспеченные среднемесячные данные о термохалинной структуре верхнего 100-метрового слоя моря (Косарев, Тужилкин, 1995), данные о спектральных характеристиках течений (Бондаренко, 1993; Терзиев и др., 1992). Наиболее достоверными, по-видимому, являются данные о суточной, внутригодовой и межгодовой изменчивости уровня моря. По данным измерений биохимических параметров были обнаружены вихревые образования вдоль южного берега Каспия (Катунин, Сапожников, 1997). Космические снимки температуры поверхности моря позволили более детально увидеть известные из данных натурных измерений пространственно-временные характеристики апвеллинга вдоль восточного берега Среднего Каспия (Sur et al, 2000). Наряду с накопленными данными традиционных океанографических наблюдений и их интерпретацией (Kosarev, Yablonskaya, 1994), исследования с применением моделей гидродинамики моря, проведенные в последние несколько лет, позволили лучше понять характер крупномасштабной циркуляции вод моря (Ибраев и др., 2002; Ибраев и др., 2001; Тужилкин и др., 1997; Trukhchev et al, 1995).
Особенности гидрофизического режима Каспийского моря и его изменчивости необходимо учитывать при разработке модели термогидродинамики моря. Во-первых, Каспийское море замкнутое, не имеет связи с Мировым океаном. Поэтому одной из наиболее чувствительных к изменчивости внешнего воздействия характеристик является водный баланс. Во-вторых, в Каспийском море, в отличие от окраинных морей, вертикальная стратификация вод определяется прежде всего термическими факторами. Как следствие - повышенная чувствительность вертикальной стратификации и течений в верхнем (-100 м) слое моря (выше и ниже пикноклина) к изменчивости (суточного и синоптического масштаба) атмосферного воздействия. В-третьих, Каспийское море - мелкое, две трети площади моря - это шельф с глубинами менее 100 м. Вышесказанное определяет основные требования к разрабатываемой модели динамики моря, которые заключаются в следующем: модель должна описывать характерные динамические особенности моря (струйные береговые течения, вихревую динамику, апвеллинг, вертикальную стратификацию), изменчивость уровня (а следовательно и массы воды) моря и работать в режиме синоптического атмосферного форсинга.
Модельные исследования циркуляции вод и термодинамики Каспийского моря до сих пор относительно немногочисленны. Геометрические характеристики мелководного Северного Каспия и глубоководных Среднего и Южного Каспия обусловили требования к моделям циркуляции для этих подбассейнов. Большинство существующих моделей можно разделить на два класса, а именно модели Северного Каспия, основанные на уравнениях мелкой воды, и трехмерные модели всего моря.
Первой задачей по моделированию всего Каспийского моря была задача диагноза течений по данным гидрологической съемки. В 1976 году Саркисяном и др. (Саркисян и др., 1976) для Каспийского моря была впервые использована диагностическая модель динамики бароклинного океана. В ней была оценена циркуляция вод Каспийского моря в зависимости от типа ветра и сезона, выявлено значительное влияние ветра на формирование поля течений. Было получено, что летом влияние термохалинности проявляется больше, чем зимой.
Сезонная изменчивость климатической циркуляции моря в рамках адаптационного подхода реконструкции течений исследовалась в работах Тужилкин и др., 1997, Trukhchev et al, 1995. Результаты расчетов демонстрируют преемственность течений в Среднем и Южном Каспии от сезона к сезону. В Среднем Каспии постоянно существует система из циклонического вихря в северо-западной части и антициклонического вихря в юго-восточной. Южный Каспий охвачен антициклоническим вихрем в северо-западной части бассейна и Гаы-т ициклоническим вихрем в юго-восточной части. Сезонная изменчивость циркуляции, по результатам этих работ, заключается в пространственной изменчивости и интенсивности указанных вихрей, составляющих основу циркуляции. Однако реалистичность воспроизведения климатической циркуляции вод моря была ограничена недостаточным пространственно-временным разрешением гидрологических,данных.
В работе Ибраев и др., 1997 на основе выполненных автором ранних диагностических расчетов была решена задача распространения пассивной примеси, поступающей в море со стоком р. Волга.
Диагноз синоптических течений моря с использованием данных гидрологической съемки, выполненной в августе 1986 г., был сделан в работе Ахвердиев, Демин, 1989. Результаты численных расчетов показали наличие сложного характера циркуляции в Южном Каспии (наличие двух антициклонических и одного циклонического круговоротов), выявили апвеллинг у восточного берега Среднего Каспия, дали возможность описать обмен вод между Средним и Южным Каспием. Было показано, что центр циклонического круговорота в Среднем Каспии расположен над самым глубоким местом этого района.
До конца 90-х годов единственной моделью всего Каспийского моря, в которой воспроизводились нестационарные процессы динамики течений, была модель из работы Бадалов, Ржеплинский, 1989. Эта модель относится к классу двухслойных моделей. В модели предполагается, что в море существуют два слоя, динамика которых различается пространственно-временными масштабами. Динамика течений верхнего слоя моря (0-30 м в летний сезон и 0-100 м в зимний сезон) описывается трехмерными нестационарными нелинейными уравнениями гидродинамики моря. Динамика течений нижнего слоя описывается диагностической моделью Саркисяна (Саркисян, 1977). В работе моделируется генерация циркуляции вод Каспийского моря под воздействием характерных ветровых ситуаций. Показано, что Южный Каспий быстрее приспосабливается к новой ветровой ситуации, чем Средний Каспий, так как над Южным Каспием расположена зона слабых ветров.
Настоящая работа является продолжением исследований изменчивости термогидродинамических процессов Каспийского моря с применением Модели гидродинамики внутреннего моря (МГВМ) (Ibrayev, 2001). Исследования крупномасштабных характеристик внутригодовои изменчивости гидродинамики моря с применением МГВМ (Ибраев, 2002; Ибраев и др., 2002) показали, что: циркуляция поверхностных вод Среднего и Южного Каспия не является циклонической в течение года; течения вдоль восточного берега Среднего Каспия носят двухслойный характер. В работе Ибраев, 2002 также была исследована сезонная изменчивость уровня моря и ее зависимость от параметров приводного слоя атмосферы и температуры поверхности моря (ТПМ), определяемой как взаимодействием атмосферы и моря, так и трехмерными термогидродинамическими процессами. Ряд процессов, связанных с мезомасштабной изменчивостью характеристик моря и ее ролью в циркуляции вод, не мог быть описан в рамках разработанной ранее модели, в частности из-за недостаточного пространственного разрешения и использования среднемесячных атмосферных условий. В настоящей работе предпринята попытка расширить спектр описываемых в модели процессов. В частности, представлены результаты моделирования гидродинамики моря с более детальным описанием топографических вихрей и фронтальных зон. Использование более детальной пространственной сетки по вертикали позволило лучше описать термогидродинамические процессы верхнего слоя моря. Особенно важно то, что в настоящей работе на верхней границе задаются синоптические атмосферные условия, что позволит воспроизвести термогидродинамические свойства моря в режиме, более близком к реальному.
Применение синоптического атмосферного форсинга приводит к генерации в модельном решении процессов, которые при осредненном за месяц форсинге не проявляются. Временная дисперсия атмосферного воздействия играет важную роль в определении потоков импульса и тепла, а следовательно - в формировании верхнего перемешанного слоя и ТПМ, которая, в свою очередь, оказывает влияние на испарение с поверхности моря и на динамику уровня моря. Влияние временного осреднения атмосферного форсинга на воспроизведение в модели гидродинамических процессов рассматривается в целом ряде работ. В частности, одномерная задача об интенсификации турбулентности в верхнем слое моря во время шквального ветра и о распаде турбулентности после выключения ветра детально рассмотрена в работе Richardson et al, 1999.
При синоптическом атмосферном форсинге в модели генерируются движения инерционного временного масштаба (ДИВМ), которые вносят значительный вклад в изменчивость решения. Под ДИВМ в настоящей работе подразумеваются гироскопические волны, инерционно-гравитационные волны (ИГВ), ИГВ при учете краевых условий (волны Пуанкаре), инерционные колебания (LeBlond, Mysak, 1978). В моделях гидродинамики морей ДИВМ часто не учитываются из-за использования гладко меняющихся атмосферных условий, либо из-за фильтрации их численной схемой. Влияние инерционных колебаний на вертикальный обмен исследовалось в работе D Asaro, 1985, где рассмотрены различные трехмерные процессы, обусловленные штормовыми и послештормовыми явлениями. Там же можно найти краткий обзор публикаций на эту тему. Как показывает анализ результатов представленной в данной работе модели, при явном разрешении ДИВМ они оказывают влияние на крупномасштабные характеристики решения, в частности - на вертикальное перемешивание, а следовательно - на ТПМ, испарение и ход уровня моря.
В настоящей работе проанализированы среднемесячные характеристики циркуляции вод, моментальные характеристики вихревых структур, мезомасштабные термогидродинамические процессы моря, изменчивость уровня моря, полученные по расчетам с помощью модифицированной МГВМ с применением синоптического атмосферного воздействия.
Целями настоящей работы являются:
1) усовершенствование МГВМ для исследования мезомасштабной изменчивости гидродинамики Каспийского моря;
2) исследование крупномасштабной циркуляции, термохалинной структуры вод и годового хода уровня Каспийского моря на основе результатов расчетов по МГВМ;
3) анализ мезомасштабных гидродинамических процессов Каспийского моря по результатам численных экспериментов и их верификация по данным наблюдений.
Работа состоит из введения, шести глав, заключения и списка литературы. Работа содержит 37 рисунков и 9 таблиц.
В первой главе приводится краткое описание физико-географических характеристик Каспийского моря и региональных климатических особенностей. Приведены морфометрические данные, важные с точки зрения описания циркуляционных и термодинамических особенностей моря. Дается краткая характеристика ключевых компонент атмосферного воздействия, в том числе ветра, солнечного излучения, осадков. Приведены количественные характеристики составляющих водного баланса (речного стока, осадков, испарения), ввиду особой важности воспроизведения годового хода уровня моря для верификации модельного результата.
Во второй главе описана МГВМ и две ее подмодели, особенно важные для целей данной работы. Модель гидродинамики внутреннего моря является совместной моделью термогидродинамики моря, взаимодействия атмосферы и моря, термодинамики морского льда. Модель основана на трехмерных полных уравнениях геофизической гидродинамики (Саркисян, 1977). Приняты приближения гидростатики, несжимаемости морской воды и приближение Буссинеска. Подвижная верхняя граница моря описывается уравнением свободной поверхности в квазилинейном приближении, позволяющем описывать как высокочастотные колебания, так и изменение среднего уровня моря вследствие ненулевого водного баланса. Математическая формулировка уравнений модели и граничных условий (раздел 2.1), их конечно-разностная аппроксимация (раздел 2.2) приведены в сжатом виде, поскольку полное описание МГВМ можно найти в работах Ибраев, 2002 и Ibrayev, 2001. Среди физических, математических и технологических усовершенствований МГВМ по сравнению с предшествующей версией модели отметим следующие. Горизонтальная вязкость описывается бигармоническим оператором; адвективный перенос тепла и соли аппроксимируется так называемой схемой коррекции потоков (flux-corrected transport scheme), которая квазимонотонна и консервативна. В разделе 2.3 рассматривается подмодель расчета потоков импульса, тепла и влаги на границе моря и приводного слоя атмосферы, основанная на теории подобия Монина-Обухова. Эта подмодель была реализована и в предыдущей версии МГВМ, но подверглась существенной переработке с целью ускорения входящих в состав ее численной реализации итерационных процедур. Входными данными подмодели являются скорость ветра, температура и влажность воздуха, а также рассчитываемая температура поверхности моря. Важность применения рассчитываемых потоков свойств на границе атмосферы и моря, против задания потоков в виде граничных условий напрямую, обусловлена необходимостью учета активной реакции обменных процессов на реализуемую в модели ТПМ.
В разделе 2.4 описана параметризация вертикальной турбулентной вязкости и диффузии, основанная на схеме Мэллора-Ямады уровня 2.5 (Mellor, Yamada, 1982). В предыдущих версиях МГВМ параметризация вертикального обмена основывалась на более простых полуэмпирических закономерностях (Pacanowski, Philander, 1981), в которых коэффициенты вязкости и диффузии зависят только от числа Ричардсона. Формулы, определяющие зависимость коэффициентов от числа Ричардсона, содержат набор констант, подлежащих настройке для конкретной задачи. В этом состоит и преимущество простой параметризации вертикального обмена, т.к. есть возможность ее тонкой настройки, и недостаток, поскольку снижается универсальность модели. С другой стороны, безразмерные константы, входящие в уравнения схемы Мэллора-Ямады, предполагаются не зависящими от исследуемого объекта и модели. Схема строится путем замыкания уравнений Рейнольдса для вторых моментов турбулентных пульсаций, масштабирования слагаемых уравнений по степени отклонения от локальной изотропии и введения ряда последовательных упрощающих предположений, каждое из которых связывается с определенным уровнем в иерархии семейства схем параметризации турбулентности. Схема уровня 2.5 в настоящее время находит наибольшее применение в геофизическом моделировании, поскольку является отражением компромисса между физической релевантностью и вычислительной трудоемкостью. В настоящей диссертационной работе схема Мэллора-Ямады уровня 2.5 реализована таким образом, что она сводится к решению двух эволюционных уравнений относительно величин: турбулентной кинетической энергии (ТКЭ) и произведения масштаба длины на ТКЭ («у2 и q2l, соответственно). Коэффициенты вертикальной турбулентной вязкости и диффузии вычисляются с помощью соотношений Прандтля KM=lqSM и K„=lqSH, где SM, SH - функции устойчивости, зависящие от градиента плотности. Таким образом, система уравнений МГВМ дополняется двумя эволюционными уравнениями для q2 и q2l. Достоинством схемы Мэллора-Ямады является теоретически обоснованный подход к выводу соотношений, определяющих коэффициенты вертикального турбулентного обмена.
В третьей главе приводится порядок постановки численных экспериментов. Описывается модельная дискретизация расчетной области, источники батиметрических данных, способ задания начальных условий, приводятся характеристики атмосферного воздействия и стока рек.
Раздел 3.3 посвящен описанию граничных условий. Для задания граничных условий на верхней границе моря и в модели взаимодействия атмосферы и моря использовались данные реанализа атмосферной циркуляции ERA-15 ECMWF (Gibson et al, 1999), выполненные для периода 1979-1993 гг. Влияние притока речных вод в МГВМ учитывается посредством задания скорости течений, температуры и солености (равной нулю) речных вод на открытых участках боковой границы.
В разделе 3.4 рассматривается вопрос о классификации гидродинамических процессов Каспийского моря в соответствии с характерными для них пространственно-временными масштабами. Анализ изменчивости гидрофизических характеристик по результатам модельного расчета выявляет, что общей особенностью для всех районов и глубин является провал в спектрах всех характеристик в частотном диапазоне, соответствующем периодам от суток до 2-3 суток, т.е. приблизительно до нижней границы атмосферного синоптического периода. Такая особенность давно известна из данных наблюдений и характерна как для океана (напр., Каменкович и др., 1987), так и для внутренних и окраинных морей. Провал в спектре на этих периодах служит обоснованием для введения разделения процессов на крупномасштабные и мезомасштабные именно на этих периодах. Соответственным образом данная диссертационная работа разбивается на главы, посвященные процессам, которые удается с той или иной степенью обоснованности отнести к крупномасштабным или мезомасштабным.
В четвертой главе описываются крупномасштабные характеристики модельного решения. Рассматриваются термохалинные поля, годовой ход уровня моря, циркуляция бассейнового и суббассейнового масштаба. Специальное внимание уделяется анализу тех крупномасштабных явлений, для которых возможна верификация по известным натурным данным (Kosarev, Yablonskaya, 1994 и другие источники) и по результатам имеющихся модельных расчетов (Ибраев, 2002).
Анализ среднемесячных характеристик циркуляции вод (раздел 4.3) показывает, что наиболее стабильными в течение года являются течения в подповерхностных (ниже 20-50 м) слоях глубоководных частей моря, которые носят циклонический характер. Наиболее вероятным двигателем крупномасштабной циклонической завихренности подповерхностных течений в глубоководных районах представляется бароклинность моря за счет притока пресных вод на шельф и летнего прогрева вод на шельфе. В сторону усиления циклонической циркуляции подповерхностных слоев работает также локальная закачка завихренности, как отмечено в" работе Тужилкин и др., 1997, особенно в зимний период, когда над Средним и Южным Каспием конвергенция ветра и интенсивный вертикальный обмен стимулируют прибрежный даунвеллинг. Поверхностная система течений формируется под влиянием изменчивости ветра и вертикального обмена импульсом с подповерхностными течениями и характеризуется высокой степенью изменчивости в течение года. Циклонический характер крупномасштабной циркуляции вод в Среднем и Южном Каспии наблюдается только в декабре и январе. В летние месяцы преобладает экмановский характер течений в верхнем слое под действием преимущественного направления ветра на юг и юго-запад, соответственно перенос вод на поверхности направлен на юго-запад и запад.
В разделе 4.3. обсуждается также ряд явлений суббассейнового масштаба. Исследуется динамический апвеллинг вдоль восточного побережья Среднего Каспия в летние месяцы - яркий пример естественной геофизической лаборатории для изучения экмановского апвеллинга. Описаны струйные течения на западном шельфе и их меандрирование; двухслойная система течений на северо-восточном участке шельфа Среднего Каспия; вихреобразование на струйных течениях Южного Каспия. Изучены течения на взморье р. Куры, измерения которых дают противоречивую информацию. Рассмотрен вопрос о возникновении аномалии теплых поверхностных вод на востоке Среднего Каспия в зимние месяцы. По результатам настоящей работы подтверждается отмеченная в работе Ибраев и др.У 2002 особенность образования такой аномалии теплых вод, которая поддерживается благодаря постоянному подповерхностному северному течению, несмотря на то, что на поверхности начиная с января течения поворачивают на юг.
В пятой главе рассматривается изменчивость гидродинамических полей в мезомасштабном временном диапазоне. Классифицируются виды мезомасштабной изменчивости, характерные для внутренних и окраинных морей и для океана. Описываются имеющиеся данные высокочастотных измерений скорости течений. Натурные измерения выявляют наличие в Каспийском море колебательных и волновых процессов различной природы, характерных для внутренних морей и для океана: инерционных, сейшевых и суточных колебаний, небольших приливов, внутренних волн, береговых захваченных волн, и др. В разделе 5.3 приводятся спектральные характеристики модельного решения. В разделах 5.4, 5.5 и 5.6 на основе модельного результата исследуются конкретные виды мезомасштабной изменчивости: инерционно-гравитационные волны, сейшевые колебания, суточные колебания.
В шестой главе приводится размерностный анализ модельного решения. Изучается относительный вклад в уравнения модели каждого из слагаемых с целью выявления процессов, имеющих преимущественную значимость при определенных условиях.
В заключении кратко сформулированы основные результаты диссертационной работы.
Результаты настоящей работы докладывались на следующих конференциях, семинарах и научных школах: NATO Advanced Study Institute "The Ocean Carbon Cycle and Climate , Анкара, Турция, 2002; Конференция молодых ученых по проблемам Черного и Азовского морей, Севастополь, Украина, 2003; Nonlinear Processes in Geophysics, Таллинн, Эстония, 2003; Physical Processes in Natural Waters, Лунд, Щвеция, 2004; на семинарах Института вычислительной математики РАН, Москва, 2001-2003. Полностью работа докладывалась на семинаре Лаборатории геофизической гидродинамики Института океанологии им. П.П. Ширшова РАН, Москва, 2004.
Автор выражает глубокую благодарность своим научным руководителям, академику А.С. Саркисяну и д.ф.-м.н. Р.А. Ибраеву, без участия которых эта работа была бы невозможна, а равно и другие работы автора.
Разностная аппроксимация уравнений и граничных условий
Многосуточные измерения течений в Каспийском море проводятся начиная с)50-х годов. Достаточно много измерений проводилось в Северном Каспии, на Апшеронском пороге, вдоль восточного берега Среднего Каспия. В начальный период измерения течений был направлены прежде всего на выявление закономерностей зависимости течений от ветра. Были построены эмпирические закономерности для стационарных ветровых ситуаций в Северном Каспии, подробно описаны преимущественные направления течений и их структура в зависимости от силы, направления и продолжительности действия ветра в разных регионах моря (см., например, обзор в Бондаренко, 1993). Однако было очевидно, что изменчивость течений, даже в мелководном Северном Каспии, определяется широким спектром вынужденной реакции моря на изменчивость атмосферы и собственной изменчивостью моря.
Натурные измерения выявляют наличие в Каспийском море колебательных и волновых процессов различной природы, характерных для внутренних морей и для океана: инерционных, сейшевых и суточных колебаний, небольших приливов, внутренних волн, береговых захваченных волн, и др. (Терзиев и др., 1992; Бондаренко, 1993).
В работе Бондаренко, 1993 описаны результаты натурных измерений течений в Северном и Среднем Каспии. В числе основных энергонесущих частот оказались частоты, близкие к инерционным (период 17-19 часов). Измерения течений в этой работе производились главным образом на шельфе. На рис. 5.2.1 приведены примеры спектральной плотности компонент скорости
Обобщая имеющиеся данные о спектральных характеристиках скорости течения (Байдин, Косарев, 1986; Терзиев и др., 1992; Бондаренко, 1993), фиксируемую наблюдениями изменчивость течений можно описать следующим образом. Периодами, на которых существуют доминирующие максимумы спектральной плотности, являются суточный и инерционный. Суточные колебания скорости характерны для Северного Каспия и мелководья, инерционные (или околоинерционные) - для глубоководных районов. Оба этих периода выделяются в спектрах практически по всем наблюдениям. В области более высоких частот также бывают энергетические всплески, но как правило они не являются статистически обеспеченными. Например, в работе Бондаренко, 1993 указывается на существование максимумов спектральной плотности компонент скорости течения в некоторых пунктах измерений на периодах 13-14 ч. Эти максимумы расцениваются как недостоверные или малодостоверные, но включаются в анализ, поскольку постоянно возникают в спектрах течений (Бондаренко, 1993). Двигаясь по спектру в сторону низких частот обнаруживаем провал в области периодов 30-40 часов, и затем максимумы на длинных периодах, соответствующих, например, атмосферной синоптической изменчивости. Считается, что в низкочастотной области спектра ввиду недостаточной продолжительности измерений спектральная плотность искажается, и исследовать долгопериодную изменчивость целесообразно другими методами. Таким образом, общепризнано, что основными энергонесущими частотами (и не только в Каспийском море) являются околоинерционная и суточная, а основные виды волновых движений связаны с инерционно-гравитационными волнами, гироскопическими волнами, их предельным случаем - инерционными колебаниями, а также шельфовыми волнами и другими типами береговых захваченных волн.
Спектральные характеристики уровня моря как правило выявляют колебания с периодами 24 и 12 ч, соответствующие приливам, а также многочисленные сейшевые колебания, с периодами от 2.7-3 ч до 2-5 суток. Инерционный максимум в спектре уровня моря также фиксируется, но выражен он слабо, поскольку инерционно-гравитационные волны значительных изменений уровня моря не вызывают. В работе Байдин, Косарев, 1986 описаны короткопериодные внутренние волны в Каспийском море, которых в модельном решении не наблюдается. Зафиксированные по данным наблюдений короткопериодные внутренние волны имели амплитуду до 3-4 метров и период около 15 минут. Длина волны не была измерена, но если предположить, что она составляла несколько километров, то окажется, что фиксация такого явления находится на грани возможностей модели. Напомним, что в модели шаг сетки по вертикали 2 м (в верхних слоях), шаг по времени 6 мин, горизонтальное пространственное разрешение около 4 км. По-видимому, процессы, порождающие такие волны, в применяемом атмосферном форсинге отсутствуют.
Граничные условия
Для задания граничных условий на верхней границе моря и в модели взаимодействия атмосферы и моря использовались данные реанализа атмосферной циркуляции ERA-15 ECMWF (Gibson et al, 1999), выполненные для периода 1979-1993 гг. Анализ используемого в данной работе атмосферного форсинга ERA-15 для различных акваторий проводится в ряде работ (Nigam et al, 2000, Huang, Shukla, 1996, Schrum et al, 2001).
Можно сделать вывод о том, что данные реанализа хорошо соответствуют наблюдаемым полям атмосферных характеристик. Сезонная изменчивость среднемесячных ветров над Каспийским морем была детально рассмотрена в работе Ибраев и др., 2002. Модель сезонной изменчивости интегрировалась в режиме ежегодно повторяющихся условий на границах. Поскольку в 1979-1993 гг. уровень моря из-за дисбаланса внешних условий претерпевал значительные изменения, был выбран год с минимальным годовым изменением уровня за этот период, для того чтобы в модельное решение вносились минимальные возмущения в начале года (Ибраев и др., 2002). Анализ уровня моря, измеренного на уровенных постах в Махачкале, Форте-Шевченко, Красноводске и Баку, показал, что наименьшее изменение среднего уровня с начала января по конец декабря было в 1982 г.
Поэтому в модели задавались внешние условия, соответствующие 1982 г. В настоящей работе тем более является обоснованным выбор атмосферного форсинга 1982 г., поскольку есть возможность проведения совместного анализа новых и предшествующих (Ибраев, 2002) результатов моделирования крупномасштабных процессов с целью верификации модели мезомасштабных процессов.
В МГВМ в качестве граничных условий на верхней границе области задаются: поток импульса, потоки явного и скрытого тепла, потоки солнечной и длинноволновой радиации, потоки массы (осадки, испарение, намерзание и таяние льда). Поток импульса, потоки явного и скрытого тепла, а также испарение рассчитываются в подмодели взаимодействия атмосферы и моря, основанной на теории подобия Монина-Обухова (раздел 2.3 настоящей работы). Исходными данными для расчета потоков служат температура и влажность воздуха на уровне 2 м, температура поверхности моря, скорость ветра на уровне 10 м, атмосферное давление на уровне моря, а также реализуемые в модели скорость течения и температура верхнего слоя моря.
Для задания условий в приводном слое атмосферы использованы данные с дискретностью по времени 6 час и пространственным разрешением около 1.125 градуса. Влияние притока речных вод в МГВМ учитывается посредством задания скорости течений, температуры и солености (равной нулю) речных вод на открытых участках боковой границы. Всего в модели учитывается приток из четырех рек (Волга, Урал, Терек и Кура), причем приток из р. Волга описывается тремя основными рукавами. В 1980 году пролив, соединяющий море с заливом Кара-Богаз-Гол, был закрыт, поэтому в модели нет бокового стока вод. На твердых боковых границах области и на дне задается условие тепло- и солеизоляции, условие непротекания и свободного скольжения. атмосферная изменчивость проиллюстрирована на рис
Во временном ходе атмосферного давления усматривается изменчивость с синоптическим периодом (5-7 суток), амплитуда изменчивости позволяет допустить, что вклад горизонтального градиента атмосферного давления в уравнения движения сравним с баротропным и бароклинным градиентами давления. Детально этот вопрос в данной работе не изучался, но определенное подтверждение приводится в главе, посвященной размерностному анализу модельного решения. Считается, что изменения атмосферного давления являются одним из энергетических источников сейшевых колебаний, речь о которых также пойдет в дальнейшем.
Временной ход скорости ветра в целом, продолжительность и частота возникновения штормов, повторяемость направлений и силы ветра над разными акваториями Каспийского моря, по данным форсинга ERA-І 5, хорошо согласуются с данными наблюдений. Верификация проводилась на основе работы Терзиев и др., 1992. среднем больше, чем летом (рис. 3.3.2), но, как показывает модельный результат, всплески кинетической энергии верхнего слоя моря, как отклик на шторм, зимой и летом сравнимы, поскольку летом верхний слой изолирован от нижележащих устойчивой стратификацией и "сдувается" белее охотно без вязкостного взаимодействия с нижележащими слоями.
Термохалинные поля
Горизонтальное распределение температуры поверхности воды в зимнее время носит квазимеридиональный характер (рис. 4.1.1а).
Значительная часть Северного Каспия покрыта льдом с декабря по март, температура воды подо льдом немногим менее О С. К югу температура воды увеличивается и в Южном Каспии в яеваре составляет 13-14 С. Особенностью зимнего распределения поверхностной температуры воды является наличие языка теплых вод вдоль восточного берега Среднего Каспия.
Летом распределение ТПМ характеризуется высокой степенью прогрева мелководного Северного Каспия и Южного Каспия (рис. 4.1.1b). С июня по сентябрь вдоль всего восточного побережья Среднего Каспия наблюдается пониженная ТПМ ввиду апвеллинга подповерхностных вод. Разница температур между прибрежным мелководьем и глубоководными районами Среднего Каспия нередко превышает 10 С. Более подробно речь об апвеллинге пойдет в разделе 4.3.
Горизонтальное распределение солености верхних слоев моря (рис. 4.1.2) характеризуется в первую очередь наличием фронта в Северном Каспии, обусловленно граспресняющим воздействием р. Волга. Прослеживается также распространение пресных вод других рек, особенно стока р. Кура на фоне повышенной солености вод Южного Каспия. Интенсивное испарение в период с июня по октябрь с широкой шельфовой зоны Среднего и Южного Каспия вызывает осолонение поверхностных вод, которые опускаются по склону шельфа и переносятся подповерхностными циклоническими круговоротами. 4.2. Годовой ход уровня моря
Изучение межгодовых вариаций уровня Каспийского моря не является задачей, решаемой с помощью МГВМ, поскольку водный баланс моря определяется водным балансом обширной территории, включающей весь водосборный бассейн р. Волга, т.е. значительную часть европейской территории России. В ряде работ изучается связь между уровнем Каспийского моря и осцилляциями планетарного масштаба, например явлением Эль-Ниньо. Для описания межгодовой изменчивости уровня Каспийского моря МГВМ может использоваться как составная часть совместной климатической модели, причем ее роль будет сводиться к воспроизведению испарения с поверхности моря. В данной работе рассматривается внутригодовой ход уровня моря и его мезомасштабная изменчивость.
Основными факторами, определяющими пространственную изменчивость реализуемого в модели уровня Каспийского моря, являются приток плавучести в результате стока рек, баланс осадков и испарения и воздействие ветра. Все перечисленные факторы поддерживают асимметричное распределение солености, приводя к повышенной плотности воды вдоль восточного берега и пониженной плотности в северной и юго-западной частях моря, что хорошо видно на картах топографии уровня моря во все сезоны (рис. 4.2.1). Диапазон пространственной изменчивости среднемесячного уровня моря минимален в июле (9 см) и максимален в декабре (17 см). Эти величины в 2-3 раза меньше амплитуды сезонной изменчивости уровня моря.
Сезонная изменчивость осредненных за месяц характеристик циркуляции вод моря является следствием сезонной изменчивости основных внешних факторов: ветра, потока тепла между атмосферой и морем и притока пресных вод с впадающими в море реками. Выделим основные режимы, характерные для глубоководных частей моря.
Отрицательные потоки тепла между атмосферой и морем (с сентября по февраль) приводят к усилению вертикального перемешивания, увеличению толщины верхнего перемешанного слоя и к образованию однородных по вертикали течений. Зимнее перемешивание достигает глубин 100 м в Среднем и 50 м в Южном Каспии.
Положительные потоки тепла (с марта по август) приводят к росту плавучести верхних слоев моря, уменьшению толщины верхнего перемешанного слоя, обострению пикноклина (в условиях Каспийского моря это определяется прежде всего термическими факторами), а, следовательно, уменьшению обменных процессов через пикноклин, и к формированию двухслойности течений по вертикали. Течения в верхнем слое открытых частей моря, вдали от берегов, формируются под влиянием ветрового воздействия (экмановские течения) и бароклинных градиентов. Ниже верхнего -20 м (летом) - 50 м (зимой) экмановского слоя моря течения определяются в основном бароклинными факторами. Анализ горизонтальной циркуляции вод можно начать с глубоководных слоев, поскольку течения здесь более структурированы и стабильны во времени, чем в поверхностных слоях. По результатам моделирования течения в подповерхностных (ниже 20-50 м) слоях глубоководных частей моря в течение всего года носят циклонический характер (рис. 4.3.1). По-видимому, причиной этого является приток распресненных вод на шельф в течение всего года, более интенсивный прогрев шельфовых вод летом, более интенсивный обмен по вертикали на шельфе. Возможной причиной глубоководной циклонической завихренности течений является закачка завихренности в глубинные слои моря в результате опускания вод. В зимний период глубинные течения более интенсивны, что, видимо, связано с усилением циклонической завихренности в результате действия ветра.
Высокочастотная изменчивость по данным наблюдений
Многосуточные измерения течений в Каспийском море проводятся начиная с)50-х годов. Достаточно много измерений проводилось в Северном Каспии, на Апшеронском пороге, вдоль восточного берега Среднего Каспия. В начальный период измерения течений был направлены прежде всего на выявление закономерностей зависимости течений от ветра. Были построены эмпирические закономерности для стационарных ветровых ситуаций в Северном Каспии, подробно описаны преимущественные направления течений и их структура в зависимости от силы, направления и продолжительности действия ветра в разных регионах моря (см., например, обзор в Бондаренко, 1993). Однако было очевидно, что изменчивость течений, даже в мелководном Северном Каспии, определяется широким спектром вынужденной реакции моря на изменчивость атмосферы и собственной изменчивостью моря.
Натурные измерения выявляют наличие в Каспийском море колебательных и волновых процессов различной природы, характерных для внутренних морей и для океана: инерционных, сейшевых и суточных колебаний, небольших приливов, внутренних волн, береговых захваченных волн, и др. (Терзиев и др., 1992; Бондаренко, 1993).
В работе Бондаренко, 1993 описаны результаты натурных измерений течений в Северном и Среднем Каспии. В числе основных энергонесущих частот оказались частоты, близкие к инерционным (период 17-19 часов). Измерения течений в этой работе производились главным образом на шельфе. На рис. 5.2.1 приведены примеры спектральной плотности компонент скорости
Обобщая имеющиеся данные о спектральных характеристиках скорости течения (Байдин, Косарев, 1986; Терзиев и др., 1992; Бондаренко, 1993), фиксируемую наблюдениями изменчивость течений можно описать следующим образом. Периодами, на которых существуют доминирующие максимумы спектральной плотности, являются суточный и инерционный. Суточные колебания скорости характерны для Северного Каспия и мелководья, инерционные (или околоинерционные) - для глубоководных районов. Оба этих периода выделяются в спектрах практически по всем наблюдениям. В области более высоких частот также бывают энергетические всплески, но как правило они не являются статистически обеспеченными. Например, в работе Бондаренко, 1993 указывается на существование максимумов спектральной плотности компонент скорости течения в некоторых пунктах измерений на периодах 13-14 ч. Эти максимумы расцениваются как недостоверные или малодостоверные, но включаются в анализ, поскольку постоянно возникают в спектрах течений (Бондаренко, 1993). Двигаясь по спектру в сторону низких частот обнаруживаем провал в области периодов 30-40 часов, и затем максимумы на длинных периодах, соответствующих, например, атмосферной синоптической изменчивости. Считается, что в низкочастотной области спектра ввиду недостаточной продолжительности измерений спектральная плотность искажается, и исследовать долгопериодную изменчивость целесообразно другими методами. Таким образом, общепризнано, что основными энергонесущими частотами (и не только в Каспийском море) являются околоинерционная и суточная, а основные виды волновых движений связаны с инерционно-гравитационными волнами, гироскопическими волнами, их предельным случаем - инерционными колебаниями, а также шельфовыми волнами и другими типами береговых захваченных волн.
Спектральные характеристики уровня моря как правило выявляют колебания с периодами 24 и 12 ч, соответствующие приливам, а также многочисленные сейшевые колебания, с периодами от 2.7-3 ч до 2-5 суток. Инерционный максимум в спектре уровня моря также фиксируется, но выражен он слабо, поскольку инерционно-гравитационные волны значительных изменений уровня моря не вызывают. В работе Байдин, Косарев, 1986 описаны короткопериодные внутренние волны в Каспийском море, которых в модельном решении не наблюдается. Зафиксированные по данным наблюдений короткопериодные внутренние волны имели амплитуду до 3-4 метров и период около 15 минут. Длина волны не была измерена, но если предположить, что она составляла несколько километров, то окажется, что фиксация такого явления находится на грани возможностей модели. Напомним, что в модели шаг сетки по вертикали 2 м (в верхних слоях), шаг по времени 6 мин, горизонтальное пространственное разрешение около 4 км. По-видимому, процессы, порождающие такие волны, в применяемом атмосферном форсинге отсутствуют.
