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

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

Рубрика Астрономия и космонавтика
Вид статья
Язык русский
Дата добавления 16.11.2018
Размер файла 134,0 K

Отправить свою хорошую работу в базу знаний просто. Используйте форму, расположенную ниже

Студенты, аспиранты, молодые ученые, использующие базу знаний в своей учебе и работе, будут вам очень благодарны.

Размещено на http://www.allbest.ru/

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

Диденко А.В., Усольцева Л.А.

Резюме

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

Введение

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

Основное содержание исследования

атмосфера геостационарный наземный космический

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

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

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

Состав сухого воздуха в отношении его основных компонентов остается постоянным до значительных высот. Его основными составляющими (более чем на 99.0%) являются молекулы азота и кислорода; кроме постоянно входящих в его состав молекул аргона и углекислого газа в количестве, немного меньшем 1%, на долю всех остальных газовых примесей приходится меньше 0,005% по объему.

В содержании азота, кислорода и аргона не отмечается никаких закономерностей и сколько-нибудь заметных изменений в зависимости от времени и географического положения мест. Некоторым изменениям подвержено содержание в атмосфере углекислого газа, в среднем принятое равным 0.033%. В дневное время его содержание меньше, чем в ночное; оно несколько различно в разное время года - больше зимой и меньше летом и осенью. Содержание СО2 в полярных странах меньше: в Антарктике, например, оно составляет около 0.02%. Оно меньше над океанами и больше над сушей. С высотой, по крайней мере до 20 км, среднее содержание СО2 изменяется мало и лишь выше оно медленно убывает. Имеющиеся данные говорят о непрерывном увеличении среднего содержания СО2 в атмосфере.

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

Содержание озона в нижних слоях атмосферы до высоты 10 км мало, после 10 км начинается его увеличение, отчетливо выраженное выше 15 км. На высоте 20-25 км наблюдается максимум содержания озона, выше его количество постепенно убывает и становится ничтожным на 60-70 км.

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

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

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

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

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

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

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

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

В Астрофизическом институте АН-МН РК задачи подобного типа были поставлены впервые академиком В.Г. Фесенковым [3-5] вслед за проблемами получения информации о земной атмосфере по данным фотометрических наблюдений полных лунных затмений. Вопросы применения результатов измерения освещенности затененного диска Луны для исследования распределения озона в верхних слоях атмосферы Земли неоднократно рассматривались в научной литературе, в том числе, выполнялись расчеты ослабления солнечной радиации при тех или иных упрощающих предположениях [6-10].

Из-за большой разряженности атмосферы Луны и оптической рефракции, увеличивающей почти в 15 раз изображение земной атмосферы, видимой на поверхности Луны, этот метод неоднократно успешно применялся, но его использование имеет ряд ограничений. Во-первых, полные лунные затмения бывают достаточно редко: максимум 2 или 3 раза в год. Во-вторых, отражательная способность поверхности Луны недостаточно известна и может меняться при переходе от одного места к другому. Кроме того, для проведения успешного эксперимента необходима непрерывная запись интенсивности радиации, отраженной от нескольких выбранных мест поверхности Луны во время затмения в течение 5 - 6 часов.

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

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

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

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

Точное выражение для яркости первичных сумерек имеет вид:

(1),

где I - общий световой поток от Солнца, соответствующий телесному углу , k () - индикатриса рассеяния в функции угла рассеяния , (h) - рассеивающая способность атмосферы на высоте h, j - функция распределения светового потока от соответствующих горизонтальных сечений солнечного диска, причем угол относится к центру Солнца, а представляет собой изменение этого угла в пределах солнечного радиуса 0, оз и - оптические толщи атмосферы для озона и аэрозолей, соответствующие лучам, которые пронизывают всю атмосферу и проходят на минимальной высоте от земной поверхности h0. Для данного состояния атмосферы оз и - функции только параметра h0. Точное рассмотрение всякой проблемы, связанной с освещением горизонтальными лучами каких-либо элементов на достаточных высотах над земной поверхностью, требует еще учета фактора f рефракционной дисперсии. Для наблюдателя, находящегося на большой высоте и наблюдающего постепенный заход Солнца, его изображение из-за рефракционной дисперсии теряет свою круговую форму и сжимается по вертикали. Это равносильно тому, что каждый освещенный элемент атмосферы и, в частности, искусственный спутник, входящий в земную тень, оказывается в поле лучей, расходящихся вследствие различия преломления на разных высотах. Для элементов атмосферы, вносящих вклад в сумеречные явления, этот фактор изменяется примерно от 0.5 до 1.0 и обычно не учитывается. Однако для искусственных спутников, обращающихся вокруг Земли на больших высотах, его величина меняется почти на целый порядок и должна приниматься во внимание.

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

. (2)

Наблюдения захода в тень будут носить относительный характер, так как альбедо поверхности спутника на момент наблюдений предполагается неизвестным. Фактически регистрируется только изменение яркости ИСЗ, начиная от момента его входа в тень и до полного исчезновения. Задача заключается в том, чтобы определить относительную яркость спутника при различных моделях распределения озона и атмосферных аэрозолей, т.е. в зависимости от хода изменения оз и . Для этого надо определить параметр h0 для каждого положения ИСЗ (рис.1). Кроме того, на каждый момент наблюдений необходимо знать радиус-вектор или высоту над земной поверхностью, угол между направлениями на центр солнечного диска и радиусом - вектором r, т.е. необходимо иметь экваториальные координаты Солнца и ИСЗ с точной привязкой по времени. Эти данные можно получить только из специальных астрометрических наблюдений ИСЗ.

Рис.1. Траектория светового луча (полярные координаты ), проходящего на минимальной высоте h0 от земной поверхности.

Из соотношений (1) и (2) следует, что монохроматическая освещенность ИСЗ при некотором геоцентрическом фазовом угле имеет вид [11]:

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

Величину I (r) можно найти по формуле Миннарта [12]:

, (4)

где I0 () - яркость центра солнечного диска в данной длине волны ; b () - константа, определяемая опытным путем для каждой ; r - радиус солнечного диска. Она удовлетворительно описывает угловой ход яркости вплоть до 0.95 радиуса Солнца.

С учетом всего сказанного соотношение (3) может быть представлено в следующем виде, наиболее удобном для расчетов на ЭВМ:

(5)

На практике наблюдения обычно производятся в некотором спектральном диапазоне, ограниченном длинами волн 1 и 2. Очевидно, что в этом случае освещенность ИСЗ B (1, 2, ) . Даже из качественных соображений следует, что интегральный коэффициент прозрачности в атмосфере для малых длин волн значительно меньше, чем для больших. Поэтому, если пучок радиации проходит через атмосферу, то в нем короткие волны ослабляются сильнее, чем длинные. В результате изменяется состав пучка и его "оптический центр" смещается в сторону более длинных волн, т.е. увеличивается прозрачность.

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

Предположим, что атмосфера разбита на j сферических поверхностей высотой h1, h2,., hn от поверхности Земли, где

, (6)

а . (7)

Слой, заключенный между поверхностями j и j+1, называется j-ым слоем. Его средняя высота , i-ый луч - это солнечный луч, проходящий через центр i-того слоя (j=i) вертикально над точкой захода Солнца . Высота этого луча - это кратчайшее расстояние, на которое он удален от поверхности Земли: . В любом j-ом слое (j1) i-ый луч делится на два отрезка по обе стороны от вертикали, проходящей через точку захода Солнца. В случае пренебрежения оптической рефракцией эти два отрезка имеют одинаковую длину, равную 1/2 aij. Тот же i-ый луч после отражения от поверхности спутника образует с j-ым слоем (1 j n) отрезок длиной . Из простых геометрических соображений следует, что:

, если i=j; (8)

, если i j;

(9).

Здесь rз - радиус Земли, Z - зенитное расстояние спутника для i-ого луча в j-ом слое. При малых зенитных расстояниях спутника () справедлива плоскопараллельная аппроксимация атмосферы, поэтому можно записать . При этом условии обозначим:

- отношение интенсивностей солнечной радиации, приходящей на верхнюю границу атмосферы в длинах волн 1 и 2;

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

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

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

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

Предположим, что в нулевом слое солнечная радиация не ослабляется. Тогда:

(10).

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

Величину с определяют при калибровке фотометра по стандартным звездам (см. ниже). Для описываемого метода значение с используется в процессе моделирования блеска ИСЗ при варьировании параметров атмосферы.

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

, (11)

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

- общее число ослабляющих частиц, содержащихся в столбе сечением 1 см2 в направлении наземный фотометр - спутник; - концентрация частиц в j-ом слое над наблюдателем. Оно задается выражением:

(12)

Из (10), (11) и (12) следует:

(13).

Рассмотрим более общий случай, когда спутник освещается солнечным лучом i, который претерпел ослабление как до, так и после отражения от поверхности спутника. Отношение интенсивностей, регистрируемое на ИСЗ и Земле , имеет вид:

(14)

и

(15).

Введем следующие обозначения:

(16).

Тогда:

(17).

При i, изменяющемся от 1 до n, для n лучей имеем n соотношений, аналогичных (17), поэтому:

(18).

где [A] - матрица n-го порядка, а {X} и {L} - столбцы векторов. Если атмосферная рефракция не учитывается, то матрица [A] определяется выражением (9).

Элементы вектора {L} получаются из наблюдений. Для наземной станции, ведущей наблюдения, выбирается произвольная, но достаточно обоснованная величина . Используя , (13), (16) и (17) можно последовательно получить Таким образом, соответственно n величинам получается n элементов

Далее, из (18) можно получить n элементов Xj вектора {x}. Согласие величин и необязательно, так как последняя величина характерна для вертикали, проходящей через точку захода Солнца, а первая характеризует ослабление в направлении наблюдатель - ИСЗ.

Задачу легко обобщить на случай, когда в атмосфере вместо одного имеется m ослабляющих компонентов. Пусть концентрация k-того компонента в j-том слое будет и разность поперечных сечений ослабления для длин волн 1 и 2 будет . Тогда вместо (17) и (18) будем иметь:

(19)

и

(20)

Особый интерес представляет случай, когда в атмосфере имеются две компоненты, по-разному ослабляющие 1 и 2. Компоненту 1 мы отождествляем с озоном (с плотностью, примерно равной от суммарной плотности частиц), а за компонент 2 примем нейтральный воздух. Для удобства обозначим где числа в скобках соответствуют озону и воздуху. Из (19) и (20) имеем:

(21)

и

(22).

Основываясь на (22), можно провести следующий эксперимент, для которого необходимо:

1) выбрать две длины волн 1 и 2, приходящие к поверхности Земли с интенсивностью, доступной для измерения, для которых

;

2) предположим, что yi является неизвестной функцией высоты, так что вектор - столбец {Y} заранее определен;

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

4) затем из (22) определить вектор - столбец {X}.

Это первая методика, принятая нами для определения хода атмосферного озона с высотой.

Рассмотрим влияние оптической рефракции. Искривление луча вследствие изменения показателя преломления при изменении плотности атмосферы приведет к тому, что суммарная длина i-ого луча будет больше значения, получаемого из (9). Согласно [15] угол полной рефракции можно представить в виде:

(23)

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

Поскольку - высота перигея луча), то , что позволяет найти tg i.

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

(24)

здесь Вычисления через , и достаточно сложны, т.к. в неявном виде входит в . Однако, задавая различные значения с достаточным шагом, можно найти численные значения правой части равенства (24), а затем использовать их для определения значений интерполяционным путем по исходным параметрам , и . Таким образом, для решения прямой задачи - расчета блеска ИСЗ и их сопоставления с результатами реальных наблюдений - есть все необходимые данные.

В настоящее время опубликовано достаточно большое число работ, посвященных зависимости всех рассматриваемых выше величин от высоты при различных метеоусловиях в разных участках спектра. Соответствующие обобщения можно найти, например, в работах [13, 14].

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

Результаты эксперимента и их анализ. Согласно [16] наиболее подходящим диапазоном длин волн для выполнения поставленной задачи является полоса Шаппюи.

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

Оптимальными парами длин волн для измерений являются 600/700 нм, 600/475 нм, 590/529 нм. В этом случае поглощение излучения другими атмосферными составляющими (кислород, водяной пар и т.д.) можно не учитывать из-за их несущественного влияния.

Искусственные спутники, используемые для наблюдений в данном эксперименте, должны удовлетворять следующим требованиям:

иметь простые геометрические формы и большую площадь отражения;

положение спутника и наклон орбиты должны удовлетворять требованиям эксперимента;

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

Координатная информация по ИСЗ была получена с помощью автоматизированной астрономической телевизионной системы (ААТС-2) [17],

Фотометрическая - с автоматизированной фотометрической системой телескопа Цейсса (D=100см). Для учета атмосферной экстинкции и калибровки фотометрической информации используются звезды каталога [18].

Синхронные наблюдения с борта ОКС "Мир" должны были проводиться на аппаратуре "Фиалка - МВ", в состав которой входит спектрофотометрический блок S2000 Ocean Optics Inc. с управляющим компьютером Toshiba CDS/1.4.

В течение 1997 года получена наземная координатная и фотометрическая информация для выбранных эталонных геостационарных ИСЗ. С 15 февраля по 15 мая 1998 года согласно плану-графику на Высокогорной Ассы-Тургеньской обсерватории проводились фотометрические наблюдения ГСС, получено 19 массивов измерений.

Аппаратура "Фиалка - МВ" на борту КС “Мир” в этот период была неисправна. Поэтому, после дополнительных согласований с представителями РКК “Энергия”, в период с 03.08. по 07.08.98 года на борту были проведены пять сеансов наблюдений; 08.08.98 года в течение часа проводилась калибровка аппаратуры по Луне. Данный период являлся неблагоприятным для наземных наблюдений, но, понимая сложившуюся ситуацию, мы были вынуждены провести синхронные наземные наблюдения двух либрационных КО с большими наклонами (i=15) для того, чтобы “привязать” наземные наблюдения весеннего периода.

В осенне-весенний период 1997-1998 гг. были проведены наземные наблюдения 12 космических аппаратов. Как правило, они проводились в трех спектральных диапазонах. По всем перечисленным КА в течение каждой ночи наблюдений получены измерения блеска для построения полной фазовой кривой (одна из них в качестве примера приведена на рис.2). Полученные результаты позволили при дальнейшей обработке определить коэффициент направленного отражения света для данного КА и использовать его при определении интегральной составляющей озона на соответствующих оптических трассах. Зарегистрированы 25 заходов - выходов КА из тени Земли.

Рис.2. Фазовая кривая блеска объекта Spacenet

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

Большинство доступных для наблюдений оптическими методами КО имеют характерные площади отражающих поверхностей от 1 до 85 метров. Их формы не отличаются особым разнообразием [19]. Как правило, корпус космического аппарата представляет собой цилиндр или куб, различным образом ориентированный в пространстве, на боковых поверхностях которого находятся элементы солнечных батарей (СБ). Они могут быть жестко связаны с корпусом или вынесены на специальных панелях за его пределы. Каждый КА имеет также одну или несколько приемно-передающих антенн различного типа. Следует отметить, что активный КА, как правило, стабилизирован и жестко удерживается в заданной подспутниковой точке, что следует из его функционального назначения.

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

, (25)

где N* - отсчет на объект, исправленный за влияние фона, темнового тока и нелинейность приемной аппаратуры; Е - распределение энергии в спектре объекта за пределами атмосферы; Т - кривая реакции системы телескоп + фотометр с конкретными фильтрами; Р - коэффициент пропускания атмосферы; M - воздушная масса; C - постоянная, определяющая нуль - пункт системы.

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

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

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

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

Определенную сложность представляет учет фона неба, т.к. при движении КO на фоне звезд не исключается возможность их попадания в поле рабочей диафрагмы электрофотометра. Например, при наблюдениях на метровом телескопе звезда 12 величины даст ошибку в 3% при времени накопления 30 сек. Вклад эффекта Форбса невелик. Ошибка, обусловленная им, в полосе V не прeвышает 0.005m, а в полосе B близка к 0.01 m.

Подробная методика влияния перечисленных параметров на блеск КО и их учет изложены в работе [17].

Ряды измерений, получаемые на фотометре, представляют собой последовательность отсчетов в квантах за указанный интервал времени, фиксируется также время начала экспозиции, индекс фильтра, отсчет на фон, количество точек в экспозиции и координатная информация для данного КО. Дальнейшая обработка проводится с помощью программы BVRGSS, созданной в Лаборатории наблюдений ИСЗ АФИФ МН-АН РК. В результате наблюдаемый блеск КО приводится в стандартную фотометрическую систему. При этом ведется учет фона и других вышеперечисленных параметров. Блеск КО вычисляется по сумме квантов за указываемый промежуток, в нашем случае чаще всего за 10 или 5 сек., что соответствует изменению высоты прохождении светового луча от центра Солнца в атмосфере Земли приблизительно на 4 и 2 км соответственно. При этом проводится предварительное сглаживание отсчетов скользящим средним по 15 - 20 точкам. На каждый момент определения блеска КО вычисляются его экваториальные топоцентрические координаты. Для расчета геометрии захода Солнца используются геоцентрические координаты, связь которых с топоцентрическими имеет вид:

(26)

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

После обработки исходной информации мы имеем все необходимые данные для дальнейших расчетов.

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

Рис.3. Выход из тени объекта Радуга 31 в фильтре V, наблюдения 08.09.1997 г.

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

, (26)

где

(27)

и

(28)

Рис.4. Заход в тень объекта Intelsat 6 F-4, наблюдения 13.03.1998 г.

При выполнении расчетов использована общепринятая модель высотного распределения озона [13,14], предполагающая, что верхняя граница озонового слоя находится на высоте 70 км. Для получения {B} необходимо знать вектор {M}. Он зависит от характера первичного рассеяния в атмосфере. При его расчете использована стандартная плотность частиц воздуха [20]. Вычислив {B}, легко найти {X}, элементы которого дают плотность озона в соответствующих слоях. Все расчеты проводились на ПЭВМ IBM Pentium 200 MMX по программе OZON.

Следует отметить, что приведенная схема расчета высотного распределения озона не дает полного представления о действительном распределении озона из-за недостаточной точности результатов, получаемых из единичных измерений. Исходя из предположения устойчивости высотного распределения озона на близких сутках, мы объединили и просуммировали результаты измерений, полученные в пределах 2-3 суток по разным КО на близких оптических трассах. При наблюдении ИСЗ весной в восточной зоне точка захода имела координаты, близкие к долготе 53 в. д. и широте 50 с. ш. При наблюдении весной в западной зоне точка захода имела координаты, близкие к долготе 80 в. д. и широте 50 с. ш. Полученные результаты с большой достоверностью характеризуют распределение озона по высоте (рис.5). Ход кривой распределения приводится в относительных единицах, о величине плотности можно судить по ее сопоставлению с теоретической кривой, приведенной на этом же рисунке. Напомним, что вычисления по методике 1 основаны на использовании формул (21) - (22), а по методике 2 - по формуле (5).

Рис.5. Вертикальное распределение озона, полученное на основе фотометрических измерений блеска ИСЗ, 1 - наблюдения 09.09.1997 г. (методика 1, усреднение по 6 объектам), 2 - наблюдения 09.09.1997 г. (методика 2, усреднение по 7 объектам), 3 - наблюдения 29.03.1998 г. (методика 2, усреднение по 7 объектам), 4 - теоретическое распределение озона, 5 - наблюдения 29.03.1998 г. (методика 1, усреднение по 4 объектам).

Моделирование блеска проводилось с помощью доработанной программы ALAN [11] с последующим сравнением вычисляемых кривых блеска ИСЗ с результатами измерений.

При расчетах использовались исходные данные, задаваемые для конкретной длины волны:

вертикальный профиль концентрации озона (зима, средние широты северного полушария - (файл Pozim. dat);

тоже для летнего периода - (файл Pozl. dat);

вертикальный профиль коэффициента аэрозольного ослабления (км_1), фоновая модель для умеренных широт северного полушария - (файл SS. dat);

вертикальный профиль показателя преломления воздуха, модель стандартной атмосферы для северного полушария - (файл SN. dat);

число молекул в 1см3 х 10-29 по модели стандартной атмосферы (файл AN. dat);

атмосферная рефракция в функции высоты (в градусах) - (файл gam. dat);

вспомогательные величины, связанные с геометрией входа и выхода из тени ИСЗ - (файл SX. dat).

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

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

Ослабление освещенности геостационарного ИСЗ, входящего или выходящего из тени Земли, для длин волн от 450 нм до 700 нм было рассчитано на PC Реntium 200. При этом взяты модели высотного распределения парциальной плотности озона [14], высотного распределения коэффициентов ослабления аэрозолей [13] и значения плотности атмосферного воздуха [20].

В качестве примера на рис.6 показана зависимость освещенности ИСЗ от длины волны при различных фазовых углах спутника.

Рис.6. Зависимость освещенности ИСЗ от длины волны для фазовых углов

1=9.03, 2=8.80, 3=8.70, 4=8.50, 5=8.40.

Видно, что по мере погружения солнечного диска в земную атмосферу ослабление становится больше в длинах волн, находящихся ближе к УФ-области спектра. В последний момент захода спутника разница в освещенности для =500 нм и =700 нм составляет уже три порядка, что обусловлено более сильным аэрозольным ослаблением вблизи 500 нм и релеевским рассеянием, величина которого пропорциональна 1/л4. Влияние озона заметно на длине волны 600 нм, это соответствует максимуму озонного поглощения в полосе Шапюи. Для полной уверенности необходимо, чтобы фотометрические кривые относились, по крайней мере, к двум участкам спектра - вне и внутри полосы Шапюи.

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

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

Рис.7. Зависимость освещенности ИСЗ от геоцентрического фазового угла для длин волн: 1 - 650 нм, 2 - 700 нм, 3 - 530 нм (фоновая модель аэрозоля), 4 - 530 нм (среднециклическая модель аэрозоля), 5 - 500 нм, 6 - 470 нм.

Плавный ход кривых нарушается в области =8.85, которая соответствует погружению центральной части солнечного диска за средние и нижние слои атмосферы, где сильно проявляется неравномерность ослабления в разных длинах волн. Из их сравнения видно, что различие в аэрозольных моделях начинает сказываться после ш = 8.6 и остается почти стабильным до ш = 8.45, составляя примерно 3-4% от величины lg B. В области углов, меньших ш = 8.45, начинают проявляться эффекты, связанные с неоднородным высотным распределением аэрозолей в тропосфере, характерным для выбранных моделей ослабления.

Рис.8. Изменение освещенности ИСЗ в зависимости от геоцентрического фазового угла, наблюдения 09.03 1998 г., объект Spacenet 2; 5 - фильтр V, 2 - фильтр R; теоретические кривые для длин волн: 1 - 700 нм, 3 - 600 нм, 4 - 450 нм.

В качестве примера на рис.8 представлены результаты наилучшего совпадения теоретической и экспериментальной кривых блеска ИСЗ для фильтра V по наблюдениям 9 марта 1998 года.

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

Как видно из рис.8, не удается получить достаточно хорошего согласования реальных и теоретических кривых блеска. Освещенность ИСЗ после входа центра Солнца в атмосферу убывает несколько быстрее, чем дает теория. Такое же разногласие отмечалось в работе В.Г. Фесенкова [2], которое, по его мнению, можно частично устранить, существенно увеличив содержание озона в слоях, лежащих ниже 25 км.

В нашем случае такой подход не дает необходимого согласия. Следует так же отметить (см. рис.6) разногласия в высотных профилях, полученных двумя способами, а именно, отсутствие верхнего "горба" на H=50 км при расчете блеска ИСЗ по второй методике (формула (5)).

Литература

1. Фесенков В.Г. О зондировании оптических свойств атмосферы при помощи искусственных спутников. // Наблюдения ИСЗ. 1966. № 5. С.97-107.

2. Фесенков В.Г. Об исследовании вхождения ИСЗ в тень Земли // АЖ. 1967. Т.44, № 1. С.3-11.

3. Фесенков В.Г. О возможности использования лунных затмений для зондирования оптических свойств атмосферы // АЖ. 1970. Т.47. № 2. С.237-245.

4. Фесенков В.Г. О зондировании оптическим методом распыленного космического вещества в высокой атмосфере // АЖ. 1968. Т.45. № 3. С.622-633.

5. Матягин В.С. Некоторые результаты и перспективы наблюдений искусственных спутников Земли при заходе их в земную тень // Наблюдения ИСЗ. 1966. № 5. С.112-117.

6. Spitzer L. Jr. Earth satellites as research vehicles. Franklin Inst. Monograph 2, 69, 1956.

7. Линк Ф. Лунные затмения.М. 1962.200 с.

8. Розенберг В.Н. Сумерки.М. 1963.380 с.

9. Гущин Г.П. Озон и аэросиноптические условия в атмосфере.Л. 1964.342 с.

10. Пятелина С.В., Диденко А.В., Павлов В.Е. Освещенность геостационарного спутника, входящего в тень Земли // Изв. АН Каз. ССР. Сер. физ. - мат. 1990. № 4. С.87-92.

11. Миннарт М. Солнце. Солнечная система. М., 1957. Т.1.700 с.

12. Зуев В.Е., Креков Г.М. Оптические модели атмосферы.Л., 1986.256 с.

13. Г.М. Гречко, А.С. Гурвич, и др. Оптические явления при выходе и заходе Солнца в космосе // Труды ГОИ., т.71. Вып. № 205.120 с.

14. Тверской П.Н. Курс метеорологии.Л., 1962.700 с.

15. Озон в земной атмосфере. Сборник статей. Под редакцией Г.П. Гущина.Л. 1966.200 с.

16. Демченко Б.И., Диденко А.В., Матягин В.С., и др. Автоматизация наблюдений подвижных космических объектов. А-Ата. 1990.160 с.

17. Диденко А.В., Синяева Н.В., Усольцева Л.А. Сводный каталог UBVRI величин стандартных звезд для наблюдений КО // Наблюдения ИНТ. М., 1990. Т.85. С.24-30.

18. Демченко Б.И., Диденко А.В., Усольцева Л.А. и др. Зональный каталог геостационарных спутникоа. А-Ата. 1996.88 с.

19. Атмосфера стандартная. Параметры. М., 1981.180 с.

20. "Космические исследования в Казахстане", Алматы, КазГосИНТИ, 2002, с.355-373

Размещено на Allbest.ru


Подобные документы

  • Ознакомление с строением Солнечной системы. Анализ научных данных и сведений по планетам земной группы. Рассмотрение особенностей Меркурия, Венеры, Земли и Марса. Изучение размеров, массы, температуры, периодов обращения вокруг оси и вокруг Солнца.

    реферат [26,8 K], добавлен 28.01.2015

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

    презентация [1,5 M], добавлен 22.10.2011

  • Планеты Земной группы: Земля и сходные с ней Меркурий, Венера и Марс. Венера - самая горячая планета группы. Планеты-гиганты: Юпитер, Сатурн, Уран и Нептун. Блеск Юпитера, кольца Сатурна. Основные характеристики планеты Уран. Нептун и его спутники.

    презентация [2,1 M], добавлен 08.04.2011

  • Изучение основных параметров планет Солнечной Системы (Венера, Нептун, Уран, Плутон, Сатурн, Солнце): радиус, масса планеты, средняя температура, среднее расстояние от Солнца, структура атмосферы, нналичие спутников. Особенности строения известных звезд.

    презентация [1,4 M], добавлен 15.06.2010

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

    презентация [10,9 M], добавлен 02.10.2013

  • Преимущества расположения телескопа "Хаббл" вне пределов земной атмосферы. Ключевые события в освоении космоса: исследование Солнца аппаратами "Уилисс", "Сохо" и Юпитера станцией "Галилео", посадка на Марс первого марсохода, исследование Сатурна.

    презентация [2,2 M], добавлен 22.01.2014

  • Венера как землеподобная планета, происхождение её имени. Современная модель внутреннего строения Венеры, состав её атмосферы и слабость магнитного поля. Основные различия Земли и Венеры (чего не хватает Венере, чтобы стать второй обитаемой "Землей"?).

    презентация [709,0 K], добавлен 29.11.2016

  • Знакомство с основами знаний о Млечном пути как скоплением звезд. Учение Галилея о спиральной галактике и современные данные. Солнечная система, ее происхождение, состав и размеры. Основные характеристики Солнца. Планеты Земной группы и их спутники.

    презентация [1,4 M], добавлен 23.02.2014

  • Естественные и искусственные космические объекты. Изучение верхней атмосферы и космического пространства с помощью экспериментов и проведения непосредственных измерений на больших высотах с помощью искусственных спутников Земли и космических ракет.

    презентация [2,4 M], добавлен 04.02.2017

  • История образования атмосферы планеты. Баланс кислорода, состав атмосферы Земли. Слои атмосферы, тропосфера, облака, стратосфера, средняя атмосфера. Метеоры, метеориты и болиды. Термосфера, полярные сияния, озоносфера. Интересные факты об атмосфере.

    презентация [399,0 K], добавлен 23.07.2016

Работы в архивах красиво оформлены согласно требованиям ВУЗов и содержат рисунки, диаграммы, формулы и т.д.
PPT, PPTX и PDF-файлы представлены только в архивах.
Рекомендуем скачать работу.