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

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

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

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

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

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

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

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

Исследования последних десятилетий показали, что температурное поле и динамика стратосферы и мезосферы тесно связаны с атмосферным озоном. Он располагается в атмосфере в виде сферического слоя, внутренняя поверхность которого совпадает с поверхностью Земли, а высота может достигать 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 - функция распределения светового потока от соответствующих горизонтальных сечений солнечного диска, причем угол относится к центру Солнца, а представляет собой изменение этого угла в пределах солнечного радиуса ₀, _оз и - оптические толщи атмосферы для озона и аэрозолей, соответствующие лучам, которые пронизывают всю атмосферу и проходят на минимальной высоте от земной поверхности h₀. Для данного состояния атмосферы _оз и - функции только параметра h₀. Точное рассмотрение всякой проблемы, связанной с освещением горизонтальными лучами каких-либо элементов на достаточных высотах над земной поверхностью, требует еще учета фактора f рефракционной дисперсии. Для наблюдателя, находящегося на большой высоте и наблюдающего постепенный заход Солнца, его изображение из-за рефракционной дисперсии теряет свою круговую форму и сжимается по вертикали. Это равносильно тому, что каждый освещенный элемент атмосферы и, в частности, искусственный спутник, входящий в земную тень, оказывается в поле лучей, расходящихся вследствие различия преломления на разных высотах. Для элементов атмосферы, вносящих вклад в сумеречные явления, этот фактор изменяется примерно от 0.5 до 1.0 и обычно не учитывается. Однако для искусственных спутников, обращающихся вокруг Земли на больших высотах, его величина меняется почти на целый порядок и должна приниматься во внимание.

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

. (2)

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

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

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

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

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

, (4)

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

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

(5)

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

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

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

, (6)

а . (7)

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

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

, если i j;

(9).

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

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

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

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

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

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

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

(10).

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

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

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

, (11)

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

- общее число ослабляющих частиц, содержащихся в столбе сечением 1 см² в направлении наземный фотометр - спутник; - концентрация частиц в 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 элементов X_j вектора {x}. Согласие величин и необязательно, так как последняя величина характерна для вертикали, проходящей через точку захода Солнца, а первая характеризует ослабление в направлении наблюдатель - ИСЗ.

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

(19)

и

(20)

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

(21)

и

(22).

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

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

;

2) предположим, что y_i является неизвестной функцией высоты, так что вектор - столбец {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]. Как правило, корпус космического аппарата представляет собой цилиндр или куб, различным образом ориентированный в пространстве, на боковых поверхностях которого находятся элементы солнечных батарей (СБ). Они могут быть жестко связаны с корпусом или вынесены на специальных панелях за его пределы. Каждый КА имеет также одну или несколько приемно-передающих антенн различного типа. Следует отметить, что активный КА, как правило, стабилизирован и жестко удерживается в заданной подспутниковой точке, что следует из его функционального назначения.

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

, (25)

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

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

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

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

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

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

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

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

(26)

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

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

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