Эффективные алгоритмы моделирования переноса излучения в атмосфере тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат технических наук Миляков, Алексей Владимирович

• Миляков, Алексей Владимирович
• кандидат технических наук
• 2003

Оглавление диссертации кандидат технических наук Миляков, Алексей Владимирович

ВВЕДЕНИЕ.

ГЛАВА 1 СОВРЕМЕННЫЕ КОМПЬЮТЕРНЫЕ МОДЕЛИ АТМОСФЕРНОЙ РАДИАЦИИ.

1.1 Расчет характеристик атмосферной радиации.

1.2 Компьютерные модели радиации в атмосфере.

1.3 Базы спектроскопической информации.

Введение диссертации (часть автореферата) на тему «Эффективные алгоритмы моделирования переноса излучения в атмосфере»

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

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

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

Рис. 1 Упрощенная модель радиационного баланса атмосферы [1]

1-а)50/4==/'' = <т7;4 0) где =1376 Вг/м - солнечная постоянная, Н - уходящий поток ИК излучения, Те- радиационная температура Земли, и сг постоянная Стефана-Больцмана. Альбедо Земли а из наблюдений [3] было определенно л а = 0.30 ±0.01. Следовательно, /<' = 241Вт/м . Пренебрегая зависимостью пропускания или альбедо от температуры, можем записать /г = еоТ*, где 1\ - средняя температура поверхности планеты (7] - 288 К ), коэффициент излучения атмосферы £- = 0.62 и, следовательно -^- = 3.3. Величина, обратная этой, обычно используется в климатологии, для выражения чувствительности средней температуры земной поверхности к изменению уходящего (или равного ему приходящего) потоков излучения.

Д7; = (2)

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

-у 'У личных оценок значение X (от 0.3 К/(Вт*м ) до 1.4 КУ(Вт*м )) сильно зависит от используемой модели общей циркуляции, в основном из-за различного учета облачности. Таким образом, незначительные возмущения потоков могут повлечь заметные изменения в климате. Иллюстрацией этого факта может служить извержение вулкана Пинатубо в 1991 г., благодаря которому в стратосферу было выброшено большое количетсво аэрозоля, что, в свою очередь, повлекло глобальное понижение температуры Земли на 0.4 К.

Учет радиационных процессов в модели общей циркуляции атмосферы дТ осуществляется через скорость радиационного выхолаживания ~ К ^ сутки которая представляет собой скорость изменения температуры воздуха со временем под действием радиационного притока. Эта величина связана с дивергенцией вектора потока I7. Для плоскопараллельной атмосферы она записывается следующим образом: дТ 1 дН = д( с р дг = 1л - Г1 где /'' и /'' - значения восходящего и нисходящего потока, соответственно, на некоторой высоте.

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

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

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

Первые попытки моделирования процесса поглощения излучения в атмосфере производились еще в позапрошлом веке. В 1860 году известный английский физик Тиндал при описании механизма парникового эффекта рассматривал поглощение радиации углекислым газом [4].

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

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

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

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

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

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

В атмосфере кроме поглощения происходят также и процессы рассеяния, которые играют существенную роль при переносе излучения в атмосфере в видимом и в УФ диапазоне. Обзор методов решения уравнения переноса приведен в монографиях [6,7,8]. Для описания процессов рассеяния излучения в атмосфере в последние годы широко используются алгоритмы на основе метода Монте-Карло [7, 9]. Однако, для задач, в которых блок расчета рас* сеяния не является основным, широко используются приближенные методы расчета радиации, которые наряду с более высокой скоростью дают и вполне приемлемую точность.

Следует отметить, что несмотря на небольшой вклад УФ излучения в радиационный баланс Земли, роль УФ радиации велика. Это связано с тем, что УФ излучение оказывает сильное воздействие на биологические и химические процессы.

Спектральный состав УФ-радиации в атмосфере определяется поглощением озоном в полосах Хартли, Хюгинса и кислородом (в короткой УФ области - (180 - 220 нм)). С точки зрения эффективности биологического воздействия принято выделять три области длин волн [10]: А (315<А„<400нм), В (280<Х<315нм) и С (А.<280нм). На прохождение УФ радиации помимо процессов поглощения сильно влияет аэрозольное рассеяние и рассеяние в облаках.

В данной работе предложена адаптация широко известного приближенного метода Эддингтона к расчетам радиации в УФ области спектра. Основа метода изложена в монографии Соболева В.В. [11], а краткое описание метода дано в работах [6, 12]. В работе [11] приведены формулы для коэффициентов пропускания, отражения, альбедо и диффузного пропускания атмосферы для случая «черной» подстилающей поверхности (альбедо поверхности равно нулю). При наличии отражающей поверхности в работе [11] приведены формулы упомянутых характеристик лишь для чистого рассеяния. В работе [6] приведены выражения альбедо и пропускания атмосферы для «черной» подстилающей поверхности.

Таким образом, полный набор формул, необходимый для моделирования интенсивности и потоков рассеянной радиации в приближении Эддинг-тона, в литературе отсутствует. В данной работе приводятся все необходимые формулы для моделирования диффузного излучения атмосферы приближенным методом Эддингтона, на основе которых было выполнено моделирование УФ радиации (области длин волн А, В). Рассмотрен также двухпо-токовый метод. Проведено сравнение с методом Монте-Карло и данными эксперимента.

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

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

1. Усовершенствовать существующий метод прямого расчета функций пропускания атмосферных трасс.

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

3. Провести исследование применимости метода Эддингтона для расчета рассеянной радиации в УФ области спектра.

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

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

Защищаемые положения

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

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

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

4. Созданный комплекс программ, позволяет проводить моделирование потоков атмосферной радиации.

Научная новизна

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

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

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

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

Практическая значимость

1. Проведено сравнение четырех алгоритмов расчета контура Фойгга (Wells - 99, Kuntz - 99, Humlicek - 82, Drayson - 76). Результаты сравнения могут использоваться для выбора алгоритма расчета контура Фойгта при создании программ полинейного счета коэффициента поглощения.

2. Рассчитанные полинейным методом коэффициенты поглощения в широком диапазоне длин волн (от 0 до 20000 см"1) для условий земной атмосферы, составили основу банков коэффициентов поглощения высокого и низкого спектрального разрешения.

3. Созданный банк коэффициентов поглощения низкого разрешения (Ли = 20см"1) и прикладное программное обеспечение, позволяют моделировать потоки и скорость радиационного выхолаживания в атмосфере.

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

5. Результаты диссертационной работы использовались при подготовке отчетов по федеральной целевой программе «Интеграция» (проекты А0061, Б0050/1038), фантов РФФИ (проекты №00-0790175, №02-07-03034, в последнем автор был руководителем), международной программы ЕС ИНКО Коперникус 2 (проект ИСИ-РЕМ, контракт ICA2 СТ- 2000 - 10024)

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

Апробация работы. Результаты работы докладывались и обсуждались на следующих конференциях: Региональная научно-техническая конференция «Радиотехнические и информационные системы и устройства» (17 мая 2000г.), 7 международный симпозиум «Оптика атмосферы и океана» (16-19 июля 2000г.), международная конференция «Моделирование, базы данных и информационные системы для атмосферных наук» (25-29 июня 2001г.), 9 международный симпозиум «Оптика атмосферы и океана» (2-5 июля 2002г.), международная конференция «Моделирование, базы данных и информационные системы для атмосферных наук» (2002г.). Результаты внедрены в институте оптического мониторинга СО РАН, в институте оптики атмосферы СО РАН. (акты внедрения приведены в приложении)

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

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

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

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

Ивану Валентиновичу Бойченко - за ценные советы по распределению времени и усилий при написании диссертации,

Коллективу кафедры АСУ за поддержку и благожелательное отношение.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Миляков, Алексей Владимирович

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

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

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

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

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

ЗАКЛЮЧЕНИЕ

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

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

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

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

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

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

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

Список литературы диссертационного исследования кандидат технических наук Миляков, Алексей Владимирович, 2003 год

1. Fouquart Y. Radiative transfer in climate models. Phisically-based Modelling and Simulation of Climate and Climatic Cange, Shlesinger M.E. Ed. Kluwer Academic Publishers, 1988, Part 1. - P. 223-283.

2. Stephens, G.L., G.G. Campbell and Т.Н. vonder haar. 'Earth radiation budgets.' // J. Geophys. Res. 1981, v. 86. - P. 9739-9760.

3. А.А. Мицель, K.M. Фирсов, Б.А. Фомин. Перенос оптического излучения в молекулярной атмосфере. Томск: STT, 2000. - 444 с.

4. Stephens G.L. The Parametrization of Radiation for Numerical Weather Prediction and Climate Models. // Mounthly weather review 1984, v. 112. - P. 826867.

5. Ленобль Ж. Перенос радиации в рассеивающих и поглощающих атмосферах. Стандартные методы расчета. Л.: Гидрометеоиздат, 1990. - 263 с.

6. Марчук Г.И., Михайлов Г.А., Назаралиев М.А., Дарбинян Н.А., Каргин Б.А., Елепов Б.П. Метод Монте-Карло в атмосферной оптике Новосибирск: Наука, 1976.-215 с.

7. Зуев В.Е., Титов Г.А. Современные проблемы атмосферной оптики. Том

8. Оптика атмосферы и климат. Из-во «Спектр» Института оптики атмосферы СО РАН, 1966. - 271 с.

9. Назаралиев М.А. Статистическое моделирование радиационных процессов в атмосфере. Новосибирск: Наука. Сиб. Отд-ние, 1990. - 227 с.

10. Безопасность России. Экологическая диагностика. М: МГФ «Знание», «Машиностроение», 2000 - глава 6. - С. 249-260

11. Соболев В.В. Рассеяние света в атмосферах планет. М.: Наука, 1972. - 331 с.

12. Joseph J.H. and Wiscombe W.J. The Delta-Eddington approximation for radiative flux transfer//J.A.S., 1976, v33. P. 2452-2459

13. Ellingson R.G., Ellis J., Fels S. The lntercomparison of Radiation Codes Used in Climate Model. Long Wave Results // J. G. R. 1991. V. 96. № D5. -P. 8929-8953.

14. Kneisys F.X. et. al. User Guide to LOWTRAN-7. 11 ARGL-TR-86-01777, ERP № 1010, AFGI, Hansom AFB, MA 01731August 1988.

15. Смеркалов В.А. Прикладная оптика атмосферы. С-Пб.: Гидрометеоиз-дат, 1997.-334 с.

16. Зуев В. Е., Белов В. В., Веретенников В В. Теория систем в оптике дисперсных сред. Томск: изд-во «Спектр» ИОА СО РАН, 1997 - 402 с.

17. Розанов Е.В., Тимофеев Ю.М., Троценко А.Н. Сравнение приближенного и эталонного методов расчета характеристик радиационного теплообмена в атмосфере // Изв. АН СССР, Физика атмосферы и океана., 1990. Т. 26, № 6. - С. 602-606.

18. Selby J.E.A., McClatchey R.A. Atmospheric Transmittance from 0.25 to 28.5 jim: Computed Code LOWTRAN 2 // AFCRL-TR-72-0745. 1972. - AD 763 721.

19. G.P.Anderson, F.X. Kneizys, J.H. Chetwynd, J. Wang, M.L. Hoke, L S. Rothman, L.M. Kimball, R.A. McClatchey. FAS-CODE/MODTRAN/LOWTRAN: Past/Present/ Future // 18th Annual Review Conference on Atmospheric Transmission Models. 6-8 June 1995

20. McClatchey R.A., R.W. Fenn, J.E.A. Selby, F.E. Votz and J.S. Garing. Optical properties of the Atmosphere // AFCRPL-70-0527. AD A715270. - 1970.

21. McClatchey R.A., R.W. Fenn, J.E.A. Selby, F.E. Votz and J.S. Garing, Optical properties of the Atmosphere (3rd edition) // AFCRPL-72-0497 AD A753075- 1972.

22. Kneizys F.X., Shettle E.P., Gallery W.O., Chetwynd J.H., Abreu L.W., Selby J.E.A., Fenn R.W., McClatchey R.A. Atmospheric Transmittance/Radiance. Computer Code LOWTRAN 5 // AFGL-TR-80-0067. Hanscom. AFB. Mass. -1980.-233 p.

23. Kneizys F.X., Shettle E.P., Gallery W.O., Chetwynd J.H., Abreu L.W., Selby J.E.A., Clough S.A., Fenn R.W. Atmospheric Transmittance/Radiance! Computer Code LOWTRAN 6 // AFGL-TR-83-0187. Hanscom. AFB. Mass. -1983.-200 p.

24. Abreu L.W., Kneizys F.X., Anderson G.P., Chetwynd J.H., Mod-tran/Lowtran: Current Status. Future Plans // Proceedings of the 14th Annual Review Conference on Atmospheric Transmission Models. Hanscom. AFB. Mass.- 1991. P. 65-72.

25. Berk A., Berustein L.S., Robertson D.C. MODTRAN. A Moderate Resolution Model for LOWTRAN 7 // L-TR-89-0122. 1989.

26. Smith H.J.P., D.J. Dube, M E. Gardner, S.A. Clough, F.X. Kneizys, L S. Rothman, FASCODE Fast Atmospheric Signature Code (Spectral Transmittance and Radiance) // AFGL-TR-78-0081. - Air Force Geophys. Lab.: Hanscom Mass, USA, 1978.

27. Clough S.A. and F.X. Kneizys, Convolution Algorithm for the Lorenz Function // Appl. Opt., 18, 2329 1979.

28. S.A. Clough, F.X. Kneizys, L.S. Rothman, W.O. Gallery, Atmospheric Spectral Transmittance and Radiance: FASCODE IB. // Proc. Of SPIE, 277. Atm. Transm. -1981.

29. Chetwynd J.G., J. Wang, G.P. Anderson, FASCODE: An Update and Application in Atmospheric Remote Sensing. // Proc. of SPIE 2266 Optical Spectroscopic Techniques for Atmospheric Research. - 1994.

30. Robertson D., Bernstein L., Duff J. et al. SHARC, an Atmospheric Radiation and Transmittance Code for Altitudes from 50 to 300 km // Ibid. P. 107-118.

31. Edwards D.P. GENLN2. A general line-by-line atmospheric transmittance and radiance model. Version 3.0 // Description and user's guide: NCAR Technical Note. 1992. -NCARЯN-367+STR.- Boulder, Colrado.

32. Edwards D.P. Recent Development of the GENLN2 Line-by-Line Model Studies in Support of the UARS Project // Ibid. P. 119-128.

33. Clough S., Kneizis F., Davies R. Line Shape and the Water Vapor Continuum // Atm. Research. 1989. - #23. - P. 229-241.

34. Мицель A.A., Руденко В.П. Пакет прикладных программ для расчета энергетических потерь оптического излучения в атмосфере LARA-1 // Препр. 57. Томск: Изд ТФ СО АН СССР, 1988. - С.55.

35. Комаров B.C., Мицель A.A., Михайлов С.А., Пономарев Ю.Н., Руденко В.В., Фирсов K.M. Программное и информационное обеспечение задач оптики атмосферы. // Оптика атмосферы. 1988. Т. 1. № 5. - С. 82-89.

36. Firsov K M., Kataev M.Yu., Mitsel A.A. Atlas of atmospheric and trace gas absorption for gas analyses radiative transfer // Proc. SPIE. 1991. V. 1811. -P. 385-388.

37. Mitsel A.A., Ponomarev Yu.N., Firsov K.M., Ptashnic I.V., and Kataev M.Yu. The computer codes LARA and AIR A for transmittance and radiance: current status // J.Q.S.R.T. 1995. V. 54. № 3. - P. 559-572.

38. Зуев B E., Комаров B.C. Статистические модели температуры и газовых компонент атмосферы. Л.: Гидрометеоиздат, 1986.

39. Fomin В.A. Effective line-by-line technique to compute radiation absorption in gases / Preprint IAE-5658. Moscow: Russian Research Center "Kurchatov Institute", 1993. - 13 p.

40. Фомин Б.А., Троценко A.H., Романов С.В. Эффективные методы расчета оптических свойств газообразных сред // Оптика атмосферы и океана. 1994. -Т.7, № 9. - С. 1457-1462.

41. F. Perrin. Review of Existing Spectroscopic Databases for Atmospheric Applications // Spectroscopy from space. Edited by Jean Demaison, Kamil Sarka and Edward A. Cohen Series II: Mathematics, Physics and Chemistry - Vol.20. - P. 235-259.

42. Goldman, A., Rinsland, C.-P., Flaud, J.-M., Perrin, A. HN03 line parameters; 1996 HITRAN update and new results. // J. Quant. Specrosc. Radiat. Transf. -1998, V. 60. P. 851-861

43. Пташник И.В. Математическое моделирование оптического абсорбционного зондирования газового состава атмосферы. Канд. дисс. физ.-мат. наук. Томск. - 1996.

44. Carter С., Handy N.C. The variational method for the calculation of ro-vibrational energy levels. // Сотр. Phys. Rep. 1986. V.5. - P. 115-172

45. Псннер С. С. Количественная молекулярная спектроскопия и излуча-тельная способность газов,- М.: Изд-во иностр. лит-ры, 1963. 492 с.

46. Tennyson J. The calculation of the vibration-rotation energies of triatomic molecules using scattering coordinates // Сотр. Phys. Rep. 1986. V.4. - P. 1-36.

47. Ben-Aryeh Y., Effect of centrifugal distortion on line intensities in asymmetric-rotor molecules//J. Opt. Soc. Amer. 1970. -V. 60, N 11 - P. 1409-1475.

48. Braslawsky J., Ben-Aryeh Y. First-order intensity perturbations for llie vibration-rotation lines of asymmetric rotor: theory and application // J. Chem. Phys. -1909 -V. 51.-P. 2233-2241.

49. Legay F. Inlensite des raies d'unc bande dc vibration-rotation // Gahiers Physique. 1958. - Vol. 12. N 99. - P. 416-436.

50. Войцеховская О. К., Макушкин Ю. С. Влияние случайных резонансов па интенсивности колебательно-вращательных; линий молекул типа асимметричного волчка // Оптика и Спектроскопия. 1976. - Т. 41, вып. 1. - С. 40-44.

51. Собельман II. И. Введение в теорию атомных спектров. М.: Физмат-из, 1993.-640 с.

52. Несмелова JI. И., Творогов С. Д., Фомин В. В. Спектроскопия крыльев линий. Новосибирск: Наука. Сиб. отд-ние, 1977 - 141 с.

53. Несмелова Л. И., Родимова О. Б., Творогов С. Д. Контур спектральной линии и межмолекулярное взаимодействие. Новосибирск: Наука. Сиб. отд-ние, 1986.-215 с.

54. Armstrong R. L. Effect of collisions narrowing on atmospheric transmit-tance II Appl. Opt. 1975. - Vol. 14, N 1. - P. 56-60.

55. R.J. Wells, Rapid approximation to the Voigt/Faddeeva function and its derivatives // JQSRT 1999. - V.62. - P.29-48.

56. M. Kuntz, M.Hopfner Efficient line-by-line calculation of absorption coefficient // JQSRT 1999. - V.63. - P.97-114.

57. J. Humlicek. //JQSRT 1982. - V.27. - P. 437.

58. S.R. Draison, Rapid computation of the Voigt profile // JQSRT 1976, V.16, N7.-P. 611-614.

59. F. Schreier, The Voigt and complex error function: a comparison of computational methods // JQSRT 1992, V.48, N 5/6. - P.743-762.

60. А.А.Мицель, И.В.Пташник, А.В.Миляков Оптимизация line-by-line алгоритма расчета молекулярного поглощения. // Оптика атмосферы и океана. 2000,13, №12. - С. 1137-1141.

61. Агеев Б.Г., Ионин А.А., Мицель А.А., Несмелова Л И. и др. Исследование коэффициентов поглощения атмосферных газов на длинах волн генерации СО-лазера // Препринт ИОА АН СССР. 1989. - № 65. - 42 с.

62. Mlawer E.J., Clough S.A., Brown P.D., Tobin D.S. Collision-indused effects and the water vapor continuum. // Proceedings of the Eight Atmospheric Radiation Measurement (ARM) Science Team Meeting. Tuscon, Arisona, 1998. - P. 503511.

63. Арефьев В Н. Молекулярное поглощение водяным паром излучения в окне относительной прозрачности атмосферы 8-13 мкм // Оптика атмосферы. 1989. -Т.2,№10. - С. 1034-1054.

64. Robert J. Selby, L. Biberman. Infrared Continuum Absorption by Atmospheric Water Vapor in the 8-12-|дт Window. // Appl. Opt. 1976, V.15, №9. - P. 2085-2090.

65. Burch, D.E. and Alt. R.L. (1984) Continuum Absorption by H20 in the 700-1200cm'1 and 2400-2800 cm"1 Windows. // AFGL-TR-84-0128 AD A147391.

66. Clough S.A., F.X. Kneizys, R.W. Davis, R.R Gamache, and R.H. Tipping, Theoretical line shape for H20 vapor: Application to the continuum, in Atmospheric Water Vapor, edited by A. Deepak, T. Wilkerson, and L.Ruhnke. Academic Press, 1980.

67. Burch D.E., Continuum absorption by H20. // Tech. Rep. AFGL-TR-81 -0300 - Air Force Geophysics Laboratory, 1981.

68. Burch D.E., Absorption by H20 in narrow windows between 3000-4200 cm"1 // Tech. Rep. AFGL-TR-85-0036 - Air Force Geophysics Laboratory, 1985.

69. Burch D.E. and R.L.Alt, Continuum absorption in 700-1200 cm"' and 24002800 cm* windows // Tech. Rep. Air Force Geophysics Laboratory, 1984.

70. Tobin D.C., L.L. Strow, W.J.Lafferty, and W.B.Olson, Experimental investigation of the self- and n2-broadening continuum within the v2 band of water vapor. // Applied Optics, 1996. N 35 - P. 4724-4734.

71. Han Y., J.A. Shaw, J.H.Churnside, P.D. Brown, and S.A. Clough, Infrared spectral radiance measurements in the tropical pacific atmosphere // Journal of Geophysical Research, 1997 N. 102 - P. 4353-4356.

72. Tobin DC., WWW home page, 1999b, http://tvler.ssec.wisc.edii/davet/liome2/H2Qcontinua/

73. Draison S.R. Rapid computation of the Voigt profile. // JQSRT. 1976. V.16. No.7. - P. 611-614.

74. Мицель A.A., Пташник И.В., Фирсов K.M., Фомин Б.А. Эффективный метод полинейного счета пропускания поглощающей атмосферы // Оптика атмосферы и океана. 1995, Т.8, N 10.-С. 1547.

75. В Mitsei A.A., Firsov K.M. A Fast Line-by-Line Method // JQSRT. 1995. V. 54. N. З.-Р. 549-557.

76. E.E.Whiting //JQSRT 1968, V.8. -P.1379.

77. R.W. Hamming, Numerical methods for scientist and engineering. MC Graw-Hill Book Company, Inc. New York, San Francisco, Toronto, London, 1962.

78. The Infrared Handbook // Editor WolfW.L., Zissis G.J., Chapters., La Rocca A.J. Atmospheric Absorption. Washington, 1978. - 132 p.

79. Гуди P. Атмосферная радиация. M.: Мир, 1966. - 417 с.

80. Пеннер С.С. Количественная молекулярная спектроскопия и излуча-тельная способность газов. М.: ПЛ, 1963. - 494 с.

81. Зуев В.Е. Распространение видимых и инфракрасных волн в атмосфере.- М.: Советское радио, 1970. 496 с.

82. Ladenberg R., Reiche F. Uber selektive Absorption // Ann. Phys. 1913. V. 42. -P. 181.

83. Godson W.L. The computation of infrared transmission by atmospheric water vapour//J. Meteorol. 1955. V. 12. - P. 272.

84. Зуев B E., Творогов С.Д. // Изв. ВУЗов. Сер. Физика. 1965. №1. -С. 185.

85. Plass G.N. Models for spectral band absorption // J. Opt. Soc. Am. 1958. -V. 48.-P. 690-703.

86. Malkmus W. Random Lorentz Band with Exponential Tailed S-1 Line - Intensity Distribution Function // J. Opt. Soc. Am. 1967. - V. 57. № 3. - P. 323-329.

87. McClatchey R.A., Fenn R.W., Selby J.E.A., Volz F.E., GaringJ.S. Optical Properties of the Atmosphere // AFCRL-71-0279. ERP N 354. Bedford. Mass. 1971.-88 p.

88. Киселева М.С., Непорент B.C., Федорова Е.О. Поглощение инфракрасной радиации при неразрешенной структуре спектра для наклонных путей в атмосфере (действие (Н20 и С02) // Изв. АН СССР. Сер. ФАО. 1967. Т. 3. № 5. - С. 640-649.

89. Киселева М.С. Определение влажности газовых примесей по инфракрасным спектрам поглощения. Закономерности поглощения инфракрасной радиации парами воды при неразрешенной структуре спектра // Оптика и спектроскопия. 1968. Т. 24. № 3. - С. 401-407.

90. Гасилевич Е.С., Федорова Е.О., Киселева М.С., Гальцев А.П., Осипов В.М. О закономерностях поглощения инфракрасной радиации атмосферной углекислотой при больших оптических толщинах // Изв. АН СССР. Сер. ФАО. 1978. Т. 14. - С. 222-225.

91. Голубицкий В.М., Москаленко Н.И. Функции спектрального пропускания в полосах паров Н20 и С02 // Изв. АН СССР. Сер. ФАО. 1968. Т. 4. №3,-С. 346-359.

92. Голубицкий В.М., Москаленко Н.И. Измерения спектрального поглощения С02 в условиях искусственной атмосферы // Изв. АН СССР. Сер. ФАО. 1968. Т. 4. № 1. - С. 85-89.

93. Филиппов В.Л., Макаров А.С. О влиянии температуры на спектральную прозрачность атмосферы (области полос водяного пара и углекислого газа) // Изв. АН СССР. Сер. ФАО. 1970. Т. 15. № 5. - С. 557-560.

94. Филлипов B.JI. Спектральное пропускание инфракрасной радиации атмосферными газами по данным лабораторных и натурных исследований // Журн. прикл. спектроскопии. 1983. - Деп. в ВИНИТИ, per. № 1045-83.

95. Козлов С.Д., Макаров А.С., Филиппов B.JI. Исследование функции пропускания паров воды в области 17,5-23,5 мкм в лабораторных условиях // Изв. АН СССР. Сер. ФАО. 1982.-Т. 18. №2.-С. 161-165.

96. Филлипов В.Л., Козлов С.Д., Румянцева Н.А., Зиатдинова Н.М., Макаров А.С. Прозрачность атмосферы в диапазоне 1-14 мкм при высокой дальности видимости // Изв. вузов. Физика. 1984. - Деп. в ВИНИТИ, per. № 2483-84. - 52 с.

97. Броунштейн A.M., Фролов А.Д. О методике расчета спектрального пропускания в ИК окнах прозрачности атмосферы для приземных горизонтальных трасс // Труды ГГО. 1985. Вып. 496. - С. 70-79.

98. Tjemkes S.A., Holmiund К., Schmetz J. Current problems of the atmospheric radiation // Proceedings of the international Radiation Simposium, IRS'96 Fairbanks, Alaska, 19-24 August 1996. - P. 463-467.

99. Pierlussi J., Tomiyama K. Numerical methods for the generation of empirical and analytical transmittance functions with applications to atmospheric trace gases // Appl. Opt. 1980. - V. 19, No. 14. - P.2298-2309.

100. Arking A., Grossman K. The influense of line shape and band structure on temperatures in planetary atmospheres // Journ. Atmosph.Sci. 1972. - V.29. - P. 937-949.

101. Творогов С.Д. Применение рядов Дирихле в атмосферной спектроскопии // Оптика атмосферы и океана, 1997. Т. 10, №4-5 - С. 403-412.

102. Domoto G.A. Frequency integration for radiative transfer problem involving homogeneous non-gray gases: the inverse transmission function // Journ. Quant. Spectrsc. and Radiat. Transf., 1974 Vol. 14. - P. 935-942.

103. Lacis A.A., Oinas V. A description of the K-distribution methods for modelling nongray gaseous absorption, thermal emission, and multiple scattering in vertically inhomogeneous atmospheres. // J. Geph.Res., 1991 V.96, No.D5 - P. 9027-9063.

104. Goody R., West R., Chen L., Crisp D. The correlated-k method for radiation calculations in nonhomogeneous atmospheres. // J. Quant. Spectrosc. Radiat .Transfer., 1989, V.42, N.6. - P. 539-550.

105. Burch D.E., Howard J.N., and Williams D. Infrared transmission in synthetic atmospheres: Absorptions laws for overlapping bands // J. Opt. Soc. Am., 1956. -N46.-P. 452-455.

106. Hoover G.M., Hathaway C.E., and Williams D. Infrared Absorption of overlapping bands of atmospheric gases // Appl. Opt. 1967. N 6. - P. 481-487.

107. Tubbs L.D., Hathaway C.E., and Williams D. Futher studies of overlapping bands of atmospheric gases // Appl. Opt. 1967. N 6. - P. 1422-1423.

108. Isaaks R.G., Wang W.C., Worsman R.D., Goldenberg S. Multiple scattering and FASCODE models. // Appl.Opt., 1987. V.26, N.7 - P.1272-1281.

109. K.M.Firsov, A.A.Mitsel, Yu.N.Ponomarev, I.V.Ptashnik, Parametrization of transmittance for application in atmospheric optics. // J.Q.S.R.T. 1998, V. 59, N. 3-5.-P. 203-213.

110. А.В. Миляков, Аппроксимация функций пропускания рядом экспонент в многокомпонентной молекулярной атмосфере Земли. // Автоматизированные системы обработки информации, управления и проектирования: Сб.науч.тр., Т.7. Томск: изд-во ТУСУР, 2002 - С. 78-85.

111. Milyakov A.V. Approximation of the transition function by exponential series in multicomponent molecular atmosphere of the earth // Оптика атмосферы и океана. Тез. 9 межд. симп. 2002г. Томск, Россия., 2-5 июля - С. 58.

112. Шифрин К.С. Рассеяние света в мутной среде. M.,JI.: Гостехиздат, 1951 -288 с.

113. В В. Белов, М.В. Винарский, И.И. Ипполитов, А.В. Миляков, А.А. Ми-цель. оделирование рассеянной УФ радиации на основе приближенных методов // Оптика атмосферы и океана. В печати.

114. W.E.Meador, W.R.Weaver. Two-Stream Approximation to Radiative Transfer in Planetary Atmospheres: A Unified Description of Existing Methods and a New Improvement. //J. Atmos. Sci. V. 37. - P. 630-643.

115. Liow K.N. A numercal experiment on Chandrasecar's discrete-ordinate method for radiative transfer: Applications for cloudly and hazy atmospheres. // J. Atmos. Sci, 1973. V. 30 - P. 1303-1326.

116. Coakley, J.A., Jr., P.Chylek. The two-stream approximation in radiative transfer: Including the angle of the incident radiation. // J. Atmos. Sci.,1975. V. 32. - P. 409-418.

117. C.M. Chu, S.W. Churchill. Numerical solution of problems in multiple scattering of electromagnetic radiation. // J. Phys. Chem., 1955 V. 59. - P. 855-863.

118. Emslie, A G., and J.R. Aronson. Spectral reflectance and emittance of particulate materials. 1. Theory. //Appl. Opt., 1973. -V. 12. P. 2563-2572.

119. Ellingson,R.G., and Y.Fouquart. The Intercomparison of Radiation Codes Used in Climate Models: An Overwiew // J.Geophys.Res., 1991 V. 96. - P. 8929-8953.

120. B.A.Fomin and Yu.V.Gershanov. Tables of the Benchmark Calculations of Atmospheric Fluxes for the ICRCCM Test Cases. Part 1 : Long-Wave Clear-Sky Results // Preprint IAE-5981/1 1996 - 49 p.

121. Белинский В.А., Гараджа М.П., Меженная Л.М., Незваль Е.И. Ультрафиолетовая радиация Солнца и неба. М: МГУ, 1968. - 228 с.

122. Molina L.T., Molina L.J. Absolute absorption cross section of ozon in the 185 to 350 nm wavelength range // J. Geoph. Res. V.91, N D13, P. 14501-14508.

123. Елисеев A.A. Ипполитов И.И. Кабанов M.B. и др. Исследования ультрафиолетовой радиации на поверхности Земли // Оптика атмосферы и океана., 1994.-Т. 7, №5-С. 8-12.

124. Винарский Максим Владимирович, Исследование изменчивости потоков ультрафиолетовой радиации на поверхности Земли. Канд. дисс. физ.-мат. наук. Томск. - 2002. - 126 с.

125. Облака и облачная атмосфера. Справочник // Под ред. И.П. Мазина, А.Х. Хргиана. Л.: Гидрометеоиздат, 1989. - 646 с.

126. Д|фсКlop ( ИопрСКо! О ПСИ I p;i K.IIIMil I ожо.кч ических исследований и порачовапи:

127. У I 1Я-,1'/1(ДЛК) Прорек юр 1УСУТ1. Г "" офессор-.11. I ордов200л I.5 J1. I. I1.li.ioiiichki 1200л I1. ЛК1 nilIJU'IlIlVI

128. АКТ ВНЕДРЕНИЯ Пакета программ Spectral Transmittance Universal' V . i>

129. Экономический эффект от использования результатов работы пакета не определялся. Настоящий акт не является основанием для финансовых претензий.1. От ИОА СО РАН

130. Ведущий научный сотрудник, д.ф.-м.н.1. К.М.Фирсов1. С'с2003 г.

131. От ТУСУР Зав. кафедрой АСУ д.т.п.

132. А.М.Кориков « ¿0 » 2003 г.4Р