Napisałem https://github.com/barrycarter/bcapps/blob/master/STACK/bc-solve-astro-13008.c, aby znaleźć historyczne przesilenia i równonoce. Pełne wyniki są dostępne pod adresem:
https://github.com/barrycarter/bcapps/blob/master/ASTRO/solstices-and-equinoxes.txt.bz2
Dokładność:
- Używam modelu precesji EARTH_IAU_1976 oraz modeli nutacji i pochylenia EARTH_IAU_1980, które są jedynymi modelami obsługiwanymi przez bibliotekę CSPICE, której używam:
http://hesperia.gsfc.nasa.gov/ssw/stereo/gen/exe/icy64/doc/frames.req
To to samo biblioteka i te same modele, których NASA używa do wykonywania własnych obliczeń, ale modele są znane z ograniczonej dokładności. Cytując „ http://aa.usno.navy.mil/faq/docs/SpringPhenom.php”:
Czasy podane w tabelach są z dokładnością do dwóch lub trzech godzin przez 25 do 5 pne i jedną lub dwie godziny przez 4 do 38 roku n.e. Niepewność w tych czasach wynika z gwałtownej, czyli nieprzewidywalnej zmiany szybkości obrotu Ziemi.
Dokładność pogarsza się im dalej wstecz LUB dalej idziesz: nie tylko nie wiemy dokładnie o rotacja Ziemi w przeszłości, ale nie możemy nawet dokładnie przewidzieć rotacji Ziemi w przyszłości.
Przykładowe dane wyjściowe:
EQU 511720269.432607 AD 2016-03-20 04: 30: 01SOL 519734120.174820 AD 2016-06-20 22: 34: 11 EQU 527826138.004142 AD 2016-09-22 14: 21: 09 SOL 535589116.137776 AD 2016-12-21 10:44:07
Format:
-
Pierwsza kolumna wskazuje, czy jest to przesilenie, czy równonoc.
-
Druga kolumna to czas efemeryd przesilenie / równonoc. Jeśli zajmujesz się poważnymi pracami astronomicznymi, to jest ta kolumna, której powinieneś użyć. Przybliżona definicja to liczba sekund od epoki J2000. Dokładna definicja stosowana w CSPICE to: https://naif.jpl.nasa.gov/pub/naif/toolkit_docs/FORTRAN/req/time.html#Ephemeris%20Time%20(ET)
-
Pozostałe kolumny wskazują czas UTC przesilenia / równonocy w formacie bardziej czytelnym dla człowieka:
- Biblioteki CSPICE zakładają, że reformacja kalendarza gregoriańskiego nastąpiła 4 października 1582 r. (co oznacza dzień po 4 października 1582 był 15 października 1582). Patrząc na te wiersze:
SOL -13191695511.794357 AD 1581-12-11 20: 07: 27 EQU -13183992131.003845 AD 1582-03-10 23 : 57: 07SOL -13175951250.170158 AD 1582-06-12 01: 31: 48 EQU -13167875920.223862 AD 1582-09-13 12: 40: 38 SOL -13160138634.917915 AD 1582-12-22 01: 55: 23 EQU -13152434815.793312 AD 1583-03- 21 05: 52: 23SOL -13144394485.035870 AD 1583-06-22 07: 17: 53EQU -13136319216.460808 AD 1583-09-23 18:25:42
widać, że daty przesilenia / równonoce skok 10 dni na reformację.
-
Biblioteki CSPICE używają kalendarza juliańskiego sprzed 4 października 1582 roku. W rzeczywistości kalendarz juliański został wprowadzony w 46 roku pne: https://en.wikipedia.org/wiki/Julian_calendar
-
Przed 46 rokiem pne w użyciu były inne systemy kalendarzy, ale biblioteki CSPICE zakładają kalendarz julian wraca na czas nieokreślony:
https://en.wikipedia.org/wiki/Proleptic_Julian_calendar
Moje obliczenia sięga 13201 roku p.n.e. (ograniczenia DE431, efemerydy, których używam) i możliwe, że ludzie nie używali wtedy nawet regularnie kalendarzy: cytując „ https://en.wikipedia.org/ wiki / History_of_calendars # Prehistory ":
Mezolityczny układ dwunastu dołów i łuku znaleziony w Warren Field w Aberdeenshire w Szkocji, datowany na około 10000 lat temu, został opisany jako kalendarz księżycowy i nazwany „najstarszym znanym kalendarzem świata” w 2013 r.
Uwagi:
-
Obrót Ziemi wokół Słońca względem gwiazd stałych zajmuje 365.256363004 dni, ale czas między równonocą wiosenną jest nieco krótszy (365,242190402 dni), ponieważ pozycja równonocy wiosennej przesuwa się
(precesje) względem gwiazd. Źródło: http://hpiers.obspm.fr/eop-pc/models/constants.html
-
Średnia długość dnia w kalendarzu gregoriańskim to 365,2425 jest znacznie bliższa 365,242190402 dni niż średnia długość dnia 365,25 kalendarza Julian, ale nadal nie jest idealna. Jak zauważono w Czy przesilenia i równonoce zmieniają się w czasie?, jeśli nadal będziemy używać kalendarza gregoriańskiego w dalekiej przyszłości, równonoce i przesilenia będą dryfować wstecz. Do 17090 (limit DE431) będą wyglądać tak:
EQU 476198945887.238159 AD 17090-02-22 08: 56: 59SOL 476207018540.040894 AD 17090-05 -26 19: 21: 11EQU 476214808218.146362 AD 17090-08-24 23: 09: 09SOL 476222655067.609985 AD 17090-11-23 18:49:59
około miesiąc za „zwykłym” razy.
MOJA WCZEŚNIEJSZA ODPOWIEDŹ CZĘŚCIOWA W CELU:
Od HORYZONÓW ( http://ssd.jpl.nasa.gov/?horizons) i SPICE ( http://naif.jpl.nasa.gov/naif/tutorials.html) może obliczyć pozycję ekliptyki i Słońca z tamtych czasów, powinno być możliwe obliczenie równonocy i przesilenia z rozsądną dokładnością. Jednak nie udało mi się znaleźć witryny, która faktycznie zawiera te daty (jestem prawie pewien, że USNO zrobiło to w pewnym momencie, ale nie mogę znaleźć ich listy). Inne potencjalnie pomocne źródła / pytania: