Open Access Library Journal
Vol.05 No.04(2018), Article ID:83592,65 pages
10.4236/oalib.1104250

The Sun’s Movement in the Sky Now and in the Past

Joseph J. Smulsky

Institute of the Earth’s Cryosphere, Tyumen Scientific Centre, Siberian Branch of Russian Academy of Sciences, Tyumen, Russia

Copyright © 2018 by author and Open Access Library Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: December 12, 2018; Accepted: April 5, 2018; Published: April 8, 2018

ABSTRACT

The paper considers the theory of the phenomena associated with the Sun’s passage in the sky when the parameters of the orbital and rotational motion of the Earth are changing. The author found that the obliquity changes from 14.8˚ to 32.1˚. These results explain the variations of paleoclimate and are consistent with several testimonies of paleoastronomy. The program SunPhenomn.mcd developed in MathCad to calculate the duration of solar day, equation of time, daylight, azimuths of the daily movement of the Sun, sunrises and sunsets, the length and azimuth of the gnomon shadow, duration of the polar days and nights, of seasons and other phenomena. The phenomena of the Sun in the contemporary epoch and in 4 extreme epochs in the interval of 50 thousand years ago are calculated. The example of the funnel on Yamal shows the application of the program SunPhenomn.mcd in aerospace studies of the relief. The work is of interest to a wide range of experts in the field of Earth Sciences and can be used by students and post-graduate students to carry out theses and dissertations.

Subject Areas:

Functional Analysis, Mathematical Analysis

Keywords:

Changing Parameters, Orbit, The Earth’s Axis, Duration, Day, Polar Night, Gnomon

1. Введение

Со временем мы все больше узнаем о прошлом человечества. Не так давно, найденные в Диринг-Юряхе, в Якутии, архаичные орудия древнего человека, по мнению археолога Ю. А. Мочанова, имеют возраст порядка 2-х - 3-х миллионов лет [1] . С другой стороны, на основании сравнения ДНК древних останков человека и современного местного населения разных континентов другие ученые пришли к выводу о том, что современный тип людей вышел из Африки несколько десятков тысяч лет назад [2] . Эти противоречивые выводы требуют дальнейшего изучения прошлого человечества.

Наиболее достоверным свидетельством о древнем человеке являются следы его деятельности, в том числе и древние календари. Использование календарей свидетельствует о высоком уровне развития древнего общества. Чтобы создать календарь должна быть не менее чем тысячелетняя база осознанных наблюдений явлений природы. Например, Геродот (484 - 425 гг. до н. э.) пишет, что увиденные им захоронения 330 верховных жрецов свидетельствуют об истории Египта на протяжении десяти тысячелетий. На основании наблюдений древних жрецов был создан Египетский календарь, в котором год состоял из 12 месяцев по 30 дней и пяти дополнительных дней.

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

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

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

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

Я рассматриваю точные результаты решения задачи двух тел: Земли и Солнца в полярной системе координат r(φ). В зависимости от радиуса r или угла φ получены время t движения по орбите при разные ее видах: эллипс, парабола, гипербола или радиальная прямая [3] . Поэтому всю теорию я вывожу по-новому, несмотря на то, что теория практически всех явлений, как я уже упоминал, по-видимому, есть. Все результаты новой теории сопоставлены с результатами традиционной. Поэтому у меня нет сомнения, что где-то могла закрасться ошибка.

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

Созданная М. Миланковичем Астрономическая теория изменения палеоклимата почти 100 лет назад основывается на решении трех проблем: 1) как изменяется орбита Земли; 2) как изменяется ось вращения Земли; 3) как изменяется освещенность (или инсоляция) Земли Солнцем в зависимости от параметров орбитального и вращательного движений Земли. Все эти три задачи я решил по-новому и по-другому [3] [4] [5] . По второй задаче получены существенные отличия от прежней теории. Они заключаются в том, что ось Земли имеет большие колебания: от 32˚ до 15˚ при современном наклоне оси Земли 23.4˚ к перпендикуляру плоскости ее орбиты. По прежним теориям этот наклон изменяется от 22.26˚ до 24.32˚. Таким образом, колебания оси Земли в новой теории в 7 - 8 раз больше по сравнению с прежними колебаниями. Во столько раз больше колебания инсоляции. Кроме того, получены другие периоды колебаний и другие моменты наступления экстремумов этих колебаний.

При сопоставлении новых результатов Астрономической теории с палеоклиматом, например, в Западной Сибири, выяснилось, что экстремумы палеоклимата совпадают с экстремумами инсоляции: потепление в Голоцене 4.2 тысяч лет назад (т.л.н.), Сартинское оледенение 15.9 т.л.н., Каргинское межледнивековье 31.3 т.л.н. и Ермаковский ледниковый период 46.4 т.л.н. [3] .

В ледниковые периоды на шельфе и побережье Карского моря образовался Ледниковый щит, который соприкасался со Скандинавским щитом. Высота ледниковых щитов достигала 1.5 - 2 км. Перекрывался сток рек Оби и Енисея в Ледовитый океан, и Западно-Сибирская равнина превращалась в пресноводное море. Уровень его повышался и через Тоболо-Тургайскую ложбину шел сток в Приаралье. Затем по южному проходу вода из Западно-Сибирского моря попадала в Каспий. При повышении уровня Каспия больше +27 м вода сибирских рек по Кумо-Манычской впадине могла проходить в Черное море.

Сток льда с Карского ледникового щита создавал ряд морен, самая южная гряда которых образовала Сибирские Увалы. Параллельными рядами остатки древних морен простираются до самого Карского моря.

В эти же эпохи происходили аналогичные события в Европе [6] и в Северной Америке [7] . Ледниковые периоды приводили к кардинальным изменениям в природе, в том числе и в судьбе человечества. В вышеупомянутых генетических исследованиях [2] миграция человека из Африки вдоль побережья Индийского океана в Полинезию и Австралию, затем вдоль берегов Тихого океана в Китай, Японию, Чукотку и через современный Берингов пролив в Америку приурочивается к ледниковым эпохам. Как уже упоминалось, можно усомниться в одноцентровом очаге зарождения человечества, но в ледниковых периодах, как причинах миграции человека, многие авторы единодушны.

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

В древних календарях также большую роль играют феномены Луны. В данной работе они не рассматриваются. Однако база для них имеется. Эволюцию орбитального движения планет и Луны я исследовал за 100 млн. лет [8] . Выяснилось, что ось орбиты Луны вращается (прецессирует) вокруг оси подвижной орбиты Земли с периодом 18.6 лет. Вращение происходит с небольшими колебаниями. На небесной сфере это отражается в том, что плоскость орбиты Луны колеблется в диапазоне ±5˚ по отношению к эклиптике ? пути Солнца на небесной сфере. И эти колебания неизменны в течение 100 млн. лет. Поэтому на приведенный алгоритм в программе SunPhenomn.mcd нужно наложить алгоритм изменения современных феноменов Луны, и он автоматически будет подстраиваться под изменяющиеся параметры орбиты Земли. Надеюсь, в дальнейшем такое дополнение будет сделано.

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

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

Как читать эту работу? При первом чтении нужно отвлечься от формул и их выводов и познакомиться с общим содержанием работы: о чем идет речь? При втором чтении нужно вникнуть в содержание рисунков и графиков: что на них изображено? При третьем чтении уже можно познакомиться с выводом формул и обоснованием программы на SunPhnmen.mcd. А тем, кто решил ее использовать, нужно познакомиться с содержанием программы и произвести расчеты разных феноменов Солнца для широты места своего нахождения в современную эпоху T = 0 и в любую другую эпоху T ≠ 0.

2. Новые решения по эволюции вращательного движения Земли

Новые решения задачи об эволюции вращательного движения Земли [9] [10] показаны линией 1 на Рис. 1 в виде эволюции угла наклона ε за прошедшие 200 тыс. лет. Угол ε―это угол между плоскостями экватора Земли и ее орбиты. Он изменяется от 14.8˚ до 32.1˚ [3] при современном его значении 23.4˚. Примерно такой же диапазон изменения ε получен при решении задачи за будущие 200 тыс. лет [11] . В то же время по прежним теориям, например [12] (см. Рис. 1, линия 2) угол ε изменялся от 22.21˚ до 24.43˚. Из приведенных значений видно, что диапазон колебания угла наклона по новым решениям в 7 - 8 раз больше, чем по прежним теориям.

Эти результаты на протяжении нескольких лет проверялись. Задача была решена еще тремя другими методами [3] и [13] . Результаты неизменились. Кроме того, полученные в этой задаче изменения параметров вращательного

Рис. 1. Эволюция угла наклона ε (в радианах) экватора Земли к плоскости ее орбиты за 200 тыс. лет в прошлое. Сравнение новых результатов 1 с результатами прежних теорий 2 на примере работы Дж. Ляскара и др. [12] . В градусах приведены максимальные и минимальные значения угла ε. T―время в тыс. лет от 30.12.1949 г.

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

Как видно из Рис. 1, новые решения 1 отличаются от прежних 2 не только амплитудами колебания, но и моментами наступления экстремумов. Эти отличия в эволюции угла наклона ε приводят к таким же существенным отличиям в эволюции инсоляции Земли [3] и [11] . При сопоставлении нового изменения инсоляции за прошедшие 50 тыс. лет с палеоклиматом было установлено [4] и [8] , что оно объясняет все известные его колебания: суровый ледниковый период 46 тыс. лет назад, сильное потепления 31 тыс. лет назад, последний ледниковый период 16 тыс. лет назад и небольшое потепления 4 тыс. лет назад.

Существуют свидетельства, что угол наклона ε изменялся более существенно, чем это следует из прежних теорий. Например, М. И. Исрапилов [14] исследовал сотни древних календарей в Дагестане. Он пришел к выводу, что угол наклона ε изменялся от 19˚ - 20˚ до 32˚ - 33˚. Возраст древних календарей у него доходит до 65 тыс. лет. Б. Г. Тилак [15] исследовал тексты Вед, Бхагавад-Гиты, Авесты и др. древних источников. В них используются характеристики полярной области, как будто древние арии обитали в ней. Это можно объяснить только большими углами наклона ε. В этом случае полярный круг проходит южнее современного и, дополнительно, на этих широтах становится теплее. Тем самым создаются необходимые для жизни человека условия в полярной области.

Ряд календарей в Дагестане М. И. Исрапилов [14] называет обсерваториями, т.к. они позволяли определять моменты наступления лунных и солнечных затмений. Таких обсерваторий найдено достаточно много, в том числе Стоунхендж в Англии [16] и Ales Stones на юге Швеции [17] . Многие из них, в том числе последние ориентированы по азимуту восхода Солнца в дни летнего и зимнего солнцестояния. Такие обсерватории и календари называют азимутальными. В зенитных обсерваториях и календарях [14] объекты ориентированы так, чтобы можно было фиксировать зенитные углы солнцестояний, равноденствий и других моментов, важных для хозяйственной деятельности человека. В календарях третьего вида, которые были распространены на Алтае и в Сибири [18] , календари выполнялись в виде жезлов, пластин, фигур животных и даже женских украшений. На них насечками или другим способом отмечены дни года и моменты наступления различных феноменов Луны и Солнца. По этим календарям можно определить продолжительность интервалов между днями равноденствий и солнцестояний, а также временные характеристики других астрономических явлений.

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

Целью настоящей работы является создание такой теории. В отличии от существующей таковой в астрономии, в ней, как уже отмечалось, используется другой алгоритм определения долгот Солнца при его годовом движении. Кроме того, все математические операции алгоритмизированы и представлены в виде программы на SunPhnmen.mcd в среде MathCad. Для ориентировки исследователей древних календарей рассчитаны солнечные феномены при современных и при отличающихся от современных параметрах орбитального и вращательного движения Земли. Они взяты для четырех экстремумов ε за последние 50 тыс. лет (см. Рис. 1): T = 2.8; 15.32; 31 и 46.44 тыс. лет назад.

3. Геометрические характеристики движения Солнца по небосводу

3.1. Движение Солнца в течение года

Рассмотрим положение Солнца S над точкой М земной поверхности (Рис. 2). Плоскость горизонта в т. М на небесной сфере 1 нанесена горизонтальным кругом HH¢. Перпендикуляр к плоскости HH¢ пересекает поверхность небесной сферы 1 в точке зенита Z. Солнце S совершает вокруг Земли годовое движение по орбите, которая проектируется на небесную сферу в виде круга эклиптики EE¢. Движение происходит против стрелки часов с началом отсчета долготы λ Солнца в точке весеннего равноденствия g. В этой точке Солнце находится в плоскости экватора AA¢, когда из южного полушария переходит в северное.

Видимое орбитальное движение Солнца вокруг Земли обусловлено неравномерным движением Земли вокруг Солнца по орбите, которая является эллипсом. Относительно неподвижных звезд период движения по орбите Psd = 365.25636042 дней. Он называется сидерическим годом. В работах [19] [20] представлен алгоритм расчета долгот Солнца λj1, где j 1 = 1 , 2 , , 365 , по дням года T d , j 1 = 1 , 2 , , 365 . Отсчет дней начинается от момента весеннего равноденствия, т.е. j1 = 0 соответствует моменту равноденствия. Этот алгоритм основан на точном решении задачи двух тел [21] [22] . Задача рассматривается в полярной системе координат (r, jo), где r ? расстояние Земли от Солнца, а jo ? угловое положение Земли на орбите относительно перигелия орбиты. Время tfp движения Земли по орбите от точки перигелия до точки ее нахождения с углом jo рассчитывается по формуле

t fp = t fp при φ o π и t f p = 2 t a t f p при π < φ o 2 × π , (1)

Рис. 2. Основные геометрические характеристики Солнца S и наблюдателя M на земной поверхности. Характеристики небесной сферы 1: AA¢―плоскость подвижного экватора; ЕЕ¢―плоскость подвижной эклиптики, а ε―угол между плоскостями AA¢ и ЕЕ¢; N―северный полюс; Характеристики плоскости горизонта НН¢: дуга NrdN = j − географическая широта точки M на плоскости горизонта; Z − зенит точки M; Nrd, Est, Sth и Wst − точки Севера, Востока, Юга и Запада, соответственно, на круге горизонта НН¢. Характеристики годового движения Солнца по кругу эклиптики ЕЕ¢: d = SB − склонение Солнца; λ = gS − долгота Солнца; точки g, Е¢, g¢ и Е―точки нахождения Солнца в дни весеннего равноденствия, летнего солнцестояния, осеннего равноденствия и зимнего солнцестояния, соответственно; SS1 − путь Солнца за одни солнечные сутки. Характеристики суточного движения Солнца по кругу SrMdSsMn: точки Sr, Md, Ss и Mn положения центра Солнца в моменты восхода, полудня, захода и полуночи, соответственно; z = ÐZMS − зенитный угол Солнца; ω = ÐZNS − часовой угол Солнца, отсчитываемый от полудня.

где

t f p = R p v p { ( α 1 + 1 ) | sin φ o | ( 2 α 1 + 1 ) [ ( α 1 + 1 ) cos φ o α 1 ] + α 1 [ arcsin α 1 + 1 α 1 cos φ o ( α 1 + 1 ) cos φ o α 1 π / 2 ] ( 2 α 1 1 ) 3 / 2 } ; (2)

Rp―радиус перигелия;

vp―скорость Земли в перигелии.

α 1 = μ 1 / ( R p υ p 2 ) ―параметр траектории;

μ1 = −G(M + m)―параметр взаимодействия;

G―гравитационная постоянная;

M―масса Солнца;

m―масса Земли;

t a = R p v p π α 1 ( 2 α 1 1 ) 3/2 ―время движения от перигелия до афелия, которое следует из (2) при jo = p. Здесь используется параметр траектории α1, так как он известен при постановке задачи, а эксцентриситет орбиты e определяется в результате решения задачи двух тел. Для гравитационного взаимодействия эти параметры однозначно связаны: e = ( 1 + α 1 ) / α 1 .

Этот алгоритм отличается от традиционного, основанного на уравнении Кеплера. Он, как уже отмечалось, базируется на точном решении задачи двух тел и является более простым в использовании. Программа Insl2bd.mcd в среде MathCad для расчета долгот λj1 и других характеристик инсоляции Земли имеется в свободном доступе: http://www.ikz.ru/~smulski//Data/Insol/. Алгоритм для определения долгот λj1 включен в пункты 3 - 5 программы на SunPhnmen.mcd расчета феноменов Солнца, которая представлена в Приложении. Текст программы на SunPhnmen.mcd также имеется в свободном доступе: http://www.ikz.ru/~smulski//Data/Insol/.

Как уже отмечалось, формулами (1)-(2) определяется время движения Земли от точки положения ее перигелия PE (см. Рис. 2). Перигелий относительно восходящего узла g определяется углом φpg. Относительно Земли Солнце движется по той же орбите, но его перигей PS сдвинут на 180˚, т.е. угол перигея Солнца ν = φpg + π. Так как долгота Солнца λ, также отсчитывается от восходящего узла g (см. Рис. 2), то полярный угол φo в формулах (1)-(2) выразится через долготу λ так:

φ o = λ ν = λ φ p γ π (3)

Итак, формулами (1)-(3) определяется время годового движения Солнца по небесной сфере, т.е. по эклиптике EE', в зависимости от его долготы λ.

Эксцентриситет орбиты e, угол перигелия φpg и угол наклона орбиты ε к экватору изменяются со временем. Они определяются в результате решения двух задач об орбитальном и вращательном движении Земли. В рассматриваемых программах в среде MathCad эти параметры в зависимости от времени T, считываются из файла, например, файла OrAl-200.prn − за 200 тысяч лет назад (т.л.н.).

Положение Солнца S (Рис. 2) на эклиптике EE' известно в любой момент времени t. А угловое расстояние Солнца SB от экватора, т.е. угол cклонения δ, определяется известным в астрономии выражением

δ = arcsin ( sin ε sin λ ) . (4)

Вывод выражения (4) имеется в работах [19] [20] . Отметим, что все используемые в настоящей работе формулы выведены нами. Многие из них известны в астрономии. Однако необходимо осуществлять их вывод по-новому, чтобы быть уверенным в возможности их использования при измененных параметрах орбитального и вращательного движения Земли.

Если в формулу (4) подставить λ = λj1, то мы получим склонение dj1 для каждого дня года. В этой формуле и в последующих, с целью упрощения изложения, индекс j1 опускаем. Координатами λ и d полностью определяется годовое движение Солнца по небосводу. В дальнейшем потребуется угловое расстояние Солнца gB (Рис. 2) на экваторе AA¢, которое известно в астрономии как прямое восхождение αS. Здесь индекс S, обозначающий Солнце, введен для отличия при обозначении значком α других параметров.

В сферическом прямоугольном треугольнике gSB известны углы: B = π/2; γ = ε и две стороны: gS = λ, BS = δ. Согласно преобразованиям сферической геометрии, например, формулы (1.1.013) в работе [23] можем записать для неизвестной стороны gB = αS:

sin α S = tg δ ctg ε .

После подстановки в это выражение δ из (4) получаем прямое восхождения Солнца в следующем виде:

α S 0 = arcsin cos ε sin λ 1 ( sin ε sin λ ) 2 ; (5)

α S = α S 0 при λ 0. ; α S = π α S 0 при 0. < λ 1.5 π ; α S = + α S 0 при λ > 1.5 π .

В формулу (5) входит долгота Солнца λ, которая изменяется от 0 до 2π. Так как функция arcsin многозначна, то на интервалах λ больших π/2 необходимо выбирать нужные значения по представленному алгоритму.

Плоскости экватора AA¢ и эклиптики EE¢ изменяются в пространстве (Рис. 2), вследствие чего точка весеннего равноденствия g по часовой стрелке перемещается по эклиптике за год на 50².25641. Как уже отмечалось, годовое движение Солнца проходит по кругу эклиптики EE¢ против часовой стрелки, что отражается изменением долготы λ, начиная от точки весеннего равноденствия g. Поэтому время прохождения Солнцем двух последовательных точек весеннего равноденствия, т.е. тропический год Ptr = 365.24219879 дней, меньше сидерического года Psd. Для того, чтобы сезоны года не смещались по датам, современный календарь основан на тропическом годе.

Солнце S при годовом движении по эклиптике EE' (см. Рис. 2) в момент весеннего равноденствия находится в т. g (λ = 0), а в день зимнего солнцестояния Солнце проходит т. E' (λ = π/2), в момент осеннего равноденствия находится в т. g' (λ = π), а в день зимнего солнцестояния проходит т. E (λ = 3π/2). В п. 8 программы на SunPhnmen.mcd вычислены моменты в днях прохождения точек лета Tdsm, осени Tdau и зимы Tdwn. Точку весеннего равноденствия, как уже отмечалось, Солнце проходит в момент Tdsp = 0. Так как ряд долгот λj1 расположен по целочисленному ряду дней Td,j1, то эти расчеты проводятся интерполяцией по соседним значением долгот. Продолжительности сезонов: весны, лета, осени и зимы рассчитываются по разностям этих моментов, соответственно:

Δ T d s p = T d s m T d s p ; Δ T d s m = T d a u T d s m ; Δ T d a u = T d w n T d a u ; Δ T d w n = P t r T d w n . (6)

Результаты расчетов моментов наступления сезонов и их длительностей, согласно (6), приведены в табл. 1 для пяти эпох. Как уже отмечалось, эпохи выбраны в моменты наступления экстремумов угла наклона ε. Как видно из табл. 1 величина ε изменяется от 0.25841 до 0.55875 в радианах, или от 14.8˚ до 32.1˚ в градусах. В современную эпоху (T = 0) наибольшую длительность имеет лето (ΔTdsm= 93.654 дня), а наименьшую―зима (ΔTdwn = 88.981 дня). В эпоху T = 2.8 т.л.н лето имеет еще большую длительность, а наименьшую―осень. В эпоху T = 15.32 т.л.н лето имеет наименьшую длительность, а наибольшую―осень. Аналогичные изменения длительности сезонов происходят в остальные эпохи. Однако эти изменения не очень большие и не превышают 6 дней.

В Таблице 1 приведены длительности сезонов для северного полушария. Эти данные применимы и для южного полушария, если весну заменить осенью, а лето―зимой. Следует отметить, что приведенные в табл. 1 длительности сезонов для современной эпохи (T = 0) совпадают с известными в астрономии.

3.2. Суточное движение Солнца

Земля совместно с наблюдателем M (см. Рис. 2), кругом горизонта HH¢ и меридианом NZE¢MdA¢ вращается вокруг оси вращения Земли MN, с угловой скоростью ω E . Вращение происходит против часовой стрелки. Поэтому Солнце относительно Земли и, в частности, относительно круга горизонта HH¢ по часовой стрелке перемещается по кругу SrMdSs параллельно экватору AA¢. В точке Sr оно восходит над горизонтом HH¢, в точке Md находится в полдень, а в точке Ss заходит за горизонт, а в точке Mn находится в полночь. Часовой угол ω Солнца S отсчитывается от меридиана NZE¢MdA¢, проходящего через точку полдня Md. Часовой угол ω равняется дуге A¢B на круге экватора AA¢. Угловое расстояние Солнца от зенита Z определяется дугой ZS = z, где z называется зенитным углом.

В сферическом треугольнике NZS известны две стороны: NZ = π/2 − φ; NS = π/2 − δ и угол N = ω между ними. Тогда по теореме косинусов [23] cos Z S = cos N Z cos N S + sin N Z sin φ N S cos N для зенитного угла получим выражение в следующем виде [19] [20] :

cos z = sin δ sin φ + cos δ cos φ cos ω . (7)

Таблица 1. Количество дней до начала сезонов в Северном полушарии и их продолжительность в эпохи T наступления экстремумов угла наклона ε: T − время в тыс. лет от 30.12.1949 г.

Зенитный угол z отсчитывается от точки зенита Z против часовой стрелки. В точке восхода Sr он равен: z = −π/2, в точке захода Ss угол z = π/2. При этих углах из выражения (7) определяется часовой угол заходов и восходов Солнца.

ω 0 = ± arcos ( tg φ tg δ ) . (8)

В представленном на Рис. 2 положении наблюдателя M и Солнца S длительность дня больше длительности ночи. При нахождении Солнца S в полдень в точке E¢ длительность дня будет наибольшая. Это точка летнего солнцестояния, по-другому, летнего солнцеворота. До этого момента Солнце каждый день приближалось к зениту Z, а в последующие дни оно будет удаляться от зенита. При нахождении Солнца S в точке E наибольшей будет длительность ночи. Это точка зимнего солнцестояния. А при нахождении Солнца S в точках g или g¢ его суточное перемещение будет происходить по кругу экватора AA¢. Этот круг пересекает круг горизонта HH¢ по его диаметру, поэтому время нахождения Солнца над горизонтом и под ним одинаково, т.е. длительность дня равна длительности ночи.

Если наблюдатель в точке M на Рис. 2 будет находиться на большей широте, т.е. дуга NrdN будет больше, то окружность SMd не пересечет круг горизонта HH¢. В этом случае для наблюдателя M наступит полярный день. При нахождении Солнца S¢ в южном полушарии вблизи т. E круг его суточного движения также не пересечет линию горизонта HH¢. В этом примере широты для наблюдателя в точке M наступит полярная ночь.

4. Длительность солнечных суток

Длительность солнечных суток Tsd определяется периодом прохождения Солнца через точку полдня Md (Рис. 2). За время Tsd Солнце S по эклиптике EE¢ переместится в точку S1 с прямым восхождением αs1 = γB1, которое определяется формулой (5). В связи с этим полдень наступит при часовом угле ωsd0 < 2π на величину Δα0 = αS1 − αS, т.е. часовой угол солнечных суток будет

ω s α 0 = 2 π Δ α 0 . (9)

Как следует из (5) величина Δα0,j1 определяется разностью долгот Солнца λj1 за два соседних дня. Величина Δλj1 изменяется по двум причинам: из-за неравномерного движения Солнца по эллиптической орбите и из-за наклона орбиты под углом ε к плоскости экватора. Из-за наклона одинаковые дуги SS1 на круге эклиптики EE¢ проектируются в неодинаковые дуги BB1 на круге экватора AA¢. Разность прямых восхождений за два соседних дня запишем в виде:

Δ α 0 , j 1 = α s , j 1 α s , j 1 1 .

Найдем среднюю за год разность прямых восхождений

Δ α m = j 1 = 1 365 Δ α 0 , j 1 365 . (10)

Тогда по отношению к средней величине Δαm разность долгот за солнечные сутки будет Δα = Δα0 − Δαm и часовой угол солнечных суток запишется

ω s d = 2 π ( Δ α 0 Δ α m ) . (11)

А длительность солнечного дня в часах будет иметь вид:

T s d j 1 = 24 ω s d = 12 ( 2 Δ α 0 , j 1 Δ α m π ) , (12)

где 24/2π―коэффициент преобразования дуги, измеряемой в радианах, в часы.

Нетрудно убедиться, что средняя за год длительность солнечных суток, согласно (12), Tsdm = 24 часа. Тогда отклонение длительности солнечных суток в минутах от средних составляет:

Δ T s d = 6 0 ( T s d 24 ) . (13)

Приведем некоторые значения для 2015 года. Отклонение солнечных суток от средних в день весеннего равноденствия j1 = 1 будет ΔTsd = 0.297 мин, наибольшее значение ΔTsd = 0.358 минуты при j1 = 181, и наименьший солнечный день при j1 = 278: ΔTsd = −0.497 минуты. В современной цивилизации счет времени m происходит по средним солнечным суткам Tsdm. Они делятся на 24 часа, 1 час состоит из 60 минут, а минута―из 60 секунд. За счет отклонения ΔTsd солнечных суток от средних Tsdm накапливается отличие η0 солнечного времени от среднего. Последовательное суммирование отклонений запишем в виде

η 0 , j 1 = η 0 , j 1 1 + Δ T s d , j 1 1 . (14)

Найдем среднее отклонение за год

η 0 m = ( j 1 = 1 365 η 0 , j 1 ) / 365 . (15)

Средняя величина отклонения η 0 m образуется при счете времени по средним солнечным суткам. Тогда отклонение солнечного времени от среднего солнечного времени будет

η j 1 = η 0 , j 1 η o m . (16)

Величина η в астрономии называется уравнением времени. Поэтому истинное солнечное время ma в часах будет выражаться через среднее солнечное время m так:

m a = m + η , (17)

где m ? отсчитывается от полуночи.

График для уравнения времени η(Td) приведен на Рис. 3. Отклонения времени η для дня весеннего равноденствия j1 = 1 равно η = −7.47 мин,

Рис. 3. Уравнение времени в современную эпоху 30.12.1949 г.: η―в мин; Td―в днях от момента весеннего равноденствия.

наименьшее значение ηmn = −14.25 мин при j1 = 329, а наибольшее ηmx = 16.43 мин при j1 = 229. Величины отклонения η не отличаются от известных в астрономии [24] . Отличие имеется в начале отсчета времени Td в днях: в астрономии дни отсчитываются с 1 января, а на Рис. 3―с момента весеннего равноденствия.

Следует отметить, что по величине отклонения времени η можно также определить отклонение истинных солнечных суток от средних

Δ T s d , j 1 = η j 1 + 1 η j 1 . (18)

Клавдий Птолемей использовал равноденственные часы [25] . Как показано выше, весенние равноденственные сутки больше среднесолнечных на ΔTsd = 0.297 мин. Если он и предшествующие ему астрономы использовали равноденственное время, оно может отличаться от среднесолнечного времени на величину порядка 10 минут. Роберт Ньютон [26] обвинил Клавдия Птолемея в искажении моментов древних наблюдений, потому что они не совпадают с современной теорией примерно на такое же количество минут. Как видим, одной из причин несовпадения может быть разная длительность использованных часов времени.

Длительность солнечных суток, согласно (12), определялась разностями долгот λj1 ? λj1−1 за средние солнечные сутки Tsdm = 24 часа. Был рассчитан скорректированный ряд долгот λс,j1 по фактической длительности суток Tsd,j1 и повторены расчеты по формулам (9)-(16). Наибольшее отличие скорректированной длительности суток ΔTsdc от ΔTsd равно 0.092 сек при j1 = 278. Это составляет относительную погрешность отклонения длительности солнечных суток 0.3%. Поэтому алгоритм расчета (9)-(16) для длительности суток Tsd и уравнения времени η можно использовать без коррекции.

На Рис. 4 изменение длительности солнечных суток ΔTsd в современную эпоху (T = 0) сопоставлено с этими изменениями в другие четыре эпохи. Как видно время экстремумов ΔTsd примерно совпадает и приходится на моменты равноденствий (Td ≈ 0 и Td ≈ 180) и солнцестояний (Td ≈ 90 и Td ≈ 280). При этом изменение длительности солнечных суток происходят в пределах от −0.5 мин до 0.35 мин от средней длительности в 24 часа.

В табл. 2 основные параметры уравнения времени η, представленного на

Рис. 4. Отклонение ΔTsd длительности солнечных суток в эпохи T от средних за год в эти же эпохи: ΔTsd―в минутах; T―время в тыс. лет от 30.12.1949 г.

Таблица 2. Изменения параметров уравнения времени по эпохам T: T―время в тыс. лет от 30.12.1949 г.

Рис. 3 для современной эпохи (T = 0), сопоставлены с таковыми для других четырех эпох. Среднее отклонение η0m, согласно (15), как видно из Таблица 2 изменяется в широких пределах: от −4.7 мин в эпоху 31 т.л.н. до 7.8 мин в современную эпоху. Максимальное отклонение ηmx и минимальное ηmn также изменяются в 2 - 3 раза. В то же время средняя разность прямых восхождений Δαm за одни сутки, согласно (10), остается неизменной Δαm = 0.0172. Это является следствием неизменности периодов орбитального и вращательного движения Земли.

На Рис. 2 мы рассматривали вращения Земли со средней угловой скоростью ω E . Угловая скорость Земли имеет свои независимые колебания [3] [10] , но их амплитуда на три порядка меньше амплитуды изменения солнечных суток. Поэтому колебания угловой скорости вращения Земли здесь не учитываются.

5. Долгота светового дня

5.1. Долгота светового дня в течение года

Длительность светового дня, или долгота дня D0 в часах определяется (см. Рис. 2) суммой часовых углов ω0 восходов и заходов центра Солнца (8):

D 0 = ( 24 / 2 π ) 2 ω 0 = ( 24 / π ) ω 0 . (19)

Так как часовой угол ω0 отсчитывается до центра Солнца, то величина D0 определяет длительность нахождения над горизонтом HH¢ центра Солнца. Световой день начинается с появления края Солнца и заканчивается с опусканием его за горизонт. Кроме того, из-за уменьшения плотности атмосферы с высотой происходит преломление (рефракция) светового луча так, что Солнце видно на угол ρ90 под горизонтом. В общем случае для наблюдаемых на небе объектов отклонение ρ90 при зенитном угле z = 90˚, или рефракция, зависит от ряда факторов, которыми для Солнца можно пренебречь. Поэтому в астрономии рефракция принята постоянной и равно ρ90 = 34' [24] .

На Рис. 2 часовой угол z показан между центром Солнца и зенитом Z наблюдателя M, находящегося на поверхности Земли. А формулой (7) угол z определяется от центра Земли. Угол p, под которым расстояние наблюдателя от центра Земли видно с Солнца, называется его параллаксом. Положение Солнца на небесной сфере будет отличаться на этот угол для двух наблюдателей, находящихся в точке зенита и в точке восхода, например. На угол p необходимо уменьшить рассчитанную по (7) величину зенитного угла. Наибольшая величина параллакса p не превышает 8''.8 [24] . Эта величина в 230 раз меньше величины рефракции ρ90, поэтому параллаксом можно пренебречь.

Имеется еще один источник искажения положения Солнца на небосводе (см. Рис. 5), который обусловлен движением наблюдателя M относительно Солнца S со скоростью v . Это явление называется аберрацией света.

При движении наблюдателя относительно источника света S происходят изменения всех характеристик света. Точное выражение для измененных характеристик света приведено в работе [27] . При скорости движения v, значительно меньшей скорости света c для угла смещения источника справедливо выражение

σ = β sin θ o b , (20)

где β = v/c

Истинное положение источника S отличается от наблюдаемого положения Sob на угол σ. Упрощенно, как это принято в астрономии [24] , можно рассмотреть движение наблюдателя M относительно Солнца за счет суточного вращения Земли vrt и за счет ее орбитального движения (vor) относительно Солнца. При суточном вращении в точках восхода и захода Солнца углы θ между ним и скоростью vrt близки к 0 и π. Поэтому, согласно

Рис. 5. Угол θob видимого положения Солнца Sob при движении наблюдателя M со скоростью v и угол θ действительного положения солнца S.

(20), угол смещения σ будет близкий к нулю. За счет орбитального движения в этих же точках угол θ между Солнцем и скоростью vor близок к π/2. В этом случае угол наблюдаемого Солнца будет на величину σ меньше, чем угол истинного. Это будет приводить к более раннему восходу Солнца и более раннему его заходу. Так что на длительность светового дня это явление не оказывает влияние. При этом угол смещения σ не превышает 20'' [24] . Эта величина в 100 раз меньше величины рефракции ρ90, а также меньше неопределенности величины ρ90. Поэтому аберрацией света можно пренебречь.

Следует отметить, что в силу пренебрежения аберрацией света такое рассмотрение этого явления по составляющим скорости vrt и vor допустимо. Однако при точном его рассмотрении, необходимо рассматривать полный вектор скорости, его три угла наклона к линии MSob (Рис. 5). Затем вычислять изменение этих углов по точным формулам в работе [27] .

Теперь рассмотрим восходы и заходы верхнего края Солнца. На Рис. 6 показано смещение края диска Солнца на угловое расстояние GK = ρ90 ниже уровня горизонта HH¢ в момент его захода. Видимый радиус Солнца KSs1 согласно [24] равен:

K S s 1 = 961 " .18 ρ = 4.6599 × 10 3 ρ , (21)

где ρ ? расстояние Солнца от Земли, отнесенное к большой полуоси a орбиты Земли. Величина ρ рассчитывается по уравнению траектории в зависимости от долготы λ [21] [22] .

С учетом (21) дуга GSs1 будет равна в радианах

G S s 1 = 10 3 ( 9.8902 + 4.6599 ρ ) . (22)

Дуга GSs1 расположена на окружности большого круга, проходящего через точку зенита Z (Рис. 2), поэтому она является приращением зенитного угла, т.е. Δz = GSs1. Тогда зенитный угол видимого края заходящего Солнца будет

z 0 a = π / 2 + G S s 1 . (23)

Так как часовой угол ω связан с зенитным углом z выражением (7), то часовой угол видимого захода края Солнца можно записать так:

cos ω 0 a = cos z 0 a sin δ sin φ cos δ cos φ . (24)

Рис. 6. Аберрационное смещение GK = ρ90 края диска Солнца в момент его захода: HH¢―линия горизонта; Ss1―положение центра Солнца во время захода его края.

Так как правая часть по модулю может превышать 1, то перепишем выражение (24) в виде:

ω 0 a = arcos ( F n ) , (25)

где функция Fn определяется в зависимости от промежуточной функции

F n 0 = cos z 0 a sin δ sin ϕ cos δ cos ϕ (26)

так:

F n = F n 0 при 1 F n 0 1 ; F n = 1 при F n 0 < 1 ; F n = 1 при F n 0 > 1 . (27)

Как следует из Рис. 2, зенитный угол видимого края восходящего Солнца будет z0ra = −π/2 − GSs1, т.е. по модулю такой же, как и z0a. Поэтому модуль часового угла восхода Солнца будет определяться также выражением (25). Тогда долгота светового дня в радианах равна 2∙ω0a, а в часах запишется так:

D = 24 ω 0 α / π . (28)

На Рис. 7 показано изменение долготы дня D в течение года на разных широтах. Например, на широте φ = 60˚ в момент весеннего равноденствия (Td = 1) D = 12.31 часа, затем она увеличивается до 18.87 часа в момент летнего солнцестояния. В момент осеннего равноденствия D = 12.26 часа, а затем уменьшается до D = 5.87 часа в день зимнего солнцестояния. С уменьшением широты φ экстремумы долготы дня D уменьшаются и приближаются к 12.1 часа на экваторе (φ = 0˚).

С увеличением широты φ > 60˚ наибольшая долгота дня в день летнего солнцестояния приближается к 24 часам. На полярном круге (φ = 90˚ − ε˚) наступает полярный день, который длится 31 сутки. С увеличением широты длительность полярного дня увеличивается до половины года на полюсе. В день зимнего солнцестояния долгота дня с увеличением широты уменьшается до D = 2.17 часа на полярном круге. С дальнейшим увеличением широты наступает полярная ночь.

Представленная на Рис. 7 картина изменения долготы дня совпадает с

Рис. 7. Продолжительность (долгота) светового дня в году на разных широтах φ Северного полушария в современную эпоху 30.12.1949 г.

таковой в астрономии [24] .

5.2. Полярные дни и ночи: моменты наступления и длительность

Рассмотрим моменты наступления полярных дней и ночей и их длительность. Для этого необходимо рассчитать долготы λ этих событий. Как уже отмечалось, полярному дню соответствует долгота (длительность) дня Dd = 24 часа, а полярной ночи Dn = 0 часов. Тогда из (28) получаем часовые углы полярного дня ω0ad =π и полярной ночи − ω0an = 0. После подстановки в (25) вместо ω0a этих значений часового угла, получаем Fnd = −1 для полярного дня и Fnn = 1 для полярной ночи. Здесь функция Fn определяется выражением (26). Перепишем его в следующем виде:

cos z 0 a sin δ sin φ = F n cos δ cos φ .

Подставив в это выражение значение sinδ, согласно (4), получим следующее уравнение для долготы λ

cos z 0 a sin ε sin λ sin φ = F n 1 sin 2 ε sin 2 λ cos φ . (29)

После возведения (29) в квадрат и преобразований с учетом того, что для рассматриваемого случая Fn2 = 1, получаем для sinλ квадратное уравнение

sin 2 λ 2 sin φ cos z 0 a sin ε sin λ cos 2 φ cos 2 z 0 a sin 2 ε = 0. (30)

Оно имеет два решения, отличающиеся знаками “+” и “−” перед квадратным корнем. Запишем это решение в виде одной формулы:

sin λ = F n 01 ; F n 02 , (31)

где функция со знаком “+” имеет вид:

F n 01 = [ ( sin φ cos z 0 a ) 2 + cos 2 φ cos 2 z 0 a + sin φ cos z 0 a ] / sin ε . (32)

Функция Fn02 перед корнем имеет знак “−”.

Чтобы исключить значения по модулю превышающие 1 вводятся функции

F n 1 = F n 01 при 1 F n 01 1 ; F n 1 = 1 при F n 01 > 1 ; (33)

F n 2 = F n 02 при 1 F n 02 1 ; F n 2 = 1 при F n 02 < 1 . (34)

Летнее полугодие характеризуется долготами 0 < λ < π, а зимнее ?π < λ < 2π. Так как функция Fn1 > 0, то она будет давать летние долготы, а отрицательная функция Fn2―зимние долготы. Тогда в соответствии с (31) летние долготы полярных дней будут λpd = arcsin(Fn1). Эта функция будет давать два решения: для первого λpd и второго квадранта λpd1. Как следует из Рис. 7 меньшее значение соответствует началу полярного дня, а большее―концу. Окончательно долготы начала и конца полярных дней запишутся в виде:

λ p d = arcsin ( F n 1 ) ; λ p d 1 = π arcsin ( F n 1 ) . (35)

Аналогично долготы начала и конца полярных ночей определятся так:

λ p n = + arcsin ( F n 2 ) ; λ p n 1 = π arcsin ( F n 2 ) . (36)

Чтобы найти моменты по времени Td, соответствующие долготам λpd, λpd1, λpn, λpn1, необходимо найти индексы j1 соседних долгот. Затем интерполяцией определяются моменты: Tdd и Tdd1―начала и конца полярного дня, а также Tdn и Tdn1―начала и конца полярной ночи. Алгоритм вычисления приведен в п. 11 программы SunFhnm.mcd в Приложении, а на Рис. 8 они представлены в виде графиков. В области I приблизительно до Td = 90 дней изображены начала полярных дней в зависимости от широты φ. Разными линиями и точками показаны графики для 5 эпох. Аналогично в области II приблизительно при 90 < Td < 180 дней этими же линиями показаны моменты Td окончания полярных дней. Аналогично приблизительно при 180 < Td < 270 и Td > 270 изображены начала (III) и окончания (IV) полярных ночей.

Например, в современную эпоху T = 0 на широте φ = 70˚ в Td = 56.88 день от момента весеннего равноденствия наступает полярный день. Он заканчивается в Td = 128.82 день. Полярная ночь наступает в Td = 250.05 день и заканчивается в Td = 302.39 день. С увеличением широты φ начала полярных дней и ночей происходит раньше, а окончание позже. С уменьшением широты φ начала полярных дней отдаляются и приближаются к Td = 92.8 дню, а начала полярных ночей приближаются к Td = 275.5 дню. Это происходит на широтах близких к широте полярного круга, которая для современной эпохи равна 66.56˚.

В другие эпохи графики начала и конца полярных дней и ночей, как видно из Рис. 8, идентичны, но широта их начала может существенно смещаться. При этом изменяется продолжительность полярных дней и ночей. Например, в эпоху 31 т. л. н. широта начала полярных дней и ночей

Рис. 8. Моменты Td наступления (I, III) и окончания (II, IV) полярных дней (I, II) и ночей (III, IV) на разных широтах Северного полушария в экстремальные эпохи за последние 50 т.л.н.

смещается до56˚ и 58˚, соответственно, а в эпоху 46.44 т. л. н.―до 74˚ и 76˚ соответственно.

Следует отметить, что древние арии [15] в эпоху 31 т. л. н. могли обитать на территории близкой по широте г. Тюмени (φ = 57.15˚) и наблюдать полярные дни и ночи. С приближением последнего ледникового периода 15.32 т. л. н. они вынуждены были мигрировать на южные территории.

Начала полярных дней и ночей и их продолжительность ΔTdd = Tdd1 − Tdd и ΔTdn = Tdn1 − Tdn соответственно, на разных широтах рассчитаны для пяти разных эпох. Они представлены в табл. 1П в Приложении работы [28] , которая имеется в свободном доступе. Для сравнения в этой таблице приведены длительности полярных дней и ночей в астрономии в современную эпоху. Они обозначены как EA по названию “Эфемеридная астрономия” источника [24] . Данные EA приведены с точностью 1 день. Округление начала и конца полярных дней и ночей может приводить к разности их длительности в 2 дня. В табл. 1П такое отличие наблюдается. А в целом представленные на Рис. 8 и табл. 1П результаты для современной эпохи совпадают с известными данными в астрономии.

5.3. Долготы светового дня в моменты солнцестояний

Рассмотрим долготу светового дня в моменты летнего и зимнего солнцестояний в неполярных широтах. Как уже отмечалось, наибольшая долгота светового дня будет в момент летнего солнцестояния λs = π/2, а наименьшая − в момент зимнего солнцестояния λw = 1.5∙π. Тогда склонение Солнца δ, согласно (4), в этих случаях будет δs = ε и δw = −ε, соответственно. Тогда после подстановки угла δs в (26) часовой угол видимых восходов и заходов края Солнца в момент летнего солнцестояния запишется аналогично (25)-(27) так:

ω 0 s a = arcos ( F n s ) ; (37)

F n 0 s = cos z 0 a 1 sin ε sin ϕ cos ε cos ϕ ; (38)

F n s = F n 0 s при 1 F n 0 s 1 ; F n s = 1 при

F n 0 s < 1 ; F n s = 1 при F n 0 s > 1 . (39)

В формуле (30) использован зенитный угол z0a1, в котором безразмерное расстояние от Земли до Солнца ρ, согласно (22), из-за малого влияния его изменения принято ρ = 1.

После подстановки угла δw = −ε в (26) часовой угол видимых восходов и заходов края Солнца в момент зимнего солнцестояния запишется аналогично (37)-(39) так:

ω 0 w a = arcos ( F n w ) ; (40)

F n 0 w = cos z 0 a 1 + sin ε sin φ cos ε cos φ ; (41)

F n w = F n 0 w при 1 F n 0 w 1 ; F n w = 1 при

F n 0 w < 1 ; F n w = 1 при F n 0 w > 1 . (42)

Тогда долгота светового дня в моменты летнего и зимнего солнцестояний запишутся, соответственно, так:

D s = 24 ω 0 s α / π ; D w = 24 ω 0 w α / π . (43)

Результаты расчетов долготы дня в моменты солнцестояний по формулам (37)-(43) на Рис. 9 даны для пяти разных эпох. Графики представляют изменение долготы дня D по широтам φ северного полушария. Для зимнего солнцестояния долгота дня D ≤ 12 часов отмечена как область I, а для летнего солнцестояния―как область II.

В современную эпоху, как видно из Рис. 9, на широте 50˚ в зимнее солнцестояние день длится 8 часов, а в летнее―16.37 часа. С уменьшением широты день в зимнее солнцестояние увеличивается и приближается к 12 часам на экваторе. А в летнее солнцестояние день уменьшается и на экваторе также приближается к 12 часам. С увеличением широты φ долгота зимнего дня уменьшается и приближается к нулю. То есть наступает полярная ночь, а долгота летнего дня увеличивается до 24 часов и наступает полярный день.

Для остальных эпох изменения долготы дня изменяются подобным образом, но с разными наклонами зависимостей D(φ). Например, для эпохи 31 т. л. н. полярные ночи и дни начинаются на меньших широтах, а в эпоху 46.44 т. л. н.―на больших широтах.

В табл. 2П работы [28] приведены долготы дней для летнего солнцестояния (Ds) и зимнего (Dw) для пяти разных эпох. Изменение широты φ дается через 2˚, начиная с экватора. Широты изменяются до величин, после которых наступает полярный день (Ds = 24 часа) и полярная ночь (Dw = 0 часов).

Рис. 9. Долгота светового дня D в моменты зимнего солнцестояния (I) и летнего солнцестояния (II) на разных широтах φ˚ Северного полушария в экстремальные эпохи T за последние 50 т.л.н.: T―время в тыс. лет от 30.12.1949 г.; D―долгота светового дня в часах.

6. Азимуты Солнца

6.1. Азимуты движения Солнца в течение суток

Азимут Солнца AS отсчитывается на круге горизонта HH' (Рис. 2) от точки севера Nrd до проекции Солнца в т. D за часовой стрелкой

A S = N r d S t h + S t h D = 2 π D N r d . (44)

В треугольнике NrdZD, две стороны которого NrdZ и ZD равны π/2, дуга DNrd равняется углу Z. В треугольнике NZS известны две стороны: NS = π/2 − δ и ZS = z, а также угол N = ω. По теореме синусов

sin Z / sin N S = sin N / sin Z S ,

находим

sin Z = sin ω cos δ / sin z , (45)

где дуга z является зенитным углом Солнца и согласно (7) определяется выражением:

z = arcos ( sin δ sin φ + cos δ cos φ cos ω ) . (46)

Тогда с учетом (45) азимут Солнца (44) AS = 2π − Z запишется в радианах так:

A S , r a d = 2 π arcsin ( sin ω cos δ / sin z ) . (47)

Так как функция arcsin в (47) неоднозначна и после достижения экстремальных значений ±π/2 имеет изломы, то проводятся следующие операции. Вводится функция

F n 03 = arcsin ( sin ω cos δ / sin z ) . (48)

Находятся максимальные значения mxFn функции Fn03 и минимальные ее значения mnFn. При изменении часового угла ω от −π до π определяются индексы часового угла in9, соответствующего mnFn, и in 10, соответствующего mxFn. В зависимости от этих индексов рассчитываются диапазоны индекса k часового угла ωk и значения функции Fn3 так:

F n 3 = F n 03 при i n d 9 k i n 10 ; F n 3 = π F n 03 при k < i n 9 ;

F n 3 = π F n 03 при k > i n 10 . (49)

Как видно из (48) при ω = 0, в полдень, функция Fn03 может быть равной 0 или π. Согласно Рис. 2 в полдень Fn03 = NrdD = π. Тогда с учетом этого и (48)-(49) азимут Солнца в градусах в соответствии с (47) запишется так:

A S = 18 0 + ( 18 0 / π ) F n 3 . (50)

Следует отметить, что минимальное и максимальные значения Fn03 могут принадлежать не центральному участку. В этом случае на азимуте AS при большом интервале между делениями ω будут всплески при переходе через 0.5・π и 1.5・π. В этом случае нужно использовать следующее условие:

Fn3 = Fn03 при ind9 < k < in10; Fn3 = −π − Fn03 при k ≤ in9; Fn3 = π − Fn03 при k ≥ in10.

Этот алгоритм реализован в п. 13 программы SunPhenmn.mcd. На Рис. 10 приведены азимуты движения Солнца по небосводу на широте Москвы в современную эпоху. Линией 1 показано изменение азимута в день весеннего равноденствия. Азимут Солнца в полночь ω = −12 часов равен нулю и увеличивается на рассвете (ω ≈ −6 часов) до 90˚. В течение дня продолжает увеличиваться: в полдень (ω = 0 часов) до 180˚ и в полночь (ω = 12 часов) достигает значения 360˚. Вертикальными линиями показаны часы восхода и захода Солнца.

В дни солнцестояний: летнего (линия 2) и зимнего (линия 3) изменение азимутов происходит аналогичным образом, но границы дня существенно изменяются. Кроме того, по сравнению с изменением азимута в дни равноденствия, в эти дни характер изменения азимутов более существенно отличается от линейного закона. Нелинейный характер изменения AS свидетельствует, что определение промежутков времени по азимуту Солнца или по тени гномона, будет характеризоваться соответствующей неравномерностью хода времени. На Рис. 10 часовые промежутки нанесены кружками. Так как в дни равноденствий азимут AS 1 изменяется наиболее близко к линейному закону, то равноденственные часы наиболее равномерны.

Условия (49) при вычислении азимута Солнца справедливы для нетропических широт. В тропических широтах Северного полушария широта φ меньше угла наклона ε, т.е. φ − ε < 0. Для зенитного угла z в полдень (ω = 0) из (7) следует

z n = φ δ . (51)

Поэтому, как уже отмечалось, в день летнего солнцестояния (δ = ε см.

Рис. 10. Азимуты AS на широте Москвы (j˚ = 55.7522˚) суточного движения Солнца в дни весеннего равноденствия (1), летнего (2) и зимнего (3) солнцестояний в современную эпоху 30.12.1949 г. Вертикальными линиями отмечены часовые углы восходов (ω < 0) и заходов Солнца (ω > 0). Кружками 4 отмечены часовые промежутки времени.

Рис. 2) зенитный угол в тропических широтах изменяется от

z n s = φ ε (52)

до

z n w = φ + ε (53)

в день зимнего солнцестояния (δ = −ε). Из (52) видно, что в день летнего солнцестояния зенитный угол zns становится отрицательным, т.е. для наблюдателя M (см. Рис. 2) Солнце находится в северной стороне неба. В этом случае в полдень (ω = 0) азимут Солнца ASn = 0, а не π, как в нетропических широтах. Следует отметить, что в этом случае суточное движение Солнца для наблюдателя M происходит на южной стороне неба справа налево, а не как на северной стороне слева направо, как это наблюдается в нетропических широтах.

Ситуации с отрицательным углом zn, согласно (51), имеют место не только для дня летнего солнцестояния, но и для других дней. Поэтому необходимо определять номера дней, для которых zn < 0. Тогда для рассмотренного случая тропических широт азимут Солнца запишется так:

ASt = (180/π)・Fn03 при zn < 0. (54)

Полностью алгоритм расчета азимута Солнца приведен в п. 13 программы SunPhenmn.mcd в Приложении. На Рис. 11 представлено изменение азимута Солнца в течение суток на широте φ = 20˚. В день весеннего равноденствия 1 и зимнего солнцестояния 3 азимут Солнца изменяется от 180˚ в полночь (ω = −12 часов) до некоторого минимального значения. Для зимнего солнцестояния 3 оно приходится на ночь, а для весеннего равноденствия―на рассвет. В полдень азимут Солнца снова равняется 180˚, т.е. Солнце находится в южной стороне неба. При движении Солнца от полдня к полночи оно проходит максимальные значения азимута и в полночь снова приходит на юг (AS = 180˚).

В день летнего солнцестояния (линия 2, правая шкала графика) Солнце в

Рис. 11. Азимуты в тропических широтах (j˚ = 20˚) суточного движения Солнца в дни весеннего равноденствия 1, летнего 2 и зимнего 3 солнцестояний в современную эпоху 30.12.1949 г. Шкала азимута дня летнего солнцестояния 2 сдвинута на 180˚ (см. справа графика). Вертикальными линиями отмечены часовые углы восходов (ω < 0) и заходов Солнца (ω > 0).

полночь находится на Севере (AS = 0˚), затем движется по северной стороне неба (AS < 90˚). Перед полднем (ω ≈ 1 час) достигает максимального отклонения к востоку и в полдень возвращается на север (AS = 0˚). После полдня Солнце переходит в западную сторону неба и достигает максимального отклонения к западу (AS ≈ 70˚). К полночи Солнце снова возвращается на север (AS = 0˚). Следует отметить, что здесь мы использовали для определения азимута шкалу не от 0 до 360˚, а от −180˚ до 180˚. Это позволяет избежать разрыва графика на интервале часового угла ω от −12 часов до +12 часов.

В нетропических широтах суточное движение Солнца происходит по полному кругу вокруг наблюдателя M, т.е. AS = 0 ÷ 2π. В тропических широтах суточный круг движения Солнца происходит в стороне от наблюдателя: летом (линия 2) в северной стороне неба, зимой, включая дни равноденствия, (линии 1 и 3) ? в южной стороне неба.

6.2. Проекция дневного пути Солнца на плоскости горизонта

Рассмотрим проекцию на плоскость горизонта дневное движение Солнца. Положение Солнца S (Рис. 2) над плоскостью горизонта HH' определяется зенитным углом z и углом азимута AS, который определяется дугой NrdSrSthD. На плоскости горизонта введем две оси координат с началом в точке M: NS ? направленная на Север Nrd и ES ? направленная на Восток Est. Тогда проекции центра Солнца на эти оси запишутся так:

N S = sin z cos A S ; E S = sin z cos A S , (55)

где AS ? азимут Солнца в радианах.

Выражением (55) также определяются координаты NS, ES точек восхода и захода верхнего края Солнца, если углами z и AS будут определены их моменты времени. На Рис. 12 проекции дневного пути Солнца показаны для трех случаев: в день весеннего равноденствия и в дни солнцестояний на широте Москвы.

В день весеннего равноденствия центр Солнца всходит на Востоке и заходит на Западе. А восход и заход верхнего края Солнца немного сдвигается на Север. Летом восход и заход Солнца находится в северной части горизонта. С увеличением широты φ траектория 2 Солнца будет приближаться к окружности. На широте, где траектория станет окружностью, наступит полярный день.

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

Следует отметить, что конечные точки на траектории Солнца в дни равноденствий 1 и в дни солнцестояний 2 и 3 играют ключевую роль в древных астрономических календарях [17] [18] . Представленный в п. 15 программы на SunPhnmen.mcd алгоритм позволяет рассчитать их для любой эпохи.

Рис. 12. Проекция дневного пути Солнца на плоскости горизонта на широте Москвы (j = 55.7522˚) в современную эпоху 30.12.1949 г.: 1―в дни весеннего равноденствия; 2―летнего и 3―зимнего солнцестояний. Ось NS направлена на Север (Nrd), а ось ES―на восток (Est). Положение наблюдателя (M) в начале координат отмечено крестиком. Точками на линиях отмечены восходы верхнего края Солнца (с правой стороны) и его заходы (с левой стороны).

6.3. Азимуты восходов и заходов Солнца

В отличие от азимутов суточного движения Солнца азимуты точек восхода и захода необходимо отсчитывать до точек появления и исчезновения верхнего края Солнца. Как видно из Рис. 2 и Рис. 6, азимут наблюдаемой точки G захода края Солнца состоит из дуг на круге горизонта HH¢

A G s = N r d E s t S t h W s t + W s t S s + S s G , (56)

где Est и Wst ? точки Востока и Запада на круге горизонта.

Дуга NrdEstSthWst равна 1.5∙π (см. Рис. 2). Через точку Ss и полюс N проведена дуга большого круга SsF, поэтому она перпендикулярна дуге экватора AA', т.е. угол F = 0.5∙π. Дуга SsF равна δ. В треугольнике SsFWst угол Wst = π/2 − φ. Поэтому дуга WstSs определяется по теореме синусов

W s t S s = arcsin ( sin δ / cos φ ) . (57)

Сторона GSs может быть найдена из треугольника GSs1Ss. (см. Рис. 6), в котором угол Ss = Wst = π/2 − φ. В виду малости треугольника, его можно рассматривать как плоский и линейный со сторонами RE・GSs1 и RE・GSs, где RE―радиус Земли. Для угла Ss можно записать: tg(π/2 − φ) = RE・GSs1/(RE∙GSs). Отсюда получаем:

G S s = G S s 1 tg φ . (58)

Подставляя в (56) составляющие получим азимут точки захода Солнца в градусах в виде:

A G s = ( 18 0 / π ) ( 1.5 π + arcsin ( sin δ / cos φ ) + G S s 1 tg φ ) . (59)

Ввиду того, что sinδ/cosφ по модулю может превышать 1, введем обозначение:

F n 04 = sin δ / cos φ (60)

и определим функцию:

F n 4 = F n 04 при 1 F n 04 1 ; F n 4 = 1

при F n 04 < 1 ; F n 4 = 1 при F n 04 > 1 . (61)

Тогда азимут точки захода Солнца будет:

A G s = ( 18 0 / π ) ( 1.5 π + arcsin ( F n 4 ) + G S s 1 tg φ ) . (62)

Чтобы рассчитать азимут точки восхода, обозначим точку Sr1 видимого восхода края Солнца на дуге MnSr (на Рис. 2 эта точка не приведена). Она находится под горизонтом HH¢ и дуга Sr1Sr равна по величине дуге Ss1Ss (см. Рис. 6), т.е. Sr1Sr = Ss1Ss. Поэтому точка наблюдаемого восхода края Солнца Gs будет ближе к точке севера Nrd на величину дуги SrG1 = SsG. Тогда азимут наблюдаемой точки Gs восхода будет определяться следующими дугами:

A G r = N r d E s t E s t S r S r G r , (63)

где точка Est находится в пересечении кругов горизонта HH¢ и экватора AA', а ее азимут NrdEst = π/2. Дуга EstSr = WstSs. Подставляя в (63) все составляющие азимут точки восхода края Солнца с учетом (60)-(61) получаем в виде:

A G r = ( 18 0 / π ) ( 0. arcsin ( F n 4 ) G S s 1 tg φ ) . (64)

Как видим из (62) и (64) сумма азимутов восходов и заходов

A G s + A G r = 36 0 . (65)

Сумма AGs + AGr может отличаться от 360˚, если горизонт в восточной части отличается от горизонта в западной части высотой окружающей местности. Приведенные формулы для азимутов не учитывают этой особенности и дают симметричные азимуты восходов и заходов относительно направления MNrd на север (Рис. 2).

На Рис. 13 показано изменение азимута точки восхода края Солнца AGr в течение года на разных широтах. Например, на широте φ = 60˚ в момент весеннего равноденствия (Td = 1) азимут AGr = 87.76˚, затем он уменьшается до 35.84˚ в день летнего солнцестояния. В момент осеннего равноденствия AGr = 88.23˚, а затем увеличивается до AGr = 141.27˚ в день зимнего солнцестояния.

С уменьшением широты φ экстремумы азимута Солнца AGr уменьшаются и на экваторе (φ = 0˚) приближаются к 66.56˚ в день летнего солнцестояния и к 113.44˚ в день зимнего солнцестояния. С увеличением широты φ > 60˚ наименьший азимут точки восхода края Солнца AGr в день летнего солнцестояния приближается к 0˚, т.е. Солнце восходит на Севере. На полярном круге (φ = 90˚ − ε˚) наступает полярный день. С увеличением широты в день зимнего солнцестояния (Td = 276) азимута Солнца AGr увеличивается до 180˚ на полярном круге (φ = 90˚ − ε˚), т.е. Солнце всходит на Юге. С дальнейшим увеличением широты наступает полярная ночь.

Рис. 13. Азимуты восходов верхнего края Солнца в современную эпоху 30.12.1949 г. в течение года на разных широтах φ˚ Северного полушария: AGr―азимуты восходов Солнца в градусах.

Следует отметить, что графики на Рис. 13 напоминают графики долготы светового дня на Рис. 7, если их перевернуть вокруг горизонтальной линии AGr = 90˚. Азимуты точки захода края Солнца AGs, как следует из (65), AGs = 360˚ − AGr, поэтому их графики будут похожи на графики долготы светового дня на Рис. 7.

6.4. Экстремальные азимуты восходов и заходов Солнца в дни солнцестояний

Наибольшие азимуты заходов и наименьшие азимуты восходов верхнего края Солнца имеют место в день летнего солнцестояния λs = π/2. В день зимнего солнцестояния λs = 1.5∙π экстремумы азимутов становятся обратными: азимуты заходов―наименьшие, а азимуты восходов―наибольшие. Приведем предельные азимуты восходов (64) в день летнего солнцестояния:

A G r s m = ( 18 0 / π ) ( 0.5 π arcsin ( sin ε / cos φ ) ( 9.8902 × 10 3 + 4.6599 × 10 3 / ρ s ) tg φ ) (66)

и в день зимнего солнцестояния:

A G r w n = ( 18 0 / π ) ( 0.5 π + arcsin ( sin ε / cos φ ) ( 9.8902 × 10 3 + 4.6599 × 10 3 / ρ w ) tg φ ) (67)

Азимуты заходов Солнца (62) в дни солнцестояний запишутся аналогичным образом. При вычислениях по формулам (66)-(67) в алгоритм необходимо ввести условия, чтобы избежать |sinε/cosφ| > 1 (см. п. 14 программы в Приложении).

Результаты расчетов азимуты восходов Солнца в дни солнцестояний по формулам (66) - (67) для пяти разных эпох даны на Рис. 14. Графики представляют изменение азимута AGr по широтам φ Северного полушария. Для летнего солнцестояния азимуты AGr ≤ 90˚ отмечены как область I, а для зимнего солнцестояния―как область II.

В современную эпоху, как видно из Рис. 14, на широте 50˚ в летнее

Рис. 14. Азимуты восходов верхнего края Солнца в моменты летнего солнцестояния (I) и зимнего солнцестояния (II) на разных широтах φ˚ Северного полушария в экстремальные эпохи T за последние 50 т.л.н.: T―время в тыс. лет от 30.12.1949 г.; AGr―азимуты восходов Солнца в градусах.

солнцестояние азимут восхода Солнца AGr равен 50.77˚, а в зимнее − 127.24˚. С уменьшением широты азимут восхода увеличивается и приближается к 66.56˚ на экваторе. А в зимнее солнцестояние азимут восхода уменьшается и на экваторе приближается к 113.44˚. С увеличением широты φ азимут восхода AGr дня летнего солнцестояния уменьшается и приближается к нулю. То есть Солнце всходит на Севере и наступает полярный день. А азимут восхода AGr дня зимнего солнцестояния с увеличением широты φ увеличивается до 180˚и наступает полярная ночь.

Для остальных эпох изменения азимута восхода Солнца AGr происходит подобным образом, но с разными наклонами зависимостей AGr(φ). Например, для эпохи 31 т. л. н. полярные ночи и дни начинаются на меньших широтах, а в эпоху 46.44 т. л. н. − на больших широтах.

Следует отметить, что графики на Рис. 14 напоминают графики долготы светового дня D в разные эпохи на Рис. 9, однако области I и II у них разные. Азимуты точки захода края Солнца, например, в день летнего солнцестояния AGss, как следует из (65), AGss = 360˚ − AGrsm. Поэтому графики AGssm и AGswn будут похожи на графики долготы светового дня на Рис. 9.

В табл. 3П работы [28] приведены азимуты восходов Солнца для летнего солнцестояния (AGrsm) и зимнего (AGrwn) для пяти разных эпох. Изменение широты φ дается через 2˚, начиная с экватора. Широты изменяются до величин, после которых наступает полярный день (AGrsm = 0˚) и полярная ночь (AGrwn = 180˚).

7. Солнечная тень гномона

7.1. Относительная длина тени гномона

Гнóмоном (др.-греч. − указатель) называется предмет, тень от которого используется в солнечных и лунных часах и календарях. На Рис. 15 показана длина тени lsh вертикального гномона длиной lg, которая возникает от верхнего края Солнца. Величина lsh зависит от зенитного угла z Солнца и его углового радиуса S. Ранее формулой (21) величина S приводилась в виде дуги KSs1 на Рис. 6. Как следует из Рис. 15, относительная длина тени гномона будет:

l s h 1 = l s h / l g = tg ( z S ) = tg ( z 4 . 6599 × 10 3 / ρ ) , (68)

где зенитный угол z Солнца определяется выражением (46).

Чтобы избежать отрицательных значений относительной длины теней lsh1 расчет проводится по следующему алгоритму:

F n 06 = z 4.6599 × 10 3 / ρ ; (69)

F n 6 = F n 06 при 0. 5 π F n 0 0. 5 π ;

F n 6 = 0. 5 π при F n 06 < 0. 5 π ;

F n 6 = 0. 5 π при F n 06 > 0. 5 π ; (70)

l s h 10 = tg ( F n 6 ) ; l s h 1 = l s h 10 при l s h 10 > 0 ; l s h 1 = 0 при l s h 10 < 0 . (71)

7.2. Изменение длины тени в течение дня

По алгоритму (69)-(71) рассчитано изменение относительной длины тени lsh1 гномона в течение суток для трех разных дней на широте г. Москвы и для одного дня на широте г. Тюмени (Рис. 16). В левой части графика вертикальными линиями показаны часовые углы восходов Солнца, а в правой―часовые углы его заходов. Так как в эти моменты длина тени стремится к бесконечности, и вертикальный размер графика приходится ограничивать, то асимптотическое приближение линий lsh1 к вертикальным линиям восходов и заходов находится вне рисунка.

Как видно из Рис. 16, в день весеннего равноденствия (линия 1) длина тени lsh1 после восхода Солнца (ω < 0) уменьшается и в полдень (ω = 0) достигает минимального значения lsh1 = 1.433. Затем она увеличивается до бесконечности в момент захода Солнца (ω > 0). В день летнего солнцестояния (линия 2) минимальная длина тени уменьшается до lsh1 = 0.626, а в день зимнего солнцестояния минимальная длина тени увеличивается до lsh1 = 5.108. Таким образом, на широте г. Москвы длина полуденной тени изменяется в пределах 5.108 ≥ lsh1 ≥ 0.626.

Рассмотрим траекторию движения тени конца гномона в координатах

Рис. 15. Длина тени lsh гномона длиной lg при освещении его солнцем S.

Рис. 16. Изменение длины солнечной тени гномона единичной длины в течение суток в дни весеннего равноденствия 1, летнего 2 и зимнего 3 солнцестояний на широте Москвы (j˚ = 55.7522˚) в современную эпоху 30.12.1949 г. 4―длина солнечной тени гномона 15 августа 2015 г. вблизи Тюмени (j˚ = 57.301575˚) по расчету и 5―по наблюдениям. Вертикальными линиями отмечены часовые углы восходов (ω < 0) и заходов Солнца (ω > 0).

Nsh, Esh, где ось Nsh направлена на Север, а ось Esh ? на Восток. Азимут тени Ash отсчитывается от оси Nsh и выражается через азимут Солнца так: Ash = AS + 180˚. Тогда проекции тени на плоскость горизонта запишутся:

N s h = l s h 1 cos A s h = l s h 1 cos A S ; E s h = l s h 1 sin A S . (72)

На Рис. 17 показана траектория дневной тени конца гномона относительной длины для трех дней на широте г. Москвы и одного дня на широте г. Тюмени. В начале координат крестиком отмечено положение гномона. На широте г. Москвы в день весеннего равноденствия (линия 1) конец тени гномона движется с Запада (ω < 0) на Восток почти по прямой линии. При этом тень находится севернее гномона. В день летнего солнцестояния тень гномона (линия 2) начинает утром движение с Юго-запада и заканчивает движение вечером на Юго-востоке. Только вблизи полудня, ±2 часа, тень находится севернее гномона. Линией 3 показано, что в день зимнего солнцестояния тень полностью движется на Севере от гномона.

7.3. Измерение длины тени гномона

С целью проверки полученных результатов на местности, недалеко от г. Тюмени, были выполнены наблюдения за перемещением тени гномона. На Рис. 18 показан нехитрый инструментарий из подсобных материалов: лист ватмана, компас, карандаш, металлический прут, рулетка и отвес, с помощью которого выполнены наблюдения. Наблюдения проводились 15.08.2015 г., что составляет Td = 148 день после дня весеннего равноденствия

Рис. 17. Суточная траектория тени конца гномона единичной длины на плоскости горизонта Nsh(Esh) в дни весеннего равноденствия 1, летнего 2 и зимнего 3 солнцестояний на широте Москвы (j˚ = 55.7522˚) в современную эпоху 30.12.1949 г. 4―траектория тени гномона 15 августа 2015 г. вблизи Тюмени (j˚ = 57.301575˚) по расчету и 5―по наблюдениям. Ось Nsh направлена на Север, а ось Esh―на восток. Положение гномона в начале координат отмечено крестиком.

Рис. 18. Наблюдение тени гномона 15 августа 2015 г. вблизи Тюмени (j˚ = 57.301575˚) и основные инструменты: лист ватмана, компас, карандаш, металлический прут диаметром 6 мм―гномон, рулетка и отвес.

21 марта. Вдоль тени карандашом проводились линии на ватмане и отмечались конец тени и время ее наблюдения. Затем были замерены азимуты AshM теней от направления от направления на Север NM, определенный по компасу.

На Рис. 19 показаны результаты наблюдений в виде зависимости относительной длины тени lsh1 от времени в г. Тюмени tm. Как видно,

Рис. 19. Результаты наблюдения относительной длины тени lsh1 гномона в зависимости от времени tm в г. Тюмени 15 августа 2015 г. на широте j˚ = 57.301575˚.

минимум lsh1 приходится на время tm отличное от 12 часов. На Рис. 16 данные наблюдения точками 5 нанесены на рассчитанную согласно (69)-(71) зависимость 4lsh1(ω). Положение полдня (ω = 0) приходится на 12 ч. 44 мин тюменского времени.

Для рассматриваемой широты места наблюдения и дня Td = 148 по формуле (50) был рассчитан азимуты AS Солнца. Аналогично Рис. 10 он был представлены в виде графика AS(ω). В этих же координатах нанесены определенные азимуты Солнца ASM по замеренным азимутам тени ASM = AshM − 180˚. Азимуты ASM оказались эквидистантно сдвинуты на ΔA = 14.47059˚ вниз по отношению к азимутам AS, т.е. ASM < AS. Это обусловлено тем, что Северный магнитный полюс сдвинут к Востоку на величину ΔA по отношению к географическому Северному полюсу.

Скорректированные данные измерения азимута тени и времени ее наблюдения использованы для расчета по формулам (72) ее траектории. На Рис. 17 она нанесена точками 5. Как видим, измеренная траектория тени гномона 5 совпала с рассчитанной 4. Некоторый разброс точек 5 вокруг линии 4 объясняется погрешностью измерений.

Итак, выполненные наблюдения и измерения тени гномона подтвердили представленный алгоритм её расчета. Кроме того, измерения позволили определить истинный полдень в 12 ч. 44 мин тюменского времени и восточное отклонение ΔA магнитной стрелки. Следует отметить, что азимут солнечного полдня и его время наступления можно непосредственно определить по данным наблюдения, представленным на Рис. 19. Но для этого наблюдения необходимо производить более часто и за больший промежуток времени.

7.4. Полуденная длина тени гномона в течение года

Как видно из Рис. 16, наименьшая тень наступает в полдень (ω = 0). В этом случае зенитный угол zn центра Солнца определяется выражением (51), тогда в соответствии с (68) единичная длина полуденной тени запишется так

l s 1 n = tg ( φ δ 4.6599 × 10 3 / ρ ) . (73)

Алгоритм расчета относительной длины тени в п. 17.3 программы на SunPhnmen.mcd приведен с учетом особенностей функции tg. Результаты расчетов изменения длины полуденной тени в течение года на разных широтах Северного полушария показаны на Рис. 20. В этих расчетах ввиду несущественного влияния относительного расстояния ρ принято ρ = 1. Как видно из Рис. 20, на широте φ = 60˚ в день весеннего равноденствия (Td = 0) длина относительной полуденной тени ls1n = 1.687. Затем она уменьшается и в день летнего солнцестояния достигает минимального значения 0.734. Далее увеличивается и принимает максимальное значение 8.351 в день зимнего солнцестояния. С увеличением широты φ минимальные и максимальные значения возрастают, а с уменьшением φ ? уменьшаются. В тропических широтах (φ < ε) зависимость ls1n(φ) имеет два максимума и два минимума. В минимумах длина тени равна нулю: в эти моменты Солнце находится в зените. А максимумы ls1n наступают в дни солнцестояний.

7.5. Экстремальные длины солнечной тени гномона в дни солнцестояний

Как показано выше, экстремальные длины теней происходят в дни солнцестояний, которые характеризуются углом δ = ±ε. Тогда из (73) с учетом ρ = 1 получаем относительные длины теней

l s 1 s = tg ( φ ε 4.6599 × 10 3 ) ; l s 1 w = tg ( φ + ε 4.6599 × 10 3 ) , (74)

для дней летнего и зимнего солнцестояний, соответственно. С учетом особенностей функции tg алгоритм расчета приведен в п. 17.3 программы на SunPhnmen.mcd.

Результаты расчетов экстремальных длин тени для пяти разных эпох в зависимости от широты φ даны на Рис. 21. Шкала для длин теней в летнее солнцестояние ls1s приведена справа, а для зимнего солнцестояния ls1w

Рис. 20. Длина полуденной солнечной тени гномона единичной длины ls1n в зависимости от дней года Td в современную эпоху 30.12.1949 г. на разных широтах φ˚ Северного полушария: (a) в обычном масштабе; (b) −при увеличении ординаты ls1n в 2 раза.

Рис. 21. Тени гномона единичной длины в полдни летнего солнцестояния (ls1s) и зимнего солнцестояния (ls1w) на разных широтах φ˚ Северного полушария в экстремальные эпохи T за последние 50 т.л.н.: T―время в тыс. лет от 30.12.1949 г.; масштабы длин ls1w и ls1s―разные и отличаются в 4 раза.

слева. В современную эпоху, как видно из Рис. 21 при T = 0, в полдень летнего солнцестояния на экваторе (φ = 0) длина тени ls1s = 0.428. Затем она уменьшается на тропике φ ≈ ε до 0, а с дальнейшим увеличением широты длина тени возрастает до ls1s = 2.277 на широте φ = 90˚.

В другие эпохи графики полуденной тени в летнее солнцестояние почти эквидистантно сдвигаются по широте φ, в зависимости от угла наклона ε. На широтах φ = ε длина тени равна нулю: ls1s = 0.

В современную эпоху (T = 0) в полдень зимнего солнцестояния на экваторе (φ = 0) длина тени ls1w = 0.428, т.е. такая же, как и в полдень летнего солнцестояния. Затем она с увеличением широты φ растет до бесконечности до широты φ = 90˚ − ε˚ полярного круга, где начинается полярная ночь. В другие эпохи полуденная тень на экваторе изменяется от 0.259 (T = 46.44 т.л.н.) до 0.618 (T = 31 т.л.н.). С увеличением широты длины теней в день зимнего солнцестояния увеличиваются во все эпохи и тем в большей степени, чем меньше широта полярного круга (90˚ ? ε˚).

В табл. 4П работы [28] приведены длины полуденных теней в дни летнего (ls1s) и зимнего (ls1w) солнцестояний для пяти разных эпох. Изменение широты φ дается через 2˚, начиная с экватора до широты 80˚. Бесконечная длина тени выражена числом 1.6e+16.

8. Основные феномены Солнца

Подведем итоги по основным феноменам Солнца. Освещение Солнцем поверхности в точке нахождения наблюдателя M (см. Рис. 2), т.е. ее инсоляция, зависит от эксцентриситета e орбиты Земли, угла наклона ε плоскости экватора Земли к плоскости ее орбиты, угла φ между перигелием орбиты Земли и восходящим углом γ и широты места φ [19] . Эти четыре параметра создают разнообразие инсоляции поверхности Земли, из которого наблюдатель может определить следующие основные 6 явлений, или феноменов (см. Таблица 3). Первым феноменом является количество дней до начала сезона Td, отсчитываемое от момента весеннего равноденствия. Для весны Tdsp = 0, для лета наблюдатель определяет количество дней до летнего солнцестояния Tdsm, для осени―до момента осеннего равноденствия Tdau и для зимы―до момента зимнего солнцестояния Tdwn. С этим феноменом также связана продолжительность сезонов ΔTd, которая, согласно (6) рассчитывается по разности величин Td. В рассмотренном выше алгоритме начала сезонов определяются долготами λ = 0; π/2; π; 3π/2. По зависимости λ(Td), определяется время Td в днях до начала каждого сезона. Этот феномен не зависит от широты φ наблюдателя, а определяется только параметрами e, δ и φ.

Второй феномен, характерный для высоких широт, ―это длительность в днях полярных дней (ΔTdd) и полярных ночей (ΔTdn). Наблюдатель на широте φ своего места может определить длительность полярных ночи и дня. А в теории полярные ночи и дни определяются долготами λpd и λpn (см. формулы (35)-(36)), по которым длительность этих феноменов рассчитывается по зависимости λ(Td). Этот феномен сильно зависит от широты φ наблюдателя. Поэтому он в наибольшей степени находит для уточнения широты палеонаблюдателя, если он находился в высоких широтах.

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

z n s p = φ , z n s m = φ ε , z n a u = φ и z n w n = φ + ε . (75)

Таблица 3. Основные феномены Солнца.

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

l = π / 2 z n . (76)

Четвертый феномен―это долгота светового дня во время летнего солнцестояния Ds и зимнего солнцестояния Dw. Для наблюдателя существует сложность в измерении этого феномена. Эти измерения требуют высокой технологии определения времени и хранения сведений о нем. Тем не менее, возможны приемы выполнения этой работы. В теории этот феномен определяется выражением (43).

Пятый феномен―азимуты восходов и заходов Солнца в дни солнцестояний и равноденствий. Они могут быть зафиксированы наблюдателем в зависимости от какого-то направления. В теории азимут восхода Солнца в день летнего солнцестояния определяется выражениями (66), а зимнего―(67). В дни равноденствий при δ =0, согласно (60) и (64) азимут восхода Солнца будет:

A G r s p = A G r a u = 90 [ 1 ( 9.8902 × 10 3 + 4.6599 × 10 3 ) tg φ / 0.5 π ] . (77)

Шестой феномен―длина теней гномона в полдни. Для дней солнцестояний: летнего (ls1s) и зимнего (ls1w) определяется согласно (74), а в дни равноденствий при δ = 0, согласно (73), длина тени будет:

l s 1 s p = l s 1 a u = tg ( φ 4.6599 × 10 3 ) . (78)

Наблюдатель может зафиксировать эти феномены по отношению к гномону, длина которого неизвестна. В этом случае могут потребоваться дополнительные данные для определения относительной длины тени. В теории относительные длины теней определяются выражениями (74) и (78).

9. Восстановление параметров наблюдателя и движения Земли по древним феноменам Солнца

Параметрами наблюдателя являются его древние широта местонахождения и направление на Север. Исследователи древних календарей отмечают изменение этих параметров наблюдателя. Например, М.И. Исрапилов пришел к выводу о смещении Северного полюса в Якутию несколько тысяч лет назад [14] . Рассмотрим как определить параметры наблюдателя по феноменам Солнца.

В дни равноденствий зенитный угол z, согласно (75), в полдень будет равняться широте места znsp = φ. А длина тени ls1sp согласно (77) также зависит от широты φ. Тогда широта места наблюдателя определяется по одному из этих феноменов:

φ = z n s p либо φ = arctg ( l s 1 s p + 4.6599 × 10 3 ) . (79)

В дни равноденствий суточное движение Солнца происходит по кругу экватора AA’ (Рис. 2). Поэтому центр Солнца восходит на Востоке, а заходит на Западе. Палеонаблюдатель отмечает восходы и заходы Солнца по его верхнему краю. По формуле (77) при известной широте φ может быть найден азимут AGrsp восхода верхнего края Солнца. Обозначим наблюдаемое значение азимута палео-Востока как AEapa = 90˚, тогда из (77) можем записать:

A E a p a = A G 1 s p + 90 ( 9.8902 × 10 3 + 4.6599 × 10 3 ) tg φ / 0.5 π . (80)

Тогда азимут палео-Севера запишется так:

A N p a = A E a p a 90 = A G 1 s p a 90 [ 1 ( 9.8902 × 10 3 + 4.6599 × 10 3 ) tg φ / 0.5 π ] (81)

где AG1spa ? наблюдаемый азимут древнего восхода Солнца во время весеннего равноденствия.

Азимут палео-Севера также можно определить через азимут восхода Солнца в момент осеннего равноденствия, а также через азимуты заходов в дни равноденствий. Четырехкратные способы определения предоставляют возможности для проверки результатов определения ANpa.

Теперь на основе феноменов Солнца определим параметры движения Земли. В полдень солнцестояний (λ = π/2, 1.5π) зенитные углы z Солнца, согласно (75), и длины тенней в полдень, согласно (74), зависят от φ ? ε и φ + ε. Поэтому при известной широте φ угол наклон ε плоскости экватора к плоскости орбиты Земли может быть определен по одной из четырех зависимостей:

ε = φ z n s m ; ε = ( φ z n s m ) ; ε = φ arctg ( l s 1 s m ) 4 . 66 × 10 3 ;

ε = [ φ arctg ( l s 1 w n ) ] 4 . 66 × 1 0 3 . (82)

Как видим, зависимости (82) позволяют разными четырьмя способами определить угол наклона ε при известной широте φ. Если широта φ неизвестна, то соотношения (82) позволяют определить и широту φ и угол ε, при этом двумя способами.

Предельные азимуты восходов в дни солнцестояний определяются выражениями (66) и (67), которые зависят от ε и φ. Вычтем из (67) выражение (66) и выразим угол наклона ε:

ε = arcsin [ sin ( π Δ A s r / 36 0 ) cos φ ] , (83)

где ΔAsr = AGrwn ? AGrsm ? разница азимутов восходов в дни солнцестояний.

Аналогичное выражение для угла наклона запишется через разницу ΔAss азимутов заходов в дни солнцестояний. Поэтому, если величины ΔAss и ΔAsr имеются, то этим предоставляется возможность проверки определения ε. Выражение (83) знаменательно тем, что в него входит разность азимутов восходов. Это исключает возможное, отличное от современного, толкование восхода и определение его момента. В этом случае на определение ε также меньше влияют географические особенности места наблюдения.

Долготы светового дня в моменты солнцестояний определяются выражениями (43). С учетом (37)-(40) найдем разность светового дня ΔD = Ds − Dw

Δ D = 24 π ( arccos ( cos z 0 a 1 sin ε sin φ cos ε cos φ ) arccos ( cos z 0 a 1 sin ε sin φ cos ε cos φ ) ) . (84)

Упростим выражение (84), рассматривая его для центра Солнца, т.е. согласно (23) z0a1 = π/2. Тогда оно приобретет вид:

Δ D = 24 π [ arccos ( t g φ t g ε ) arccos ( t g φ t g ε ) ] , (85)

После упрощения (85) выразим угол наклона ε:

ε = arctg { [ sin ( 0.5 π Δ D u d ) ] / tg φ } , (86)

где

Δ D u d = Δ D / 24 (87)

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

Выражение (86) позволяет определить угол наклона ε при известной широте места φ. Оно не требует знания способа деления суток на части, так как в (86) входит доля от суток, на которую солнечный день в летнее солнцестояние больше дня в зимнее солнцестояние. После нахождения ε из (86), по выражению (84) можно уточнить ΔD, а затем по (86) определить уточненное значение угла наклона ε.

Определенные по представленным формулам широта φ и угол наклона ε в дни равноденствий и солнцестояний содержат ошибки, связанные с определением палеонаблюдателем азимута, длины тени, продолжительности разности дня ΔD. Кроме того в формулах содержится ошибка, связанная с тем, что принятые моменты равноденствий и солнцестояний на заходе и восходе Солнца или в полдень могут не совпадать с моментом прохождения Солнца через соответствующие точки γ, γ', E и E' на эклиптике EE' (Рис. 2). Разность моментов может достигать 0.5 дня. Он дает погрешность пропорциональной отношению 0.5 дня к периоду прецессии оси Земли, равному 25,780 тыс. лет. То есть, эта погрешность несущественна. Погрешность наблюдений палеонаблюдателя также может быть оценена. Для этого необходимо современному наблюдателю повторить наблюдения в месте нахождения палеонаблюдателя. Обработка их результатов позволит определить погрешность палеонаблюдателя.

При известных палеопараметрах φ и ε остальные параметры: эксцентриситет e орбиты Земли и долгота φ перигелия могут быть определены по количеству дней Tdsm, Tdau, Tdwn до начала сезона и по длительности полярных дней ΔTdd или ночей ΔTdn. Эти параметры зависят от долготы λ, которая рассчитывается численно в прилагаемой программе на SunPhnmen.mcd. Методом последовательных приближений с ее помощью параметры e и φ могут быть найдены. Они могут быть определены разными способами, так как для нахождения двух неизвестных e и φ существует 5 уравнений для величин Tdsm, Tdau, Tdwn, ΔTdd и ΔTdn.

Изменение параметров ε, e и φ во времени известно в результате решения задач об орбитальном движении Земли и о ее вращательном движении. В файле OrAl-5kyr.prn величины ε, e и φ даны с интервалом 1 год за 5 т.л.н., а в файле OrAl-200ky.prn ? с интервалом 40 лет за 200 т.л.н. Поэтому, по трем параметрам ε, e и φ должен однозначно определиться возраст эпохи палеонаблюдателя. Если из археологических исследований примерный возраст эпохи наблюдения известен, то по одной из зависимостей ε, e или φ от времени этот возраст может быть уточнен.

Рассмотрим схематически алгоритмы определения параметров e и φ в двух случаях: 1) при известном примерном возрасте эпохи наблюдения и 2) при неизвестном. Пусть и в первом и во втором случае определены два каких-либо параметра палеоклимата, например, Tdsm,p и ΔTdd,p. Индекс p, здесь и далее добавленный к обозначению параметра, обозначает, что этот параметр определен по древнему (палео-) календарю.

В первом случае пусть T1 и T2―время начальной и конечной эпох, в пределах которого находится предполагаемый возраст древнего календаря. С помощью программы на SunPhnmen.mcd на интервале от T1 до T2 рассчитываются Tdsm и ΔTdd и строятся графики Tdsm) и ΔTdd(e). По этим графикам по величинам Tdsm,p и ΔTdd,p определяются угол перигелия φ,p и эксцентриситет ep, которые были в эпоху древнего календаря. Время эпохи Tp определяется из файла данных OrAl-5kyr.prn или OrAl-200ky.prn.

Во втором случае, когда возраст календаря неизвестен, при заданных параметрах φ и ε с помощью программы на SunPhnmen.mcd производится первая серия расчетов при неизменном угле перигелия φ варьируется эксцентриситет e. Во второй серии при неизменном эксцентриситете e варьируется угол перигелия φ. Затем строится первая серия графиков Tdsm) при разных e и вторая серия графиков ΔTdd(e) при разных φ. Затем на графиках Tdsm) проводится горизонтальная линия Tdsm,p = const, а на графиках ΔTdd(e) ? ΔTdd,p = const. По пересечению этих линий с линиями графиков Tdsm) и ΔTdd(e) определяются параметры: угол φ,p и эксцентриситет ep, которые были в эпоху древнего календаря. Время эпохи Tp определяется из файла данных OrAl-5kyr.prn или OrAl-200ky.prn по трем параметрам: ε, φ,p и ep.

10. Использование результатов в аэрокосмических исследованиях рельефа

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

Осенью 2013 г. оленеводами на Ямале, в 30 км Бованенковского месторождения (70˚21'44'' с.ш. и 68˚26'46'' в.д.) [29] обнаружен провал в земле, который еще называют воронкой. Он имел диаметр 25 м и глубину больше 50 м [30] . Вверху провал завершается коническим раструбом диаметром около 60 м. По космическому снимку от 09.06.2013 г. (см. Рис. 22(a)) было установлено [30] , что на месте провала находился бугор 1. Цифрой 2 показана тень от бугра.

Солнце находится в направлении противоположном тени, и, как следует из Рис. 22(a), азимут Солнца равен AShl = 326.5˚. Широта бугра φ = 70.3622˚ , а дата его наблюдения соответствует Td = 81 день от дня равноденствия 21 марта. По зависимости азимута Солнца AS от часового угла ω (50) время Td и широта бугра φ позволяют рассчитать часовой угол Солнца ω , который обозначим как ωhS. Эти расчеты выполнены в в п. 18 программы на SunPhnmen.mcd.

Ряд азимутов Солнца AS,k по формуле (50) рассчитан для ряда часовых углов ωh,k . По значению азимута AShl из ряда азимутов AS,k выбирается индекс in12 , при котором азимут AS,in12 находится вблизи значения AShl , но меньше его. Тогда часовой угол Солнца находится интерполяцией по следующей формуле.

Рис. 22. Бугор (a) от 9 июня 2013 г. вблизи Бованенково на Ямале [30] перед появлением осенью 2013 г. воронки и профиль бугра в метрах (b), рассчитанный по его тени: точки―по результатам измерения тени, линия―апроксимационная зависимость (90).

ω h S = ω h , i n 12 + ( ω h , i n 12 + 1 ω h , i n 12 ) ( A S , i n 12 + 1 A S h l ) A S , i n 12 + 1 A S , i n 12 . (88)

По часовому углу Солнца ωhS с помощью формулы (71) рассчитывается единичная длина гномона lsh1hl в месте нахождения бугра. Так как формула (71) и входящие в нее выражения являются сложными нелинейными зависимостями, то вычисления в п. 18 программы на SunPhnmen.mcd производятся аналогично определению часового угла по формуле (88).

Наблюдаемую длину тени бугра lOsh можно определить по космоснимку (Рис. 22(a)) с учетом приведенного на нем масштаба длины, равному 70 м (см. Рис. 22). Тогда высота бугра определяется так:

H = l O s h l s h 1 h l . (89)

Такие расчеты выполнены для ряда точек тени и получен профиль бугра. На Рис. 22(b) точками 1 представлен профиль бугра, определенный по замерам длины тени. Линией 2 приведена аппроксимация профиля бугра зависимостью

H c s = H m x cos [ 0.5 π ( r R h ) 2 ] , (90)

где Hmx = 8.76 м ? высота бугра в центральной точке;

Rh = 21.2 м ? радиус бугра.

В работе [30] по стереопаре снимков оценили диаметр бугра в основании: 45 - 48 м и его высоту: 5 - 6 м. Как видим, полученные по тени размеры бугра согласуются с этими оценками: диаметр бугра по тени меньше на 9.7%, а высота − больше на 37%. В работе [30] сообщается, что в месте нахождения воронки наблюдается уклон местности. Возможно, наличие уклона увеличивает высоту холма.

Рассчитаем тень от гномона, когда она падает не на горизонтальную поверхность, а имеющую уклон β к горизонтальной поверхности HH' (см. Рис. 15). Обозначим длину тени гномона на наклонной поверхности lshβ = OB . В треугольнике AOB угол B = 90 ( z S ) β . Тогда по теореме синусов l g / sin B = O B / sin ( z S ) находим длину тени гномона на наклонной поверхности OB

l s h β = l g sin ( z S ) cos ( z S + β ) . (91)

После преобразования (91) длину тени гномона единичной длины на поверхности с углом уклона β получаем в виде:

l s h 1 β = l s h 1 cos β l s h 1 sin β , (92)

где lsh1―длина тени гномона единичной длины на горизонтальной поверхности;

β―угол наклона поверхности в направлении тени.

При наличии уклона в месте нахождения бугра длина тени согласно (92) увеличится. Это приведет, в соответствии с (89), к уменьшению высоты H бугра.

Представляет интерес полученный профиль (90) бугра. Если бугор образован движением льдистого материала, то его профиль даже при движении с неизменным ускорением будет такой же, как и профиль скорости. При ламинарном течении вязкой жидкости в трубе профилем скорости является парабола:

v = v max [ 1 ( r R h ) 2 ] . (93)

Тогда профиль бугра будет иметь вид:

H q = H m x [ 1 ( r R h ) 2 ] . (94)

Профиль скорости, аналогичный профилю бугра (90), получен для вертикальной скорости в сходящемся к центру течении жидкости, например, в пористой трубе с одним заглушенным торцом (см. стр. 59 [31] ). Профили (90) и (94) для лучшего сравнения на Рис. 23 нанесены на профиль тени. При этом в выражениях (90) и (94) вместо Hmx использована максимальная длина тени lOshmx = 70.4 м. Как видно из графиков, профиль (90) для сходящегося течения значительно лучше согласуется с профилем тени, нежели профиль ламинарного течения (94). Имеющееся в правой

Рис. 23. Сравнение профиля тени с апроксимационными зависимостями: 1―для сходящегося течения (90); 2―для ламинарного течения (94); точки―по результатам измерения тени; lsh и r―длина тени бугра и его радиус в м.

части отличие наблюдаемого профиля тени от зависимости (90) вызвано влиянием снега на изображение правой стороны бугра (Рис. 22(a)).

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

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

11. Заключение

Большие колебания оси вращения Земли объясняют колебания палеоклимата и имеют косвенное подтверждение в некоторых древних календарях и в исторических свидетельствах. Расшифровка древних календарей затруднена сложностью процессов, от которых зависят солнечные явления, фиксируемые календарями. Разработана теория феноменов Солнца, которая алгоритмизирована и реализована в виде программы в среде MathCad. В отличие от традиционной теории, основанной на уравнении Кеплера, в новой теории долгота λ годового движения Солнца в зависимости от номера дня года определяется по результатам точного решения задачи двух тел.

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

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

Cite this paper

Smulsky, J.J. (2018) The Sun’s Movement in the Sky Now and in the Past. Open Access Library Journal, 5: e4250. https://doi.org/10.4236/oalib.1104250

References

  1. 1. Дирингская культура (2016) Википедия.
    https://ru.wikipedia.org/wiki/

  2. 2. Оппенгеймер, С. (2004) Изгнание из Эдема. М.: Из-во Эксмо. - 640 с.

  3. 3. Смульский, И.И. (2014) Основные положения и новые результаты астрономической теории изменения климата. Институт криосферы Земли СО РАН. Тюмень,. - 30 с.: ил: 16. - Библиогр.: 44 назв. - Рус. Деп. в ВИНИТИ РАН 30.09.2014, No. 258-В2014.
    http://www.ikz.ru/~smulski/Papers/OsPoATLP3.pdf

  4. 4. Smulsky, J.J. (2016) New Results on the Earth Insolation and Their Correlation with the Late Pleistocene Paleoclimate of West Siberia. Russian Geology and Geophysics, 57, 1099-1110.
    https://doi.org/10.1016/j.rgg.2016.06.009

  5. 5. Смульский, И.И. (2016) Эволюция оси Земли и палеоклимата за 200 тысяч лет. Saarbrucken, Germany: “LAP Lambert Academic Publishing”, 228 с.
    http://www.ikz.ru/~smulski/Papers/InfEvEAPC02MEn.pdf

  6. 6. Svendsen, J.I., Astakhov, V.I., Bolshiyanov, D.Y., Demidov, I., Dowdeswell, J.A., Gataflin, V., Hjort, C., Hubberten, H.W., Larsen, E., Mangerud, J., Melles, M., Moller, P., Saarnisto, M. and Siegert, M.J. (1999) Maximum Extent of the Eurasian Ice Sheets in the Barents and Kara Sea Region during the Weichselian. Boreas, 28, 234-242.
    https://doi.org/10.1111/j.1502-3885.1999.tb00217.x

  7. 7. Флинт, Р.Ф. (1978) История Земли. М.:"Прогресс", 358 с.

  8. 8. Смульский, И.И. (2015) Новые инсоляционные периоды и последние похолодания в плиоцене/В сб. Арктика, Субарктика: мозаичность, контрастность, вариативность криосферы: Труды международной конференции/Под ред. В.П.Мельникова и Д.С.Дроздова. Тюмень: Изд-во Эпоха, С. 360-363.
    http://www.ikz.ru/~smulski/Papers/smulsky_J_J2015_03_15c1.pdf

  9. 9. Смульский И.И. and Сеченов К.Е. (2007) Уравнения вращательного движения Земли и их решения при воздействии Солнца и планет / Институт криосферы Земли СО РАН. Тюмень, 35 с. - ил.: 7. Библиогр.: 19 назв. - Рус. - Деп. в ВИНИТИ 02.05.07 г. No. 492-В2007.
    http://www.ikz.ru/~smulski/Papers/UVrVzSPc.pdf

  10. 10. Smulsky J.J. (2011) The Influence of the Planets, Sun and Moon on the Evolution of the Earth’s Axis. International Journal of Astronomy and Astrophysics, 1, 117-134.
    http://www.scirp.org/journal/ijaa/ https://doi.org/10.4236/ijaa.2011.13017

  11. 11. Smul’skii, I.I. (2013) Analyzing the Lessons of the Development of the Orbital Theory of the Paleoclimate. Herald of the Russian Academy of Sciences, 83, 46-54.
    https://link.springer.com/article/10.1134%2FS1019331613010073
    https://doi.org/10.1134/S1019331613010073

  12. 12. Laskar, J., Robutel, P., Joutel, F., Gastineau, M., Correia, A.C.M. and Levrard, B. (2004) A Long-Term Numerical Solution for the Earth. Icarus, 170, 343-364.
    https://doi.org/10.1016/j.icarus.2004.04.005

  13. 13. Smulsky, J.J. (2016) Fundamental Principles and Results of a New Astronomic Theory of Climate Change. Advances in Astrophysics, 1, 1-21.
    http://www.isaacpub.org
    http://www.isaacpub.org/Journal/AdAp https://doi.org/10.22606/adap.2016.11001

  14. 14. Исрапилов, М.И. (2003) Наскальные рисунки Дагестана и колебания полюсов и наклона оси Земли в голоцене. Махачкала: Издательство “Юпитер”. 432 с.

  15. 15. Тилак, Б.Г. (2002) Арктическая родина в Ведах. М.: ФАИР-ПРЕСС., 528 с.

  16. 16. Chippindale, Ch. (1983) Stonehenge Complete: Thames & Hudson and Ithaca, London.

  17. 17. Morner, N.-A. and Lind, B.G. (2012) Stonehenge Has Got a Younger Sister Ales Stones in Sweden Decoded. International Journal of Astronomy and Astrophysics, 2, 23-27
    https://doi.org/10.4236/ijaa.2012.21004

  18. 18. Ларичев, В.Е. (1989) Мудрость змеи: Первобытный человек, Луна и Солнце. Новосибирск: “Наука”, Сибирское Отделение, 272 с.

  19. 19. Смульский, И.И. and Кротов, О.И. (2013) Новый алгоритм расчета инсоляции Земли/Институт криосферы Земли СО РАН. - Тюмень, 38 с. - Деп. в ВИНИТИ 08.04.2013 No. 103-В2013.
    http://www.ikz.ru/~smulski/Papers/NwAlClI2c.pdf

  20. 20. Smulsky, J.J. and Krotov, O.I. (2014) New Computing Algorithm of the Earths’ Insolation. Applied Physics Research, 6, 56-82.
    https://doi.org/10.5539/apr.v6n4p56

  21. 21. Смульский, И.И. (1999) Теория взаимодействия. - Новосибирск: Из-во Новосиб. ун-та, НИЦ ОИГГМ СО РАН, 294 с.
    http://www.ikz.ru/~smulski/TVfulA5_2.pdf

  22. 22. Смульский, И.И. (2007) Математическая модель Солнечной системы/В сб. Теоретические и прикладные задачи нелинейного анализа. Российская Академия Наук: ВЦ им. А.А. Дородницына. М.: ВЦ РАН А.А. Дородницына. С. 119-138.
    http://www.ikz.ru/~smulski/Papers/MatMdSS5.pdf

  23. 23. Справочное руководство по небесной механике и астродинамике (1976) /Под ред. Дубошина Г.Н. Изд. 2-е, доп. и перераб. М., Наука, 862 с.

  24. 24. Труды ИПА РАН. Вып. 10. Эфемеридная астрономия (2004) Санкт-Петербург: ИПА РАН, 488 с.

  25. 25. Птолемей, К. (1998) Альмагест. Математическое сочинение в тринадцати книгах. М.: Наука. Физматлит, 672 с.

  26. 26. Ньютон, Р.Р. (1985) Преступление Клавдия Птолемея: пер. с англ. М.: Наука, 384 с.

  27. 27. Smulsky, J.J. (2014) Exact Equations for the Light Doppler Effect. Journal of Modern Physics, 5, 1602-1607.
    https://doi.org/10.4236/jmp.2014.516161

  28. 28. Смульский, И.И. (2016) Феномены Солнца в исторической перспективе/Институт криосферы Земли СО РАН. - Тюмень, 66 с. - Илл.: 23.- Библиогр.: 24 назв. - Рус. Деп. в ВИНИТИ РАН 11.01.2016, No. 9-В2016.
    http://www.ikz.ru/~smulski/Papers/SunPhnmen.pdf

  29. 29. Богоявленский, В.И. and Гарагаш, И.А. (2015) Обоснование процесса образования кратеров газового выброса в Арктике математическим моделированием. Арктика: экология и экономика, No. 3 (19), с. 12-17.

  30. 30. Кизяков, А.И., Сонюшкин, А.В., Лейбман, М.О., Зимин, М.В. and Хомутов, А.В. (2013) Геоморфологические условия образования воронки газового выброса и динамика этой формы на центральном Ямале. Криосфера Земли. Т. XVII. No. 4. С. 36-47.

  31. 31. Смульский, И.И. (1992) Аэродинамика и процессы в вихревых камерах. Новосибирск: ВО "Наука", 301 с.
    http://www.ikz.ru/~smulski/Aerpro/aerpro.djvu