Положение Солнца
Положение Солнца на небе является функцией как времени , так и географического места наблюдения на . поверхности Земли Когда Земля вращается вокруг Солнца называемой в течение года , кажется, что Солнце движется относительно неподвижных звезд на небесной сфере по круговой траектории, эклиптикой .
Вращение Земли вокруг своей оси вызывает суточное движение , так что кажется, что Солнце движется по небу по траектории Солнца наблюдателя , которая зависит от географической широты . Время, когда Солнце проходит наблюдателя, меридиан зависит от географической долготы .
Поэтому, чтобы найти положение Солнца в данном месте в данный момент времени, можно выполнить следующие три шага: [1] [2]
- вычислить положение Солнца в эклиптической системе координат ,
- перевести в экваториальную систему координат и
- преобразовать в горизонтальную систему координат для местного времени и местоположения наблюдателя. Это система координат, обычно используемая для расчета положения Солнца с точки зрения угла зенита Солнца и угла азимута Солнца , и эти два параметра могут использоваться для изображения траектории Солнца . [3]
Этот расчет полезен в астрономии , навигации , геодезии , метеорологии , климатологии , солнечной энергетике и проектировании солнечных часов .
Примерное положение
[ редактировать ]Эклиптические координаты
[ редактировать ]Эти уравнения из Астрономического Альманаха : [4] [5] может использоваться для расчета видимых координат Солнца , среднего равноденствия и эклиптики даты с точностью около 0°,01 (36 дюймов) для дат между 1950 и 2050 годами. Подобные уравнения закодированы в программе Fortran 90 на языке Ссылка. [3] и используются для расчета угла зенита Солнца и угла азимута Солнца , наблюдаемого с поверхности Земли.
Начните с расчета n , количества дней (положительных или отрицательных, включая дробные дни) с полудня по Гринвичу по земному времени 1 января 2000 года ( J2000.0 ). Если известна юлианская дата искомого времени, то
Средняя долгота Солнца с поправкой на аберрацию света равна:
Средняя аномалия Солнца (на самом деле Земли на ее орбите вокруг Солнца, но удобно представить, что Солнце вращается вокруг Земли) равна:
Помещать и в диапазоне от 0° до 360°, прибавляя или вычитая кратные 360° по мере необходимости и действительно подлежат оценке ( мод. 360).
Наконец, эклиптическая долгота Солнца равна:
Эклиптическая широта Солнца составляет около:
- ,
поскольку эклиптическая широта Солнца никогда не превышает 0,00033° (чуть больше 1 дюйма), [6] а расстояние Солнца от Земли в астрономических единицах равно:
- .
Наклон эклиптики
[ редактировать ]Если наклон эклиптики не получен где-либо еще, его можно аппроксимировать:
Экваториальные координаты
[ редактировать ], и образуют полное положение Солнца в эклиптической системе координат . Ее можно преобразовать в экваториальную систему координат , вычислив наклон эклиптики : и продолжая:
- , где находится в том же квадранте, что и ,
Чтобы получить RA в правом квадранте в компьютерных программах, используйте функцию Arctan с двойным аргументом, например ATAN2(y,x).
и склонение ,
- .
Прямоугольные экваториальные координаты
[ редактировать ]Правые прямоугольные экваториальные координаты в астрономических единицах :
- Где ось находится в направлении мартовского равноденствия , ось к июньскому солнцестоянию , а оси к Северному полюсу мира . [7]
Горизонтальные координаты
[ редактировать ]Склонение Солнца, вид с Земли
[ редактировать ]Кажется, что Солнце движется на север во время северной весны , пересекая небесный экватор в мартовское равноденствие . Его склонение достигает максимума, равного углу наклона оси Земли (23,44° или 23°26'). [8] [9] в день июньского солнцестояния , затем уменьшается до достижения минимума (-23,44° или -23°26') в день декабрьского солнцестояния , когда его значение является отрицательным по отношению к наклону оси. Эта вариация создает времена года .
Линейный график склонения Солнца в течение года напоминает синусоидальную волну с амплитудой 23,44°, но, помимо прочих отличий, один лепесток волны на несколько дней длиннее другого.
Следующие явления произошли бы, если бы Земля была идеальной сферой , находящейся на круговой орбите вокруг Солнца, и если бы ее ось была наклонена на 90 °, так что сама ось находилась в плоскости орбиты (аналогично Урану ). В какой-то день в году Солнце будет находиться прямо над Северным полюсом , поэтому его склонение будет +90°. В течение следующих нескольких месяцев подсолнечная точка будет двигаться к Южному полюсу с постоянной скоростью, пересекая круги широты с постоянной скоростью, так что склонение Солнца будет линейно уменьшаться со временем. В конце концов, Солнце окажется прямо над Южным полюсом со склонением -90 °; тогда он начнет двигаться на север с постоянной скоростью. Таким образом, график солнечного склонения, если смотреть с этой сильно наклоненной Земли, будет напоминать треугольную волну, а не синусоидальную волну, зигзагообразную между плюс и минус 90°, с линейными сегментами между максимумами и минимумами.
Если уменьшить осевой наклон на 90°, то абсолютное максимальное и минимальное значения наклона уменьшятся, чтобы сравняться с осевым наклоном. Кроме того, формы максимумов и минимумов на графике станут менее острыми и будут изогнуты, напоминая максимумы и минимумы синусоидальной волны. Однако даже когда осевой наклон равен наклону реальной Земли, максимумы и минимумы остаются более острыми, чем у синусоидальной волны.
На самом деле Земли эллиптическая . орбита [примечание 1] Земля движется вокруг Солнца быстрее вблизи перигелия в начале января, чем вблизи афелия в начале июля. Из-за этого такие процессы, как изменение склонения Солнца, происходят в январе быстрее, чем в июле. На графике это делает минимумы более острыми, чем максимумы. Кроме того, поскольку перигелий и афелий не совпадают с датами солнцестояний, максимумы и минимумы немного асимметричны. Темпы изменений до и после не совсем равны.
Таким образом, график кажущегося солнечного склонения во многом отличается от синусоидальной волны. Его точный расчет требует некоторой сложности, как показано ниже.
Расчеты
[ редактировать ]Склонение Солнца , δ ☉ , представляет собой угол между лучами Солнца и плоскостью экватора Земли. Земли Наклон оси (называемый астрономами наклоном эклиптики ) — это угол между осью Земли и линией, перпендикулярной орбите Земли. Наклон оси Земли меняется медленно в течение тысяч лет, но его нынешнее значение около ε = 23 ° 44 'почти постоянно, поэтому изменение склонения Солнца в течение одного года почти такое же, как и в течение следующего года.
В дни солнцестояний угол между лучами Солнца и плоскостью земного экватора достигает максимального значения 23°44'. Следовательно, δ ☉ = +23°44’ в день северного летнего солнцестояния и δ ☉ = –23°44’ во время южного летнего солнцестояния.
В момент каждого равноденствия кажется, что центр Солнца проходит через небесный экватор , а δ ☉ составляет 0°.
Склонение Солнца в любой момент рассчитывается по формуле:
где EL — долгота по эклиптике (по сути, положение Земли на ее орбите). Земли Поскольку эксцентриситет орбиты невелик, ее орбиту можно аппроксимировать кругом, что приводит к ошибке до 1 °. Приближение круга означает, что EL будет на 90° впереди солнцестояний на орбите Земли (в точках равноденствия), так что sin(EL) можно записать как sin(90+NDS)=cos(NDS), где NDS — число дней после декабрьского солнцестояния. Используя также приближение, согласно которому arcsin[sin(d)·cos(NDS)] близко к d·cos(NDS), получается следующая часто используемая формула:
где N — день года, начинающийся с N=0 в полночь по всемирному времени (UT) в начале 1 января (т. е. дни, составляющие порядковую дату -1). Число 10 в (N+10) представляет собой приблизительное количество дней после декабрьского солнцестояния до 1 января. Это уравнение завышает склонение вблизи сентябрьского равноденствия на величину до +1,5°. Аппроксимация синусоидальной функции сама по себе приводит к ошибке до 0,26°, и ее не рекомендуется использовать в приложениях солнечной энергетики. [2] Формула Спенсера 1971 года [10] (на основе ряда Фурье ) также не рекомендуется из-за погрешности до 0,28°. [11] Дополнительная ошибка до 0,5 ° может возникнуть во всех уравнениях, касающихся равноденствий, если не использовать десятичный знак при выборе N для корректировки времени после полуночи UT для начала этого дня. Таким образом, приведенное выше уравнение может иметь погрешность до 2,0°, что примерно в четыре раза больше угловой ширины Солнца, в зависимости от того, как оно используется.
Склонение можно более точно рассчитать, не делая двух приближений, используя параметры орбиты Земли для более точной оценки EL: [12]
который можно упростить, вычислив константы до:
N — это количество дней, прошедших с полуночи по всемирному стандартному времени с начала 1 января (т. е. количество дней в порядковой дате -1), и оно может включать десятичные дроби для корректировки местного времени позже или раньше в тот же день. Земли Число 2 в (N-2) — это приблизительное количество дней после 1 января до перигелия . Число 0,0167 — это текущее значение эксцентриситета орбиты Земли. Эксцентриситет меняется со временем очень медленно, но для дат, довольно близких к настоящему, его можно считать постоянным. Самые большие ошибки в этом уравнении составляют менее ± 0,2 °, но меньше ± 0,03 ° для данного года, если число 10 корректируется в большую или меньшую сторону с дробными днями, что определяется тем, насколько далеко до или после произошло декабрьское солнцестояние в предыдущем году. полдень 22 декабря. Эти точности сравниваются с расширенными расчетами NOAA. [13] [14] которые основаны на алгоритме Жана Меуса 1999 года с точностью до 0,01°. [15]
(Приведенная выше формула связана с достаточно простым и точным расчетом Уравнения Времени , которое описано здесь .)
Более сложные алгоритмы [16] [17] Исправьте изменения эклиптической долготы, используя термины в дополнение к поправке на эксцентриситет 1-го порядка, указанной выше. Они также корректируют наклон 23,44°, который со временем меняется очень незначительно. Поправки могут также включать влияние Луны на смещение положения Земли от центра орбиты пары вокруг Солнца. После получения склонения относительно центра Земли применяется дополнительная поправка на параллакс , которая зависит от расстояния наблюдателя от центра Земли. Эта поправка составляет менее 0,0025°. Ошибка расчета положения центра Солнца может составлять менее 0,00015°. Для сравнения, ширина Солнца составляет около 0,5°.
Атмосферная рефракция
[ редактировать ]Описанные выше расчеты склонения не учитывают эффекты преломления света в атмосфере, из-за чего видимый угол подъема Солнца, видимый наблюдателем, оказывается выше фактического угла подъема, особенно при низких высотах Солнца. [2] Например, когда Солнце находится на высоте 10°, кажется, что оно находится на высоте 10,1°. Склонение Солнца можно использовать вместе с его прямым восхождением для расчета его азимута, а также истинной высоты, которую затем можно скорректировать с учетом рефракции, чтобы получить его видимое положение. [2] [14] [18]
Уравнение времени
[ редактировать ]Помимо годовых колебаний видимого положения Солнца с севера на юг, соответствующих описанному выше изменению его склонения, существует также меньшее, но более сложное колебание в направлении восток-запад. Это вызвано наклоном оси Земли, а также изменением скорости ее орбитального движения вокруг Солнца, вызванным эллиптической формой орбиты. [2] Основными эффектами этих колебаний с востока на запад являются изменения во времени таких событий, как восход и закат, а также в показаниях солнечных часов по сравнению с часами, показывающими местное среднее время . Как показывает график, солнечные часы могут идти быстрее или медленнее примерно на 16 минут по сравнению с часами. Поскольку Земля вращается относительно Солнца со средней скоростью один градус каждые четыре минуты, это 16-минутное смещение соответствует сдвигу видимого положения Солнца на восток или запад примерно на четыре градуса по сравнению с его средним положением. Сдвиг на запад приводит к тому, что солнечные часы опережают часы.
Поскольку основной эффект этого колебания касается времени, его называют уравнением времени , используя слово «уравнение» в несколько архаичном смысле, означающем «коррекцию». Колебания измеряются в единицах времени, минутах и секундах, что соответствует тому, на сколько солнечные часы опережают часы. Уравнение времени может быть положительным или отрицательным.
Аналемма
[ редактировать ]Аналемма — это диаграмма , показывающая годовое изменение положения Солнца на небесной сфере относительно его среднего положения, если смотреть из фиксированного места на Земле. (Слово аналемма также иногда, но редко, используется в других контекстах.) Его можно рассматривать как изображение видимого движения Солнца в течение года , которое напоминает восьмерку. Аналемму можно изобразить, наложив фотографии, сделанные в одно и то же время суток с интервалом в несколько дней в течение года .
Аналемму также можно рассматривать как график склонения Солнца , обычно построенный вертикально, в зависимости от уравнения времени , построенного горизонтально. Обычно масштабы выбираются так, чтобы равные расстояния на схеме представляли собой равные углы в обоих направлениях на небесной сфере. Таким образом, 4 минуты (точнее 3 минуты 56 секунд) в уравнении времени представлены тем же расстоянием, что и 1° по склонению , поскольку Земля вращается со средней скоростью 1° каждые 4 минуты относительно Солнца. .
Аналемма рисуется так, как ее видел бы на небе наблюдатель, смотрящий вверх. Если север вверху показан запад то справа . , Обычно это делается даже тогда, когда аналемма отмечена на географическом глобусе , на котором континенты и т. д. показаны западом влево.
Некоторые аналеммы отмечены, чтобы показать положение Солнца на графике в разные даты с интервалом в несколько дней в течение года. Это позволяет использовать аналемму для простых аналоговых вычислений таких величин, как время и азимуты восхода солнца и захода . Аналеммы без отметок даты используются для корректировки времени, указываемого солнечными часами . [19]
Световые эффекты
[ редактировать ]Мы видим свет Солнца примерно в 20 угловых секундах от того места, где находится Солнце, когда мы видим свет. См. Солнечную годовую аберрацию .
См. также
[ редактировать ]Примечания
[ редактировать ]- ^ В небесной механике орбита Земли будет называться орбитой Кеплера с эксцентриситетом орбиты менее 1.
Ссылки
[ редактировать ]- ^ Меус, Жан (1991). «Глава 12: Преобразование координат». Астрономические алгоритмы . Ричмонд, Вирджиния: ISBN Willmann Bell, Inc. 0-943396-35-2 .
- ^ Jump up to: а б с д и Дженкинс, Алехандро (2013). «Положение Солнца на небе». Европейский журнал физики . 34 (3): 633–652. arXiv : 1208.1043 . Бибкод : 2013EJPh...34..633J . дои : 10.1088/0143-0807/34/3/633 . S2CID 119282288 .
- ^ Jump up to: а б Чжан Т., Стэкхаус П.В., Макферсон Б. и Миковиц Дж.К., 2021. Формула солнечного азимута, которая делает ненужным косвенное рассмотрение без ущерба для математической строгости: математическая установка, применение и расширение формулы, основанной на подсолнечной точке и функция atan2. Возобновляемая энергия , 172, 1333–1340. DOI: https://doi.org/10.1016/j.renene.2021.03.047.
- ^ Военно-морская обсерватория США; Гидрографическое управление Великобритании, Управление морского альманаха Ее Величества (2008 г.). Астрономический альманах на 2010 год . Правительство США. Типография. п. С5. ISBN 978-0-7077-4082-9 .
- ^ Практически тот же набор уравнений, охватывающий период с 1800 по 2200 год, можно найти в разделе «Приблизительные солнечные координаты» на веб-сайте Военно-морской обсерватории США . Архивировано 31 января 2016 г. на Wayback Machine . графики погрешности этих уравнений по сравнению с точными эфемеридами . Также можно просмотреть
- ^ Меус (1991), с. 152
- ^ Морской альманах Военно-морской обсерватории США (1992). П. Кеннет Зайдельманн (ред.). Пояснительное приложение к Астрономическому альманаху . Университетские научные книги, Милл-Вэлли, Калифорния. п. 12. ISBN 0-935702-68-7 .
- ^ «Избранные астрономические константы, 2015 г. (PDF)» (PDF) . Военно-морская обсерватория США. 2014. с. К6–К7.
- ^ «Избранные астрономические константы, 2015 (TXT)» . Военно-морская обсерватория США. 2014. с. К6–К7.
- ^ Дж. В. Спенсер (1971). «Представление положения Солнца в ряд Фурье» .
{{cite journal}}
: Для цитирования журнала требуется|journal=
( помощь ) - ^ Спроул, Алистер Б. (2007). «Вывод солнечных геометрических соотношений с использованием векторного анализа». Возобновляемая энергия . 32 (7): 1187–1205. doi : 10.1016/j.renene.2006.05.001 .
- ^ «СанАлайн» . Архивировано из оригинала 9 марта 2012 года . Проверено 28 февраля 2012 г.
- ^ «Солнечный калькулятор NOAA» . Лаборатории исследования системы Земли . Проверено 28 февраля 2012 г.
- ^ Jump up to: а б «Детали солнечного расчета» . Лаборатории исследования системы Земли . Проверено 28 февраля 2012 г.
- ^ «Астрономические алгоритмы» . Проверено 28 февраля 2012 г.
- ^ Бланко-Мюриэль, Мануэль; Аларкон-Падилья, Диего С; Лопес-Мораталья, Теодоро; Лара-Койра, Мартин (2001). «Вычисление солнечного вектора» (PDF) . Солнечная энергия . 70 (5): 431–441. Бибкод : 2001SoEn...70..431B . дои : 10.1016/s0038-092x(00)00156-0 .
- ^ Ибрагим Реда и Афшин Андреас. «Алгоритм положения Солнца для приложений, связанных с солнечным излучением» (PDF) . Проверено 28 февраля 2012 г.
- ^ «Приближение атмосферной рефракции» . Национальное управление океанических и атмосферных исследований . Проверено 28 февраля 2012 г.
- ^ Солнечные часы#Отметки полудня
Внешние ссылки
[ редактировать ]- Алгоритм положения Солнца Национальной лаборатории возобновляемых источников энергии на веб-сайте Центра данных по возобновляемым ресурсам .
- Калькулятор положения Солнца на pveducation.org . Интерактивный калькулятор, показывающий путь Солнца по небу.
- Солнечный калькулятор NOAA NOAA Earth System Research Laboratories на веб-сайте отдела глобального мониторинга .
- Калькулятор склонения и положения солнца NOAA
- Система HORIZONS на сайте JPL . Очень точные положения объектов Солнечной системы на основе эфемерид серии JPL DE .
- Общие эфемериды тел Солнечной системы , на сайте IMCCE . Положения объектов Солнечной системы по эфемеридам серии INPOP.
- Солнечное положение в пакете R.Insol.