Солнечно-земная физика, 2016, том 2, № 3
Бесплатно
Основная коллекция
Издательство:
Институт солнечно-земной физики СО РАН
Наименование: Солнечно-земная физика
Год издания: 2016
Кол-во страниц: 82
Количество статей: 9
Дополнительно
Вид издания:
Журнал
Уровень образования:
Дополнительное профессиональное образование
Артикул: 349900.0007.01
Тематика:
ББК:
УДК:
ГРНТИ:
Скопировать запись
Фрагмент текстового слоя документа размещен для индексирующих роботов.
Для полноценной работы с документом, пожалуйста, перейдите в
ридер.
СОЛНЕЧНО-ЗЕМНАЯ ФИЗИКА Свидетельство о регистрации средства массовой информации от 29 января 2014 г. ПИ № ФС77-56768 Издается с 1963 года ISSN 2412-4737 DOI: 10.12737/issn. 2412-4737 Том 2, № 3, 2016. 81 с. Выходит 4 раза в год Учредитель: Федеральное государственное бюджетное учреждение науки Институт солнечно-земной физики Сибирского отделения Российской академии наук SOLAR-TERRESTRIAL PHYSICS Certificate of registration of mass media from January 29, 2014. ПИ № ФС77-56768 The edition has been published since 1963 ISSN 2412-4737 DOI: 10.12737/issn. 2412-4737 Vol. 2, Iss. 3, 2016. 81 p. Quarterly Founder: Institute of Solar-Terrestrial Physics of Siberian Branch of Russian Academy of Sciences Состав редколлегии журнала Editorial Board Жеребцов Г.А., академик — главный редактор, ИСЗФ СО РАН Zherebtsov G.A., Academician, Editor-in-Chief, ISTP SB RAS Степанов А.В., чл.-корр. РАН — заместитель главного редактора, ГАО РАН Stepanov A.V., Corr. Member of RAS, Deputy Editor-in-Chief, GAO RAS Потапов А.С., д-р физ.-мат. наук — заместитель главного редактора, ИСЗФ СО РАН Potapov A.S., D.Sc. (Phys.&Math), Deputy Editor-in-Chief, ISTP SB RAS Члены редколлегии Members of the Editorial Board Алтынцев А.Т., д-р физ.-мат. наук, ИСЗФ СО РАН Altyntsev A.T., D.Sc. (Phys.&Math.), ISTP SB RAS Мордвинов А.В., д-р физ.-мат. наук, ИСЗФ СО РАН Mordvinov A.V., D.Sc. (Phys.&Math.), ISTP SB RAS Максимов В.П., д-р физ.-мат. наук, ИСЗФ СО РАН Maksimov V.P., D.Sc. (Phys.&Math.), ISTP SB RAS Леонович А.С., д-р физ.-мат. наук, ИСЗФ СО РАН Leonovich A.S., D.Sc. (Phys.&Math.), ISTP SB RAS Тащилин А.В., д-р физ.-мат. наук, ИСЗФ СО РАН Tashchilin A.V., D.Sc. (Phys.&Math.), ISTP SB RAS Перевалова Н.П., д-р физ.-мат. наук, ИСЗФ СО РАН Perevalova N.P., D.Sc. (Phys.&Math.), ISTP SB RAS Белан Б.Д., д-р физ.-мат. наук, ИОА СО РАН Belan B.D., D.Sc. (Phys.&Math.), IAO SB RAS Муллаяров В.А., канд. физ.-мат. наук, ИКФИА СО РАН Mullayarov V.A., C.Sc. (Phys.&Math.), IKFIA SB RAS Деминов М.Г., д-р физ.-мат. наук, ИЗМИРАН Deminov M.G., D.Sc. (Phys.&Math.), IZMIRAN Обридко В.Н., д-р физ.-мат. наук, ИЗМИРАН Obridko V.N., D.Sc. (Phys.&Math.), IZMIRAN Гульельми А.В., д-р физ.-мат. наук, ИФЗ РАН Guglielmi A.V., D.Sc. (Phys.&Math.), IPE RAS Трошичев О.А., д-р физ.-мат. наук, ААНИИ Troshichev O.A., D.Sc. (Phys.&Math.), AARI Ермолаев Ю.И., д-р физ.-мат. наук, ИКИ РАН Yermolaev Yu.I., D.Sc. (Phys.&Math.), IKI RAS Мареев Е.А., чл.-корр. РАН, ИПФ РАН Mareev E.A., Corr. Member of RAS, IAP RAS Сафаргалеев В.В., д-р физ.-мат. наук, ПГИ РАН Safargaleev V.V., D.Sc. (Phys.&Math.), PGI KSC RAS Лазутин Л.Л., д-р физ.-мат. наук, НИИЯФ МГУ Lazutin L.L., D.Sc. (Phys.&Math.), SINP MSU Стожков Ю.И., д-р физ.-мат. наук, ФИАН Stozhkov Yu.I., D.Sc. (Phys.&Math.), LPI RAS Сомов Б.В., д-р физ.-мат. наук, ГАИШ МГУ Somov B.V., D.Sc. (Phys.&Math.), SAI MSU Лестер М., проф., Университет Лестера, Великобритания Lester M., Prof., University of Leicester, UK Йихуа Йан, проф., Национальные астрономические обсерватории Китая, КАН, Пекин, Китай Yan Yihua, Prof., National Astronomical Observatories, Beijing, China Панчева Дора, проф., Национальный институт геодезии, геофизики и географии БАН, София, Болгария Pancheva D., Prof., Geophysical Institute, Bulgarian Academy of Sciences, Sofia, Bulgaria Салахутдинова И.И., канд. физ.-мат. наук, ученый секретарь, ИСЗФ СО РАН Salakhutdinova I.I., C.Sc. (Phys.&Math.), Scientific Secretary, ISTP SB RAS Полюшкина Н.А., ответственный секретарь редакции, ИСЗФ СО РАН Polyushkina N.A., Executive Secretary of Editorial Board, ISTP SB RAS
СОДЕРЖАНИЕ Потехин А.П., Сетов А.Г., Лебедев В.П., Медведев А.В., Кушнарев Д.С. Перспективный радар НР-МСТ: потенциал и диагностические возможности ……………………….......................... 3–16 Клунко Е.В., Еселевич М.В., Тергоев В.И. Наблюдения ТГК «Прогресс» на оптическом телескопе АЗТ-33ИК …………………………………………………………………………………………………... 17–23 Еселевич М.В., Хахинов В.В., Клунко Е.В. Параметры оптических сигналов на телескопе АЗТ-33ИК, зарегистрированных в активном космическом эксперименте «Радар–Прогресс» …... 24–32 Липко Ю.В., Пашинин А.Ю., Рахматулин Р.А., Хахинов В.В. Эффекты в магнитном поле Земли от работы двигателей космического корабля …………………………………………………………… 33–40 Пилипенко В.А., Козырева О.В., Федоров Е.Н., Успенский М.Г., Кауристи К. Широтная амплитудно-фазовая структура МГД-волн: моделирование наблюдения на радаре STARE и магнитометрах сети IMAGE ………………………………………………………………………………….. 41–51 Ларюнин О.А. Численный синтез ионограмм в горизонтально-неоднородной ионосфере на основе модели комбинированного параболического слоя ………………………………………………... 52–58 Гуляева Т.Л. Модификация индексов солнечной активности в международных справочных моде лях ионосферы IRI и IRI-Plas в связи с пересмотром ряда чисел солнечных пятен ………………….. 59–68 Васильев П.А., Карпов И.В., Кшевецкий С.П. Моделирование воздействия внутренних гравитаци онных волн на состояние верхней атмосферы в период внезапного стратосферного потепления …. 69–73 Михалев А.В., Подлесный С.В., Стоева П.В. Свечение ночной атмосферы в RGB цветовом представлении ……………………………………………………………………………...……………... 74–80 CONTENTS Potekhin A.P., Setov A.G., Lebedev V.P., Medvedev A.V., Kushnarev D.S. Prospective IS-MST radar. Potential and diagnostic capabilities ………………………………………………………………... 3–16 Klunko E.V., Eselevich M.V., Tergoev V.I. Progress cargo spacecraft observed with the AZT-33IK optical telescope …………................................................................................................................................... 17–23 Eselevich M.V., Khakhinov V.V., Klunko E.V. Parameters of optical signals registered with the AZT-33IK telescope in active Radar–Progress space experiment ……………………………………… 24–32 Lipko Yu.V., Pashinin A.Yu., Rakhmatulin R.A., Khakhinov V.V. Geomagnetic effects caused by exhaust streams of spacecraft engines ……………………………………………………………………… 33–40 Pilipenko V.А., Kozyreva O.V., Fedorov E.N., Uspensky M.G., Kauristie K. Latitudinal amplitude-phase structure of MHD waves: STARE radar and image magnetometer observations and modeling ………………………………………………………………………………………………………… 41–51 Laryunin O.A. Numerical synthesis of ionograms in horizontally inhomogeneous ionosphere on the basis of compound parabolic layer model ……………………………………………………………………….. 52–58 Gulyaeva T.L. Modification of the solar activity indices in the International Reference Ionosphere IRI and IRI-Plas models due to recent revision of sunspot number time series ………………………….. 59–68 Vasilyev P.A., Karpov I.V., Kshevetskii S.P. Modeling of the effect of internal gravity waves on upper atmospheric conditions during sudden stratospheric warming ……............................................................. 69–73 Mikhalev А.V., Podlesny S.V., Stoeva P.V. Night airglow in RGB mode …………………………….. 74–80
Солнечно-земная физика. 2016. Т. 2, № 3 Solar-Terrestrial Physics. 2016. Vol. 2. Iss. 3 3 УДК 621.396.6, 551.551.5, 550.388.2 Поступила в редакцию 06. 05.2016 DOI: 10.12737/19444 Принята к публикации 15.07.2016 ПЕРСПЕКТИВНЫЙ РАДАР НР-МСТ: ПОТЕНЦИАЛ И ДИАГНОСТИЧЕСКИЕ ВОЗМОЖНОСТИ PROSPECTIVE IS-MST RADAR. POTENTIAL AND DIAGNOSTIC CAPABILITIES А.П. Потехин Институт солнечно-земной физики СО РАН, Иркутск, Россия, potekhin@iszf.irk.ru А.Г. Сетов Институт солнечно-земной физики СО РАН, Иркутск, Россия, setov@iszf.irk.ru В.П. Лебедев Институт солнечно-земной физики СО РАН, Иркутск, Россия, lebedev@iszf.irk.ru А.В. Медведев Институт солнечно-земной физики СО РАН, Иркутск, Россия, medvedev@iszf.irk.ru Д.С. Кушнарев Институт солнечно-земной физики СО РАН, Иркутск, Россия, ds_k@iszf.irk.ru A.P. Potekhin Institute of Solar-Terrestrial Physics SB RAS, Irkutsk, Russia, potekhin@iszf.irk.ru A.G. Setov Institute of Solar-Terrestrial Physics SB RAS, Irkutsk, Russia, setov@iszf.irk.ru V.P. Lebedev Institute of Solar-Terrestrial Physics SB RAS, Irkutsk, Russia, lebedev@iszf.irk.ru A.V. Medvedev Institute of Solar-Terrestrial Physics SB RAS, Irkutsk, Russia, medvedev@iszf.irk.ru D.S. Kushnarev Institute of Solar-Terrestrial Physics SB RAS, Irkutsk, Russia, ds_k@iszf.irk.ru Аннотация. В ближайшие годы в районе Иркутска планируется строительство радара, сочетающего в себе возможности радаров некогерентного рассеяния (НР) и мезосферно-стратосферно-тропосферных радаров (МСТ) [Жеребцов и др., 2011]. Радар НР-МСТ представляет собой фазированную решетку, состоящую из двух разнесенных антенных полотен, с цифровой многоканальной приемной системой, позволяющей проводить детальную пространственновременную обработку сигнала обратного рассеяния. Показаны характеристики, конфигурация и возможности антенной и приемопередающей систем радара. Оценивается потенциал радара в основных режимах работы, включающих исследование ионосферы методом НР на высотах более ~100 км и изучение атмосферы с помощью сигналов, рассеянных на флуктуациях показателя преломления, вызванных атмосферными турбулентными неоднородностями, на высотах до 100 км. Моделирование показало, что радар позволит проводить регулярные измерения параметров нейтральной атмосферы на высотах до 26 км, а также наблюдать мезосферное летнее эхо на высотах ~85 км в присутствии заряженных ледяных частиц (при увеличении числа Шмидта) и мезосферное зимнее эхо на высотах ~65 км при повышении фоновой электронной концентрации. Оценка возможностей радара при работе в режиме НР в высотных диапазонах 100–600 и 600–2000 км показала, что в дневных условиях при накоплении 10 мин верхняя граница определения электронной концентрации и температуры ионосферной плазмы составляет ~1500 и ~1300 км соответственно при стандартном отклонении не более 10 %. Верхний предел определения скорости дрейфа составляет ~1100 км при стандартном отклонении 45 м/с. Оценка интерферометрических возможностей Abstract. In the next few years, a new radar is planned to be built near Irkutsk. It should have capabilities of incoherent scatter (IS) radars and mesospherestratosphere-troposphere (MST) radars [Zherebtsov et al., 2011]. The IS-MST radar is a phased array of two separated antenna panels with a multichannel digital receiving system, which allows detailed space-time processing of backscattered signal. This paper describes characteristics, configuration, and capabilities of the antenna and transceiver systems of this radar. We estimate its potential in basic operating modes to study the ionosphere by the IS method at heights above 100 km and the atmosphere with the use of signals scattered from refractive index fluctuations, caused by turbulent mixing at heights below 100 km. The modeling shows that the radar will allow us to regularly measure neutral atmosphere parameters at heights up to 26 km as well as to observe mesosphere summer echoes at heights near 85 km in the presence of charged ice particles (an increase in Schmidt number) and mesosphere winter echoes at heights near 65 km with increasing background electron density. Evaluation of radar resources at the IS mode in two height ranges 100–600 and 600–2000 km demonstrates that in the daytime and with the accumulation time of 10 min, the upper boundaries of electron density and ionospheric plasma temperature are ~1500 and ~1300 km respectively, with the standard deviation of no more than 10 %. The upper boundary of plasma drift velocity is ~1100 km with the standard deviation of 45 m/s. The estimation of interferometric capabilities of the MST radar shows that it has a high sensitivity to objects of angular size near 7.5 arc min, and its potential accuracy in determining target angles can reach 40 arc sec.
А.П. Потехин, А.Г. Сетов, В.П. Лебедев, A.P. Potekhin, A.G. Setov, V.P. Lebedev, А.В. Медведев, Д.С. Кушнарев A.V. Medvedev, D.S. Kushnarev 4 радара показала, что он обладает высокой чувствительностью к угловым размерам объекта порядка 7.5 угл. мин, а потенциальная точность определения углов цели может достигать 40 угл. сек. Ключевые слова: радар некогерентного рассеяния, МСТ-радар, фазированная решетка, потенциал атмосферных радаров, ионосфера, атмосфера. \ Keywords: incoherent scatter radar, MST radar, phased array, atmospheric radar potential, ionosphere, atmosphere. ВВЕДЕНИЕ В рамках проекта Национального гелиогеофизического комплекса планируется создание радара для комплексного исследования атмосферы. Высокая мощность излучения и большая апертура нового радара позволят изучать ионосферу методом некогерентного рассеяния (НР) и нейтральную атмосферу на высотах мезосферы — стратосферы — тропосферы (МСТ). Разработка радара осуществляется Радиотехническим институтом им. акад. Минца (РТИ), имеющим опыт проектирования радиолокационных установок. Проект радара основан на конструкции существующих установок, разработанных РТИ, и модифицируется в соответствии с задачами исследования атмосферы. Необходимо изучить конструктивные особенности нового радара НР-МСТ и оценить его научный потенциал. Основной задачей радара НР-МСТ является исследование структуры и динамики нейтральной атмосферы и ионосферы, а также взаимодействий между ними. Наблюдения нейтральной атмосферы включают измерение полного вектора скорости ветра, параметров турбулентности и атмосферных слоев на высотах тропосферы, стратосферы и мезосферы. Особый интерес представляют мезосферные эхосигналы [Rapp, Lubken, 2004], которые появляются в летнее и зимнее время, иногда сопровождаясь видимыми серебристыми облаками [Romejko et al., 2003]. МСТ-измерения включают также исследования метеоритов на высотах мезосферы и нижней термосферы, рассеяния на тепловых флуктуациях плазмы, атмосферных гравитационных волн и отражающих атмосферных слоев [Hocking, 2011]. Метод НР позволяет получать параметры ионосферной плазмы: электронную концентрацию, температуры ионов и электронов, ионный состав, скорость дрейфа плазмы — на высотах 100–2000 км и наблюдать ионосферные неоднородности, когерентные эхосигналы, D-слой ионосферы. Таким образом, радар позволит исследовать физические процессы и связи в системе «нейтральная атмосфера — ионосфера — плазмосфера». В задачи радара входят также интерферометрические наблюдения тонкой структуры атмосферных неоднородностей и сосредоточенных целей (спутников, космического мусора и метеоритов), пассивные радиоастрономические наблюдения Солнца и космических радиоисточников. Национальный гелиогеофизический комплекс включает ряд оптических инструментов, которые будут размещаться рядом с радаром НР-МСТ и позволят расширить его диагностические возможности. Так, наблюдения нижней и средней атмосферы совместно с лидарными измерениями позволяют исследовать процессы переноса [Bertin et al., 2001]. Всего в мире насчитывается около десяти радаров НР и двадцать крупных МСТ-радаров. За последние пять лет введены в строй такие инструменты, как МСТ-радар MAASRY в Норвегии [Latteck et al., 2012] и МСТ/НР-радар PANSY в Антарктиде [Sato et al., 2014]. Планируется создание многофункционального радара EISCAT-3D [McCrea et al., 2015]. По своим техническим характеристикам радар НРМСТ не должен уступать современным установкам мирового уровня. В настоящее время планируется разместить радар вблизи о. Байкал, что позволит исследовать атмосферу и ионосферу над Восточной Сибирью — регионом, где раньше не было диагностического инструмента со столь широким высотным охватом, и дополнить мировую сеть радаров. Кроме того, это даст возможность осуществлять экологический мониторинг атмосферы вблизи озера. В статье рассматриваются конфигурация антенны и характеристики приемопередающей системы радара, оцениваются рабочий диапазон высот, потенциал при работе в режимах МСТ- и НР-наблюдений и возможности интерферометрических и радиоастрономических исследований. КОНФИГУРАЦИЯ И ХАРАКТЕРИСТИКИ РАДАРА Антенная система радара НР-МСТ состоит из двух фазированных решеток (антенных полотен), расположенных вдоль магнитного меридиана на удалении 100 м друг от друга и наклоненных в северном и южном направлении на 20° (рис. 1). Каждое антенное полотно имеет размеры 40×40 м и включает 1536 крест-вибраторов, расположенных в узлах треугольной сетки. Сектор обзора каждого антенного полотна составляет ±45°, а общий сектор обзора радара равен ±65° в направлении север — юг и ±45° в направлении восток — запад относительно магнитного меридиана. Следует учесть, что при качании луча фазированной решетки происходит расширение луча и уменьшение коэффициента усиления, поэтому для данной конфигурации радара потенциал повышен в направлении 20° относительно горизонта и понижен при вертикальном направлении луча. В приемной системе обрабатываются сигналы двух поляризаций. В режиме излучения сигнал имеет линейную поляризацию, что позволит измерять абсолютную электронную концентрацию на основе метода фарадеевских замираний [Shpynev, 2004].
Перспективный радар НР-МСТ: потенциал и… Prospective IS-MST radar. Potential and diagnostic capabilities 5 Рис. 1. Радар НР-МСТ (вверху) и диаграмма направленности решетки. Наклон решетки скомпенсирован качанием луча на 20° для обеспечения вертикальных измерений Местоположение радара выбиралось таким образом, чтобы уменьшить радиогоризонт — расстояние до ближайших препятствий, способных создать местный помеховый отраженный сигнал. Это необходимо, чтобы уменьшить минимальное расстояние, на котором возможно измерение характеристик среды. Остаточные местные помехи будут компенсироваться специальной обработкой сигналов. В табл. 1 представлены основные характеристики радара НР-МСТ. Частота повторения импульсов PRF, скважность q и длительность импульса τи определяются выбранными режимами работы: МСТ — в высотном интервале 5–100 км, НР-1 — 100–600 или НР-2 — 600–2000 км. Данные высотные интервалы приведены в качестве примера, а реальные будут задаваться в соответствии с научной задачей. Так, при проведении измерений в нижней и средней атмосфере нет необходимости в высокой частоте повторения импульсов, поскольку ширина доплеровского спектра МСТ-сигналов не превышает 50 Гц, что позволяет работать по нескольким направлениям одновременно и проводить дополнительное когерентное накопление сигнала для увеличения соотношения сигнал/шум. Каждое антенное полотно разбито на 12 фрагментов, состоящих из 128 крест-вибраторов и подключенных к отдельным приемопередающим модулям. В модуле происходит частичное диаграммообразование для сигналов, принятых с 128 антенн, в результате чего Таблица 1 Характеристики радара НР-МСТ Характеристика Обозначение Значение диапазон рабочих частот f 154–162.7 МГц конфигурация радара – 2 квадратных решетки 40×40 м тип антенны – Крест-вибратор число излучателей Nt 3072 пиковая мощность Pt 2 МВт эффективная площадь антенны (2 решетки) Ae 2500 м2 средний уровень боковых лепестков – –35 дБ шумовая температура Te ~500 К ширина луча θ0.5 2.5×2.5° (1 решетка) ширина полосы приемника ΔB 625 кГц длина дискрета широкополосного импульса τд 3.2 мкс минимальное разрешение по высоте Δr 480 м частота повторения импульсов (режим МСТ/НР-1/НР-2) PRF 1333.3/250/75 Гц длительность импульса (режим МСТ/НР-1/НР-2) τи 51.2/700/6500 мкс скважность (режим МСТ/НР-1/НР-2) q 6.8/17.5/37.5 % произведение средней мощности на эффективную площадь Pср · Ae 3.4·108/8.75·108/18.75·108 Вт·м2 сектор обзора – ±45° получаются четыре выходных приемных канала (по две поляризации и квадратуры), доступных для дальнейшей обработки. Таким образом, каждая антенная решетка радара позволит сформировать один узкий луч 2.5×2.5° на излучение и 12 широких независимых приемных лучей (в 12 раз шире). Приемные лучи формируются фазовыми центрами, разнесенными по поверхности антенного полотна, а сами полотна находятся на расстоянии 100 м друг от друга. Это даст возможность комбинировать сигналы с отдельных приемных каналов для повышения отношения сигнал/шум, получения пучка узких приемных лучей и проведения интерферометрических наблюдений. Излучаемый сигнал формируется в цифровом синтезаторе, переносится на промежуточную частоту и переводится в аналоговую форму, после чего разветвляется на все передатчики, где фазируется в соответ
А.П. Потехин, А.Г. Сетов, В.П. Лебедев, A.P. Potekhin, A.G. Setov, V.P. Lebedev, А.В. Медведев, Д.С. Кушнарев A.V. Medvedev, D.S. Kushnarev 6 ствии с требуемым направлением. Фазирование в передающей системе осуществляется 4-битовыми аналоговыми фазовращателями с дискретом фазы 22.5°. Выбранный шаг квантования фазы аналогового фазовращателя приводит к снижению коэффициента усиления на 0.06 дБ и повышению уровня боковых лепестков. Сигналы синхронизации также формируются единым гетеродином и разветвляются на все приемопередающие модули, поэтому радар когерентен по фазе относительно передаваемых и принимаемых радиоимпульсов. Приемная система, в отличие от передающей, является полностью цифровой: сигналы каждой поляризации с антенн усиливаются в малошумящем усилителе, фильтруются и оцифровываются. Это обеспечивает широкий динамический диапазон (менее 83 дБ), позволяющий принимать как слабые сигналы некогерентного рассеяния, так и мощные когерентные сигналы от нейтральной и ионизированной атмосферы. Выбор динамического диапазона определяет также минимальную рабочую высоту, поскольку мощные местные помехи приводят к насыщению в приемном тракте, что препятствует их фильтрации. Принятые цифровые сигналы разделяются на квадратуры, фазируются и частично суммируются, образуя 12 каналов с каждой решетки в двух поляризациях и двух квадратурах (всего 96 цифровых каналов). Кроме того, приемопередающая система позволит задавать произвольное амплитудно-фазовое распределение в антенном полотне для сигналов с подрешеток из 16 крест-вибраторов, поэтому возможно получение диаграмм направленности (ДН) с малым уровнем боковых лепестков или нулевым в выборочном направлении излучения и приема. На рис. 1 представлена ДН для одного антенного полотна (1536 элементов) в плоскости земли. ДН построена по классическому методу, в основе которого лежит представление всего антенного поля одним фазовым центром, что позволяет уменьшить время вычисления. При этом не учитывается взаимное влияние между антеннами, которое приводит к искажению ДН отдельных элементов решетки. Для компенсации физического наклона решетки элементы полотна сфазированы таким образом, чтобы главный лепесток был направлен вертикально. ДН обладает несимметричной изогнутой формой из-за треугольного расположения элементов и преобразования координат к плоскости земли. В дальнейшем будет проводиться более точное моделирование ДН, учитывающее поле отдельных антенн и взаимодействие между ними. Характеристики современных НР-радаров можно найти, например, в [Медведев, 2014]. В большинстве этих радаров, построенных еще до 2000 г., используются параболические антенны с коэффициентами усиления больше 40 дБ и ламповые передатчики мощностью 1–5 МВт. Некоторые МСТ-радары перечислены в [Hocking, 1997]. К ним стоит добавить такие современные радары, как EAR, MAARSY и PANSY. В отличие от НР-радаров, практически все они представляют собой многолучевые фазированные решетки, что, скорее всего, обусловлено менее жесткими требованиями к мощности излучаемого сигнала. Благодаря развитию и удешевлению цифровых интегральных устройств, преимущество фазированных решеток, позволяющих быстро изменять направление луча, исследовать несколько направлений одновременно и менять форму ДН различных амплитудно-фазовых распределений, становится очевидным. Радар НР-МСТ, имеющий пиковую мощность 2 МВт, апертуру 80×40 м и цифровую антенную систему, позволит проводить наблюдение атмосферы на уровне мировых установок. В настоящее время для исследования ионосферы в Институте солнечно-земной физики СО РАН используется Иркутский радар НР [Жеребцов и др., 2002; Потехин и др., 2008], основные недостатки которого — невозможность проведения измерений на высотах до 150–200 км и фиксированный период повторения импульсов (24.4 Гц). Место расположения нового радара выбиралось с учетом требования малого радиогоризонта, поэтому он не будет обладать подобными недостатками. Уникальность установки также состоит в разнесении двух антенных полотен на расстояние 100 м и наклоне вдоль магнитного меридиана. Такое положение решеток позволит проводить интерферометрические наблюдения и исследовать вытянутые вдоль магнитного поля неоднородности. Стоит отметить, что в настоящее время в мире существуют всего две установки, работающие на частоте 158 МГц: Харьковский и Иркутский радары, построенные по схожей схеме. ПОТЕНЦИАЛ РАДАРА ПРИ ИССЛЕДОВАНИИ НИЖНЕЙ, СРЕДНЕЙ И ВЕРХНЕЙ АТМОСФЕРЫ Приоритетными задачами радара являются диагностика динамических процессов в верхней, средней и нижней атмосфере, исследование вариаций параметров ионосферной плазмы и изучение взаимодействий в системе «атмосфера — ионосфера — плазмосфера». Следовательно, необходимо, чтобы радар охватывал широкий высотный диапазон. При этом параметры среды по мере роста высоты значительно изменяются и для диагностики атмосферы должны использоваться различные методы. Основные источники эхо-сигналов в этом диапазоне — некогерентное рассеяние на тепловых флуктуациях ионосферной плазмы и брегговское когерентное рассеяние на турбулентных флуктуациях в области мезосферы — стратосферы — тропосферы, иногда также называемое когерентным рассеянием в нейтральной атмосфере. Режимы работы радара НР-МСТ показаны на рис. 2. Разделение радара на два разнесенных антенных полотна позволяет ускорить процесс сканирования атмосферы в разных направлениях в случаях, когда мощность сигналов обратного рассеяния велика (например, при зондировании высотного интервала 100–600 км в режиме НР-1). Выше 600 км в ионосфере (режим НР-2) требуется одновременная работа двух решеток в одном направлении для повы-
Перспективный радар НР-МСТ: потенциал и… Prospective IS-MST radar. Potential and diagnostic capabilities 7 Рис. 2. Схематичное изображение режимов работы радара НР-МСТ. Сверху вниз: режимы МСТ, НР-1, НР-2 вышения потенциала радара. Радар НР-МСТ имеет линейную поляризацию излучаемого сигнала, что позволит получать профиль электронной концентрации в ионосфере как на основе фарадеевского вращения [Shpynev, 2004], так и по абсолютным измерениям рассеянной мощности. Режимы МСТ и НР-1 позволят исследовать несколько направлений с временным разрешением, достаточным для определения полного вектора скорости движения атмосферы, ионосферного дрейфа и нейтрального ветра, а также для детального изучения атмосферных неоднородностей.
А.П. Потехин, А.Г. Сетов, В.П. Лебедев, A.P. Potekhin, A.G. Setov, V.P. Lebedev, А.В. Медведев, Д.С. Кушнарев A.V. Medvedev, D.S. Kushnarev 8 В нижней и средней атмосфере для каждого излученного импульса в каждом высотном интервале, равном разрешению по высоте, получают одиночную комплексную выборку для дальнейшей обработки. Сигналы на этих высотах обладают временем когерентности ~0.1 с, поэтому частота следования импульсов может быть небольшой (десятки и сотни герц), но длительность измерений должна быть достаточной для обеспечения необходимого разрешения по частоте (десятки секунд). В целом с ростом высоты ширина спектра сигналов обратного рассеяния увеличивается и сигналы, приходящие от разных импульсов, теряют когерентность. В связи с этим для полноценной обработки данных необходимо чаще получать выборки сигналов, приходящих с одного и того же высотного интервала (повысить частоту следования импульсов). Кроме того, чем дальше находится исследуемая область, тем реже можно посылать импульсы, чтобы не возникало неопределенностей по дальности. При этом разрешение по дальности необходимо поддерживать малым, т. е. использовать короткие импульсы либо применять сложные широкополосные сигналы вплоть до высоты ~100 км, где атмосферные процессы характеризуются меньшими масштабами. Эти масштабы, в свою очередь, зависят от скорости изменения параметров атмосферы с высотой (шкалы высот) и протяженности атмосферных неоднородностей и волновых возмущений в плоскости, перпендикулярной лучу радара. Таким образом, на высотах слоев D и E необходимо использовать специальные многоимпульсные режимы передачи [Farley, 1972] из-за жестких требований к разрешению по частоте и дальности. Начиная с высот 100–200 км требования к ширине импульса смягчаются, сигнал становится более некогерентным и для каждого высотного интервала возможно получение необходимого количества выборок за время следования одного импульса. Сложности, связанные с обработкой сигналов, могут усугубляться состоянием самой среды. На некоторых высотах возмущенность атмосферы становится недостаточной для получения эхо-сигналов и оценки параметров среды, что усложняет задачу исследования связей между разными атмосферными слоями и делает ее зависимой от вариаций среды. Большинство МСТ-радаров работает в диапазоне частот 45–55 МГц [Fukao et al., 1990; Hocking, 1997; Latteck et al., 2012; Sato at el., 2014]. Диапазон исследуемых высот 2–30 км соответствует регулярным наблюдениям тропосферы и нижней стратосферы, а в диапазоне 60–95 км ведутся наблюдения мезосферы, в том числе мезосферных эхо-сигналов. Из-за затухания малых турбулентных неоднородностей в мезосфере мощность рассеянного сигнала падает с ростом рабочей частоты радара, что может повлиять на частоту наблюдения мезосферных эхо-сигналов. Тем не менее мезосферные эхо-сигналы, хоть и значительно ослабленные, наблюдались даже на радарах с частотой 930 МГц [Rapp et al., 2008]. Рабочая частота радара НРМСТ составляет ~158 МГц, поэтому необходимо оценить зависимость мощности сигналов, рассеянных на турбулентных неоднородностях, от частоты. Методика оценки потенциала Диапазон высот, которые позволяет исследовать радар, можно оценить, сравнивая мощность сигнала обратного рассеяния с чувствительностью радара. Удобнее использовать удельную эффективную площадь рассеяния (ЭПР) единичного объема среды σср [1/м], которая зависит только от параметров атмосферы. Тогда чувствительность радара определяется минимальной удельной ЭПР σмин, которой должна обладать атмосфера для получения заданной статистической достоверности оценки параметров среды. Если σср/σмин>1, то проведение измерений на данной высоте возможно. Параметр σмин находится из основного радиолокационного уравнения для случая объемного рассеяния [Bowles et al., 1962]: 2 2 t r 3 2 , , sin 128 PL c P G d d R (1) где Pt — пиковая излучаемая мощность, L — потери в системе радара и при распространении сигнала в атмосфере, τ — длительность зондирующего импульса, R — расстояние до центра рассеивающего объема, G(θ, φ) — коэффициент усиления антенны в направлении (θ, φ). Квадрату коэффициента усиления G2, учитывающему прием и излучение одной и той же антенной системой, соответствует двунаправленная ДН шириной 1.8×1.8° для одного сегмента радара. В данном случае σ определяет долю изотропно рассеиваемой мощности на единицу рассеивающего объема. Принятую мощность Pr выразим через отношение сигнал/шум S/N: r / , P kT BS N (2) где k — постоянная Больцмана, T — шумовая температура радара, ΔB — полоса рабочих частот приемника. Предполагается, что в приемнике присутствует только белый тепловой шум радара и отсутствуют местные помехи (считается, что они исключены после обработки данных). С учетом того, что шумовые характеристики известны, отношение S/N определяет уровень сигнала, требуемый для достижения статистической достоверности измерений, например для ограничения дисперсии оценки параметров среды в заданных пределах. Интеграл в (1), учитывающий форму ДН, можно аппроксимировать гауссовой кривой для случая антенной решетки с квадратной апертурой [Hagen, Baumgartner, 1996]: 2 e 2 2 . 4 , sin 0.44 A G d d (3) Подставляя (2) и (3) в (1), получим 2 мин 3 t e 18.2 , cos S N R kT B PL cA (4) где cos3 θ учитывает потери при качании луча (cos θ из-за уменьшения проекции апертуры и cos2 θ как оценка ДН одного крест-вибратора). На данном этапе создания радара НР-МСТ потери L не известны. Обычно они уточняются на этапе калибровки радара. В дальнейших вычислениях примем L=0.25.
Перспективный радар НР-МСТ: потенциал и… Prospective IS-MST radar. Potential and diagnostic capabilities 9 В случае рассеяния на турбулентных неоднородностях σср вычисляется по формуле [Lübken, 2014]. Для расчета атмосферных параметров используются модель нейтральной атмосферы NRLMSIS-00 и модель ионосферы IRI-2012. Интенсивность турбулентности, определяемая скоростью диссипации турбулентной энергии ε, задана равной 10–3 м2с–3 в тропосфере, 10–5 м2с–3 в стратосфере и 10–1 м2с–3 в мезосфере, что соответствует умеренным возмущениям [Hocking, 1985; Watkins et al., 1988; Kantha, Hocking, 2011; Gavrilov, 2013]. Предполагается, что параметры нейтральной атмосферы будут оцениваться по задержкам автокорреляционной функции. Основным измеряемым параметром в нижней и средней атмосфере является скорость ветра, поэтому отношение сигнал/шум в (4) задается таким образом, чтобы стандартное отклонение при определении скорости не превышало 1 м/с в стратосфере — тропосфере и 2 м/с в мезосфере. Выражения для дисперсии оценки скорости по фазе автокорреляционной функции взяты из [Zrnic, 1979]. Потенциал при исследовании тропосферы — стратосферы На МСТ-радарах проводятся регулярные измерения параметров сигналов, рассеянных на турбулентных флуктуациях среды [Hocking, 2011]. Турбулентное перемешивание нейтральной составляющей в стратосфере и тропосфере приводит к появлению градиента показателя преломления. Рассеяние происходит на флуктуациях показателя преломления с характерными масштабами, равными половине длины падающей волны (условие Брегга). Удельная ЭПР σср флуктуаций показателя преломления, вызванных турбулентными неоднородностями в атмосфере, зависит от рабочей частоты радара, формы спектра турбулентных флуктуаций и параметров атмосферы [Hocking, 1985]. Энергетический спектр турбулентностей делится на длинноволновый интервал, где происходит основная генерация энергии, интервал инерции, где энергия медленно убывает, и интервал вязкости, где турбулентности резко затухают, а энергия преобразуется в тепло. Высота перехода между интервалами инерции и вязкости зависит от частоты. Так, для радаров метрового диапазона переход наступает на высотах более 40 км [Hocking, 1985]. В средней и верхней стратосфере, где параметры атмосферы меняются медленно и интенсивность турбулентностей мала, мощность рассеянного сигнала падает, что приводит к исчезновению эхо-сигналов на высотах 30–60 км. Важная особенность измерений в нижней и средней атмосфере в вертикальном направлении — наличие мощных когерентных отражений от горизонтальных протяженных слоев. Эти эхо-сигналы обладают узким спектром шириной менее 1 Гц и чувствительностью к углу наклона сканирующего луча. Далее будет оценен потенциал радара при рассеянии исключительно на турбулентных флуктуациях. На рис. 3 показано отношение σср/σмин в тропосфере и стратосфере для трех рабочих частот: 50 МГц (стандартная частота МСТ-радаров), 158 МГц (частота радара НР-МСТ) и 233 МГц (частота радара EISCAT VHF). Последняя частота выбрана в связи с тем, что на радаре EISCAT VHF, который обладает сравнительно низкой для НР-радаров рабочей частотой (233 МГц), также проводились МСТ-наблюдения [Strelnikova, Rapp, 2010]. Поскольку для режима работы в нижней атмосфере нет особых требований к мощности излучения, на рис. 3 указано также отношение σср/σмин при работе с одним антенным полотном, что позволит улучшить временное разрешение. При моделировании были заданы следующие характеристики: τи=7 мкс, PRF=1000 Гц, длительность накопления 20 с, ширина доплеровского спектра 5 Гц. При моделировании изменялись только частота радара и эффективная температура, но не размер апертуры. Видно, что при оценке скорости ветра в тропосфере и нижней стратосфере повышение частоты позволяет увеличивать потенциал (прежде всего — разрешение по времени) до тех пор, пока рассеяние происходит на турбулентных флуктуациях в инерционном интервале спектра. Это связано с тем, что при повышении частоты увеличивается доплеровский сдвиг и его проще различить на фоне помех. После перехода в интервал вязкости (с ростом высоты), который происходит раньше для меньших длин волн, интенсивность рассеяния резко падает. Чтобы увеличить потолок измерения скорости ветра, необходимо существенно улучшить характеристики радара и перейти на более низкие частоты [Hocking, 1997]. В заданных условиях радар НР-МСТ (158 МГц) позволит измерять скорости ветра до ~26 км со стандартным отклонением не более 1 м/с. Потенциал при исследовании мезосферы В мезосфере на высоте D-слоя ионосферы сигнал дополнительно усиливается за счет повышения электронной концентрации. В то же время с ростом высоты кинематическая вязкость атмосферы увеличивается, что приводит к затуханию турбулентностей малого размера в мезосфере и переходу в интервал вязкости и, следовательно, к ослаблению рассеянного сигнала [Hocking, 1985]. Тем не менее на высотах ~85 км в летнее время и ~70 км в зимнее время наблюдаются мезосферные летние (MSE, mesospheric Рис. 3. Отношение σср/σмин в стратосфере для разных частот и размеров используемого антенного полотна
А.П. Потехин, А.Г. Сетов, В.П. Лебедев, A.P. Potekhin, A.G. Setov, V.P. Lebedev, А.В. Медведев, Д.С. Кушнарев A.V. Medvedev, D.S. Kushnarev 10 summer echo) и зимние (MWE, mesospheric winter echo) эхо-сигналы. Эхо регистрируется в широком диапазоне от средних до сверхвысокочастотных волн [Терещенко и др., 2002; Rapp et al., 2008]. Несмотря на то, что мезосферные эхо-сигналы в основном наблюдались в полярном регионе, регулярные наблюдения их проводились также и в средних широтах. В работе [Zeller et al., 2006] показано, что для среднеширотного радара частота появления эхо-сигналов составила 7 и 0.3 % от общего времени наблюдения для MSE и MWE соответственно. МСТ-радары позволяют также исследовать рассеяние на следах метеоритов на высотах 80–100 км, при этом можно определить температуру и скорость нейтрального ветра. В интервале вязкости происходит затухание турбулентностей, однако существует механизм, который уменьшает ослабление рассеянного сигнала и, как полагают, приводит к появлению MSE [Lübken et al., 1998]. Летом на высотах ~85 км образуются ионизированные ледяные частицы, которые уменьшают диффузию электронов, что приводит к усилению флуктуаций электронной плотности на малых пространственных масштабах [Rapp, Lubken, 2004]. Степень ослабления диффузии характеризуется числом Шмидта Sc, которое может достигать нескольких тысяч единиц [Strelnikov et al., 2009]. Часто вместе с MSE на высотах ~82 км возникают серебристые облака [von Zahn, Bremer, 1999], видимые невооруженным глазом [Терещенко и др., 2010]. На высоте образования MWE (~70 км) кинематическая вязкость среды еще достаточно низкая и образование эха не требует наличия ледяных частиц. Считается, что MWE появляются, когда электронная концентрация и турбулентность среды повышены, что подтверждается комплексными измерениями атмосферных параметров [Lübken et al., 2007]. Использованная ранее формула для определения удельной ЭПР σср [Lübken, 2014] позволяет варьировать число Шмидта Sc. Для моделирования MSE задан слой с повышенным значением Sc как необходимым условием появления эха. На рис. 4 показано отношение σср/σмин для трех частот и значений числа Шмидта. Стандартное отклонение оценки для мезосферы задано 2 м/с в связи с тем, что наблюдаемые скорости ветра в мезосфере выше, чем в тропосфере—стратосфере. При моделировании выбрано время суток с высокой электронной концентрацией в ионосфере, так как она влияет на мощность эхо-сигналов [Lübken, 2014]. Число Шмидта плавно возрастает от 1 на высоте 82 км до заданного значения на высоте 85 км. В отличие от тропосферы— стратосферы, в мезосфере затухание турбулентностей малых размеров становится значительным и радары с меньшими рабочими частотами обладают преимуществом независимо от оцениваемого параметра. Видно, что радар НР-МСТ позволит наблюдать MSE при повышенном значении Sc. Выше 87 км затухание становится значительным и уровень сигнала резко убывает. Появление MWE не сопровождается увеличением Sc, однако, как правило, имеет место повышенная Рис. 4. Отношение σср/σмин при наблюдении летнего мезосферного эха для разных частот и чисел Шмидта Рис. 5. Отношение σср/σмин при наблюдении зимнего мезосферного эха для разных частот и уровней электронной концентрации. Вертикальной линией отмечен порог чувствительности радара электронная концентрация. На рис. 5 показано отношение σср/σмин при наблюдении MWE и кратном увеличении электронной концентрации относительно модельного значения IRI. Модель IRI не дает значения ниже 65 км, что отражено на рисунке. Черной вертикальной линией показана граница превышения сигналом чувствительности радара (σср/σмин=1). Видно, что незначительное эхо на частоте 158 МГц возможно при повышении на порядок концентрации в D-слое ионосферы. Частота 233 МГц не позволяет измерить скорость нейтрального ветра в слое MWE в заданных условиях. Некогерентное рассеяние Метод НР является самым информативным средством анализа верхней атмосферы Земли. Он позволяет измерять параметры ионосферной плазмы: электронную концентрацию, температуру ионов и электронов, скорость дрейфа плазмы, ионный состав и частоту столкновения с нейтральной составляющей в D-слое. Удельная эффективная площадь рассеяния для НР-измерений имеет вид e e НР 2 2 e i , 1 1 N T T (5) где σe≈10–28 м2 — эффективная площадь рассеяния электрона, α=4πD/λ, D — радиус Дебая, Te — тем
Перспективный радар НР-МСТ: потенциал и… Prospective IS-MST radar. Potential and diagnostic capabilities 11 пература электронов, Ti — температура ионов, Ne — электронная концентрация. Параметры ионосферной плазмы оцениваются с помощью многопараметрической подгонки автокорреляционной функции или спектра принятого сигнала. Для анализа дисперсии оценок получаемых параметров удобно воспользоваться эмпирическим приближением, представленным в [Murdin, 1979] для электронной концентрации, температуры и скорости дрейфа плазмы. Параметры ионосферы значительно меняются с ростом высоты, что существенно влияет на спектр сигнала обратного рассеяния и характерные масштабы протекающих процессов (шкалу высот, длину приливных волн). Соответственно должен меняться режим работы радара при исследовании разных высот. Существуют два базовых режима для наблюдения ионосферы на высотах от 100 км: НР-1 и НР-2 (см. рис. 2). Первый режим предназначен для работы на высотах 100–600 км и характеризуется следующими параметрами: длительность излучения τи=700 мкс, PRF=250 Гц. Кроме того, режим НР-1 предполагает использование только одного сегмента радара. В режиме НР-2 для высот от 600 км используются более протяженные импульсы τи=6500 мкс, а частота их повторения снижается до PRF=75 Гц. На рис. 6 показано отношение σср/σмин для режимов НР-1 и НР-2 при оценке электронной концентрации и температуры с относительным стандартным отклонением, не превышающим 10 %, и скорости дрейфа плазмы со стандартным отклонением не более 45 м/с. Длительность накопления выбрана равной 10 мин. Моделирование осуществлялось для летнего полудня, когда электронная концентрация достаточно высока. Видно, что оценить скорость дрейфа плазмы, даже с высоким статистическим отклонением, становится невозможно на меньшей высоте, чем концентрацию и температуру. На границе смены режима работы (600 км) длительность импульса повышается, но снижается частота повторения импульсов, поэтому отношение σср/σмин повышается незначительно. Выше 600 км отношение σср/σмин мало даже в предельном режиме работы, поэтому если реальный потенциал радара окажется меньше из-за неучтенных потерь в трактах или системе обработки, потолок определения параметров плазмы может опуститься до 800–1000 км. Тем не менее в возмущенных условиях возможно снижение времени накопления и, следовательно, улучшение временного разрешения, что позволит детально исследовать процессы в ионосфере. Для сравнительной оценки потенциала на рис. 7 показано отношение σср/σмин, построенное по той же методике для ИРНР, который используется для исследования ионосферы более десяти лет [Жеребцов и др., 2002; Медведев, 2014]. ИРНР имеет частотный принцип сканирования и, по сравнению с радаром НР-МСТ, меньшую эффективную площадь антенны, принимает сигналы обратного рассеяния только в одной поляризации, излучает импульсы длиной не более 900 мкс с частотой повторения импульсов 24.4 Гц. Данные ограничения не позволяют эффективно изучать ионосферу на высотах более 600 км. Более того, окружающая помеховая обста- Рис. 6. Отношение σср/σмин при оценке электронной концентрации, температуры и скорости дрейфа плазмы. Число антенн режимов НР-1 и НР-2 определяет работу одним (1536 антенн) или двумя (3072 антенны) сегментами радара соответственно Рис. 7. Отношение σср/σмин для ИРНР новка ограничивает нижний диапазон исследуемых высот, в результате чего невозможно наблюдение ионосферы и атмосферы на высотах до 100 км. Видно, что граница измерения скоростей находится намного ниже, чем показано на рис. 6, что обусловлено малой частотой повторения импульсов ИРНР и, соответственно, малой длиной выборки при одном и том же времени накопления (10 мин). При сравнении рис. 6 и 7 следует учитывать, что график на рис. 6 (режим