Бесплатный автореферат и диссертация по наукам о земле на тему
Глобальное распределение составляющих радиационного баланса Земли по данным ИСЗ России и США
ВАК РФ 25.00.30, Метеорология, климатология, агрометеорология

Автореферат диссертации по теме "Глобальное распределение составляющих радиационного баланса Земли по данным ИСЗ России и США"

0034Ы г ои-> На "Р3*3* рукописи

Фомина Надежда Владимировна

ГЛОБАЛЬНОЕ РАСПРЕДЕЛЕНИЕ СОСТАВЛЯЮЩИХ РАДИАЦИОННОГО БАЛАНСА ЗЕМЛИ ПО ДАННЫМ ИСЗ РОССИИ И

США

Специальность 25.00.30 - метеорология, климатология, агрометеорология

Автореферат

диссертации на соискание ученой степени кандидата географических наук

1 О ДЕК 2009

Казань-2009

003487603

Работа выполнена в Саратовском государственном университете им. Н.Г. Чернышевского на кафедре метеорологии и климатологии

Научный руководитель: Официальные оппоненты:

Ведущая организация:

доктор технических наук, профессор Скляров Юрий Андреевич

доктор физико-математических наук, профессор Вельтищев Николай Федорович доктор географических наук, профессор Френкель Марат Ошерович

ГОУ ВПО Пермский государственный университет

Защита состоится 17 декабря 2009 г. в 15 часов на заседании диссертационного совета Д.212.081.20 Казанского государственного университета по адресу: 420008, г. Казань, ул. Кремлевская, 18, корп. 2, ауд. 1512

С диссертацией можно ознакомиться в библиотеке им. Н.И. Лобачевского Казанского государственного университета.

Автореферат разослан 12 ноября 2009 г.

Ученый секретарь диссертационного совета кандидат географических наук, доцент

Ю.Г. Хабутдинов

Общая характеристика работы.

Актуальность проблемы. В Саратовском государственном университете ведутся работы как по созданию наземной, аэростатной и спутниковой аппаратуры для измерений составляющих РБЗ так и по обработке результатов наблюдений с аэростатов и ИСЗ. Настоящая диссертация является естественным продолжением двух предшествующих (А.И. Котумы и Н.В. Семеновой), в которых были отработаны основные вопросы обработки спутниковых наблюдений от «сырых» данных на орбите, до построения глобальных карт распределений. Однако остались нерешенными ряд важных вопросов. Кроме того, в связи с ожидаемыми новыми запусками ИСЗ с модернизированной аппаратурой ИКОР-2М, возникла необходимость детального анализа этих вопросов, связанных с обработкой наблюдений как с ИСЗ «Ресурс-01» №¡4 так и в более широком плане. Было решено также выполнить анализ материалов четырех атласов NASA США, обеспечить сравнимость наших материалов с американскими и подготовить программное обеспечение для обработки предстоящих спутниковых наблюдений.

Цель настоящей работы - оценка временной и пространственной изменчивости составляющих РБЗ по материалам наблюдений с ИСЗ «Ресурс-01» №4 (Россия) и ИСЗ «Nimbus-б» и «Nimbus-7» (США).

Для достижения данной цели были поставлены следующие задачи:

1. Разработать метод и выполнить расчеты среднесуточных и среднемесячных величин альбедо по данным от ИСЗ «Ресурс-01» №4 с перспективой использования в дальнейших исследованиях.

2. Завершить работы по оценке наблюдаемости земной поверхности с гелиосинхронных ИСЗ при любой ориентации плоскости орбиты по отношению к плоскости меридиана Солнца.

3. Выполнить детальный анализ материалов 4 атласов альбедо, поглощенной радиации и УДР, изданных NASA США, для учета имеющегося опыта и обеспечения сравнимости материалов с подготовленным нами атласом УКР, альбедо и поглощенной радиации с ИСЗ «Ресурс-01» №4.

4. Создать и представить атлас наблюдений радиометром ИКОР с ИСЗ «Ресурс-01» №4. Выполнить предварительный анализ и оценить качество материалов атласа.

Использованные материалы.

1. Данные наблюдений УКР аппаратурой ИКОР с гелиосинхронного ИСЗ «Ресурс-01 »№4 за 1998-99 гг.

2. Данные наблюдений УКР аппаратурой ИКОР с ИСЗ «Метеор-3» №7 за 1994-95 гг. и их обработки.

3. Атласы NASA карт распределения альбедо и поглощенной солнечной радиации за 1975-1978 гт. и 1985-1987 гг.

4. Атласы NASA уходящей длинноволновой радиации за 1979-1987 гг.

Научная новизна работы.

1. Разработан метод и алгоритм получения среднесуточных и среднемесячных величин альбедо по измерениям аппаратурой ИКОР с ИСЗ «Ресурс-01» №4.

2. Получены количественные оценки наблюдаемости земной поверхности с гелиосинхронных ИСЗ при любой ориентации плоскости орбиты.

3. Выполнен детальный анализ материалов атласов альбедо, поглощенной радиации и УДР, изданных NASA США.

4. Обработан весь массив наблюдений с ИСЗ «Ресурс-01 »№4. Создан макет атласа наблюдений спутниковым радиометром ИКОР с ИСЗ «Ресурс-01» №4 и выполнен предварительный анализ полученных материалов.

На защиту выносятся:

1. Методики и алгоритмы расчета среднесуточных и среднемесячных - величин альбедо.

2. Анализ и результаты оценки наблюдаемости ..земной поверхности с гелиосинхронных ИСЗ для радиометров, ориентированных в надир.

3. Особенности широтного распределения и временных вариаций составляющих радиационного баланса Земли по данным атласов NASA. Методика получения глобальных величин альбедо и поглощенной радиации. Результаты анализа полученного материала.

4. Материалы атласа наблюдений по данным с ИСЗ «Ресурс-01» №4 и оценка его качества.

Научно-практическая ценность работы.

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

2. Разработана программа расчета среднесуточных величин альбедо. Доработана программа визуализации «ИКОР-М», которая обеспечивает:

получение широтного распределения УКР, альбедо и поглощенной радиации для любого витка;

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

получение глобального распределения среднемесячных величин УКР, альбедо и поглощенной радиации на сетке с ячейками 5°х5° дуги большого круга Земли;

получение среднемесячных зональных величин УКР, альбедо, поглощенной радиации.

3. На наиболее крупные программы: «Программа расчета среднесуточного альбедо Альбедометр» и «Программа анализа данных спутникового измерителя ИКОР-М» получены свидетельства о государственной регистрации.

4. Создан макет атласа составляющих радиационного баланса Земли по данным ИСЗ «Ресурс-01 №4, который готовится к опубликованию.

Апробация работы.

Основные положения и материалы работы докладывались на следующих конференциях:

1. Всероссийская конференция «Современные проблемы дистанционного зондирования Земли из космоса» (Россия, Москва, ИКИ РАН 2003).

2. Третья всероссийская открытая научная конференция «Современные проблемы дистанционного зондирования Земли из космоса» (Россия, Москва, ИКИ РАН 2005).

3. Пятая Юбилейная Огкрытая Всероссийская конференция «Современные проблемы дистанционного зондирования Земли из космоса» (Россия, Москва, ИКИ РАН 2007).

4. Всероссийская научная конференция с международным участием «Окружающая среда и устойчивое развитие регионов: новые методы и технологии исследований» (Татарстан, Казань, 2009)

на отчетных научных конференциях географического факультета СГУ и на объединенных семинарах кафедры метеорологии и климатологии СГУ и лаборатории

астрономических и геофизических исследований НИИ ЕН отделения механики и физики СГУ. '

Структура и объем работы. Диссертация состоит из введения, пяти глав, заключения. Список использованных источников содержит 84 наименования. В диссертации имеется всего 147 страниц машинописного текста, включая 50 рисунков и 15 таблиц.

Публикации. По теме диссертации опубликовано 13 работ. В том числе 2 работы [2, 3] опубликованы в журнале, издаваемом РАН и входящем в список ВАК.

СОДЕРЖАНИЕ РАБОТЫ

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

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

Альбедо, полученное по величине УКР, дает его мгновенное среднее значение по ПЗ. Если в течение светового дня надирным радиометром будет измеряться альбедо одной и той же сцены, то будут получены различающиеся величины альбедо. Это зависит от ряда причин, в том числе и от условий освещенности сцены. Поскольку на ВГА единственным источником освещенности является Солнце, то альбедо сцены в каждый момент есть функция угла падения лучей Солнца, т.е. его зенитного угла Z0: а=а(Z0). Следовательно, по мгновенному измеренному альбедо конкретного участка поверхности при конкретном Z0 необходимо оценить среднесуточное значение альбедо

Конкретно использовались модели, разработанные в проекте США ERBE. Эти модели охватывают 12 различных типов сцен. Каждый тип имеет свою особенную зависимость а от Z0 При обработке спутниковых измерений надо оценить тип сцены и дальнейшие определения среднесуточных величин альбедо для данного региона, попадающего в ПЗ радиометра, вести с соответствующей моделью.

Для описания изменений потоков радиации в зависимости от зенитного угла Солнца значения этого угла (диапазон изменений ZQ от 0° до 90°) разделены на промежутки, названные «бииами», и модель представлена в средних оценках для каждого бина. Бины выбирались из условия, что cosZQ от бина к бину изменяется на 0,1. Альбедо в пределах бина не меняется.

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

(при /=1,2,... 10)

Бин ¿=1 соответствует зенитному значению альбедо для данного типа модели.

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

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

Из определения нормированной функции альбедо следует, что для измеренного мгновенного значения «(20) при фактическом зенитном угле в момент измерения можно получить соответствующее зенитное значение альбедо, т.е. 0):

О)

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

Очевидно, что для получения среднесуточного значения аСут нужно найти коэффициент пропорциональности п'}, который рассчитывается как средневзвешенное

по инсоляции значение п^ для тех бинов, в которых для измеряемой сцены может в

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

в данный день. Коэффициент Щ зависит также от времени суток (¿^ в момент измерения). Учитывая (1), получим:

асут =иу-а(г, = 0) = яу , ° (2)

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

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

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

бина находится Ъ& в момент измерения. По географическим координатам подспутниковой точки, к которой привязывается а(20), определяется номер ячейки принятой сетки. Следовательно, сразу определяется базовая модель альбедо. Базовых моделей в таблице моделей направленного альбедо всего пять: океан, суша, снег, пустыня, смесь. Остальные созданы на основе базовых при соответствующих градациях облачности. Принадлежность мгновенного измеренного значения альбедо к

соответствующему табличному типу определяется следующим образом. По

б

полученному бину ЗУ С, отнекивается табличная величина альбедо в), наиболее

близкая к измеренному значению а(Za) среди типов моделей для данной подстилающей поверхности. Базовая модель отыскивается по типу подстилающей поверхности, а конкретная - с учетом облачности. Следовательно, к a(Zs) можно применять определенное моделями направленного альбедо значение нормированной

направленной функции n'j, рассчитанной для ЗУС в пределах конкретного бина. Тогда

для измеренной величины альбедо o(Z0) сразу находится конкретное значение «зенитного» альбедо, т.е. альбедо при ZQ=0:

a(zs = 0)=^ (3)

Таким образом, определение конкретного значения зенитного альбедо по a(Z0) в случае, когда углы Z@ заключены в пределах одного бина, производится с одним и тем

же значением n'j.

Реально над данной ячейкой сетки на ВГА Солнце в течение светового дня проходит бины от 10-го (после восхода) до некоторого минимального (поддень) и затем в обратном порядке до захода. Поскольку суточные параллели Солнца симметричны относительно небесного меридиана, можно рассматривать, например, только первую половину дня.

Выражение (3) можно использовать для получения среднесуточного значения альбедо «. Из этого выражения следует: )-a(ZS)-0)n'h Если учесть, что Солнце за день проходит количество бинов, до некоторого бина к, и получить усредненное за световой день значение ") для выбранной модели, то тогда ä=a(Ze> = 0)n'J в соответствии с принципами построения направленных моделей альбедо. Очевидно, что вклад каждого п) в среднесуточное значение альбедо пропорционален промежутку времени A'i пребывания Солнца в соответствующем бине.

А сумма At'b для бинов от восхода Солнца до меридиана равна часовому углу восхода ts выраженному во времени. Тогда среднесуточное (т.е. за световой день, т.к. ночью альбедо не измеряется) значение нормированной величины n'j, очевидно, равно:

Е Л '¿и /

ТГ\ = -,

' t,

где г= 10, 9...к, при этом к — наименьший номер бина, соответствующего ^э' для выбранной ячейки сетки. В наших обозначениях номер бина совпадает с номером часового угла Солнца при пересечении им нижней границы бина. В зависимости от широты места ячейки сетки и времени года Солнце может проходить разное число бинов. Для того, чтобы определить наименьший номер бина, находится минимальное зенитное расстояние Солнца, ^э" для любой точки на ВГА по значениям ее широты и склонения Солнца в момент его верхней кульминации.

Часовые углы Солнца в моменты восхода и захода ts определяются по формуле:

cos ts=-tg<ptgS (4)

и выражаются в единицах времени. Для определения часовых углов моментов входа Солнца в бин применяется соотношение:

.- еовг.-Бшювии?

сое(I =-----'

соырсо %0

Здесь 7-'ь - значения ЗУС для нижних границ бинов, <ь - часовые углы Солнца при пересечении им нижних границ бина.

Применяя (5) к нижним границам бинов получим: ~ 'м.

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

относительно малых Z@. А при получении среднего значения ^учитываются все бины до горизонта (включая бин 10). Но у бинов с большими значениями 2'@ нормализованные функции п'} имеют высокие значения (кроме сцены «снег»), за счет чего среднесуточные величины альбедо для этого типа орбиты будут больше мгновенных в низких широтах. В высоких широтах, когда Солнце может в течение дня двигаться в пределах 1-2 бинов у горизонта, среднесуточное альбедо будет близко к мгновенному.

Однако, в высоких и особенно приполярных широтах в различные времена года из-за влияния годовых изменений склонения Солнца в пределах от -23°27' до +23°27', восходы и заходы могут отсутствовать. Далее в главе рассматриваются в общем виде особенности суточного движения Солнца на приполярной широте.

Показано, что условием для незаходящего Солнца будет следующее: 8>90°-ф, а для заходящего и восходящего: 8<90°-ф. Так как предельное положительное склонение Солнца в день солнцестояния в нашу эпоху равно +23°27', то можно легко найти широту, начиная с которой появляются полярные дни. Очевидно, что это будет условие 90°-фп=23°27' откуда фп=90о-23°27'=66°33'. Это и есть широта северного полярного круга фп. Именно на этой широте теоретически один раз в году 22 июня Солнце не заходит, а лишь коснется горизонта в точке севера N. Таким образом, широта полярного круга является пределом, до которого в течение всего года применима разработанная методика, т.к. до этой широты всегда существуют восходы-заходы Солнца, и действует формула (4).

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

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

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

При приближении к полюсу убывает величина 90°-ф и растет число дней полярного дня. Если ф=90°, т.е. площадка находится на полюсе, 90°-ф = 0, т.е. полюс мира Р совпадает с зенитом, а небесный экватор совпадет с математическим горизонтом. При любых 8>0 Солнце не будет заходить, а при 5<0 будет кружить под горизонтом. Это тот случай, когда 21 марта Солнце пересечет экватор, поднимаясь в

северное полушарие. Это и будет день восхода. Через полгода Солнце уйдет в южное полушарие. Один восход, одни заход, полгода полярный день. После 21 марта Солнце медленно поднимаясь почти по окружностям последовательно проходит 10-й, 9-й... и пересекает нижнюю границу 6-го бина (ее зенитное расстояние 66°42' а ЗУС 22 июня 66°33'). После этого бины проходятся Солнцем в обратном порядке. Нет надобности искать среднесуточное альбедо. Просто нужно определить по 50 в каком бине происходит суточное движение Солнца в данный день. Измеренное значение альбедо равно среднесуточному. А во всех других случаях, когда ф>срп, но меньше 90°, усреднение по разработанной методике производится только для дней, когда б<90°-ср, т.е. когда еще имеются восходы-заходы Солнца.

В случае незаходящего Солнца, при 5>90°-ф, усреднение можно выполнить следующим образом. Рассмотрим рисунок 1. На нем изображена ситуация, когда 5>90°-ф, восхода-захода нет, Солнце в нижней кульминации проходит над точкой севера (Ы).

Г;

/ ! /

Рис. 1. К расчету среднесуточных альбедо в условиях «полярного дня»

Для оценки ЗУС воспользуемся формулой зенитного угла Солнца в нижней кульминации:

z™=m-<p-s (6)

Эта формула легко находится из рисунка 1: Z®*' - это ZZO® (над точкой севера N). Он состоит из 2-х углов: ZZOP = 90°-ф и ZPO© = 90°-8. Их сумма дает формулу (6). Находя по формуле (6) Zq1' и сопоставляя его с табличными величинами нижних границ бинов, получим место Солнца в соответствующем бине в момент нижней кульминации. Изображен случай, когда Солнце кульминирует над нижней границей 9-го бина. Далее, вместо определения часового угла восхода Солнца сразу принимаем часовой угол в нижней кульминации 12 часов. Затем, подставляя в формулу (5) нижнюю границу 8-го бина, получим и время пребывания Солнца в бине 9: (12 час -

и далее, как и в предложенной методике. Таким образом, в данном примере бин 10-й исключается из осреднения. Аналогично и для других приполярных широт.

Для получения среднесуточного значения поглощенной радиации на какой-либо площадке необходимо знать среднесуточную величину приходящего солнечного излучения Е [Вт/м2] и среднесуточную величину альбедо этой же площадки Gf. Тогда суточное среднее поглощенной радиации Q можно выразить как:

ß=£(l-ä) (7)

... Величину Е можно получить, пользуясь известной формулой Миланковича для расчета дневной суммы солнечной радиации Е2 на верхней границе атмосферы (ВГА):

гт

^^^-(¡^трзтЗ+созрсовЗыи!,), (8)

г л

где Е0 - величина солнечной постоянной, То - продолжительность суток (То=86400 с); - выраженный в радианах часовой угол Солнца в моменты восхода ^ либо захода и, г-текущее расстояние Земля-Солнце в астрономических единицах, <р -широта места, 8 - склонение Солнца. Для ВГА (нет рефракции) величина находится по формуле (4), откуда +15 = V. -Ь| = 1в-

Выражение (8) дает полную энергию Е^ [Дж/м2] на квадратный метр за световой день. Очевидно, что эта же энерпия получается площадкой и за сутки. Тогда для получения Е необходимо выражение (8) разделить на Т0 и расчетная формула для

Е принимает вид:

_ ^

£' = £ +соевом (Ып/,) (9)

г к

Следует отметить, что поскольку в течение суток расстояние г меняется незначительно, то можно при расчетах Е=Е0/г2 пользоваться единым значением для текущих суток.

В приполярных широтах формула приобретает вид (^=71=12 ч.): г

На полюсе (<р=90°): £ = -^соэ г&

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

Глава 1 Описывается методика оценки наблюдаемости земной поверхности с гелиосинхронных орбит радиометрами, ориентированными в надир. Раннее, в работах ЮА Склярова с сотрудниками была решена конкретная задача оценки наблюдаемости земной поверхности с ИСЗ «Ресурс-01» №4. В настоящей работе метод обобщен. Созданная на основе полученных соотношений программа позволяет оценить годовой ход наблюдаемости земной поверхности для орбит любой ориентации.

Практически у всех гелиосинхронных орбит угол наклонения ¡' задается в пределах 98-99 градусов. Следовательно, величиной, которая может изменяться в пределах 0°...+90° и 0°...-90° для каждого из узлов орбиты, является угол ориентации со, т.е. угол в плоскости экватора между заданным узлом и меридианом Солнца. Если 67=0°, то спутник является полуденным (он пересекает экватор в 12 часов местного солнечного времени). При со > 0°, при величинах 20°-45° говорят о дополуденном ИСЗ (время пересечения экватора 10,5-9 часов). При со близком к 90° спутник является утренним. Наоборот, при £»<0° спутники являются послеполуденными и вечерними для выбранного узла. При одном и том же значении со, но с различными узлами (восходящим, либо нисходящим на освещенной стороне Земли), картина наблюдаемости изменяется для северного и южного полушарий. Рассмотрен ряд примеров анализа.

Спутник с восходящим узлом на освещенной стороне Земли (см. рис. 2). Солнце находится в ЮП, склонение его 5 отрицательно, а в СП зима. Линия

ю

терминатора TET' теперь ограничивает освещенную часть поверхности Земли справа со стороны Солнца. Северный полюс Р находится в зоне полярной ночи. Восходящий узел К на солнечной стороне. Спутник движется из ЮП в СП и линия терминатора в СП для него является линией захода за нее. Условия наблюдаемости изменились, т.к. изменился порядок прохождения точек it\ (щ -момент времени, когда терминатор открывает точку перегиба М, выводя ее на освещенное полушарие) и тг2 (п2 - аналогично, момент закрытия точки M терминатором) по отношению к точкам равноденствий Y и -■ По мере движения Солнца из ЮП в СП (Солнце приближается к экватору), линия терминатора сначала проходит точку Q при ¿Ю, это день весеннего равноденствия, И лишь через интервал времени Д1, пройдет точку перегиба M (момент 7ti), выводя ее на освещенную часть параллели +81°.

захода подспутниковой точки за терминатор.

В СП с этого момента начинается период двукратных измерений до и после прохождения спутником точки М для любой параллели с ¡р < +81° вплоть до минимального значения <рт1„ 22 июня. Затем, при обратном ходе Солнца к экватору из СП в ЮП линия терминатора сначала пересечет точку перегиба М (момент п2) и лишь через отрезок времени Д1 пройдет точку 0 (день осеннего равноденствия). Ситуация изображена на рисунке 3 (А). Из него следует, что в СП период, благоприятный для наблюдений высоких широт (от точки щ до я2), теперь меньше, чем для случая нисходящего узла на освещенной стороне. Для рассматриваемого случая промежуток щп2 составляет 164 дня (186 - 2Д1), щщ 201 день (179 + 2Д(:). Аналогично изменится наблюдаемость в ЮП, где эти цифры поменяются местами. Из-за эллиптичности земной орбиты и неравномерного движения Земли по ней, промежутки времени от дня весеннего равноденствия до дня осеннего ГУ1-*£Ь и наоборот - не равны. В настоящую эпоху они составляют: <У1-*'£Ь равен 186 дней (Земля проходит афелий орбиты около 3 июля, т.е. в весенне-летний период движется медленнее); промежуток времени й-»'У равен 179 дней (перигелий Земля проходит около 4 января).

Таким.образом, гелиосинхронная орбита с восходящим узлом на солнечной стороне для дополуденного ИСЗ лучше приспособлена для наблюдения высоких широт в 1рП.

ИСЗ полуденный, со = 0. Из рисунка 2 следует, что АМР()=0, точка перегиба М совпадает с точкой () на большом круге Р<ЗЕР', перпендикулярном меридиану Солнца. Поскольку только что рассмотрена ситуация для наблюдаемости с восходящим узлом на освещенной стороне, обратимся к рисунку 2. При ®=0° восходящий спутник пересекает экватор при местном солнечном времени 12 часов (поэтому ИСЗ полуденный). Линия терминатора в СП является линией захода спутника за него, перехода ИСЗ на неосвещенную половину Земли.

А)

т s. 71, а

О 60 120 180 240 300 360

Дни года

Б)

О 60 1 20 180 240 300 360

Дни года

Рис. 3. Границы наблюдаемых и не наблюдаемых частей поверхности Земли

А) для случая, приведенного на рис. 2.

Б) для случая, для полуденных гелиосинхронных ИСЗ (а = 0, // = ±9°).

Из рисунка 2 видно, что ИСЗ в точке перегиба орбиты М выйдет на освещенную часть в день весеннего равноденствия (при 5=0°), когда линия терминатора совпадет с линией PQEP'. Ясно также, что точка М будет оставаться на освещенной части параллели +81° в течение всего весенне-летнего периода и вновь попадет в теневую часть параллели при S=0° в день осеннего равноденствия и т.д.

Результаты расчетов для этого случая (/=99°, <у=0°) приведены на рисунке 3 (Б). Видно, что кривые захода за терминатор в СП и выхода из-за терминатора в ЮП симметричны относительно точек и £ fa совпадает с rY\ а я2 - с Вновь напомним, что продолжительность периодов fY,~*— и iWY разная за счет эллиптичности орбиты Земли. В данном случае они различаются на 7 дней (186 суток первый и 179 второй).

Не трудно понять, что в случае нисходящего полуденного спутника в точке q (рисунок 2) находится нисходящий узел, общая картина не изменится. Точка перегиба М при этом располагается за плоскостью солнечного меридиана PTqP' на противоположной части круга PEP'. Годовая кривая пересечения подспутниковой точкой терминатора будет кривой выхода из-за него в СП и захода за терминатор в ЮП. Остальные обстоятельства наблюдаемости поверхности одинаковы для восходящего и

нисходящего полуденных ИСЗ. Отметим, что наиболее близким реализованным полуденным ИСЗ являлся спутник проекта ERB «Nimbus-7» с восходящим узлом.

ИСЗ послеполуденный, а <0." Для случая восходящего спутника вновь воспользуемся рисунком 2. Из него следует, что при со <0 восходящий узел будет восточнее Солнца (на экваторе за плоскостью рисунка). При прохождении спутником узла местное время в этот момент будет больше 12 часов (обычно для послеполуденных ИСЗ 14 час.30 мин. - 15 часов). В этом случае точка перегиба М окажется правее точки Q (для условий рисунка 2). Тогда, при движении Солнца к северу, терминатор сначала пересечет точку М, выводя ее на освещенную часть земного шара. И только через некоторый интервал времени At терминатор совпадает с кругом PQEP', при этом 0=0, день весеннего равноденствия. Возвращаясь, при уходе Солнца из СП, терминатор вначале пройдет точку Q (день осеннего равноденствия), и лишь через тот же интервал времени At пройдет точку перегиба М, переводя ее на затененную часть параллели +81° в СП. Ситуация в точности совпадает с той, которая рассмотрена для дополуденного (с нисходящим узлом) ИСЗ «Ресурс-01»№4. Разница в том, что теперь кривая годового хода точек пересечения терминатора в СП является кривой захода подспутниковой точки за него. Таким образом, при равных / и со, но смене положения узлов на 180° и знака со на -а условия наблюдаемости земной поверхности одинаковы (мы уже не анализируем ситуацию в ЮП, она аналогична).

ИСЗ утренний или вечерний, <у=±90°. В такой постановке представляется, что имеется 4 варианта: утренний (со = +90°) восходящий и нисходящий; вечерний (а = -90°), также два спутника. Однако в силу особенностей данных конфигураций фактически условия наблюдаемости определяются двумя вариантами.

Из рисунка 2 можно видеть, что при совпадении восходящего узла К с точкой Е экватора (перпендикулярна точке q солнечного меридиана), линия узлов совпадает с линией EW. Вокруг этой линии, при движении Солнца по широте (от -23,5° до +23,5°) в течение года, поворачивается плоскость терминатора. Теперь терминатор пересекается спутником только в узлах орбиты (т.е. в точках Е и W) при переходе из одного полушария в другое. Нет годовой линии хода широт пересечения подспутниковой точкой терминатора. В течение года плоскость терминатора дважды совпадает с плоскостью орбиты ИСЗ: при движении Солнца к северу (до 22 июня) и при его обратном движении из СП в ЮП (до 23 декабря). В такие дни подспутниковая точка весь виток совершает над терминатором. Затем, в зависимости от направления движения Солнца, либо северная половина орбиты оказывается целиком на Солнце, а южная - целиком в тени, либо наоборот. Таким образом, одномоментно происходит резкая смена условий наблюдения земной поверхности. Оценить наступление дней таких переходов несложно.

Рассмотрим случай орбиты с восходящим узлом К при аг=+90°, спутник утренний. Из рисунка 4 следует, что в СП точка перегиба орбиты М теперь расположена в плоскости, проходящей через Солнце на меридиане РМаР', противоположном Солнцу, за полюсом Р. Очевидно из предыдущего, что дуга PM=//=9° (рассматриваем орбиту с наклонением 99°). ТЕТ' - линия терминатора, © -подсолнечная точка на поверхности Земли. Склонение Солнца 6 отрицательное, Солнце - в ЮП. М' - точка орбиты, противоположная М. М'КМ - восходящая половина орбиты спутника. Для выбранного момента (до дня весеннего равноденствия) из рисунка 4 видно, что южная часть орбиты (дуга М'К, все равно что М'Е) полностью находится на солнечной стороне, а северная КМ - на теневой. Наблюдениями охвачены широты 0°...-81°, причем на всем этом интервале широт делается по два наблюдения (на

восходящей и нисходящей ветвях 9рбиты) для каждой широты. Северное полушарие полностью не наблюдаемо.

При движении Солнца к северу в момент 3= 0°, день весеннего равноденствия, терминатор совпадет с кругом РЕР', перпендикулярным меридиану Солнца. Далее, через интервал времени Д1' терминатор совпадает с треком орбиты МКМ'. Поскольку дуга РМ=9°, то это произойдет также при <50=9°, что соответствует дате 14 апреля, т.е. на 24 дня позже прохождения точки весны. Это предельное значение даты запаздывания (ранее ДЦ. Далее вся северная половина орбиты оказывается на Солнце, а южная так же одномоментно погружается в тень. При обратном (после 22 июня) движении терминатора, за 24 дня до дня осеннего равноденствия, происходит обратная смена картины. Таким образом, СП может наблюдаться с 14 апреля до 30 августа, а ЮП для этого типа орбиты с 30 августа до 14 апреля следующего года. Соответствующий график с результатами вычислений показан на рисунке 5 (А). Кажется, что покрытия наблюдениями в течение всего светлого периода от 0° до ±81° благоприятно для сбора данных. Однако, одно обстоятельство сводит на нет это преимущество. Дело в том, что спутник даже в самом лучшем случае не отходит далеко от терминатора. На терминаторе же Солнце везде находится на местном горизонте.

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

Из рисунка 4 видно, что наиболее значимая в энергетическом плане тропическая зона проходится спутником при высотах Солнца от 0° на экваторе (точка Е) до нескольких градусов. Результаты расчетов для ИСЗ с со = +90° и нисходящим узлом представлены на рисунке 5 (Б). Видно, что теперь большая часть времени наблюдений приходится на СП. Отметим еще, что при «=±90° фактически остается 2 варианта.

Рис. 4. Схема расположения основных элементов для ИСЗ утреннего (а> - 90°) с восходящим узлом

В самом деле условия наблюдений с утреннего ИСЗ (со - 90°) с восходящим узлом тождественны с условиями для вечернего ИСЗ (<и = -90°) с нисходящим узлом, т.к. можно брать такую же орбиту на ее нисходящей половине. Не рассматривая всех вариантов, отметим, что в силу указанного выше недостатка подобные орбиты на

4-Х-X

. практике не реализуются, скорее представляют теоретический интерес. Утренними и вечерними ИСЗ принято на практике называть спутники с временами пересечения ' экватора, близкими к 7 час.ЗО мин. и 19 час. 30 мин соответственно: А)

чг а

Дни гада

Б)

т а

Дни гада

Рис. 5. Границы наблюдаемых и не наблюдаемых частей поверхности Земли А) для случая, приведенного на рисунке 4 Б) для ИСЗ утреннего а = 90° с нисходящим узлом

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

Остается проблема учета переходных процессов у линии терминатора. Поле зрения СПЗ радиометров обычно имеет угловой радиус около 5° дуги большого круга на поверхности Земли. Поэтому «чистые» измерения будут наблюдаться, начиная с такого углового расстояния от терминатора (ближе - в ПЗ появляется ночная часть поверхности Земли). Для СПЗ радиометров это все же действительно не очень большая проблема. Иное дело - ШПЗ радиометры, работавшие и работающие на американских ИСЗ. Радиус их ПЗ ~30° дуги большого круга. То есть, уже на таком угловом расстоянии захватывается терминатор и часть поверхности за ним.

Понятно, что в таком случае «чистые» измерения фактически ограничены сравнительно небольшой зоной земной поверхности (приблизительно минус 30° от северной и столько же от южной частей терминатора).

Глава 3 посвящена анализу материалов атласов NASA, изданных по результатам наблюдений с американских спутников «Nimbus-б» и «Nimbus-7». На этих ИСЗ работали широко и среднеугольные радиометры. Орбиты гелиосинхронные с параметрами, близкими к параметрам орбиты российского ИСЗ «Ресурс-01» №4, так что можно сравнивать результаты наблюдений с нашим радиометром ИКОР. Всего в нашем распоряжении оказалось 4 атласа, два атласа альбедо и поглощенной солнечной

радиации, представляют в основном карты .глобальных распределений обоих параметров и сводные таблицы среднемесячных альбедо для широтных зон с шагом 5° по широте. Два атласа уходящей длинноволновой радиации (УДР) представляют карты глобальных распределений УДР на каждый месяц. Эти данные сопровождаются таблицами сферических гармонических коэффициентов разложения на каждый месяц. Одной из задач данного исследования явилась возможность проверки на материалах измерений условия равенства нулю РБЗ на верхней границе атмосферы (ВГА).

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

С увеличением широты в обоих полушариях альбедо растут, причем рост альбедо отражает сезонные особенности и соответствует увеличению зенитного угла Солнца. В тропических широтах (грубо от -25° до +25°) в любой сезон величины альбедо минимальны, имеют небольшой разброс и группируются около 0,20 для ИСЗ Nimbus-6 и около 0,25 для ИСЗ Nimbus-7. Это указывает на незначительные изменения количества облачности в среднемесячных значениях в экваториальной зоне. На всех графиках отчетливо прослеживается небольшой максимум альбедо в зоне +5°...+20° в июне и сентябре, минимум в марте и частично в декабре. Аналогично наблюдается вторичный максимум в зоне -5°... -15° в декабре-марте месяцах.

Такие изменения обычно связывают с облачностью внутритропической зоны конвергенции (ВЗК). Во все годы и сезоны величины альбедо над Антарктидой (-75°...-90°) выше, чем соответствующие величины в северных полярных широтах. В северном полушарии (СП) величины альбедо в зимне-весенний период (декабрь-март) и летне-осенний (июнь-сентябрь) различаются значительно больше, чем соответствующие по сезону величины в южном полушарии.

Очевидно, что это следствие различия в характере поверхности (материки и океаны). Наступление и сход снежного покрова в СП вносят большой вклад в эту разницу. Абсолютно наибольшие величины зонального альбедо отмечены на широтах 80° и 85° южного полушария (ЮП), Где они составляют 0,87 и 0,89, что существенно

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

НСЗ «Nimbus-7»

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

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

1 0,-175°

агл=~^ 2] <а,в; для глобального альбедо, (10)

* в,-о»

1 »/-85°

гл ~ ^ ш'"> для северного полушария, аналогично для южного, при 0/ от 90° до 175°.

Величины о, являются, весовыми коэффициентами для расчетов альбедо и находятся по формуле:

г» "„Г

со,= | ^¡пШШЛ = 0, -СОБ0^) (Ц)

где 9, - отсчитывается от северного полюса Земли, при этом 6,- - кратно 5°, т.е. О/ = 0°; 5°; 10°...85° для северного полушария и 0°; 5°; 10°... 175° для Земли в целом. Зональные значения а, выбираются из среднемесячных величин альбедо соответствующих широтных зон таблиц атласов.

Анализ полученных данных показывает следующее.

1. Глобальное альбедо имеет практически правильный полугодовой ход: минимальное значение в январе, максимальное в марте, снова минимум в июне, затем подъем до сентября - октября и минимум в декабре.

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

2. Нельзя не заметить значительной разницы в величинах глобального альбедо между данными от ИСЗ «Нимбус-6» и «Нимбус-7». Среднегодовые величины глобальных альбедо между этими спутниками отличаются практически точно на 5 единиц альбедо: 0,281 для «Нимбус-6» и 0,330 для «Нимбус-7», хотя среднегодовые величины для каждого спутника очень близки. Разница в глобальных величинах АаГЛ = 0,05 существенно превышает погрешности определения глобальных альбедо. Наш анализ показывает, что считать этот сводный ряд однородным невозможно, имеется систематическая разница в данных.

Поглощенная солнечная радиация является приходной частью РБЗ. Выше мы уже обсуждали методику получения среднесуточных величин поглощенной радиации (см. формулы 7-9).

Для получения среднемесячных значений поглощенной рациации в формулу (7) нужно подставить среднемесячные значения приходящей радиации на ВГА и альбедо. В работе была рассчитана приходящая солнечная радиация двумя методами: 1) по формуле Миланковича, для верхней и нижней параллелей (границ зон), и взято их среднее; 2) по формуле Миланковича со средним значением широты данной зоны. Приходящая радиация, рассчитанная по широтам границ зон, отличается от величин приходящей радиации, рассчитанной по средним значениям широты, не более чем на 0,1%. Поэтому в дальнейших расчетах используются данные, рассчитанные по средним значениям широты.

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

1

£ Ш1ЕР-«/> для глобального значения, <?¡=o°

S <¡>,£,0-«,) для северного полушария, и аналогично, для южного.

Здесь обозначения те же, что и в (10), а Е,- - зональные величины приходящей радиации.

Глобальное значение поглощенной радиации в течение года меняется не очень сильно. Имеется два небольших максимума в марте и сентябре, и минимум в июне-июле. На глобальное распределение поглощенной радиации большое влияние оказывают полушарные распределения. В северном полушарии минимум значений приходится на зимние месяцы (декабрь-январь), а максимум на летние (июнь-июль). В южном полушарии также минимум значений приходится на зимние месяцы, а максимум на летние. Причем в южном полушарии значение максимума больше чем в северном полушарии, а минимума меньше.Атласы УДР содержат данные за 1979 -1984 гг. и 1985 -1987 гг. Анализ среднеглобальных среднемесячных значений УДР выявил, что в холодное полугодие наблюдаются минимальные среднеглобальные значения УДР около 232 Вт/м2. Максимальные величины около 238 Вт/м2 наблюдаются в июне -августе. Это можно объяснить преобладанием материков в северном полушарии, более нагретых летом и преобладанием океанов в южном полушарии.

Одной из задач настоящего исследования было оценить вклад в УДР каждого из полушарий в отдельности. Эту задачу можно решить, используя имеющиеся в атласах данные о среднемесячных глобальных величинах УДР (Lcp) и о разности среднемесячных полушарных значений (LN-LS =ДЬ). Мы применили следующую методику.

IV = 2LCP - Ls и - Д£ + Ls; Отсюда Ls=LCP-~

По этим формулам рассчитывались значения УДР для северного и южного полушарий.

Анализ данных показал, что амплитуда изменений УДР в СП достигает 16-17 Вт/м2, а в южном около 6 Вт/м2. Очевидно, это следствие различия в поверхностях

полушарий. На материках СП большую роль играет наступление и сход снежного покрова. Поверхность суши быстрее и сильнее прогревается, чем океан. Максимум УДР в СП приходится не на июнь, a fia июль-август месяцы. Этот сдвиг почти в точности отражает уменьшение глобального альбедо СП, которое продолжается до августа месяца. А это приводит к росту поглощенной радиации и, соответственно, к росту УДР. Также прослеживается снижение в СП уровня УДР с 247 Вт/м2 в 1979 г. до 242 и менее Вт/м2 в 1986-87 гг. При этом в ЮП уровень УДР практически остается одинаковым. В обоих полушариях проявилось крупное изменение УДР в октябре 1983 г. Это свидетельствует о том, что масштабы явления Эль-Ниньо 1983-84 гг. таковы, что затронули оба полушария синхронно.

После завершения всех изложенных выше исследований можно оценить выполнение условия равенства нулю среднеглобального среднегодового радиационного баланса Земли на ВГА. К сожалению, выбор данных весьма ограничен. Есть всего два года где имеются парные данные: 1986 и 1987. Однако в 1986 г. нет данных за май месяц, а в 1987 г. за 2 месяца, ноябрь и декабрь. Но за 2 предыдущих года показания за эти месяцы были одинаковы с разницей менее 0,1%, и мы их использовали в 1987 г. С такими замечаниями результаты выглядят следующим образом.

Среднегодовое значение поглощенной радиации с формулой Миланковича получается больше, чем методом радиационного воздействия почти точно на 7 Вт/м2. В 1987 г. УДР по данным атласа в среднем за год составляет 232,48 Вт/м2, а рассчитанное значение методом «свертки» поглощенной радиации 236,2 Вт/м2. Разница всего 1,6%. С правдоподобными коррекциями в 1986 г. величина разницы тоже близка к этому. Можно считать, что для ИСЗ «Nimbus-7» РБЗ за эти годы близок к нулю.

А для ИСЗ «Nimbus-б» картина иная. При любой возможной экстраполяции данных на 1975-77 годы УДР оказывалась меньше поглощенной радиации на 6-8 процентов. По-видимому, это следствие отмеченного ранее различия в альбедо.

В главе 4 приведены характеристика наблюдательного материала с ИСЗ «Ресурс-01» №4 и анализ материала наблюдений. Одной из задач диссертации являлось построение атласа наблюдений с помощью аппаратуры ИКОР с гелиосинхронного ИСЗ «Ресурс-01» №4. Спутник «Ресурс-01»№4 за время своего активного существования с июля 1998 г. по апрель 1999 г. работал в режиме отладки. Радиометрическая аппаратура включалась нерегулярно. Измерения проводились всего в течение 109 наблюдательных суток, распределенных по 8 месяцам. В распоряжении были данные за отдельные месяцы неполного года работы ИСЗ (после чего он прекратил выдачу данных).

Отсчеты производились ежеминутно. Всего за рабочий период получено 86034 отдельных измерения УКР (соответственно альбедо). В данном эксперименте в большинстве месяцев обеспечивалось глобальное покрытие наблюдениями земной поверхности. А это главное условие возможности построения карт.

Была произведена обработка всех данных. Преследовались две цели: отработать создаваемые программы счета на данном материале и сравнить полученные результаты с аналогичными результатами обработки атласов США. В результате обработки первичных материалов измерений радиометром ИКОР создана база данных (БД) мгновенных значений УКР и альбедо и дополнительной орбитальной информации. Получение мгновенных величин УКР и альбедо велось отработанным методом «фактора формы».

Для обработки результатов была создана программа визуализации «ИКОР-М», которая позволяет визуально просматривать информацию, полученную со спутника: выводить на экран данные об УКР, альбедо в течение дневной части витка, а также выводить дополнительную информацию; производить построение графиков временного

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

Для построения карт глобальных распределений всех полученных величин был использован образец карты в цилиндрической проекции. Месячное осреднение как и ранее проводилось по системе ячеек равной площади 5°х5°. Таким образом, весь земной шар покрывается сеткой, состоящей из 1654 ячеек. Для каждой широтной зоны можно рассчитать координаты ячеек. Для каждой ячейки с известными координатами производится выборка значений из базы данных, содержащей измерения с ИСЗ «Ресурс-01 »№4. Зная координаты ячеек, легко определить тип базовой подстилающей поверхности (без облачности). Учет влияния облачности производился по методике (см. главу 1).

Получаемые из измерений мгновенные средние по ПЗ величины УКР и альбедо находят применение при построении широтных мгновенных распределений вдоль витка, при слежении за развитием различных образований и решении других задач. Однако для климатических исследований первичный интерес представляют суточные средние УКР и альбедо. А современные атласы публикуются, как правило, с использованием среднемесячных величин, получаемых из среднесуточных. Кроме того, эти же величины используются при расчетах суточных и среднемесячных значений поглощенной радиации, как одного из основных компонентов радиационного баланса Земли (РБЗ).

По рассчитанным среднесуточным значениям альбедо были получены среднемесячные их значения, осредненные по ячейкам упоминавшейся глобальной сетки. Построены карты среднемесячных значений по мгновенным и среднесуточным величинам альбедо за весь период наблюдений ИСЗ «Ресурс-01» №4. Пример карты приведен на рисунке 7.

Рисунок 7 Карта распределений альбедо по среднесуточным значениям за декабрь 1998г. с ИСЗ «Ресурс-01» №4.

20

Наличие в ряде месяцев за 1998-99 гг. глобального покрытия наблюдениями в течение нескольких суток (от 3 до 8) позволило решить задачу построения карт изолиний. Была разработана и применена методика построения карт изолиний альбедо. Наилучшие совпадения отмечены для рядов разложения, включая гармоники степени и порядка 13.

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

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

I

I

I I

I

Рисунок 8. Карта распределения поглощенной радиации за декабрь 1998 г. по данным | ИСЗ «Ресурс-01»№4

I

I

, Рассчитаны среднемесячные зональные величины для пятиградусных широтных

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

Были построены графики широтного распределения альбедо. Анализ показал, что в тропической зоне (широты от -25° до +25°) в любой сезон величины альбедо имеют | минимальные значения, и группируются около 0,25 - 0,30. Исключение составляет февраль 1999 г. В данном месяце значения альбедо отмечались в пределах 0,25. Возможно, это связано с наблюдавшимся в этот период явлением Ла-Ниньо. При приближении к полюсам в обоих полушариях альбедо растет, причем рост альбедо отражает сезонные особенности и соответствует увеличению зенитного угла Солнца.

Пользуясь разработанным методом «свертки» можно по среднемесячным I зональным величинам получить глобальные и полушарные. Результаты приведены в таблице 1.

I 2)

. Таблица 1 Глобальное распределение альбедо

Год/месяц 1998 1999

август ноябрь декабрь январь февраль март апрель

СП 0,331 0,319 0,321 0,330 0,305 0,348 0,362

ЮП 0,317 0,344 0,342 0,327 0,304 0,358 0,348

Глоб. 0,324 0,331 0,332 0,328 0,305 0,353 0,355

Аналогично были получены данные о поглощенной радиации. Результаты представлены в таблице 2.

Таблица 2 Глобальное распределение поглощенной радиации, Вт/м2

Год/месяц 1998 1999

август ноябрь декабрь январь февраль март апрель

СП 275,285 167,613 149,714 165,247 202,452 232,187. 244,453

ЮП 175,227 297,801 314,105 311,084 296,597 230,007 206,346

Глоб. 225,256 232,707 231,909 238,165 249,524 231,097 225,399

В заключении приводятся основные результаты и выводы работы:

1. Разработан "метод и проведены расчеты получения среднесуточных "и среднемесячных величин альбедо. Детально рассмотрена задача получения среднесуточных величин альбедо в приполярных широтах. Разработан алгоритм, составлены: «Программа расчетов среднесуточных альбедо Альбедометр» и «Программа анализа данных спутникового измерителя ИКОР-М», зарегистрированные установленным порядком.

2. Проведена детальная оценка условий наблюдений земной поверхности радиометром ИКОР с орбиты гелиосинхронного ИСЗ «Ресурс-01» №4 и для надирных радиометров с различных ИСЗ. Для оценки условий наблюдаемости земной поверхности была создана программа расчетов. Разработанный нами метод дает возможность рассчитать условия наблюдаемости земной поверхности для подспутниковой точки (центра ПЗ) радиометра при любой ориентации орбиты, относительно меридиана Солнца. В диссертации приведены примеры оценок для всех основных возможных положений орбиты: спутник с восходящим узлом на освещенной стороне Земли, ИСЗ полуденный, ИСЗ послеполуденный, ИСЗ утренний или вечерний.

3. Выполнен детальный анализ атласов NASA США. Отметим основные особенности в широтном распределении альбедо:

С увеличением широты в обоих полушариях альбедо растут, что в значительной степени определяется увеличением зенитных углов Солнца с ростом широты, большим вкладом облачности в средних и высоких широтах и наличием снежного (ледяного) покрова в высоких широтах вблизи полюсов. В обоих полушариях кривые роста альбедо отражают также и сезонные особенности, связанные с условиями облучения земной поверхности солнечной радиацией. Во все годы и сезоны величины альбедо над Антарктидой (-75°...-90°) выше, чем соответствующие величины в северных полярных широтах. В северном полушарии (СП) величины альбедо в зимне-весенний период (декабрь-март) и летне-осенний (июнь-сентябрь) различаются значительно больше, чем соответствующие по сезону величины в южном полушарии. В тропических широтах (от -25° до +25°) в любой сезон величины альбедо минимальны, имеют небольшой разброс и группируются около 0,20 для ИСЗ Nimbus-б и около 0,25 для ИСЗ Nimbus-7.

Это^указывает на незначительные изменения количества облачности в среднемесячных значениях в экваториальной зоне.

4. Разработана методика и получены среднемесячные и среднегодовые глобальные величины альбедо и поглощенной радиации по зональным данным атласов. Результаты анализа показали, что глобальное альбедо имеет практически правильный полугодовой ход: минимальное значение в декабре-январе, максимальное в марте, снова минимум в июне, затем подъем до сентября- октября и минимум в декабре. По полушариям картина иная. Почти точный полугодовой ход в южном полушарии с достаточно хорошей привязкой к ключевым датам - равноденствиям и солнцестояниям и значительная асимметрия в северном полушарии. Так, фаза уменьшения глобального альбедо в СП длится 6 месяцев (март - август). Такое распределение глобального альбедо объясняется значительной разницей в характере подстилающей поверхности СПиЮП.

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

6. Проведено сравнение полученных результатов расчетов поглощенной радиации, с использованием формулы Миланковича и методом радиационного воздействия. Отметим, что среднегодовое значение поглощенной радиации методом «свертки» (с формулой Миланковича) получается больше, чем методом радиационного воздействия. Разница составляет около 7 Вт/м2. Устойчивость этого значения для разных лет и различной спутниковой аппаратуры позволяет считать, что эта разница 7 Вт/м2 и есть та поправка к приближенному методу «радиационного воздействия», используемому для оценок поглощенной радиации. Отметим, что значения поглощенной радиации, полученные методом «свертки», более близки к глобальным значениям УДР из атласов.

7. Были рассчитаны значения УДР для северного и южного полушарий. По полученным значениям УДР можно сделать ряд заключений:

Амплитуда изменений УДР в СП значительно больше, чем ЮП. В северном она достигает 16-17 Вт/м2, а в южном около 6 Вт/м2. Очевидно, это следствие различия в поверхностях полушарий. Отчетливо видно, что максимум УДР в СП приходится не на июнь, а на июль-август месяцы. Этот сдвиг почти в точности отражает поведение глобального альбедо СП. Прослеживается снижение в СП уровня УДР с 247 Вт/м2 в 1979 г. до 242 и менее Вт/м2 в 1986-87 гг. При этом, в ЮП уровень УДР практически остается одинаковым. В обоих полушариях проявилось крупное изменение УДР в октябре 1983г. Это свидетельствует о том, что масштабы явления Эль-Ниньо 1983-84 гг. таковы, что затронули оба полушария синхронно.

8. Создана «Программа анализа данных спутникового измерителя «ИКОР-М». Рассчитаны значения альбедо, поглощенной радиации. Создан макет атласа составляющих радиационного баланса Земли по данным ИСЗ «Ресурс-01» №4. Предварительный анализ показывает, что с учетом неполного покрытия наблюдениями и краткости в целом периода работы ИСЗ, основные характеристики, полученные по данным атласа, достаточно близко соответствуют результатам, полученным на основе анализа атласов NASA США.

По теме диссертации опубликованы следующие работы:

1. Ю.А.Скляров, Ю.И.Бричков, А.И.Котума, Н.В.Фомина. Метод получения среднесуточных значений альбедо с использованием спутниковых направленных

моделей 11 Материалы Всероссийской конференции , «Современные проблемы . дистанционного зондирования Земли из космоса» ИКИ РАН 2003. С. 60

2. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Котума А.И., Семенова Н.В. ' Определение среднесуточных величин альбедо с использованием спутниковых направленных моделей. И Исследование Земли из космоса. 2005. №3. С. 13-21.

3. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Котума А.И., Семенова Н.В. Оценка наблюдаемости земной поверхности для надирных радиометорв с различных ИСЗ // Исследование Земли из космоса. 2006. №2. С. 1-8.

4. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Котума А.И., Семенова Н.В.Определение поглощенной солнечной радиации на верхней границе атмосферы по спутниковым измерениям. // Тезисы докладов Третьей открытой Всероссийской конференции «Современные проблемы дистанционного зондирования Земли из космоса» М„ ИКИ РАН. 2005. С.37.

5. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Котума А.И., Семенова Н.В. О наблюдаемости поверхности Земли с гелиосинхронных ИСЗ радиометрами, ориентированными в надир. // Тезисы докладов на Четвертой Всероссийской открытой конференции «Современные проблемы дистанционного зондирования Земли го космоса». М.: ИКИ РАН, 2006. С.43.

6. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Семенова Н.В., Кудряшов Д.А. О согласованности величин глобального альбедо по данным атласов наблюдений с ИСЗ «Нимбус-6» и «Нимбус-7». // Тезисы докладов на Четвертой Всероссийской открытой конференции «Современные проблемы дистанционного зондирования Земли из космоса». М.: ИКИ РАН, 2006. С.44.

7. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Семенова Н.В.О разрешающей способности среднеугольных радиометров. // Тезисы докладов на Пятой Юбилейной Всероссийской открытой конференции «Современные проблемы дистанционного зондирования Земли из космоса». М.: ИКИ РАН, 2007. С.111.

8. Скляров Ю.А., Фомина Н.В., Котума А.И., Семенова Н.В. Альбедо, поглощенная солнечная радиация и уходящая длинноволновая радиация по материалам атласов NASA США // Известия Саратовского университета. Серия Науки о Земле. 2009. Вып. 1. Т. 9. С. 44-55,

9. Скляров Ю.А., Фомина Н.В., Котума А.И., Семенова Н.В. Об особенностях определения среднесуточных величин альбедо в высоких широтах // Известия Саратовского университета Серия Науки о Земле. 2009. Вып. 1. Т. 9. С. 56-60.

10. Скляров Ю.А., Фомина Н.В. Составляющие радиационного баланса Земли по материалам атласов NASA США И Труды Всероссийской научной конференции с международным участием «Окружающая среда и устойчивое развитие регионов: новые методы и технологии исследований» Казань: Изд-во «Отечество». 2009. Т. 2. С. 200-204.

11. Скляров Ю.А., Фомина Н.В. К вопросу об оценке среднесуточных величин альбедо в высоких широтах // Труды Всероссийской научной конференции с международным участием «Окружающая среда и устойчивое развитие регионов: новые методы и технологии исследований» Казань: Изд-во «Отечество». 2009. Т. 2. С. 195-199.

12. Котума А.И., Скляров Ю.А., Фомина Н.В. Программа анализа данных спутникового измерителя ИКОР-М. Свидетельство № 2009612383 от 12 мая 2009 г.

13. Котума А.И., Скляров Ю.А., Фомина Н.В. Программа расчета среднесуточного альбедо АЛЬБЕДОМЕТР. Свидетельство № 2009612384 от 12 мая 2009 г.

Подписано в печать 06.11.2009 г. Объем - 1 печ.л. Тираж 100. Заказ № 555.

Отпечатано в типографии ООО «Техно-Декор» по адресу: г. Саратов, Московская, 160

Содержание диссертации, кандидата географических наук, Фомина, Надежда Владимировна

Глава 1 Метод и результаты расчетов среднесуточных и среднемесячных величин альбедо по данным ИСЗ «Ресурс-01» №

Глава 2 Оценка наблюдаемости земной поверхности для надирных радиометров с различных ИСЗ

2.1 Постановка и методика решения задачи

2.2 Анализ наблюдаемости земной поверхности при различных параметрах гелиосинхронных орбит

Глава 3 Анализ атласов альбедо, поглощенной радиации за период 1975-1978, 1985-1987 гг., и уходящей длинноволновой радиации за период 1979-1987 гг., изданных NASA США

3.1 Широтное распределение альбедо

3.2 Расчет глобальных величин альбедо и их вариации

3.3 Расчет приходящей и поглощенной радиации

3.4 Уходящая длинноволновая радиация и ее вариации

3.5 Радиационный баланс Земли по материалам атласов

Глава 4 О создании атласа наблюдений радиометром ИКОР с ИСЗ «Ресурс-01» №

4.1 Характеристика наблюдательного материала с ИСЗ «Ресурс-01» №

4.2 Анализ материала наблюдений с ИСЗ «Ресурс-01» №

Введение Диссертация по наукам о земле, на тему "Глобальное распределение составляющих радиационного баланса Земли по данным ИСЗ России и США"

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

В качестве основного источника и стока энергии фундаментальную роль играет радиация. Именно по этой причине определение компонентов радиационного баланса Земли (РБЗ) на верхней границе атмосферы (ВГА) и радиационного баланса на поверхности (РБП), является приоритетной задачей Всемирной программы Исследований Климата [9].

Все измерения составляющих радиационного баланса с искусственных спутников земли (ИСЗ) приводят к ВГА [28]. За ВГА условно принимается поверхность высотой 30 км над поверхностью Земли. Это дает возможность оценить, как происходит обмен радиацией с космосом любого региона планеты.

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

Компоненты РБЗ, т.е. приходящие и уходящие потоки на ВГА, определяют с одной стороны влияние радиации на климатическую систему, а с другой, они являются мерой общего отклика на это влияние. В настоящее время измерения составляющих радиационного баланса Земли составляют обязательную основу для системы мониторинга климата [5, 7].

Возрастающая актуальность данных наблюдений РБЗ и РБП определяется прежде всего их важной ролью в решении ряда ключевых задач: мониторинг пространственно-временной изменчивости климата; оценка роли радиационных факторов в формировании климата и его изменений, исследование ярко выраженных аномалий распределения облачности в тропических широтах Тихого океана (явления Эль-Ниньо и Jla-Нинья), а также соответствующих обратных связей (особенно облачно-радиационной); проверка надежности результатов численного моделирования и других. При этом особо важное значение имеет совместная интерпретация данных наблюдений РБЗ и РБП, позволяющая, в частности, определить радиационные притоки тепла ко всей толще атмосферы по данным наблюдений [5, 10-15, 19-21, 72, 87].

Большая часть инструментальных измерений составляющих РБЗ выполнена за рубежом [65, 66, 71, 73-75, 79-81, 85, 86, 88, 89, 99].

Известно, что наиболее точные и подробные данные об уходящей коротковолновой радиации (УКР), альбедо и поглощенной солнечной радиации (ПСР) получены при помощи узкоугольных сканирующих радиометров. Пространственное разрешение этих приборов 20-50 км в надире. Это обеспечивает возможность измерений чистого неба, то есть уходящей радиации без облаков.

На спутнике «Nimbus-б» сканирующая аппаратура проработала около 6 месяцев. Сам спутник и широкоугольный радиометр (ШПЗ - широкого поля зрения) проработали на ИСЗ с июля 1975 г. по октябрь 1978 г. и дали наиболее важные результаты по изучению РБЗ. На спутнике «Nimbus-7», запущенном в ноябре 1978, сканирующая аппаратура проработала около 19 месяцев, в то время как с несканирующей информация поступала до 1993 года [64, 65, 67-69, 71, 77-81, 86, 90, 99].

В восьмидесятых годах была разработана программа ERBE -Эксперимент Радиационного Баланса Земли, предназначенная для детального изучения РБЗ и его составляющих в рамках мониторинга и моделирования климата [100].

Это была трехспутниковая программа: два гелиосинхронных ИСЗ NOAA-9 и NOAA-IO на околополярных орбитах и один спутник специально для программы оценки РБЗ - ERBS (Earth radiation budget satellite). В этой программе сканирующие радиометры проработали около 6 лет, в то время как несканирующие - до конца времени жизни спутников.

Еще один проект - российско-французский ScaRab (scanner for radiation budget) [82-84]. В составе спутниковой аппаратуры работали два среднеугольных радиометра (СПЗ) Саратовского университета ИКОР и ИКОР-2.

Основная цель этого проекта - наблюдения РБЗ у ВГА. Вывод на орбиты должен был осуществляться при помощи российских спутников, запущенных на приполярные орбиты. Первый запуск (радиометры ScaRab и ИКОР) состоялся в январе 1994 года на борту российского спутника «Метеор-3» №7. Приборы работали до начала марта 1995 года, когда ИСЗ прекратил выдачу данных. С учетом месячного перерыва общее время наблюдений составило 11 месяцев. Второй запуск состоялся на ИСЗ «Ресурс-01» №4 в 1998 году, но по ряду причин, зависящих от неустойчивой работы ИСЗ, получен ограниченный материал наблюдений за 8 неполных месяцев (вновь отказала спутниковая система телеметрии).

Во всех этих экспериментах, наряду со сканирующей узкоугольной аппаратурой, использовались ШПЗ и СПЗ радиометры, которые успешно работают в течение всего срока существования соответствующего спутника. Эти радиометры позволяют решать задачи оценки зональных, глобальных величин УКР, альбедо и ПСР [88], а также регистрировать их вариации.

На основе анализа и обработки материалов наблюдений с ШПЗ радиометров были опубликованы атласы [67, 68, 101, 102 и др.], которые содержат карты среднемесячных распределений уходящей длинноволновой радиации (УДР), альбедо и поглощенной солнечной радиации. Также в этих атласах представлена информация о зональных средних альбедо за каждый месяц по широтным зонам в 5°.

В Саратовском государственном университете ведутся работы как по созданию наземной, аэростатной и спутниковой аппаратуры для измерений составляющих РБЗ [4, 6, 8, 23, 24, 36, 38-50, 91-98] так и по обработке результатов наблюдений с аэростатов и ИСЗ. Настоящая диссертация является естественным продолжением двух предшествующих (А.И. Котумы [16] и Н.В. Семеновой [25]), в которых были отработаны основные вопросы обработки спутниковых наблюдений от «сырых» данных на орбите, до построения глобальных карт распределений. Однако ряд важных вопросов остался нерешенным. Кроме того, в связи с ожидаемыми новыми запусками 1 I

ИСЗ с модернизированной аппаратурой ИКОР-2М [45, 48, 63], возникла , необходимость детального анализа вопросов, связанных с обработкой ' i наблюдений как с ИСЗ «Ресурс-01» №4 так и в более широком плане. Было решено также выполнить анализ материалов упомянутых выше четырех атласов NASA США, обеспечить сравнимость наших материалов с американскими и подготовить программное обеспечение для обработки t предстоящих спутниковых наблюдений.

Цель настоящей работы - оценка временной и пространственной изменчивости составляющих РБЗ по материалам наблюдений с ИСЗ < «Nimbus-б» и «Nimbus-7» (США) и ИСЗ «Ресурс-01» №4 (Россия).

Для достижения данной цели были поставлены следующие задачи:

1. Разработать метод и выполнить расчеты среднесуточных и , среднемесячных величин альбедо по данным от ИСЗ «Ресурс-01» №4.

2. Завершить работы по оценке наблюдаемости земной поверхности с гелиосинхронных ИСЗ при любой ориентации плоскости орбиты по отношению к плоскости меридиана Солнца. 1

3. Выполнить детальный анализ материалов 4 атласов альбедо, i поглощенной радиации и УДР, изданных NASA США, для учета I имеющегося опыта и обеспечения сравнимости материалов с , подготовленным нами атласом УКР, альбедо и поглощенной радиации с ИСЗ «Ресурс-01» №4.

4. Создать и представить атлас наблюдений радиометром ИКОР с ИСЗ «Ресурс-01» №4. Выполнить предварительный анализ и оценить качество материалов атласа. I

Использованные материалы.

1. Данные наблюдений УКР аппаратурой ИКОР с гелиосинхронного , ИСЗ «Ресурс-01» №4 за 1998-99 гг.

2. Данные наблюдений УКР аппаратурой ИКОР с ИСЗ «Метеор-3» №7 за 1994-95 гг. и результаты их обработки.

3. Атласы NASA карт распределения альбедо и поглощенной солнечной радиации за 1975-1978 гг. и 1985-1987 гг.

4. Атласы NASA уходящей длинноволновой радиации за 1979-1987 гг.

Научная новизна работы.

1. Разработан метод и алгоритм получения среднесуточных и среднемесячных величин альбедо по измерениям аппаратурой ИКОР с ИСЗ «Ресурс-01» №4.

2. Получены количественные оценки наблюдаемости земной поверхности с гелиосинхронных ИСЗ при любой ориентации плоскости орбиты.

3. Выполнен детальный анализ материалов атласов альбедо, поглощенной радиации и УДР, изданных NASA США.

4. Обработан весь массив наблюдений с ИСЗ «Ресурс-01» №4. Создан макет атласа наблюдений спутниковым радиометром ИКОР с ИСЗ «Ресурс-01» №4 и выполнен предварительный анализ полученных материалов.

На защиту выносятся следующие основные результаты:

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

2. Анализ и результаты оценки наблюдаемости земной поверхности с гелиосинхронных ИСЗ для радиометров, ориентированных в надир.

3. Особенности широтного распределения и временных вариаций, составляющих радиационного баланса Земли по данным атласов NASA. Методика получения глобальных величин альбедо и поглощенной радиации. Результаты анализа полученного материала.

4. Материалы атласа наблюдений по данным с ИСЗ «Ресурс-01» №4 и оценка его качества.

Научно-практическая ценность работы.

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

2. Разработана программа расчета среднесуточных величин альбедо. Доработана программа визуализации «ИКОР-М», которая обеспечивает: получение широтного распределения УКР, альбедо и поглощенной радиации для любого витка; возможность визуального просмотра на картах соответствующих проекций пространственно-временной изменчивости УКР и альбедо для любого витка и любого отсчета с выдачей координат подспутниковой точки, значений местного времени, зенитных углов Солнца и других сопутствующих величин в момент измерений; получение глобального распределения среднемесячных величин УКР, альбедо и поглощенной радиации на сетке с ячейками 5°х5° дуги большого круга Земли; получение среднемесячных зональных величин УКР, альбедо, поглощенной радиации.

3. На наиболее крупные программы: «Программа расчета среднесуточного альбедо Альбедометр» и «Программа анализа данных спутникового измерителя ИКОР-М» получены свидетельства о государственной регистрации [17,18].

4. Макет атласа составляющих радиационного баланса Земли по данным ИСЗ «Ресурс-01» №4 готовится к опубликованию.

Апробация работы.

Основные положения и материалы работы докладывались на следующих конференциях:

1. Всероссийская конференция «Современные проблемы дистанционного зондирования Земли из космоса» (Россия, Москва, ШСИ РАН 2003).

2. Третья всероссийская открытая научная конференция «Современные проблемы дистанционного зондирования Земли из космоса» (Россия, Москва, ИКИ РАН 2005).

3. Пятая Юбилейная Открытая Всероссийская конференция , «Современные проблемы дистанционного зондирования Земли из космоса» (Россия, Москва, ИКИ РАН 2007).

4. Всероссийская научная конференция с международным участием «Окружающая среда и устойчивое развитие регионов: новые методы и технологии исследований» (Татарстан, Казань, 2009) на отчетных научных конференциях географического факультета СГУ и на объединенных семинарах кафедры метеорологии и климатологии

СГУ и лаборатории астрономических и геофизических исследований

НИИ механики и физики СГУ.

Структура и объем работы. Диссертация состоит из введения, , четырех глав, заключения. Список использованных источников содержит 105 ' наименований. В диссертации имеется всего 150 страниц машинописного , текста, включая 50 рисунков и 15 таблиц.

Заключение Диссертация по теме "Метеорология, климатология, агрометеорология", Фомина, Надежда Владимировна

Заключение

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

1. Разработан метод и проведены расчеты получения среднесуточных и среднемесячных величин альбедо. Детально рассмотрена задача получения среднесуточных величин альбедо в приполярных широтах. Разработан алгоритм и составлена «Программа расчетов среднесуточных альбедо Альбедометр» и «Программа анализа данных спутникового измерителя ИКОР-М», зарегистрированные установленным порядком [17, 18].

2. Проведена детальная оценка условий наблюдений земной поверхности радиометром ИКОР с орбиты гелиосинхронного ИСЗ «Ресурс-01» №4 и для надирных радиометров с различных ИСЗ. Для оценки условий наблюдаемости земной поверхности была создана программа расчетов. Разработанный нами метод дает возможность рассчитать условия наблюдаемости земной поверхности для подспутниковой точки (центра ПЗ) радиометра при любой ориентации орбиты относительно меридиана Солнца. В диссертации приведены примеры оценок для всех основных возможных положений орбиты.

3. Выполнен детальный анализ атласов NASA США. Отметим основные особенности в широтном распределении альбедо: а) с увеличением широты в обоих полушариях альбедо растут, что в значительной степени определяется увеличением зенитных углов Солнца с ростом широты, большим вкладом облачности в средних и высоких широтах и наличием снежного (ледяного) покрова в высоких широтах вблизи полюсов; б) в обоих полушариях кривые роста альбедо отражают также и сезонные особенности, связанные с условиями облучения земной поверхности солнечной радиацией; в) во все годы и сезоны величины альбедо над Антарктидой (-75°. -90°) выше, чем соответствующие величины в северных полярных широтах; г) в северном полушарии (СП) величины альбедо в зимне-весенний период (декабрь-март) и летне-осенний (июнь-сентябрь) различаются значительно больше, чем соответствующие по сезону величины в южном полушарии. д) в тропических широтах (от -25° до +25°) в любой сезон величины альбедо минимальны, имеют небольшой разброс и группируются около 0,20 для ИСЗ Nimbus-б и около 0,25 для ИСЗ Nimbus-7. Это указывает на незначительные изменения количества облачности в среднемесячных значениях в экваториальной зоне;

4. Разработана методика и получены среднемесячные и среднегодовые глобальные величины альбедо и поглощенной радиации по зональным данным атласов. Результаты анализа показали, что глобальное альбедо имеет практически правильный полугодовой ход: минимальное значение в декабре-январе, максимальное в марте, снова минимум в июне, затем подъем до сентября- октября и минимум в декабре. По полушариям картина иная. Почти точный полугодовой ход в южном полушарии с достаточно хорошей привязкой к ключевым датам — равноденствиям и солнцестояниям и значительная асимметрия в северном полушарии. Так, фаза уменьшения глобального альбедо в СП длится б месяцев (март - август), что приводит к тому, что на вторичный максимум (ноябрь) и минимум (декабрь — январь) в сумме также приходится 6 месяцев. Такое распределение глобального альбедо объясняется значительной разницей в характере подстилающей поверхности СП и ЮП. Преимущественно материковая поверхность в СП и преимущественно открытый океан в ЮП.

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

6. Проведено сравнение полученных результатов расчетов поглощенной радиации, с использованием формулы Миланковича и методом радиационного воздействия. Отметим, что среднегодовое значение поглощенной радиации методом «свертки» (с формулой Миланковича) получается больше, чем методом радиационного воздействия. Разница составляет около 7 Вт/м~. Устойчивость этого значения для разных лет и различной спутниковой аппаратуры позволяет считать, что эта разница 7

•у

Вт/м~ и есть та поправка к приближенному методу «радиационного воздействия», используемому для оценок поглощенной радиации. Отметим, что значения поглощенной радиации, полученные методом «свертки», более близки к глобальным значениям УДР из атласов.

7. Были рассчитаны значения УДР для северного и южного полушарий. По полученным значениям УДР можно сделать ряд заключений: а) амплитуда изменений УДР в СП значительно больше, чем ЮП. В

9 О северном она достигает 16-17 Вт/м~, а в южном около 6 Вт/м". Очевидно, это следствие различия в поверхностях полушарий. На материках СП большую роль играет наступление и сход снежного покрова, чего нет в ЮП. Поверхность суши быстрее и сильнее прогревается, чем океан. б) отчетливо видно, что максимум УДР в СП приходится не на июнь, а на июль-август месяцы. Этот сдвиг почти в точности отражает поведение глобального альбедо СП, уменьшение которого в СП продолжается до августа месяца. А уменьшение альбедо приводит к росту поглощенной радиации и, соответственно, к росту УДР. в) отчетливо прослеживается снижение в СП уровня УДР с 247 Вт/м" в 1979 г. до 242 и менее Вт/м2 в 1986-87 гг. При этом в ЮП уровень УДР практически остается одинаковым. Таким образом, отмеченная ранее тенденция к снижению глобального значения УДР вызывается главным образом изменениями в СП. г) в обоих полушариях проявилось крупное изменение УДР в октябре 1983г. Это свидетельствует о том, что масштабы явления Эль-Ниньо 1983-84 гг. таковы, что затронули оба полушария синхронно. Это подтверждается и другими независимыми исследованиями (см. например [64, 78]).

Показано, что условие равенства нулю среднегодового РБЗ на ВГА выполняется с точностью до 1,5% для ИСЗ «Nimbus-7». Получено также, что среднегодовые значения альбедо для ИСЗ «Nimbus-б» отличаются от данных ИСЗ «Nimbus-7» на 5 единиц альбедо, и что сводный ряд для этих двух ИСЗ нельзя считать однородным.

8. Создана «Программа анализа данных спутникового измерителя «ИКОР-М». Рассчитаны значения альбедо, поглощенной радиации. Создан макет атласа составляющих радиационного баланса Земли по данным ИСЗ «Ресурс-01 №4. Предварительный анализ показывает, что с учетом неполного покрытия наблюдениями и краткости в целом периода работы ИСЗ, основные характеристики, полученные по данным атласа, достаточно близко соответствуют результатам, полученным на основе анализа атласов NASA США.

Библиография Диссертация по наукам о земле, кандидата географических наук, Фомина, Надежда Владимировна, Саратов

1. Астрономический календарь на 1990 г./ М.: Наука. Гл.ред. физ.-мат. Лит. 1989.336 с.

2. Астрономический календарь на 1993 г. / Под ред. Д.Н. Пономарева.- М.: Наука. Гл. ред. Физ.-мат. лит. 1992. 288 с.

3. Бакулин П.И., Кононович Э.В., Мороз В.И., Курс общей астрономии. М. 1970. 536 С.

4. Борисенков Е.П., Кмито А.А. Скляров Ю.А. и др. Измерения солнечной постоянной // Метеорология и гидрология. 1986. №2. С. 5-11.

5. Винников К.Я. Чувствительность климата. JL: Гидрометеоиздат, 1986. 224 с.

6. Вьюшков П.В., Скляров Ю.А. Болометрический пиргелиометр как эталонный прибор для абсолютных измерений прямой солнечной радиации// Астрономический журнал. 1964. Т. 41. Вып. 3. С. 555-558.

7. Кислов JI.B. Климат в прошлом, настоящем и будущем М.: МАИК Наука/Интерпериодика, 2001. 360 с.

8. Кмито А.А., Скляров Ю.А. Пиргелиометрия. JL: Гидрометеоиздат, 1981.-232 с.

9. Кондратьев К.Я. Всемирная исследовательская климатическая пр ограмма: Состояние, перспективы и роль космических средств наблюдений. М.: ВИНИТИ, 1972. 276 с. (Итоги науки и техники. Метеорология и климатология. Т.8).

10. Кондратьев К.Я. Глобальный климат и его изменение. JL: Наука, Ленинградское отделение, 1987. 232 с.

11. Кондратьев К.Я. Радиационные факторы современных изменений глобального климата. Л.: Гидрометеоиздат, 1980. 288 с.

12. Кондратьев К.Я. Радиационный баланс Земли как планеты // Метеорология и климатология. 1962. №1. С.28-34.

13. Кондратьев К.Я. Радиационный баланс Земли, аэрозоль и облака. М.: ВИНИТИ, 1983. 316 с. (Итоги науки и техники. Метеорология и климатология. Т. 10).

14. Кондратьев К.Я., Биненко Н.И. Влияние облачности на радиацию и климат. JL: Гидрометеоиздат, 1984. 240 с.

15. Кондратьев К .Я., Дьяченко JI.H., Козодеров В.В. Радиационный баланс Земли. JL: Гидрометеоиздат, 1988. 352 с.

16. Котума А.И. Уходящая коротковолновая радиация и альбедо системы Земля-атмосфера по наблюдениям с ИСЗ «Метеор-3» №7. Канд. дисс. 1998. 179 с.

17. Котума А.И., Скляров Ю.А., Фомина Н.В. Программа расчета среднесуточного альбедо АЛЬБЕДОМЕТР. Свидетельство № 2009612384 от 12 мая 2009 г.

18. Котума А.И., Скляров Ю.А., Фомина Н.В. Программа анализа данных спутникового измерителя ИКОР-М. Свидетельство № 2009612383 от 12 мая 2009 г.

19. Марчук Г.И., Кондратьев К.Я., Авасте О.А. и др. Межгодовая изменчивость компонентов радиационного баланса Земли по данным спутниковых измерений // Докл. АН. СССР. 1985. Т.280. №1. С.65-70.

20. Марчук Г.И., Кондратьев К.Я., Дымников В.П. Некоторые проблемы теории климата. М.: ВИНИТИ, 1981. 104 с. (Итоги науки и техники. Метеорология и климатология. Т. 7).

21. Марчук Г.И., Кондратьев К.Я., Козодеров В.В. Радиационный баланс Земли: ключевые аспекты. М.: Наука, 1988. 224 с.

22. Миланкович М. Математическая климатология и астрономическая теория колебаний климата. M.-JI. ГОНТИ. 1939. 208 с.

23. Предтеченский А.В., Скляров Ю.А. Измерения потоков прямой солнечной радиации в стратосфере // В кн.: Тезисы докладов XI

24. Всесоюзного совещания по актинометрии. Ч. И. Приборы и методы наблюдений. Таллин. Изд-во АН ЭССР. 1980. С. 21-24.

25. Предтеченский А.В., Скляров Ю.А. Болометрический пиргелиометр с автоматической компенсацией // Тр. ГГО. 1976. Вып.370. С. 3-11.

26. Семенова Н.В. Уходящая коротковолновая радиация и альбедо на верхней границе атмосферы по наблюдениям с гелиосинхронного ИСЗ «Ресурс-01» №4. Канд. дисс. 2003. 158 с.

27. Скляров Ю.А. Радиационные воздействия на климатическую систему и их роль в формировании климата / Матер. Всерос. научн. конф., посвященной 200-летию Казанского университета. Казань. 2004. С. 381383 (см. здесь также список литературы).

28. Скляров Ю.А. Измерение составляющих радиационного баланса Земли с летательных аппаратов // Проблемы оптической физики: Материалы молодежной научной школы. Саратов: Изд-во Сарат. ун-та, 1997. С. 1419.

29. Скляров Ю.А., Семенова Н.В., Котума А.И. Оценка условий наблюдений земной поверхности среднеугольным радиометром ИКОР с орбиты гелиосинхронного ИСЗ «Ресурс-01 »№4. // Исслед. Земли из космоса. 2002. №4 С. 14-20

30. Скляров Ю.А., Котума А.И. Семенова Н.В., Бричков Ю.И., Скляров В.П. О построении карт глобального распределения альбедо по наблюдениям радиометра ИКОР с гелиосинхронного ИСЗ «Ресурс-01»№4. // Исследование Земли из космоса. 2003. №3. С. 14-21.

31. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Котума А.И., Фомина Н.В. Модернизированный измеритель солнечной постоянной 2-го поколения ИСП-2 на ИСЗ "Ресурс-01"№4.// Изв. Вузов "Прикладная нелинейная динамика" Саратов. 2000. Т.8, №5. С.97-103.

32. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Котума А.И., Семенова Н.В. Определение среднесуточных величин альбедо с использованием спутниковых направленных моделей. // Исследование Земли из космоса. 2005. №3. С. 13-21.

33. Скляров Ю.А., Бричков Ю.И., Фомина Н.В., Котума А.И., Семенова Н.В. Оценка наблюдаемости земной поверхности для надирных радиометорв с различных ИСЗ // Исследование Земли из космоса. 2006. №2, С. 1-8.

34. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Котума А.И. Измерения уходящего коротковолнового излучения и альбедо радиометром ИКОР с ИСЗ «Метеор-3» № 7 // Исследование Земли из космоса. 1999. № 2. С.15-26.

35. Скляров Ю.А., Бричков Ю.И., Котума А.И., Фомина Н.В., Фейгин В.М., Липовецкий В.А. Радиометрические измерения с ИСЗ "Ресурс-01" №4. // Исследование Земли из космоса. 2000. №3. С.58-62.

36. Скляров Ю.А. О новой шкале абсолютных измерений прямой солнечной радиации // В кн.: Радиационные процессы в атмосфере и на земной поверхноти: Материалы X Всесоюзного совещания по актинометрии. JL: Гидрометеоиздат. 1979. С.64-67.

37. Скляров Ю.А. Астрономические методы в географии и метеорологии. Саратов: Изд-во Сарат. Ун-та. 1990. 36 с.

38. Скляров Ю.А. Первая сессия Международной Научной Рабочей группы проекта SCARAB (4-6 октября 1995 г., Париж) // Исследование Земли из космоса. 1996. №3. С. 121-123.

39. Скляров Ю.А., Бричков Ю.И. Актинометрические исследования в Саратовском государственном университете // География в ВУЗах России. Межвузовский сборник научных трудов. 1994. С. 80-85.

40. Скляров Ю.А., Бричков Ю.И., Воробьев В.А. и др. Измеритель коротковолновой отраженной радиации // Сб. «V Совещание по атмосферной оптике. Тезисы докладов». Томск. Томский научный центр СО АН СССР. 1991. С. 110.

41. Скляров Ю.А., Бричков Ю.И., Воробьев В.А. О двух методах спутниковых измерений уходящих радиационных потоков // Исслед. Земли из космоса. 1993. №6. С.3-11.

42. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Котума А.И. Об обработке данных спутниковых измерений уходящих радиационных потоков широкоугольными радиометрами // Исследование Земли из космоса. 1996. №3. С. 48-56.

43. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Котума А.И. Спутниковый эксперимент «Измеритель солнечной постоянной» // Письма в Астрономический журнал. 1996. Т. 22. №4. С. 318-320.

44. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Котума А.И., Пахомов Л.А., Фейгин В.М. Измеритель солнечной постоянной второго поколения на спутнике «Метеор-3» №7 // Исследование Земли из космоса. 1995. №4. С. 17-23.

45. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Попова Е.П., Сазонов Л.Б. Аэростатные измерения уходящей коротковолновой радиации // Исследование Земли из космоса. 1994. №1. С. 11-20.

46. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Предтеченский А.В. О некоторых задачах аэростатных радиометрических измерений // Сб. Краткие сообщения по физике «Исследования на высотных аэростатах». М.: ФИАН СССР. 1989. С. 49-50.

47. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Сазонов Л.Б. Некоторые результаты измерений уходящей KB радиации с аэростата // Сб. «V Совещание по атмосферной оптике. Тезисы докладов». Томск. Томский научный центр СО АН СССР. 1991. С. 74.

48. Скляров Ю.А., Бричков Ю.И., Воробьев В.А., Сазонов Л.Б. Спутниковый измеритель коротковолновой отраженной радиации и его исследования // Метеорология и гидрология. 1992. №6. С. 106-110.

49. Скляров Ю.А., Двинских В.А., Бричков Ю.И., Воробьев В.А., Котума А.И. Поиск осцилляции солнечного потока по наблюдениям со спутника «Метеор-3» №7 // Письма в Астрономический журнал. 1997. Т. 23. №10. С. 771-778.

50. Скляров Ю. А., Двинских В.А., Бричков Ю.И., Воробьев В.А., Котума А.И. Солнечная активность и осцилляции солнечного потока // Известия вузов. Сер. Прикладная нелинейная динамика. 1997. Т. 5. №5. С. 69-74.

51. Скляров Ю.А., Фомина Н.В., Котума А.И., Семенова Н.В. Об особенностях определения среднесуточных величин альбедо в высоких широтах // Известия Саратовского университета. Серия Науки о Земле. 2009. Вып. 1.Т. 9. С. 56-60.

52. Скляров Ю.А., Фомина Н.В., Котума А.И., Семенова Н.В. Альбедо, поглощенная солнечная радиация и уходящая длинноволновая радиация по материалам атласов NASA США // Известия Саратовского университета. Серия Науки о Земле. 2009. Вып. 1. Т. 9. С. 44-55.

53. Тараканов Г.Г. Тропическая метеорология. Л.: Гидрометеоиздат, 1980. 175 с.

54. Тропические муссоны / Под ред. Петросянца М.А., Белова П.Н. JL: Гидрометеоиздат, 1988. 338 с.

55. Федеральное космическое агентство. Лицензия №622К от 19 марта 2007 г. На осуществление космической деятельности. Предоставлена ГОУ ВПО Саратовский университет. Руководитель А.Н. Перминов. Per. № ООО 299.

56. Хромов С.П., Петросянц М.А. Метеорология и климатология. М.: Изд-во МГУ, 1994. 520 с.

57. Частное техническое задание на ОКР «Разработка, создание и испытания измерителя коротковолновой отраженной радиации ИКОР-М комплекса ГГАК-М для КА «Метеор-М»». М. 2003. 50 с.

58. Arking A., Vemury S. The Nimbus-7 ERB data set: A critical analysis // J.Geophys. Res. 1984. Vol. 89. №D4. P.5089-5098.

59. Barkstrom B.R., Smith G.L. The earth radiation budget experiment: Science and implementation, Rev. Geophys. 1986. V.24. P.379-390.

60. Bess T.D., Smith G.L. Atlas of wide-field-of-view outgoing longwave radiation derived from Nimbus 7 Earth radiation budget data set November 1978 to October 1985. NASA Ref. Publ.- 1186. Aug. 1987. 174 p.

61. Bess T.D., Smith G.L. Atlas of wide-field-of-view outgoing longwave radiation derived from Nimbus 7 Earth radiation budget data set November 1985 to October 1987. NASA Ref. Publ.- 1261. June 1991. 52 p.

62. Bess T.D., Green R.N., Smith G.L. Deconvolution of wide field of view radiometer measurements of Earth emitted radiation. Part 2. Analysis of first year of nimbus 6 ERB data // J. of the atmospheric sciences. 1981. Vol. 38. №3. P. 474-488.

63. Brooks D.R., Harrison E.F., Minnis P. et. al. Development of algorithms for understanding the temporal and spatial variability of the Earth's radiation balance // Rev. of Gephys. 1986. Vol. 24. №2. P. 422-438.

64. Campbell G.G., Vonder Haar Т.Н. An analysis of two years of Nimbus-6 earth radiation budget observations // Rep. CSU-ATSP-320.Colo.State Univ., Fort Collins. 1980.

65. Clouds and Earth's Radiation Energy System Experiment // Bull. Amer. Meteorol. Soc. 1996. Vol. 77. №5, May.

66. Ellis J.S., Vonder Haar Т.Н., Levitus S., Oort A.H. The annual variation in the global heat balance of the earth // J. Geophys. Res., 1978. 83, 1958-1962.

67. Hartmann D.L., Ramanathan V., Berrior A., Hunt G.E. Earth Radiation Budget Data and Climate Research // Rev. of Geophys. May 1986. V. 24. №2. P.439-468.

68. House F.B., Gruber A., Hunt G.E., Mecherikunnel A.T. History of satellite missions and measurements of the Earth radiation budget (1957-1984) // Rev. of Geophys. 1986. Vol.24. №2. P.357-377.

69. House F.B., Gruber A., Hunt G.E., Mecherikunnel A.T. History of satellite missions and measurements of the Earth radiation budget (1957-1984) // Rev. of Geophys. 1986. Vol.24. №2. P.357-377.

70. House F.B., Jaffolla J.C. One dimensional technique for enhancing Earth radiation budget observations from Nimbus 7 satellite. Presented at 1980 Int. Radiation Symp., Fort Collins. P.392-394.

71. Hucek R.R., Kyle H.L., Ardanuy P.E. «Nimbus-7» Earth radiation budget wide field of view climate data set improvement. Part I. The earth albedo from deconvolution of short wave measurements // J. Geophys. Res. 1987. V.92. № D4. P. 4107-4123.

72. Jacobowitz H., Smith W.L., Howell H.B., Nagle F.W. The first 18 months of planetary radiation budget measurements from the Nimbus-6 ERB experiment //J. Atmos. Sci. 1979. Vol.36. P.501-507.

73. Jacobowitz H., Soule H.V., Kyle H.L. et. al. The Earth radiation budget (ERB) experiment: an overview // J. Geophys. Res. 1984. Vol. 89. № D4. P.5021-5038.

74. Jacobowitz H., Tighe R.J. and Nimbus 7 experiment team. The Earth radiation budget derived from the Nimbus 7 ERB experiment // J. Geophys. Res. 1984. Vol. 89. №D4. P.501-507.

75. Kandel R.S., Monge J.L., Viollier M. et. al. The ScaRab Project: Earth Radiation Budget Observations from the METEOR Satellites // Adv. Space Research. 1994. V.14. P. 147-154.

76. Kandel R.S., Viollier M., Pakhomov L.A. et. al. The SCARAB Earth Radiation Badget Dataset // Proceeding of the Fourth International ScaRaB Science Working Group Meeting. Brussels, Belgium. June 9-11, 1997. Р.4.1.-4.27.

77. Kandel R.S., Viollier M., Raberanto P. et. al. The ScaRaB Earth Radiation Budget Dataset // Bull. Amer. Meteorol. Soc. 1998. Vol.79. №5. P.765-783.

78. Kopia L.P. The Earth Radiation Budget Experiment Scanner Instrument // Rev. of Geophys. 1986. Vol.24. №2. P.400-406.

79. Kyle H.L., Ardanuy P.E., Hurley E.J. The status of the Nimbus-7 Earth-radiation-budget data set // Bull. Amer. Meteorol. Soc. 1985. Vol.66. №11. P.1378-1388.

80. Lee III R.B., Barkstrom B.R., Smith G.L. et. al. The clouds and the Earth's radiant energy system (CERES) sensors and preflight calibration plans // J. Atmosp. Ocean Techn. 1996. V.13. №2. P.300-313.

81. Luther M.R., Cooper J.E., Taylor G.R. The Earth Radiation Budget Experiment Nonscanner Instrument // Rev. of Geophys. 1986. Vol. 24. №2. P.391-399.

82. Ohring G., Gruber A. Satellite radiation observations and climate theory // Adv. Geophys. 1983. Vol. 25. P.237-304.

83. Sellers W.D. Physical climatology. University of Chicago Press, 1969. 272 p.

84. Sklyarov Yu.A., Brichkov Yu.I., Kotuma A.I. Some Results of Measurements Received by IKOR Radiometer from "Meteor-3" №7 Satellite // Proceedings of the Fifth International ScaRaB Science Working Group Meeting. France, Paris. June 2-4, 1998. Annex 3.

85. Sklyarov Yu.A., Brichkov Yu.I., Vorobyov V.A. & Kotuma A.I. New «Solar Constant Monitor» programme // Proceedings of the Fourth International ScaRaB Science Working Group Meeting. Brussels, Belgium. June 9-11, 1997. P. 26.1-26.3.

86. Sklyarov Yu.A., Brichkov Yu.I., Vorobyov V.A. Albedo Ballon Measurements With A Satellite Monitor// Turkish Journal of Physics. 1996. Vol. 20. №4. P. 376-379.

87. Sklyarov Yu.A., Brichkov Yu.I., Vorobyov V.A. The development of the solar constant measurement's program// In: New developments and applications in optical radiometry. WRC, PMOD. Davos, Switzerland. 1990. P. 46.

88. Sklyarov Yu.A., Brichkov Yu.I., Vorobyov V.A., Bryantsev I.I. Development of a solar constant measurement programme // Metrologia 158 -1991. Vol. 28. P. 275-279

89. Sklyarov Yu.A., Dvinskikh V.A., Brichkov Yu.I., Vorobyov V.A., Kotuma A.I. Search of the Total Solar Flux Oscillations from the «Meteor-З» №7

90. Satellite Observations // Minutes of the Third International ScaRaB Science Working Group Meeting. Budapest, Hungary. November 6-8, 1996. Annex 9.4. P. 1-9.

91. Smith W.L., Hickey J., Howell H.B., Jacobowitz H., Hilleary D.T., Drummond A.J. Nimbus 6 Earth Radiathion Budget experiment. 1977. Appl. Opt. V.16. P.306-318.

92. Smith G.L., Green R.N., Raschke E. et. al. Inversion methods for satellite studies of the Earth's radiation budget: development of algorithms for ERBE mission//Rev. of Gephys. 1986. Vol. 24. №2. P. 407 421.

93. Smith G.L., Rutan D., Bess T.D. Atlas of albedo and absorbed solar radiation derived from Nimbus 6 Earth radiation budget data set July 1975 to May 1978. NASARef. Publ. - 1230. 1990. 120 p.

94. Smith G.L., Rutan D., Bess T.D. Atlas of albedo and absorbed solar radiation derived from Nimbus 7 Earth radiation budget data set November 1985 to October 1987. NASA Ref. Publ. - 1281. 1992. 58 p.

95. Smith G.L., Rutan D. Deconvolution of wide-field-of-view measurements of reflected solar radiation // J. Appl. Meteorol. 1990. V.29. №2. P. 109-122.

96. Smith G.L.,Rutan D. Observability of albedo by shortwave wide field-of-view radiometers in various orbits / Internat. Workshop on remote sensing retrieval methods.1987. Williamsburg. Virg.llp.

97. Suttles J.T., Green R.N., Minnis P. et al. Angular radiation models for Earth-atmosphere system. NASA Ref. Publ. 1184. 1988. 145 p.