Дистанционное лазерное определение параметров атмосферного аэрозоля и облаков с использованием методов многоволнового зондирования и явлений многократного рассеяния излучения тема диссертации и автореферата по ВАК РФ 05.11.07, кандидат технических наук Коренский, Михаил Юрьевич

• Коренский, Михаил Юрьевич
• кандидат технических наук
• 2007

Оглавление диссертации кандидат технических наук Коренский, Михаил Юрьевич

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

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

1.1. Информационное содержание данных, полученных с использованием лидара с переменным полем зрения.

1.2. Решение обратной задачи зондирования для лидара с переменным полем зрения.

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

1.4. Выводы к главе.

Глава 2. Основные параметры лидарной системы, разработанной для проведения исследований.

2.1. Технический облик системы.

2.2. Многоволновые измерения.

2.3. Флуоресцентные измерения.

2.4. Выводы к главе.

Глава 3. Результаты экспериментальных исследований.

3.1. Определение параметров аэрозоля методом многоволнового зондирования.

3.2. Определение параметров частиц по данным лидара с переменным полем зрения.

3.3. Измерение спектров флуоресценции.

3.4. Выводы к главе.

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

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

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

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

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

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

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

Параметры относительно мелкого аэрозоля (от 0,05 до 2 мкм) могут быть определены с использованием многоволнового лидара, данный подход получил широкое распространение в последнее десятилетие [1621]. Восстановление микрофизических параметров аэрозоля обычно разбивается на две независимые задачи: вычисление коэффициентов обратного рассеяния (3 и общего ослабления а (экстинкции) из результатов лидарных измерений и, затем, восстановление по ним микрофизических параметров аэрозоля, таких как: размер, концентрация и показатель преломления. Стабильность решения соответствующей обратной задачи сильно зависит от точности вычисления входных оптических данных (Р и а), поэтому уменьшение погрешности их определения является принципиальным моментом в многоволновых лидарных измерениях.

Сложность решения лидарного уравнения обусловлена тем, что оно содержит два неизвестных (коэффициент обратного рассеяния р и экстинкцию а) и для его решения необходимо задаваться предположениями о соотношении между ними. На этом основаны широко используемые методы Клетта, Ферналда и метод асимптотического сигнала [22-24, 4]. Как правило, подобная информация не доступна в процессе лидарных измерений и коэффициенты р и а вычисляются со значительной погрешностью.

Точность измерения р и а повышается при использовании лидара комбинационного рассеяния [25,26]. Следует оговориться, что в современных публикациях лидары комбинационного рассеяния повсеместно называются рамановскими, поэтому в процессе изложения будет использоваться именно этот термин. Рамановский метод позволяет независимо вычислять коэффициенты (3 и а за счет одновременной регистрации сигналов упругого рассеяния и сигналов на длине волны, соответствующей рамановскому сдвигу молекулы кислорода или азота. При этом предполагается, что высотное распределение плотности воздуха известно.

Обратная задача нахождения распределения аэрозоля по размерам из-набора входных оптических данных аир для различных длин волн формулируется в виде интегрального уравнения Фредгольма второго рода и относится к классу некорректных задач. Для ее решения обычно используется метод регуляризации Тихонова [27-29]. Впервые возможность восстановления оптических параметров подобным методом была продемонстрирована в работе [30], при этом точность определения параметров частиц была невысокой. Поэтому применение метода регуляризации в многоволновом лидарном зондировании стало возможным лишь относительно недавно, с появлением мощных рамановских лидаров, обеспечивающих достаточную точность измерений. Наиболее успешной реализацией такого подхода можно считать результаты, полученные в конце 90-х годов в Институте Тропосферных Исследований (ИТИ, Лейпциг), где была создана многоволновая система, использующая сигналы упругого и рамановского рассеяния для вычисления шести коэффициентов обратного рассеяния и двух коэффициентов экстинкции [17-20].

Гораздо более привлекательной является упрощенная версия лидара на основе Nd:YAG лазера с генератором третьей гармоники. Такой лидар использует упругое рассеяние (355, 532, 1064 нм) и сигналы рамановского рассеяния азота (387, 607 нм). Соответственно, могут быть определены три коэффициента обратного рассеяния и два коэффициента экстинкции. Результаты математического моделирования, проведенного в работах [3134], демонстрируют, что даже в такой упрощенной конфигурации лидар способен обеспечить достаточную точность оценки параметров аэрозоля.

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

Для получения информации о более крупных частицах (свыше 2 мкм) методы многоволнового зондирования оказываются непригодными, следовательно, необходимо использовать другие физические процессы, например, эффекты многократного рассеяния излучения в атмосфере. Многократное рассеяние излучения является источником погрешности при определении оптических параметров аэрозолей и для её минимизации измерения проводятся с минимально возможным полем зрения, либо же используют модели, учитывающие вклад многократного рассеяния в лидарный сигнал [35-39]. Вместе с тем, угловой спектр многократно рассеянного излучения содержит информацию о размерах частиц [40-44]. Возможность использования подобного метода впервые была продемонстрирована в работах Roy и Bissonnette [40-42], где было показано, что использование лидарных сигналов, соответствующих различным полям зрения, позволяет оценивать размеры частиц. Такой лидар получил название «лидар с переменным полем зрения». Несмотря на результаты, полученные в работах [40-42], многие вопросы, связанные с использованием лидара с переменным полем зрения, оставались невыясненными. В том числе не ясно какое количество полей зрения должно быть использовано в приемной системе, в каком диапазоне размеров частиц может быть использован данный метод и как этот диапазон зависит от геометрии рассеяния и погрешности измерения. Получение ответов на эти вопросы является одной из целей диссертации.

В традиционном лидаре с переменным полем зрения используется интегральная по азимутальному углу мощность рассеянного излучения. Однако хорошо известно, что для поляризованного лазерного источника мощность многократно рассеянного излучения может характеризоваться сильной азимутальной зависимостью. При этом ее вид определятся ориентацией поляризатора в приемной системе. Данная азимутальная зависимость содержит дополнительную информацию о параметрах частиц [45], которая может быть использована для решения обратной задачи зондирования. Азимутальная зависимость мощности рассеянного излучения при проведении экспериментов в аэрозольной камере, а также при зондировании атмосферных облаков, впервые наблюдалась авторами работ [46,47]. Зависимость параметров азимутальной структуры от размера частиц, их формы и оптической толщины была экспериментально исследована в аэрозольной камере в работе [48]. В частности, было показано, что азимутальные структуры характерны для процессов двукратного рассеяния излучения, что позволило авторам с их помощью оценить вклад высших порядков рассеяния. Вместе с тем, ряд принципиальных вопросов, связанных с возможностью использования азимутальных структур в лидарном зондировании, оставались открытыми. В том числе не вполне ясно, как радиальное и азимутальное распределение мощности в фокальной плоскости приемника излучения зависит от параметров частиц. Для получения этой информации необходимо разработать адекватную математическую модель, позволяющую вычислять элементы вектора Стокса рассеянного излучения для заданного поля зрения приемной системы и азимутального угла. В модель также должны быть введены такие параметры геометрии рассеяния, как высота основания облака и глубина зондирования.

Относительно простая модель, основанная на рассмотрении двукратного рассеяния излучения на сферических частицах и достаточно хорошо описывающая' экспериментальные данные, была разработана Rakovic и Kattawar [49-50]. К сожалению, эта модель не может быть использована для моделирования процесса лидарного зондирования, поскольку (а) она рассматривает сигнал, проинтегрированный по высоте, соответственно параметры азимутального распределения не зависят от геометрии рассеяния излучения; (б) предполагается, что излучение, рассеянное назад, распространяется параллельно, исходному. Таким образом, данный подход должен быть модифицирован применительно к существующей задаче лидарного зондирования. Подобная модификация была предложена в работе [51], однако в данном случае не рассматривалась азимутальная зависимость.

Помимо размера частиц в ряде приложений желательно иметь информацию об их составе и происхождении. Одним из способов получения подобной информации является измерение спектров лазерно-индуцировапиой флуоресценции. Данный подход широко используется, например, при изучении биоаэрозолей [52-56]. Флуоресцентная методика используется также во многих приложениях, таких как: контроль нефтяных пленок на поверхности океана [2-5,12,57-58], обнаружение опасных загрязняющих веществ на различных подстилающих поверхностях, контроль состояния растительного покрова [59]. Таким образом, возможность флуоресцентных измерений значительно расширяет потенциал лидарной системы. В результате, можно сделать вывод, что совокупное использование различных лидарных методов позволяет получать значительную часть информации о параметрах частиц.

Исследования в данном направлении ведутся в ряде научно-исследовательских институтов, в частности, в Институте Оптики Атмосферы им. В.Е. Зуева СО РАН, в Институте Физики Атмосферы им. A.M. Обухова РАН, в Институте Океанологии им. П.П. Ширшова и некоторых других. Результаты исследований отражены в публикациях: В.Е. Зуева, И.Э. Нааца, И.В. Самохвалова, Б.В. Кауля, С.В Самойловой и других российских ученых.

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

Результаты, представленные в работе, были получены в МГТУ им. Н.Э. Баумана при тесном сотрудничестве с Центром физического приборостроения ИОФ РАН им. А. М. Прохорова и Центром космических полетов (Годдард, США)

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

Данная цель включает следующие задачи исследования:

1) Оптимизация существующих и разработка новых алгоритмов определения коэффициентов обратного рассеяния и общего ослабления аэрозоля по данным многоволнового рамановского зондирования.

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

3) Разработка модели и проведение расчетов азимутальной зависимости мощности рассеянного излучения в фокальной плоскости премного устройства в приближении двукратно рассеяния поляризованного излучения.

4) Экспериментальное исследование возможности определения параметров частиц с использованием методов мультиспектрального и многократного рассеяния излучения.

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

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

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

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

Многоволновое зондирование.

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

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

Многократное рассеяние.

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

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

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

Флуоресцентные измерения.

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

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

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

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

Az к высоте основания облака т7=—, при этом азимутальное

Za распределение от rj не зависит.

3. Лидар с шестью полями зрения в диапазоне [0.25-5 мрад] при длине волны зондирования 1064 нм позволяет определять размеры частиц в интервале 1.5-25 мкм погрешностью -30%, при условии, что г] > 0.02.

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

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

Практическая ценность заключается в возможности использования полученных результатов:

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

- при моделировании изменений радиационного баланса планеты, обусловленных вариациями содержания аэрозоля естественного и антропогенного происхождения;

- при мониторинге загрязняющих веществ на различных подстилающих поверхностях естественного и искусственного происхождения.

Апробация результатов

Выносимые на защиту результаты опубликованы в ведущих международных научных журналах и представлены на конференциях: «23st International Laser Radar Conference» в июле 2006 г (Нара - Япония) и «Необратимые процессы в природе и технике» в МГТУ им Н.Э. Баумана в январе 2007 г. А также были проведены приемосдаточные испытания под контролем Государственной технической комиссии при Президенте России.

Публикации

По теме диссертации были опубликованы 8 печатных работ [45, 6970, 88-92]. Из них тезисы докладов конференций - 1, статьи в рецензируемых журналах из списка ВАК - 3, сборники трудов конференций - 4.

Структура и объём работы Диссертация состоит из введения, трех глав, заключения и библиографического списка из 92 наименований, изложена на 165 страницах, включая 5 таблиц и 95 рисунков.

Заключение диссертации по теме «Оптические и оптико-электронные приборы и комплексы», Коренский, Михаил Юрьевич

Выводы и заключение

В заключении сформулируем основные результаты, полученные в работе

1) Разработан метод анализа информации, содержащейся в данных лидарных измерений с переменным полем зрения, основанный на рассмотрении собственных чисел матриц ковариации, угловых распределений мощности двукратно рассеянного излучения. Показано, что при зондировании на длине волны 1064 нм и использовании полей зрения в диапазоне [0.25-5 мрад] такой лидар позволяет определять размер частиц в диапазоне радиусов 1,5<г<25 мкм. Установлено также, что при 10% погрешности измерений вариации мощности двукратно рассеянного излучения для типичных атмосферных аэрозолей содержат не более 6 независимых компонент, следовательно, использование большего количества полей зрения внутри выбранного диапазона не приводит к существенному улучшению точности решения обратной задачи.

2) На основе метода регуляризации Тихонова разработан алгоритм определения распределения частиц по размерам по данным лидара с переменным полем зрения. Достоверность работы алгоритма была подтверждена результатами численного моделирования и результатами обработки экспериментальных данных.

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

4) Предложен и реализован в виде программы алгоритм расчета коэффициентов обратного рассеяния и экстинкции аэрозоля по результатам многоволнового рамановского зондирования. Использование данного алгоритма позволяет проводить вычисления оптических коэффициентов с погрешностью менее 10%.

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

6) Исследована высотная вариация параметров аэрозоля в пограничном слое атмосферы методом многоволнового лидарного зондирования. Продемонстрировано, что высотные слои с относительной влажностью превышающей 80% характеризуются увеличением коэффициентов экстинкции и обратного рассеяния и размера частиц, что является следствием их гигроскопического роста. Проведено сравнение с данными солнечного радиометра основных параметров аэрозоля, таких как: средний и эффективный радиус, объемная концентрация и комплексный показатель преломления. Показано, что данные, полученные с использованием этих двух приборов, согласуются.

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

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

Список литературы диссертационного исследования кандидат технических наук Коренский, Михаил Юрьевич, 2007 год

1. Contribution of Working Group 1.to the Third Assessment Report of the Intergovernmental Panel on Climate Change (IPCC). / J.T. Houghton, Y. Ding, D.J. Griggs et al. //Cambridge University, England, 2001. 944 pp.

2. Межерис P. Лазерное дистанционное зондирование: Пер. с англ.- М: Мир, 1987. 510 с.

3. Инаба X. Лазерный контроль атмосферы. М.: Мир, 1979. 417 с.

4. Оптико-электронные системы экологического мониторинга природной среды. /В.И. Козинцев, В.М. Орлов, М.Л. Белов и др.; Под редакцией В.Н. Рождествина. М.: МГТУ им. Н.Э.Баумана, 2002. 528 с.

5. Зуев В.Е. Лазерное зондирование тропосферы и поверхности Земли. -Новосибирск: Наука, 1987. 321 с.

6. Зуев В.Е., Банах А.В., Покасов В.В. Современные проблемы оптики атмосферы. Т. 5. Оптика турбулентной атмосферы. М.: Гидрометеоиздат, 1988. 271 с.

7. Захаров В.М., Костко O.K., Хмелевцов С.С. Лидары и исследование климата. Ленинград: Гидрометеоиздат, 1990. 216 с.

8. Зуев В.Е., Зуев В.В. Дистанционное оптическое зондирование атмосферы-Ленинград: Гидрометеоиздат, 1992. 232 с.

9. Зуев В.Е., Титов Г.А. Современные проблемы оптики атмосферы Т. 9. Оптика атмосферы и климата. М.: изд. Спектр ИОА СО РАН, 1996.272 с.

10. Борн М., Вольф Э. Основы оптики. 2-е. изд. - М.: Наука, 1973. 856 с.

11. Ивлев J1.C., Андреев С.Д. Оптические свойства атмосферных аэрозолей-Ленинград: ЛГУ, 1986. 320 с.

12. Бункин А.Ф., Власов Д.В., Миркамилов Д.М. Физические основы лазерного аэрозондирования поверхности земли. Ташкент: издательство Фан, 1987. 272 с.

13. In-line type micropulse lidar with an annular beam: experiment. / T. Shiina, K. Yoshida, M. Ito, and Y. Okamura. //Appl. Opt. 2005. -Vol. 44, № 34. - P. 7407-7413.

14. Spinhirne J.D. Micro pulse lidar. // IEEE Trans. Geosci. Remote Sens. -1993.-Vol. 31.-P. 48-55.

15. Spinhirne J.D. Micro pulse lidar systems and applications. // Proceedings of 17th International Laser Radar Conference. 1994. - P. 162-165.

16. Scanning 6-wavelength 11-channel aerosol lidar. /D. Althausen, D. Miiller, A. Ansmann, U. Wandinger, H. Hube, E. Clauder, S. Zorner. // J. Atmos. and Oceanic Technol. 2000. - Vol. 17. - P. 1469 - 1482.

17. Miiller D., Wandinger U., and Ansmann A. Microphysical particle parameters from extinction and backscatter lidar data by inversion with regularization: theory. //Appl. Opt.- 1999. Vol. 38. - P. 2346-2357.

18. Miiller D., Wandinger U., and Ansmann A. Microphysical particle parameters from extinction and backscatter lidar data by inversion with regularization: simulation. //Appl. Opt- 1999. Vol. 38. -P. 2358-2368.

19. Miiller D., Wandinger U., and Ansmann A. Microphysical particle parameters from extinction and backscatter lidar data by inversion with regularization: experiment. // Appl. Opt- 2000. Vol. 39. -P. 879-1892.

20. Raman lidar measurements of aerosol extinction-to-backscatter ratio over southern Great Plains. / R.A. Ferrare, D. Turner, L. Brasseur, W. Feltz, O. Dubovik, T. Tooman. //J. Geophys. Res. 2001. - Vol. 106. -P. 20333-20348.

21. Raman lidar and sunphotometric measurements of aerosol optical properties over Thessaloniki, Greece during a biomass burning episode. /D.S. Balis, V. Amiridis, C. Zerefos, et al. //Atmospheric Environment. 2003. - Vol. 37. - P. 4529-4538.

22. Kovalev V.E., Eichinger W.E. Elastic Lidar: Theory, Practice, and Analysis Methods. New Jersey: John Wiley & Sons, inc., 2004. 620 pp.

23. Klett D. Stable analytical inversion solution for processing lidar returns. // Appl. Opt.- 1981. Vol. 20. - P. 211-220.

24. Fernald F.G. Analysis of atmospheric lidar observations: some comments. // Appl. Opt.- 1984. Vol. 23. - P. 652-653.

25. Тихонов A.H., Арсенин В.Я. Методы решения некорректных задач. -М.: Наука 1974. 285 с.

26. Twomey S. Introduction to the Mathematics of Inversion in Remote Sensing and Direct Measurements. New York: Elsevier, 1977. 245 c.

27. Зуев В.Е., Наац И.Э. Обратные задачи оптики атмосферы. М.: Гидрометеоиздат, 1989. 242 с.

28. Зуев В.Е Наац И.Э. Современные проблемы оптики атмосферы. Т. 7. Обратные задачи лазерного зондирования атмосферы. -Новосибирск: Наука, 1982. 287 с.

29. Inversion with regularization for the retrieval of tropospheric aerosol parameters from multiwavelength lidar sounding. / I. Veselovskii, A. Kolgotin, V. Griaznov, D. Muller, U. Wandinger, D. Whiteman. //Appl. Opt. 2002. - Vol. 41. - P. 3685-3699.

30. Retrieval of tropospheric aerosol parameters from multiwavelength lidar sounding. / I. Veselovskii, A. Kolgotin, V. Griaznov, D. Muller, U. Wandinger, C. Bockmann, D. Whiteman. //Proceedings of ILRC 21 Quebec City. July 2002. - P. 573-576.

31. Inversion of multiwavelength Raman lidar data for retrieval of bimodal aerosol size distribution. / I. Veselovskii, A. Kolgotin, V. Griaznov, D. Muller, K. Franke, D.N. Whiteman. //Appl. Opt. 2004. - Vol. 43. - P. 1180-1195.

32. Veselovskii I., Kolgotin A., Muller D. Retrieval of bimodal aerosol size distribution with multiwavelength Mie-Raman lidar. // 6-th International Symposium on Troposphere Profiling, Leipzig, Germany. 14-20 September, 2003. - P. 363-365.

33. Кауль Б.В., Самохвалов И.В. Уравнение лазерной локации атмосферы с учетом двухкратного рассеяния. // Изв. вузов, Физика, -1975.-№8.-С. 109-113.

34. Samokhvalov I.V. Double scattering approximation of lidar equation for inhomogeneous atmosphere. // Opt. Lett. 1979. - № 5. - P. 12-14.

35. Самойлова С.В. Теория кратного рассеяния и ее приложения к задачам лазерного зондирования аэрозолей. // Оптика атмосф. и океана.-2001.-Т. 14, №3,-С. 180-186.

36. Use of a polarization lidar for the retrieval of cloud optical parameters. / S.V. Samoilova, Y.S. Balin, G.G. Matvienko, M.V. Panchenko. // Seventh International Symposium on Atmospheric and Ocean Optics; Proc. SPIE. Dec. 2000. - Vol. 4341P. 488-496.

37. Самойлова C.B., Балин Ю.С. Использование поляризационного лидара для определения из космоса оптических параметров облачности. // Известия РАН, Физика атмосферы и океана. -2001. -Т. 37, №2.-С. 201-212.

38. Bissonnette L.R., Hutt D.L. Multiply scattered aerosol lidar returns: inversion method and comparison with in situ measurements. // Appl. Opt. 1995. - Vol. 34. - P. 6959-6975.

39. Estimation of cloud droplet size density distribution from multiple- fleld-of-view lidar returns. / G. Roy, L.C. Bissonnette, C. Bastille, G. Vallee, // Opt. Eng. 1997. - Vol. 36. - P. 3404-3415.

40. Retrieval of droplet-size density distribution from multiple-field-of-view cross-polarized lidar signals: theory and experimental validation. /G. Roy, L. Bissonnette, C. Bastille, G. Vallee. //Appl. Opt. 1999. - Vol. 38. -P. 5202- 5211.

41. Bissonnette L.R., Roy G., Roy N. Multiple-scattering-based lidar retrieval: method and results of cloud probing. // Appl. Opt. 2005. -Vol. 44.-P. 5565-5581.

42. THOR cloud thickness from offbeams lidar returns. / R.F. Cahalan, M. McGill, J. Kolasinski, T. Varnai, K. Yetzer. //J. of Atmospheric and Ocean Tech. - 2005. - Vol. 22. - P. 605-627.

43. Information content of data measured with a multiple-field-of-view lidar. /1. Veselovskii, M. Korenskii, V. Griaznov, D.N. Whiteman, M. McGill, G. Roy and L. Bissonnette. // Appl. Opt. 2006. - Vol. 45. - P. 68396848.

44. Carswell A.I., Pal S.R. Polarization anisotropy in lidar multiple scattering from clouds. //Appl. Opt. 1980. - Vol. 19. - P. 4123-4126.

45. Pal S.R., Carswell A.I. Polarization anisotropy in lidar multiple scattering from atmospheric clouds. // Appl. Opt. 1985. - Vol. 24. - P. 3464-3471.

46. Measurement of the azimuthal dependence of cross-polarized lidar returns and its relation to optical depth. / N. Roy, G. Roy, L.R. Bissonnette and J-R. Simard. //Appl. Opt. 2004. - Vol. 43. - P. 2777-2785.

47. Rakovic M.J. and Kattawar G.W. Theoretical analysis of polarization patterns from incoherent backscattering of light. //Appl. Opt. 1998. -Vol. 37.-P. 3333-3338.

48. Light backscattering polarization patterns from turbid media: theory and experiment. / M.J. Rakovic, G.W. Kattawar, M. Mehrubeoglu, B.D. Cameron, L.V. Wang, S. Rastegar, and G.L. Cote. // Appl. Opt. 1999. -Vol. 38, № 15. - P. 3399-3408.

49. КаиГ В.V., Samokhvalov I.V. Double scattering approximation of the atmospheric laser location equation taking polarization effects into account. // Izvestia Vuzov: Fizika. 1976. - № 1. - P. 80-85.

50. Ben-David A., Ren H. Detection, identification, and estimation of biological aerosols and vapors with a Fourier-trans form infrared spectrometer. // Appl. Opt. 2003. - Vol. 42, № 24. - P. 4887^1900.

51. Emission wavelength dependence of fluorescence lifetimes of bacteriological spores and pollens. / A. Thomas, D. Sands, D. Baum, L. To, and G.O. Rubel. //Appl. Opt. 2006. - Vol. 45, № 25. - P. 66346639.

52. Flanigan D.F. Detection of organic vapors with active and passive sensors: a comparison. // Appl. Opt. 1986. - Vol. 25, № 23. - P. 42534260.

53. Oil spill fluorosensing lidar for inclined onshore or shipboard operation. /R. Karpicz, A. Dementjev, Z. Kuprionis, S. Pakalnis, R. Westphal, R. Reuter, and V. Gulbinas. // Appl. Opt. 2006. - Vol. 45, № 25. -P. 6620-6625.

54. Laser spectroscopy of mineral oils on the water surface. / S. Patsayeva, V. Yuzhakov, V. Varlamov, R. Barbini, R. Fantoni, C. Frassanito and A. Palucci. // Proceedings of EARSeL-SIG-Workshop LIDAR. Dresden, FRG.-2000.-June 16- 17.-P. 106-114.

55. Tan S., Narayanan RM. Design and performance of a multiwavelength airborne polarimetric lidar for vegetation remote sensing. // Appl. Opt. -2006. Vol. 45, № 25. - P. 2360-2368.

56. Information content of multiwavelength lidar data with respect to microphysical particle properties derived from eigenvalue analysis. / I. Veselovskii, A. Kolgotin, D. Miiller, D.N. Whiteman. //Appl. Opt. -2005. Vol. 44. - P. 5292-5303.

57. Bruscaglioni P., Ismaelli A., Zaccanti G. Monte Carlo calculations of lidar returns: procedure ans results of Florence group. // Appl. Phys. B: Lasers and optics. 1995. - Vol. 60. - P. 325-331.

58. Efficient technique to determine backscattered light power for various atmospheric and oceanic sounding and imaging systems. / I.L. Katsev, E.P. Zege, A.S. Prikhach, and I.N. Polonsky. // J. Opt. Soc. Am. A . -1997.-Vol. 14.-P. 1338-1346.

59. Bissonnette L.R. Lidar and Multiple Scattering. //Lidar : Range-Resolved Optical Remote Sensing of the Atmosphere. / Editor Clauss Weikamp-Springer, 2005. 455 pp.

60. Eloranta E.W. Practical model for the calculation of multiply scattered lidar returns. // Appl. Opt. 1998. - Vol. 37. - P. 2464-2472.

61. Sassen K. and Zhao H. Lidar multiple scattering in water droplet clouds: toward an improved treatment. //Opt. Rev. 1995. - Vol. 2. -P. 394-400.

62. Golub G.H., Heath M., and Wahba G. Generalized cross validation as a method for choosing a good ridge parameter. // Technometrics. 1979. -Vol. 21.-P. 215-223.

63. O'Sullivan F. A statistical perspective on ill-posed inverse problems. // Statistical Science 1986. - Vol. 1. - P. 502-527.

64. Sabatier P.C. Basic concepts and methods of inverse problems in Basic Methods of Tomography and Inverse Problems Bristol and Philadelphia : Adam Hilger, 1987. - 669 pp.

65. Polonskii I.N., Zege E., and Katsev I.L. Lidar Sounding of Warm Clouds and Determination of Their Microstructure Parameters. // Izvestiya, Atmospheric and Oceanic Physics. -2001. Vol. 37, № 5. - P. 624-632.

66. Bohren F. B. and Huffman D. R. Absorption and scattering of light by small particles. New York: John Wiley, 1983. 530 pp.

67. Whiteman David N. Examination of the traditional Raman1 lidar technique. I. Evaluating the temperature-dependent lidar equations. // Appl. Opt. 2003. - Vol. 42. - P. 2571-2592.

68. Whiteman David N. Examination of the traditional Raman lidar technique. II. Evaluating the ratios for water vapor and aerosols. // Appl. Opt. 2003. - Vol. 42. - P. 2593-2608.

69. Holben B.N., et al. AERONET-A federated instrument networkand data archive for aerosol characterization. / Remote Sens. Environ. -1998.-Vol. 66.-P. 1-16.

70. Variability of Absorption and Optical Properties of Key Aerosol Types Observed in Worldwide Locations. / O. Dubovik, B.N. Holben, T.F. Eck, A. Smirnov, Y.J. Kaufman, M.D. King, D. Tanre, and I. Slutsker. //J. Atmos. Sci. 2002. - Vol. 59. - P. 590-608.

71. Tang I.N. and Munkelwitz H. Composition and temperature dependence of the deliquescence properties of hygroscopic aerosols. // Atmospheric Environment. 1993. - Vol. 27A. - P. 467-473.

72. Behrendt and Nakamura T. Calculation of the calibration constant of polarization lidar and its dependency on atmospheric temperature. //Optics Express.-2002.-Vol. 10, №. 16.-P. 805-817.

73. Heaton H. I. Principal-components analysis of fluorescence cross-section spectra from pathogenic and simulant bacteria. // Appl.Opt. 2005. -Vol. 44. - P. 6486-6495.

74. Interpretation of Raman spectra of nitro-containing explosive materials. Part I: group frequency and structural class membership. / I.R. Lewis, N.W. Daniel, JR., and P.R. Griffiths. // Appl. Spectr. 1997 - Vol. 51. -P. 1854-1867.

75. Ершов А.Д., БалинЮ.С., Самойлова C.B. Обращение лидарных данных при исследовании оптических характеристик слабозамутненной атмосферы. // Оптика атмосф. и океана. 2002. -Т. 15, № 10.-С. 894-899.

76. Самойлова С.В., Балин Ю.С., Ершов А.Д. Устойчивый .метод восстановления оптических характеристик аэрозольных полей по данным комбинационного лидарного зондирования. // Известия РАН, Физика атмосферы и океана. 2003.- Т. 39, № 3. - С. 384-393.

77. Специальный принцип отбора решений при использовании регуляризирующего алгоритма Тихонова для обработки многоволновых лидарных данных. / А.В. Колготин, В.И. Алехнович, М.Ю. Коренский, К.Н. Камша, // Измерительная техника. 2005. -№ 10.-С. 14-19.

78. Veselovskii I., Korenskii M., Vartapetov S. Remote detection of narcotics and explosives by fluorescence lidar. // Proceedings of 23rd International Laser Radar Conference, Nara, Japan. 2006. - 24-28 July. - P. 95-99.

79. Епифанов A.C., Корейский М.Ю. Быстрая гладкая сплайн-аппроксимация. // Необратимые процессы в природе и технике -Конференция в МГТУ им. Н.Э. Баумана, Сб. труд, часть 1. Москва ФИАН. - 2007 - С. 260.