Одной из главных задач атмосферной оптики является корректная интерпретация характеристик атмосферы и различных атмосферных образований [1–3] как экспериментально, так и теоретически и попытка их обобщения для дальнейшего использования. С развитием лазерных технологий все большее внимание стало уделяться возможности зондирования атмосферы лидарными методами с использованием лазеров, генерирующих мощное тераваттное и мультитераваттное излучение [4–9]. На первых этапах численного моделирования рассеяния света рассматривается излучение на длине волны 0,55 мм, например [10–14], и в дальнейшем в зависимости от поставленных задач включаются, например, различные длины волн, поляризация излучения, турбулентность атмосферы и т.д. Отдельной задачей стоит зондирование окружающей среды предельно короткими (по отношению к длине волны зондирования) импульсами. Создаются системы зондирования на основе лазеров, излучающих импульсы длительностью от нескольких десятков фемтосекунд, проводятся эксперименты по исследованию окружающей среды (атмосферы, водной и земной поверхности, растительного покрова, облачности) с целью регистрации сигналов флуоресценции, определения состава атмосферы на трассах до десятков километров, высвечиваемых излучением «суперконтинуума», возникающего в результате спектрального и временного уширения зондирующего импульса [4–9]. За счет эффектов фемтосекундной нелинейной оптики спектр трансформированного излучения расширяется на несколько спектральных октав – от ультрафиолетовой до средней и даже дальней инфракрасных областей. Параметры излучения СК остаются в пределах параметров лазерного излучения. И, как следствие, излучение СК признано перспективным для лазерного зондирования атмосферы.
Материалы и методы исследования
Основная цель данной работы – изучение возможности применения сверхкоротких лазерных импульсов в задачах зондирования. Методической основой проекта является комплексность исследований, сочетающая в себе натурные, лабораторные, численные эксперименты, использование современной материальной базы.
Основной метод, используемый в экспериментальной части данной работы, – метод лазерного зондирования окружающей среды и различных трасс распространения (включая более плотные, чем воздух, образования различной природы) мощными фемтосекундными импульсами (регулируемая энергия импульса от 10 до 100 мДж длительностью несколько десятков фемтосекунд). Применение данного активного метода зондирования по изучению распространения мощного тераваттного фемтосекундного импульса базируется на использовании фемтосекундной лазерной системы Института оптики атмосферы СО РАН [5, 6]. Основные подходы в экспериментальной части работы при использовании мощного тераваттного излучения фемтосекундного Ti-Sa лазера будут следующие:
1) использование различных схем лазерного зондирования (излучатель и приемник располагаются либо в одном месте, либо разнесены в одной горизонтальной плоскости);
2) в качестве зондирующего излучения в режиме генерации суперконтинуума будут использоваться длины волн 785–815 нм (длительность 20 фс, энергия 10–100 мДж, мощность 0,5–5 тВт);
3) использование второй и третьей гармоники.
Еще одной составляющей при выполнении данного проекта предполагается использование метода статистического и численного моделирования. В частности, использование разработанной программы в ИОА СО РАН при непосредственном участии членов научного коллектива [10–14] на основе метода геометрической оптики. Данный программный продукт позволяет решить задачу по рассеянию света на ледяных кристаллах перистой облачности, что дает возможность объяснить получаемые экспериментальные данные при дистанционном зондировании атмосферных трасс с учетом наличия на них ледяных кристаллических частиц, а также перистой облачностью. Комплексный подход по исследованию характеристик распространения мощного фемтосекундного лазерного излучения, например оптоакустических, геометрических, физических, колориметрических, позволяет получить более точную информацию об особенности распространения тераваттного излучения в атмосфере. Численное решение уравнения лазерного зондирования с различными начальными характеристиками делает возможной интерпретацию полученных результатов.
Результаты исследования и их обсуждение
В данной работе зондирование слоя облачности проводилось под углом 30 ° относительно горизонта (рис. 1). Для выделения излучения обратного рассеяния в спектре СК был установлен фильтр СЗС 23.
Спектральный диапазон его работы практически совпадет с полосой пропускания используемого ФЭУ, что позволило регистрировать обратное рассеяние от СК на длинах волн 0,35–0,65 мкм и 1–1,1 мкм. При проведении зондирования установленный фильтр толщиной 2 мм полностью не пропускал рабочую (лазерную) длину волны. На рис. 1 представлен спектр обратного рассеяния от слоя облачности.
Рис. 1. Запись лидарного сигнала от слоя облачности при зондировании на наклонной трассе: черный – излучение обратного рассеяния на сумме длин волн лазерная + СК; темно-серый – обратное рассеяние на лазерной длине волны 0,8 мкм; светло-серый – обратное рассеяние на длинах волн СК
Рис. 2. Сфера направлений рассеяния света при длине волны 0,55 мкм на ледяных кристаллических пластинчатых кристаллах перистой облачности при угле падения 30 ° относительно зенита
По данным зондирования был оценен показатель ослабления лазерного излучения вдоль трассы. Судя по его высотному распределению и протяженности слоя облачности можно предположить, что на участке 450–600 м находилась слоистая облачность, высота основания которой обычно составляет 0,1–0,7 км, а толщина слоя – от 0,2 до 0,8 км. На участке 1150–1250 м – слоисто-кучевая облачность. Высота основания составляет 0,5–1,5 км, толщина слоя – от 0,2–0,8 км.
Вторая часть работы заключалась в численном моделировании рассеяния света на длине волны 0,55 мкм на примере наиболее распространенных ледяных кристаллов перистой облачности, гексагональных пластинкок и столбиков (рис. 2).
Отметим, что актуальность проводимых исследований подтверждается тем, что, во-первых, натурные измерения влияния перистой облачности показывают, что вклад в парниковый эффект от ее наличия в атмосфере сопоставим с вкладом от парниковых газов, образующихся от жизнедеятельности человечества. Во-вторых, фемтосекундную систему ИОА СО РАН можно использовать для изучения перистых облаков как 0,55 мкм, так и с помощью волнового пакета. Дальнейшая модификация численного алгоритма для расчета на других длинах волн позволит установить основные закономерности изменений в матрице Мюллера.
Несмотря на то, что данная задача была частично рассмотрена [10–14], следует уделить отдельное внимание задаче рассеяния света на горизонтально и преимущественно ориентированных гексагональных пластинках при различных начальных параметрах. Во-первых, потому что, как уже было отмечено, формы данных кристаллических частиц являются наиболее распространенными. Во-вторых, ледяные гексагональные пластинки вносят главенствующий вклад в образование наиболее часто встречаемых гало, например, таких как паргелический круг (Parhelic circle), субпаргелический круг (Subparhelic circle), ложное солнце (Sundog) паргелий 120º (Parhelion 120 °). В-третьих, последние исследования по лидарному зондированию перистой облачности показали [15], что для определения наличия ориентированных гексагональных кристаллов предпочтительным является использование сканирующего лидара и выбора наиболее информативного элемента по поляризационным характеристикам. А от ориентации кристаллов, ее частоты существования зависит зеркальное отражение солнечного излучения, которое дает вклад в парниковый эффект.
Исходя из анализа геометрии пучков в кристаллах, были сделаны следующие выводы. В соответствии с законом Снеллиуса (q = arcsin(n2/n1)), если присутствует полное внутреннее отражение, то энергия пучка без потерь остается в пучке в данной среде распространения. При этом линейная поляризация света остается линейной с возможным наклоном. Круговая поляризация образуется лишь, когда присутствует полное внутреннее отражение. Как было показано ранее, если углы падения света на шестигранники пластинчатых кристаллов лежат промежутке π/2-θ* < θ0 < θ* , где θ* ≈ 58 ° [11–13],, то можно утверждать, что присутствует эффект полного внутреннего отражения. Вследствие геометрии наиболее распространенных кристаллов перистых облаков получим при параметре формы 0,4 с учетом, что угол между вертикалью и направлением распространения излучения будет примерно 34 градуса, будет присутствовать траектория с сохранением энергии излучения внутри кристалла. При уменьшении относительных размеров частицы рассмотренная траектория будет видоизменяться путем добавления дополнительных траекторий с дополнительными полными отражениями от шестиугольных граней. И, как следствие, будет происходить постепенное увеличение доли рассеянной компоненты M44.
Таким образом, было показано, что на основе анализа поведения рассеянного на кристаллах излучения можно сделать вывод об относительных размерах частиц. Дальнейшая модификация алгоритма предполагает как лабораторные исследования процесса распространения излучения, например, фемтосекундной длительности уже в различных кристаллах, так и натурные эксперименты на основной длине волны Ti-Sapphire лазера и конической эмиссии суперконтинуума. Совместные исследования атмосферы широким диапазоном спектра КЭ СК и TomSky позволят дать более детальную информацию о состоянии и динамике атмосферы.
Заключение
Результаты проведенных экспериментов позволили сделать вывод, что зондирование можно проводить как на основной длине волны лазера, так и на второй, третьей и даже четвертой гармонике. Для реализации гармоник, отличных от лазерной, могут быть использованы преобразователи частоты в виде нелинейных кристаллов или перестраиваемый лазер.
В экспериментальной части работы использовалась уникальная фемтосекундная система, характеристики которой кратко описаны выше, в совокупности с телескопом конструкции Ньютона с диаметром главного зеркала 30 см и фокусным расстоянием 80 см. Были получены результаты по распределению энергии тераваттного излучения между приосевой частью лазерного излучения и конической эмиссии суперконтинуума.
В численной части на основе метода геометрической оптики было исследовано поведение элементов матрицы рассеяния и определен наиболее чувствительный элемент для идентификации размеров частиц на длине волны 0,55 мкм. С одной стороны, полученные данные могут быть использованы в интерпретации получаемых обратных сигналов. С другой стороны, модификация численного алгоритма позволит сделать расчеты как для других длин волн, в частности для основной длины волны Ti-Sapphire лазера 0,8 мкм, так и для других кристаллических структур.
На основе полученных данных был сделан вывод о перспективности применения фемтосекундного излучения для широкого круга исследований вследствие отсутствия необходимости перестройки лазера на одной длине волны и возможности получения экспериментальных данных сразу на совокупности длин волн, что повышает оперативность получения информации в широком спектральном диапазоне и позволяет более точно определять разнообразные характеристики зондируемой атмосферы.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 16-35-00173 мол_а и РНФ № 16-17-10128.