Восстановление метеорологических полей по данным наблюдений тема диссертации и автореферата по ВАК РФ 05.13.18, доктор физико-математических наук Климова, Екатерина Георгиевна

• Климова, Екатерина Георгиевна
• доктор физико-математических наук
• 2005

Оглавление диссертации доктор физико-математических наук Климова, Екатерина Георгиевна

Введение

1 Схема численного анализа данных метеорологических наблюдений

1.1 Введение.

1.2 Метод интерполяции.

1.2.1 Система координат.

1.2.2 Обработка данных наблюдений.

1.2.3 Основные уравнения, используемые при интерполяции

1.2.4 Объединение близко расположенных данных в одно обобщенное наблюдение.

1.2.5 Контроль данных наблюдений с помощью статистической интерполяции

1.2.6 Уравнение анализа по данным, оставшимся после контроля

1.2.7 Моделирование статистических характеристик ошибок прогноза

1.2.8 Корреляционная функция ошибок прогнозы высоты и относительной топографии.

1.2.9 Корреляционная функция ошибок прогноза скорости ветра

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

1.2.11 Статистические характеристики ошибок наблюдений

1.3 Обработка данных наблюдений в схеме объективного анализа

1.3.1 Предварительные замечания.

1.3.2 Организация вычислений.

1.3.3 Контроль данных в схеме объективного анализа, основанной на методе трехмерной многоэлементной статистической интерполяции.

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

1.4.1 Предварительные замечания.

1.4.2 Подготовка данных наблюдений.

1.4.3 Преданализ.

1.4.4 Контроль данных с привлечением метода статистической интерполяции

1.4.5 Список исключенных данных.

1.4.6 Анализ.

1.4.7 Постанализ.

1.5 Численные эксперименты.

1.6 Основные результаты.

2 Система усвоения данных метеорологических наблюдений

2.1 Введение.

2.2 Численные эксперименты по усвоению данных с бароклинной квазигеострофической моделью атмосферы.

2.2.1 Квазигеострофическая бароклинная полусферная модель

2.2.2 Численные эксперименты по усвоению моделируемых данных с квазигеострофической моделью атмосферы

2.2.3 Численные эксперименты по усвоению реальных данных (данных ПГЭП) с квазигеострофической моделью атмосферы

2.2.4 Нормальные моды и численные эксперименты по усвоению данных.

2.3 Система усвоения метеорологических данных для Сибирского региона

2.3.1 Введение.

2.3.2 Анализ данных наблюдений.

2.3.3 Региональная модель краткосрочного прогноза погоды

2.3.4 Нелинейная инициализация для региональной модели атмосферы

2.4 Численные эксперименты с системой усвоения.

2.5 Основные результаты

3 Субоптимальные алгоритмы, основанные на фильтре Калма-на

3.1 Введение.

3.2 Фильтр Калмана.

3.3 Задача оптимальной фильтрации и субоптимальные алгоритмы

3.4 Упрощенные модели для расчета матриц ковариаций ошибок прогноза.••

3.5 Численные эксперименты по сравнению свойств упрощенных моделей

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

3.6.1 Методика расчета ковариаций, основанная на расщеплении по физическим процессам

3.6.2 Численные эксперименты.

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

3.7.1 Вывод уравнений для ковариаций ошибок прогноза

3.7.2 Расчет ковариаций ошибок прогноза в однородном изотропном случае

3.7.3 Численные эксперименты.

3.8 Основные результаты

4 Асимптотическое поведение алгоритма фильтра Калмана

4.1 Введение.

4.2 Задача оптимальной фильтрации (непрерывная постановка

4.3 Решение матричного дифференциального уравнения Риккати

4.4 Асимптотическое поведение фильтра Калмана.

4.5 Субоптимальные алгоритмы

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

4.7 Основные результаты.

5 Численные эксперименты по усвоению метеорологических данных

Оглавление б

5.1 Введение.

5.2 Упрощенные модели, используемые в численных экспериментах.

5.3 Численные эксперименты по усвоению моделируемых метеорологических данных с региональной бароклинной моделью атмосферы

5.4 Адаптивный алгоритм усвоения, основанный на фильтре Калмана

5.4.1 Введение.

5.4.2 Адаптивный алгоритм фильтра Калмана.

5.4.3 Численные эксперименты.

5.5 Применение методики усвоения данных в задачах моделирования циркуляции в океане и охраны окружающей среды.

5.5.1 Введение.

5.5.2 Численные эксперименты по усвоению данных о температуре поверхности океана в модели циркуляции Японского моря .•.

5.5.3 Применение субоптимального алгоритма, основанного на фильтре Калмана, в задаче усвоения данных наблюдений о пассивной примеси

5.6 Основные результаты

Введение диссертации (часть автореферата) на тему «Восстановление метеорологических полей по данным наблюдений»

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

В мире в настоящее время имеется обширная сеть наблюдательных станций. Наземные аэрологические станции и станции на кораблях передают данные радиозондирования. Эти данные передаются через каждые 6 часов в основные синоптические сроки: 00 час. СГВ, Об час. СГВ, 12 час. СГВ, 18 час. СГВ (СГВ - среднегринвичское время). Данные радиозондирования содержат информацию на изобарических поверхностях, до 10 мбар включительно. Синоптические данные передают наземную информацию. Эта информация также поступает в основные синоптические сроки наблюдений. Большой объем информации поступает со спутников. Спутниковые данные поступают практически непрерывно и поэтому их называют асиноптическими. Кроме основной сети наблюдений метеорологические данные передают также самолеты, дрейфующие буи. Для всех типов данных наблюдений характерно их неравномерное распределение по пространству и наличие случайных ошибок измерений.

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

Пусть вектор ж", заданный на сетке в момент времени tn (либо состоящий из коэффициентов разложения по некоторому базису), имеет размерность N. Задача прогноза может быть символически записана в виде задачи Коши: zf1 = А(х]), = д.

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

Определение 1.

Задачей численного (объективного) анализа будем называть задачу получения по данным наблюдений "наилучшего "в некотором смысле значения х) [Ц, 123, ИЗ, 112].

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

- климатические данные;

- "persistence"(анализ за предыдущий момент времени);

- прогноз по данным анализа за предыдущий момент времени.

Такую априорную информацию принято называть первым приближением.

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

- эмпирические методы;

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

- вариационный анализ (3DVAR - 3-Dimensional Variational);

- статистическая (оптимальная) интерполяция (ОИ).

Первые методы анализа данных были эмпирическими. В работах [88, 95] был предложен метод последовательных коррекций для проведения объективного анализа метеоданных. Математическая постановка задачи анализа была предложена в работе Л.С.Гандина [113] и в работе Сасаки [142]. Л.С.Гандиным был предложен метод оптимальной интерполяции, состоящий в минимизации средней по множеству реализаций ошибки анализа, при этом учитываются статистические характеристики ошибок наблюдений и некоторого поля, называемого "первым приближением". Сасаки был предложен вариационный подход к решению задачи анализа. В [124] было показано, что все известные методы численного (объективного) анализа метеорологических данных эквивалентны между собой. Что касается эмпирического метода последовательных коррекций, то в [92] показана сходимость итераций метода последовательных коррекций к методу оптимальной интерполяции при соответствующем выборе весовых функций.

С конца 70-х годов прошлого века во всех мировых прогностических цетрах мира лидирующее место занял предложенный Л.С.Гандиным метод оптимальной интерполяции [13, 15]. Лучшим по качеству был признан бокс-вариант метода трехмерной многоэлементной оптимальной интерполяции, разработанный в Европейском Центре Среднесрочных Прогнозов Погоды

ECMWF, Reading) А.Лоренком [123]. В Национальном Метеорологическом Центре США (NMC, в настоящее время NCEP - Национальный Центр Прогнозов Окружающей Среды) был разработан спектральный метод статистической интерполяции [129]. Оперативной схемой анализа в Гидрометцентре России до настоящего времени является схема, основанная на трехмерной многоэлементной оптимальной интерполяции [4, 6, 7]. В настоящее время в NCEP, ECMWF, Meteo-France (Франция) и метеослужбе UK (Великобритания) разработан и применяется оперативно вариационный метод анализа с привлечением статистической информации об ошибках данных наблюдений и первого приближения [104].

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

Задача анализа данных позволяет восстанавливать метеорологические поля по данным наблюдений за один момент времени. С 70-х годов 20 века появились наблюдательные средства, делающие замеры непрерывно (асиноп-тические данные наблюдений, например, спутниковые). Кроме того, данных наблюдений за конкретный срок может быть недостаточно для получения адекватной картины распределения метеополей. В связи со всем вышесказанным в настоящее время принято проводить пространственно-временную обработку данных с привлечением прогностических моделей для описания распределения метеополей в пространстве и во времени [112, 113].

Определение 2.

Задачу совместного учета данных наблюдений и прогностической модели для наиболее точного описания пространственно-временного распределения метеополей принято называть задачей усвоения данных (в зарубежной литературе "data assimilation").

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

В мировых прогностических центрах, таких как Европейский Центр Среднесрочных Прогнозов Погоды (ECMWF, Рэдинг) и Национальный Центр Прогнозов Окружающей Среды (NCEP, США) с 80-х годов 20 века стали действовать системы усвоения данных, которые представляли из себя 6-часовой цикл прогноз-анализ-инициализация [127, 129]. В этом цикле шаг "анализ"осуществляется в основные синоптические сроки каждые 6 часов по данным наблюдений из "окна"± 3 часа от основного синоптического срока наблюдения, для подавления ложного роста инерционно- гравитационных мод производится инициализация проанализированных данных [128], а затем прогноз на 6 часов. На этапе "анализ"при этом чаще всего используется метод трехмерной многоэлементной оптимальной интерполяции.

В последнее время такие системы усвоения, основанные на цикле прогноз-анализ-инициализация, усовершенствуются с привлечением общей вариационной постановки задачи и использованием сопряженной модели для минимизации соответствующего функционала, основываясь на известных результатах теории оптимального управления [62, 101, 114, 123]. В то же время в ряде работ рассматривается постановка задачи усвоения данных в терминах теории оценивания с привлечением алгоритма фильтра Калмана [107, 114]. Следует заметить, что задачи оптимального управления и оптимальной фильтрации связывает принцип двойственности и, кроме того, существуют алгоритмы совместного оценивания и управления процессом [10, 113, 107].

Вариационная постановка в задаче усвоения в настоящее время очень популярна в мире. Следует подчеркнуть, что впервые такая постановка задачи усвоения была предложена в 1976 г. в работе В.В.Пененко, Н.Н.Образцова [65]. Популярный в настоящее время алгоритм фильтра Калмана также исследовался в 70-х годах прошлого века нашими отечественными авторами (Д.М.Сонечкин, [77]). В работе Д.М. Сонечкина [78] показана эквивалентность постановок задачи оптимальной фильтрации и вариационного усвоения при ряде условий (в частности, при нормальном распределении ошибок наблюдений и прогноза). В работах В.В.Пененко [63, 135] показано, что алгоритм вариационного усвоения алгебраически эквивалентен процедуре обобщенного фильтра Калмана. В настоящее время теоретическое исследование проблем усвоения (в вариационной постановке) рассматривается в работах [145, 146].

Из зарубежных авторов особого внимания заслуживает А.Лоренк (A.Lorenc, Англия), посвятивший целый ряд работ теоретическим проблемам задачи усвоения данных. Так, в работе [125] обобщенный метод анализа был распространен автором на пространственно-временной случай. Полученные им формулы при определенных предположениях (некоррелируемость ошибок наблюдений по времени, гауссовские функции распределения для ошибок наблюдений и прогноза) могут быть сведены к вариационной постановке задачи усвоения, а также к динамико-стохастической. Поскольку алгоритм фильтра Калмана требует обращения матриц очень высокого порядка, авторы предложили свой подход к решению задачи - итерационный метод минимизации функционала, обобщающий метод последовательных коррекций [126].

Таким образом, в настоящее время существуют два подхода к задаче усвоения данных: вариационный(4DVAR - 4-Dimensional Variational) и дина-мико-стохастический (фильтр Калмана). Современные системы усвоения, основанные на вариационном подходе (4DVAR) более легко реализуемы на ЭВМ, но они

1). не учитывают изменчивость со временем ковариаций ошибок прогноза;

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

Как отмечается в [113, 125, 137], все задачи усвоения в принципе можно свести к задаче статистического оценивания. В этих работах задачу усвоения связывают с минимизацией квадратичного функционала. При этом, если ищется начальное значение, прогноз по которому доставляет минимум функционалу, имеем вариационное усвоение. Если осуществляется последовательное оценивание по времени за необходимый период, имеем фильтр Калмана. В случае линейной точной (без ошибок) модели оба алгоритма приведут к одному и тому же результату в конце периода усвоения. Если модель линейная, но имеет ошибки, тогда также можно установить соответствие между двумя подходами, если уравнение модели ввести в функционал в качестве слабого условия. Однако в этом случае 4DVAR будет так же сложно реализуем, как и фильтр Калмана.

В настоящее время задача усвоения данных исследуется отечественными авторами в основном в вариационной постановке [64, 132]. Следует отметить монографию [76], в которой рассматривается применимость алгоритма фильтра Калмана для проведения усвоения в океане. В Гидрометцентре России разработана коллективом автором система усвоения, основанная на цикле прогноз-анализ-инициализация [82, 106] и исследуются вопросы, связанные с ее развитием в виде применения обобщенной постановки задачи усвоения [149, 150].

В ECMWF метод 4DVAR используется в настоящее время оперативно [152], кроме того, рассматривается упрощенный вариант фильтра Калмана RRKF (Reduced Rank Kalman Filter), состоящий в аппроксимации ковариаций ошибок прогноза в 4DVAR с помощью сингулярных векторов [104]. В NCEP оперативно работает система усвоения типа 3DVAR со статистическим методом анализа PSAS (Physical-Space Statistical Analysis System) [151].

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

Система усвоения метеорологической службы Франции (Meteo-France) представляет собой реализацию глобального многоэлементного 4-мерного вариационного анализа (4DVAR), аналогичную системе усвоения ECMWF. В метеорологической службе Германии (NMC Offenbach) схема усвоения представляет собой 3-часовой цикл прогноз-анализ-инициализация со схемой анализа, основанной на методе 3-мерной многоэлементной оптимальной интерполяции. В метеорологической службе Канады в настоящее время работает глобальная система усвоения, представляющая собой 6-часовой цикл прогноз-анализ-инициализация с методом анализа 3DVAR. При этом в Канаде действует система прогнозирования ансамблей (16 member Ensemble Prediction System - EPS). В этой системе планируется заменить метод анализа, основанный на оптимальной интерполяции, на ансамблевый фильтр Калмана.

Данная работа посвящена во многом исследованию применимости алгоритма фильтра Калмана к задаче усвоения данных метеорологических наблюдений. Алгоритм фильтра Калмана был предложен Р.Е.Калманом в 1960 г. для дискретного случая и Р.Е.Калманом и Р.С.Бьюси в 1961 г. для непрерывного по времени случая для линейной динамической системы со случайным шумом. Разработаны также обобщения фильтра Калмана на нелинейный случай [107].

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

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

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

Большой обзор по применимости алгоритма фильтра Калмана в задачах усвоения данных метеорологических и океанических наблюдений приведен в - работах [112, 113].

При применении алгоритма фильтра Калмана к задаче усвоения метеорологических данных основными проблемами являются:

- высокий порядок матрицы ковариаций ошибок прогноза;

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

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

В настоящее время существует несколько подходов к преодолению проблемы большой размерности ковариационной матрицы ошибок прогноза. Один из них состоит в аналитическом описании ковариаций. Так, в работе [94] рассматриваются способы получения дифференциальных уравнений для локальных ковариаций ошибок прогноза. Второй подход состоит в использовании для расчета ковариаций упрощенной модели (меньшей размерности). Такой алгоритм предложен в [99] и называется субоптимальным алгоритмом, основанным на фильтре Калмана. Варианты субоптимальных алгоритмов рассмотрены в [99, 148]. Третий подход, приобретающий все большую популярность, состоит в применении метода Монте-Карло для расчета поведения ковариаций ошибок прогноза со временем. Это так называемый ансамблевый алгоритм фильтра Калмана. Такой подход предложен в [111]. Ансамблевый подход позволяет рассчитывать ковариации ошибки прогноза для нелинейных моделей. Однако он тоже содержит ряд сложностей. В частности, размерность ансамбля прогнозов для получения хорошо обусловленной ковариационной матрицы должна быть порядка размерности этой матрицы. Отметим здесь работу [132], в которой по существу также предлагается вариант ансамблевого фильтра Калмана. В работе [121] рассматривается комбинированный алгоритм, в котором используется метод ансамблей и при этом рассматривается оператор меньшей размерности (разложение по естественно-ортогональному базису).

Современные схемы, основанные на вариационном подходе (4DVAR), в настоящее время модифицируются, при этом в качестве весовых матриц используются ковариации ошибок прогноза, меняющиеся со временем. Для вычисления ковариационных матриц в этих схемах используется либо ансамб-левыц подход [117], либо разложение по естественно - ортогональному базису сингулярные векторы) [104].

Для успешной реализации алгоритма фильтра Калмана требуется задание матрицы ковариаций "шумов"прогностической модели. Точное значение элементов этой матрицы неизвестно. Если считать ее нулевой, элементы матрицы ковариаций ошибок прогноза в алгоритме фильтра Калмана убывают со временем, что приводит к эффекту, называемому "расходимостью "фильтра Калмана [107, 10]. В этом случае данные наблюдений входят в процедуру анализа со все меньшими весовыми коэффициентами и фактически перестают "усваиваться". Матрица "шумов"прогностической модели может быть оценена с помощью адаптивного алгоритма по данным наблюдений. Все предлагаемые к настоящему времени адаптивные алгоритмы основаны на идее, предложенной в [108]. Алгоритм основан на использовании свойств вектора "невязок"(разности между наблюденным значением и прогнозом).

Адаптивные алгоритмы усвоения, основанные на использовании вектора "невязок"и его статистических свойств предлагаются во многих работах [90, 100, 130]. В [77] предлагается полуэмпирический алгортм усвоения, также основанный на использовании векторов "невязок". Следует отметить работу [149], посвященную проблеме оценки ошибок модели и исследованию природы этих ошибок.

В [100] был предложен алгоритм адаптивного фильтра Калмана для усвоения метеорологических данных, в котором для "шумов модели"и ошибок наблюдений рассматривается дополнительный алгоритм оценивания. В работе [130] предложен адаптивный алгоритм, основанный на ансамблевом подходе. В [90] рассмотрен ряд адаптивных алгоритмов и исследована их применимость для задачи усвоения в океане.

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

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

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

- разработка системы усвоения метеорологических данных для Сибирского региона на основе цикла прогноз-анализ-инициализация, исследование ее свойств на реальных данных (данные OA ГМЦ России); создание технологической линии, реализующей эту систему;

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

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

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

Научная новизна результатов диссертационной работы.

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

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

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

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

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

Задача улучшения качества анализа метеорологических данных имеет большое практическое значение. Использование более качественных схем анализа и систем усвоения позволит получать более точное описание пространственно-временного распределения метеополей.

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

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

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

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

Материалы, вошедшие в диссертационную работу, обсуждались на семинарах ИВТ СО РАН, ИВМиМГ СО РАН, Гидрометцентра России, Главной Геофизической обсерватории. Они докладывались на Междуведомственном семинаре по научным результатам ПГЭП (Москва, 1983), Всесоюзной конференции "Статистическая интерпретация гидродинамических прогнозов с целью прогноза элементов и явлений погоды"(Одесса, 1991), Всесоюзном совещании по природным и антропогенным катастрофам (Томск, 1991), Международной конференции АМСА-95 (Новосибирск, 1995), Втором Сибирском конгрессе по прикладной и индустриальной математике ИНПРИМ-96 (Новосибирск, 1996), Международной конференции "Математические модели и методы их использования"(Красноярск. 1997), Третьем Сибирском конгрессе по прикладной и индустриальной математике ИНПРИМ-98 (Новосибирск, 1998), Первой Международной конференции по Реанализу (Silver Spring, Maryland, USA, 1997), Третьем Сибирском конгрессе по прикладной и индустриальной математике ИНПРИМ-98 (Новосибирск, 1998), Третьем международном симпозиуме по усвоению данных наблюдений в метеорологии и океанологии (Quebec City, Canada, 1999), Четвертом Сибирском конгрессе по прикладной и индустриальной математике ИНПРИМ-2000 (Новосибирск, 2000), Международной конференции "Моделирование, базы данных и информационные системы для атмосферных HayK"MODAS (Иркутск, 2001), VIII Международном Симпозиуме "Оптика атмосферы и океана. Физика атмосферы" (Иркутск, 2001), Рабочем совещании NATO "Air pollution processes in regional scale"(Kallithea, Halkidiki, Greece, 2003), Международной конференции по измерениям, моделированию и информационным системам для изучения окружающей среды ENVIROMIS-2002 (Томск, 2002), Международной конференции "Вычислительно-информационные технологии для наук об окружающей среде"Cites-2003 (Томск, 2003), Международной конференции "Математические методы в геофизике"(Новосибирск, 2003), Международной конференции по измерениям, моделированию и информационным системам для изучения окружающей среды ENVIROMIS -2004 (Томск, 2004).

Основные результаты опубликованы в 24 статьях: 13 статьях в ведущих отечественных журналах [21] - [26], [28], [29, 31, 34, 35, 39, 40], 4 докладах в трудах международного симпозиума и международного рабочего совещания [27, 45, 119, 120], 7 статьях в трудах институтов Академии Наук и Госкомгид-ромета [32] - [33], [36] - [38], [71, 72].

Диссертация состоит из введения, пяти глав, заключения и списка литературы, содержит 29 рисунков и 13 таблиц. Каждая глава разбита на разделы.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Климова, Екатерина Георгиевна

5.6 Основные результаты

Ниже перечислены основные результаты и выводы, полученные в V главе.

• Предложен алгоритм усвоения данных метеорологических наблюдений для бароклинной региональной модели атмосферы. Усвоение осуществляется с помощью субоптимального алгоритма, основанного на фильтре Калмана, предложенного в [21, 23, 24, 25]. Показано, что усвоение данных наблюдений позволяет существенно улучшить качество прогноза, при этом предлагаемый алгоритм усвоения имеет преимущество по сравнению с традиционным усвоением.

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

• Предлагаемый алгоритм адаптивного фильтра Калмана

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

- позволяет компенсировать отсутствие информации о матрице ковариаций "шумов модели "Qk-i при вычислении матрицы ковариаций ошибок прогноза Р[\

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

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

Заключение

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

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

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

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

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

2. Проведено аналитическое исследование свойств алгоритма усвоения метеорологических данных, основанного на фильтре Калмана.

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

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

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

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

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

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

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

Список литературы диссертационного исследования доктор физико-математических наук Климова, Екатерина Георгиевна, 2005 год

1. Алдухов О.А. Комплексный контроль аэрологических данных ПГЭП уровня 1.. - Метеорология и гидрология, 1983, N12, с.94-103.

2. Анцыпович В.А. Комплексный контроль геопотенциала и температуры на стандартных изобарических поверхностях. Труды Гидрометцентра СССР, 1980, вып.217, с.67-82.

3. Аоки М. Оптимизация стохастических систем. Москва, Наука, 1971.

4. Багров А.Н. Оперативная схема объективного анализа аэрологической информации для Северного полушария. труды Гидрометцентра СССР, 1978, вып. 196, с.3-30.

5. Багров А.Н., Локтионова Е.А. Восполнение недостающих данных геопотенциала и температуры на уровнях 300 и 10 мб. Труды Гидрометео-центра СССР, 1978, вып.212, с.42-46.

6. Багров А.Н., Локтионова Е.А. Оперативная схема объективного анализа барической топографии по данным аэрологического и спутникового зондирования атмосферы. Метеорология и гидрология, 1982, N9, с.25-33.

7. Багров А.Н., Гордин В.А., Цырульников М.Д. Оперативная схема объективного анализа в тропосфере и стратосфере. Метерология и гидрология, 1990, N8, с.37-45.

8. Балакришнан А.В. Теория фильтрации Калмана. Москва, Мир, 1988, 168 с.

9. Беллман Р. Введение в теорию матриц. Москва, Наука, 1969, 368 с.

10. Браммер К., Зиффлинг Г. Фильтр Калмана-Бьюси. Москва, Наука, 1982, 257 с.

11. Вейль И.Г., Кордзахия Г.И., Машкович С.А., Сонечкин Д.М. Численные эксперименты по четырехмерному анализу на основе динамико-статистического подхода. Метеорология и гидрология, 1975, №7, с.11-19.

12. Гандин J1.C. Объективный анализ метеорологических полей. J1., Гид-рометеоидат, 1963, 287 с.

13. Гандин JI.C. Четырехмерный анализ метеорологических полей. JI.: Гид-рометеоиздат, 1976. - 61 с.

14. Гордин В.А. Математические задачи гидродинамического прогноза погоды. Вычислительные аспекты. J1.: Гидрометеоиздат, 1987. - 264 с.

15. Дымников В.П., Контарев Г.Р. Объективный анализ поля геопотенциала с помощью собственных функций динамического оператора. В кн.: Статистические методы в метеорологии. 4.1, Новосибирск, ВЦ СО РАН СССР, 1969, с.123-138.

16. Дымников В.П. О связи естественных ортогональных составляющих полей метеоэлементов с собственными функциями динамических операторов. Известия АН СССР, сер.ФАО, 1988, том 24, №7, с.675-679.

17. Каленкович Е.Е., Пененко В.В. Численная схема прогноза полей метеоэлементов для Северного полушария. В кн.: Численные методы решения задач прогноза погоды и общей циркуляции атмосферы. Новосибирск, 1970, с.7-42.

18. Камке Э. Справочник по обыкновенным дифференциальным уравнениям. Москва, Наука, 1976, 576 с.

19. Ким В.Ф. Некоторые результаты численного эксперимента по квазигео-строфической модели прогноза геопотенциала на полусфере.- Изв. АН СССР, Физика атмосферы и океана, 1968, т.4, N2, с. 115-124.

20. Климова Е.Г. Методика усвоения данных метеонаблюдений на основе обобщенного субоптимального фильтра Калмана. Метеорология и гидрология, 1997, №11, с.55-65.

21. Климова Е.Г. Асимптотическое поведение схемы усвоения метеорологических данных, основанной на алгоритме фильтра Калмана. Метеорология и гидрология, 1999, №8, с.55-65.

22. Климова Е.Г. Упрощенные модели для расчета ковариационных матриц в алгоритме фильтра Калмана. Метеорология и гидрология, 2000, №6, с.18-30.

23. Климова Е.Г. Модель для расчета ковариаций однородных изотропных случайных полей ошибок прогноза. Метеорология и гидрология, 2001, №10, с.24-33.

24. Климова Е.Г. Модель для расчета ковариаций ошибок прогноза в алгоритме фильтра Калмана, основанная на полных уравнениях. Метеорология и гидрология, 2001, №11, с.11-21.

25. Климова Е.Г. Численные эксперименты по усвоению метеорологических данных с помощью субоптимального фильтра Калмана. Метеорология и гидрология, 2003, №10, с.54-67.

26. Климова Е.Г. Методика усвоения данных наблюдений на основе адаптивного субоптимального фильтра Калмана. Математические методы в геофизике. Труды Международной конференции, часть II, Новосибирск, 2003, с.400-404.

27. Климова Е.Г. Алгоритм усвоения данных наблюдений на основе адаптивного субоптимального фильтра Калмана. Метеорология и гидрология, 2005, №3, с.24-35.

28. Климова Е.Г., Воронина П.В. Модификация модельных корреляционных функций в схеме многоэлементного объективного анализа для Сибирского региона. Метеорология и гидрология, 1994, №1, с. 18-24.

29. Климова Е.Г., Ривин Г.С. О решении уравнения Булеева-Марчука. Известия АН СССР, Физика атмосферы и океана, 1978, т.15, N4, с.368-383.

30. Климова Е.Г., Ривин Г.С. О краткосрочном прогнозе погоды на полусфере по квазигеострофической модели. Труды ЗСРНИГМИ, 1979, вып.45, с.34-36.

31. Климова Е.Г., Ривин Г.С. Численные эксперименты по четырехмерному анализу с квазигеострофической моделью атмосферы. Труды ЗапСиб-НИИ, 1981, вып.53, с.53-59.

32. Климова Е.Г., Ривин Г.С. Уточнение первого приближения для численного анализа метеорологических полей. Метеорология и гидрология, 1985, №10, с.29-36.

33. Климова Е.Г., Ривин Г.С. Схема численного трехмерного многоэлементного анализа для Сибирского региона. Метеорология и гидрология, 1992, №3, с.16-23.

34. Климова Е.Г., Ривин Г.С. Численный анализ метеоданных. Изд. ИВТ СО РАН, Новосибирск. В кн.: "Вычислительные технологии", т.1, №3, 1992, с.235-240.

35. Климова Е.Г., Ривин Г.С. Восстановление полей метеоэлементов по данным наблюдений. Изд. ИВТ СО РАН, Новосибирск. В кн.: "Вычислительные технологии", т.2, №4, 1993, с.111-116.

36. Климова Е.Г., Ривин Г.С. Объективный анализ метеоданных на полусфере. Изд. ИВТ СО РАН, Новосибирск. В кн.: "Вычислительные технологии", т.2, №7, 1993, с.156-164.

37. Климова Е.Г., Ривин Г.С. Система усвоения метеоданных для Сибирского региона: численные эксперименты с реальными данными. Метеорология и гидрология, 1996, №12, с. 19-26.

38. Климова Е.Г., Ривин Г.С. Задача усвоения данных метеорологических и аэрозольных наблюдений. Химия в интересах устойчивого развития, 2002, №10, с.601-608.

39. Кондратьев К.Я. Влияние использования данных зондирования атмосферы со спутников на оправдываемость краткосрочного прогноза погоды. В кн.: Итоги науки и техники. Метеорология и климатология,т.5, М., ВИНИТИ, 1979.

40. Крупчатников В.Н., Крылова А.И. Численное моделирование распределения метана по данным наблюдений на поверхности Земли. Оптика атмосферы и океана, 2000, 13, №6-7, с.622-626.

41. Курбаткин Г.П., Зулунов С.М. Исследование проблемы геострофической адаптации метеорологических полей на сфере. Новосибирск, 1980. -34с. (Препринт /ВЦ СО АН СССР: 221).

42. Курбаткин Г.П., Зулунов С.М. Геострофическое согласование метеорологических полей с помощью собственных решений линеаризованной системы полных уравнений. Численные методы механики сплошной среды. Новосибирск, 1982. т.13, N2, с.82-97.

43. Ландау Л.Д., Лифшиц Е.М. Гидродинамика. Москва, Наука, 1986, 736 с.

44. Либерман Ю.М. Параметры статистического контроля вертикальных профилей геопотенциала и температуры. Труды ЗапСибНИГМИ, 1978, вып.296 с.62-67.

45. Марчук Г.И. Численные методы в прогнозе погоды, Гидрометеоиздат, Ленинград, 1967, 356 с.

46. Марчук Г.И. Математическое моделирование в проблеме окружающей среды. М.: Наука, 1982. - 319 с.

47. Марчук Г.И. Методы вычислительной математики. М.: Наука, 1989. -535 с.

48. Марчук Г.И., Курбаткин Г.П. и др. Оперативная квазигеострофическая схема краткосрочного прогноза погоды для пяти уровней атмосферы. -Изв. АН СССР, Физика атмосферы и океана, 1965, т.1, N1, с.119-135.

49. Марчук Г.И., Контарев Г.Р., Ривин Г.С. Краткосрочный прогноз погоды по полным уравнениям для ограниченной территории. Известия АН СССР, Физика атмосферы и океана, 1967, т.З, N11, с.1166-1178.

50. Машкович С.А. О четырехмерном анализе данных наблюдений. Метеорология и гидрология, 1971, N11, с.97-102.

51. Машкович С.А. О многоэлементном объективном анализе метеорологических элементов. Метеорология и гидрология, 1980, N5, с.5-14.

52. Машкович С.А. Многоэлементный трехмерный объективный анализ метеорологических величин. Метеорология и гидрология, 1988, N12, с. 1424.

53. Медведев С.Б. Инициализация региональной модели атмосферы с использованием вариационного метода. Метеорология и гидрология, 1998, №6, с.35-44.

54. Монин А.С. Прогноз погоды как задача физики. Москва, Наука, 1969, 184 с.

55. Монин А.С. Теоретические основы геофизической гидродинамики. Ленинград, Гидрометеоиздат, 1988, 424 с.

56. Монин А.С., Яглом A.M. Статистическая гидромеханика. Часть 2. Москва, Наука, 1967, 720 с.

57. Образцов Н.Н. Процедура быстрого преобразования Фурье, (инв. номер П-001147).- Инф. бюлл.ВНТИЦ Алгоритмы и программы, 1975, N2, с.12.

58. Окиншевич Р.И., Ривин Г.С., Уразалина З.К. Об оценке результатов численных экспериментов. Труды ЗСРНИГМИ, 1978, вып. 29, с.85-91.

59. Пененко В.В. Методы численного моделирования атмосферных процессов. Ленинград, Гидрометеоиздат, 1981, 351 с.

60. Пененко В.В. Системная организация математических моделей для задач физики атмосферы, океана и охраны окружающей среды. Новосибирск, 1985, 43 с. (Препринт/АН СССР, Сиб. отделение, ВЦ; 619).

61. Пененко В.В. Теоретические основы совместного использования данных наблюдений и моделей для исследования процессов гидротермодинамики и переноса примесей в атмосфере. Оптика атмосферы и океана, 1999, т.12, N5, с.485-462.

62. Пененко В.В., Образцов Н.Н. Вариационный метод согласования полей метеорологических элементов Метеорология и гидрология, 1976, N11, с.3-16.

63. Пененко В.В., Образцов Н.Н. Вариационно-разностный метод объективного анализа Метеорология и гидрология, 1978, N6, с. 15-25.

64. Пененко В.В., Ривин Г.С. Полная проблема собственных значений яко-биевой матрицы. В кн.: Вычислительные методы линейной алгебры. Новосибирск, ВЦ СО РАН СССР, 1969, с. 62-68.

65. Ривин Г.С. Численное моделирование фоновых атмосферных процессов и проблема переноса аэрозолей в Сибирском регионе. Оптика атмосферы, 19966 т.9, т, с.780-785.

66. Ривин Г.С., Бузова З.С., Смирнова А.И. Оперативная схема численного анализа для Сибирского региона. Метеорология и гидрология, 1990, N4, с.42-49.

67. Ривин Г.С., Климова Е.Г., Слуднов А.В. Изучение усвоения дополнительной метеорологической информации с помощью численных экспериментов типа близнецов и кузенов. В кн.: Математические модели атмосферных движений. Новосибирск, ВЦ СО РАН, 1980, с.100-113.

68. Ривин Г.С., Климова Е.Г., Медведев С.Б., Фомин В.М., Воронина П.В., Куликов А.И. Математическое моделирование процессов в атмосфере для проведения экологических экспертиз. В кн: Математические проблемы экологии. Новосибирск, ИМ СО РАН, 1994, с.90-95.

69. Ривин Г.С., Медведев С.Б. Гидродинамическая модель атмосферы для региона с применением метода расщепления. -Метеорология и гидрология, 1995, №5, с.13-22.

70. Ривин Г.С., Слуднов А.В. Модули доступа к данным DST-6 и ПГЭП на ЭВМ ЕС. Препринт ОВМ АН СССР, Москва, 1982, N29, 23 с.

71. Ривин Г.С., Сакулин В.Н., Слуднов А.В. Банк данных DST-6 и доступ к нему на ЭВМ серии ЕС. Новосибирск, 1980. - 26 с. (Препринт/ВЦ СО АН СССР: 240).

72. Саркисян А.С. Моделирование динамики океана. Санкт-Петербург, Гидрометеоиздат, 1991, 259 с.

73. Сонечкин Д.М. Обоснование четырехмерного (непрерывного) усвоения данных метеорологических наблюдений на основе динамико-стохастического подхода. Метеорология и гидрология, 1973, Ng4, с. 13-20.

74. Сонечкин Д.М. Динамико-стохастический подход к задаче объективного анализа данных разнородных метеорологических наблюдений. Труды Гидрометцентра СССР, 1976, вып. 181. с.54-76.

75. Толстых М.А. Полулагражева модель атмосферы с высоким простран-ственым разрешением для численного прогноза погоды. Метеорология и гидрология, 2001, №4, с.5-16.

76. Уилкинсон Д., Райнш К. Справочник алгоритмов на языке Алгол. Линейная алгебра. М., Машиностроение, 1976, 389 с.

77. Филатов А.Н., Шершков В.В. Асимптотические методы в атмосферных моделях. Ленинград, Гидрометеоиздат, 1988, 270 с.

78. Цырульников М.Д., Толстых М.А., Багров А.Н. Зарипов Р.Б. Развитие глобальной системы усвоения данных с переменным разрешением. Метеорология и гидрология, 2003, №4, с.5-24.

79. Шокин Ю.И., Ривин Г.С, Хакимзянов Г.С., Чубаров Л.Б. Вычислительный эксперимент как инструмент для исследования природных явлений.- Вычислительные технологии, 1992, т. 1, № 3, с. 12 33.

80. Bengtsson L. Four-dimensional asimilation of meteorological observations.-GARP Publ. Ser., 1975, N15, 76 p.

81. Bengtsson L., Kanamitsu M., Kallberg P., Uppala S. FGGE 4-dimensional data assimilation at ECMWF. Bui. of Amer.Met.Soc., 1982, v.63, N1, p.29-43.

82. Bergman K.H. Multivariate analysis of temperatures and winds using optimum interpolations. Monthly Weather Review, 1979, v. 107, N11, p. 1423-1444.

83. Bergman K.H. Role of observational errors in optimum interpolation analysis.- BAMS, 1978, v.59, N12, p.1603-1611.

84. Bergtorsson P., Doos B.R. Numerical weather map analysis. Tellus, 1955, 7, p.329-340.

85. Bermejo R., Staniforth A. The conversion of semi-lagrangian advection scheme to quasi-monotone scheme. Monthly Weather Review, 1992, v.120, p.2622-2632.

86. Blanchet I., Frankignoul C. A comparison of adaptive Kalman filters for a Tropical Pacific Ocean Model. Short-term dynamics of forecast error covariances. - Monthly Weather Review, 1997, v. 125.

87. Blumen W. Geostrophic ajustment. Rev. of Geoph. and Space Phys. - 1972. -Vol.10, /No2, p.485-528.

88. Bratseth A.M. Statistical interpolation by means of successive corrections. -Tell us, 1986, 38A, p.439-447.

89. Catz G.J. A method for solving a system of linear equations efficiently in order to optimize the analysis code: operational implementation. ECMWF Technical Memorandum, 1981, N33.

90. Cohn S.E. Short-term dynamics of forecast error covariances. In: Proceedings of ECMWF, Workshop "Variational assimilation, with special emphasis on three-dimensional aspects". - Nov. 9-12, 1992, p.157-170.

91. Cressman G.P. An operational objective analysis system. Monthly Weather Review, 1959, 87, p.367-374.

92. Courtier P., Talagrand 0. Variational assimilation of meteorological observations with the direct and adjoint shallow-water equation. Tellus, 1990, JVM2A, p.531-549.

93. Daley R., Menard R. Spectral characteristics of Kalman filter systems with application to the equatorial beta-plane. Proceedings of a Workshop

94. Variational assimilation, with special emphasis on three-dimensional aspects", ECMWF, June 1993, p.207-235.

95. Data asimilation methods. Seminar 1980. ECMWF, 1981, 454 p.

96. Dee D.P. Simplification of the Kalman filter for meteorological data assimilation. Quarterly Journal of the Royal Meteorological Society, v. 117, 1991. P.365-384.

97. Dee D.P., Cohn S.E., Dalcher A., Ghil M. An efficient algorithm for estimating noise covariances in distributed systems. IEEE Transactions on automatic control, 1985, v. AC-30, No.ll, p. 1057-1065.

98. Le Dimet F.-X., Talagrand O. Variational algorithms for analysis and assimilation of meteorological observations. Tellus, №38A, 1986. P. 97-110.

99. ECMWF data assimilation. Research manual, 1983, N1, 103 p.

100. Ehrendorfer M., Tribbia J. Optimal prediction of forecast error covariances throuh Singular Vectors. Journal of the Atmospheric Sciences, 1997, v. 54, p.286-313.

101. Fisher M., Andersson E. Development in 4D-Var and Kalman Filtering. -Technical Memorandum No.357, ECMWF, Reading, England, 2001, 36 p.

102. Forester C.K. Higher oder monotonic convective difference schemes. J. of Computational Physics, 1977, vol.23, p.1-22.

103. Jazwinski А.Н. Stochastic processes and filtering theory. Academic Press, New York, 1970, 377 p.

104. Jazwinski A.H. Adaptive filtering. Automatica, 1969, vol.5, p.475-485.

105. Heemink A.W., Segers A.J. Modeling and prediction of environmental data in space and time using Kalman filterin. Stochastic Environmental Research and Risk Assessment, 2003, 16, p.225-240.

106. Hollingsworth A. Objective analysis for numerical weather prediction. In: Short- and medium-range numerical weather prediction. Collection of papers presented at the WMO/IUGG NWP Symposium, Tokio, 4-8 August 1986. Tokyo, 1988, p.11-61.

107. Houtekamer P.L., Mitchell H.L. Data assimilation using an ensemble Kalman Filter Technique. Monthly Weather Review, 1998, v. 126, p.796-811.

108. Ghil M. Meteorological data assimilation for oceanographers. Part I: description and theoretical framework. Dynamics of Atmospheres and Oceans, 1989, v.13, N 3-4, p.171-218.

109. Ghil M., Malanotte-Rizzolli P. Data assimilation in meteorology and oceanography. In: Advances in Geophysics, v. 33, Academic Press, 1991, p.141-266.

110. Ghil M., Todling R. Tracking atmospheric instabilities with the Kalman filter. Part II: Two-layer results. Monthly Weather Review, v.122, 1994, p.2340-2352.

111. Kalnay E., Park S.K., Pu Z., Gao J. Application of quasi-inverse method to data assimilation. Monthly Weather Review, v. 128, 2000, p.864-878.

112. Kalnay Е. Atmospheric modelling, data assimilation and predictability. -Cambridge, University Press, 2003, 341 p.

113. Kalnay E., Hunt B.R., Kostelich E.J., Ott E., Patil D.J., Sauer Т., Szunyogh I., Yorke J.A., Zimin A.V. Four-dimensional ensemble Kalman filtering. -Tellus, 2004, 56A, pp.273-277.

114. Keppenne C.L., Rienecker M.R. Initial testing of a massive parallel ensemble Kalman filter with the Poseidon Isopycnal Ocean General Circulation Model.- Monthly Weather Review, v. 130, 2002, p.2951-2965.

115. Klimova E.G. Estimation of a state of an environment by the use of the observations and mathematical models. In: Air Pollution Processes in Regional Scale, NATO Science Series, IV. Earth and Enviromental Sciences- Vol.30,2003, p.185-190.

116. Lermusiaux P.E.J., Robinson A.R. Data assimilation via error subspace statistical estimation. Part I: Theory and schemes. Monthly Weather Review, 1999, v. 127, p.1385-1407.

117. Li Y. A note on the uniqueness problem of variational adjustment approach to four-dimensional data assimilation. Journal of Meteorological Society of Japan, 1991, v.69, N5, p.581-585.

118. Lorenc А.С. A global three-dimensional multivariate statistical interpolation scheme. Monthly Weather Review, 1981, v. 109, N4, p.701-721.

119. Lorenc A.C. Analysis methods for numerical weather prediction. Quarterly Journal of the Royal Meteorological Society, 1986, 112, p.1177-1194.

120. Lorenc A.C. Optimal nonlinear objective analyses. Quarterly Journal of the Royal Meteorological Society, 1988, v. 114, p.205-240.

121. Lorenc A.C., Bell R.S., Macphersson B. The Meteorologocal office analysis correction data assimilation scheme. Optimal nonlinear objective analyses. - Quarterly Journal of the Royal Meteorological Society, 1991, v.117, p.59-89.

122. Lorenc A., Rutherford I., Larsen G. The ECMEF analysis and data assimilation scheme analysis of mass and wind fields. - Technocal Report N6, 1976, ECMWF. -45p.

123. Manabe S., Holloway I.L., Spellman M.I. GFDL global 9-level atmospheric model. In: Modelling for the first GARP global experiment- GARP Publ.Ser., 1974, N14, p.7-26.

124. McPherson R.D., Bergman K.H., Kistler R.E., Rash G.E., Gordon D.S. The NMC operational global data assimilation system. Monthly Weather Revier, 1979, v.107, N11, p.1445-1461.

125. Mitchell H.L., Houtekamer P.L. An adaptive ensemble Kalman filter. -Monthly Weather Review, 2000, v. 128, p.416-433.

126. North G.R. Empirical ortogonal function and normal modes, Journal of the Atmospheric Sciences, 1984, v.41, N5, p.879-887.

127. Ogorodnikov V.A., Protasov A.V. Dynamic probabilistic model of atmospheric processes and variational methods of data assimilation. -Russian J. Numer. Anal. Math. Modelling, 1997, Vol.12, N5, p.461-479.

128. Phillips N.A. On the completness of multivariate optimum interpolation for large-scale meteorological analysis. Monthly Weather Review, 1982, v. 110, N 10, p.1329-1334.

129. Phillips N.A. The spatial statistics of random geostrophic modes and first-guess errors. Tellus, 1986, vol.38A, p.314-332.

130. Penenko V.V. Some aspects of mathematical modeling using the models together with observational data. Bull.Nov.Сотр.Center, Num.Model, in Atmosph., etc., 1996, N4, p.31-52.

131. Platov G.A., Golubeva E.N. The Japan Sea circulation modelling with data assimilation and analysis of deep water strucrure. Математические методы в геофизике. Труды Международной конференции, часть II, Новосибирск, 2003, с.374-379.

132. Rabier F., J-N Thpaut, Courtier P. Extended assimilation and forecast experiments with a four-dimensional variational assimilation system, Quarterly Journal of the Royal Meteorological Society, 1998, 124, p. 18611887.

133. Research Activities in Atmospheric and Oceaning Modelling fed. A. Staniforth. 1995, Rep. 21, WMO/TD-No.665.

134. Rivin G.S. Modern computational technologies for the estimaton of atmospheric processes on regional scales. In: Air Pollution Processes in Regional Scale, NATO Science Series, IV. Earth and Enviromental Sciences - Vol.30,2003, p.241-248.

135. Rutherford I. An operational multi-variate statistical objective analysis scheme. In: The GARP Programme on Numerical Experimentation. Rep. N11, p.98-121.

136. Sasaki Y. An objective analysis based on the variational method. -J.Meteor.Soc.Japan, 1958, 36, p.77-88.

137. Shlatter T.W. Some experiments with a multivariate statistical optimum interpolation. Monthly Weather Revie, 1975, v.103, N3, p.246-257.

138. Schlatter T.W., Branstator G.W. Estimation of errors in Nimbus-6 temperature profiles and their spatial correlation. Monthly Weather Review, 1979, v.107, N10, p.1402-1413.

139. Shutyaev V. Data assimilation for the time-dependent transport problem.-Evolution equations: existence, regularity and singularities, Banach center publications, vol.52, Institute of mathematics Polish Academy of Sciences, Warshawa, 2000, p.213-220.

140. Shutyaev V. Some properties of the control operator in the problem of data assimilation and iterative algorithms. Russ. J. Numer. Anal. Math. Modelling, 1995, N10, p.357-371.1. Литература 233

141. Talagrand О., Courtier P. Variational assimilation of meteorological observations with the adjoint vorticity equation. I: Theory. Quarterly Journal of the Royal Meteorological Society, 1987, 113, p.1311-1328.

142. Todling R., Cohn S. Suboptimal shemes for atmospheric data assimilation based on the Kalman filter. Monthly Weather Review, 1996, v. 124, p.2530-2557.

143. Tsyroulnicov M.D. Model-error models: identification and estimation. -EGS-AGU-EGU Joint Assembly, Nice, France, 6-11 April 2003.

144. Tsyroulnicov M.D. An isotropy-like property of meteorological fields. -Quarterly Journal of the Royal Meteorological Society, 2000, 126, p. 1-10.151. http://gmao.gsfc.nasa.gov152. http://www.ecmwf.int