Astronomia Sferyczna
Astronomia Sferyczna
Astronomia Sferyczna
z
z1
θ x1
O
ψ
φ y
N
x
6 Przemiana współrz˛ednych 49
6.1 Macierzowe transformacje współrz˛ednych astronomicznych . . . . . . . . . . . . . . . . 49
6.1.1 Układ ekliptyczny i układ równikowy . . . . . . . . . . . . . . . . . . . . . . . . 49
6.1.2 Układ horyzontalny i godzinny . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
6.1.3 Układ godzinny i równikowy . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
6.1.4 Układ równikowy i galaktyczny . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
6.2 Dodatek A. Nastawianie teleskopów . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52
6.3 Dodatek B. Zadania . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
8 Refrakcja 72
8.1 Wst˛ep – miejsca topocentryczne ciał niebieskich . . . . . . . . . . . . . . . . . . . . . . 74
8.2 Refrakcja — model płaskiej atmosfery . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
8.3 Wpływ refrakcji na współrz˛edne równikowe obiektu . . . . . . . . . . . . . . . . . . . . 78
8.4 Refrakcja — model atmosfery radialnie symetrycznej . . . . . . . . . . . . . . . . . . . . 79
8.5 Stałe refrakcji. Tablice refrakcji . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 85
8.6 Dodatek A. Zadania . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
9 Współrz˛edne geocentryczne 88
9.1 Współrz˛edne geocentryczne obserwatora . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
9.2 Paralaksa geocentryczna . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 92
9.3 Dygresja. Jednostka astronomiczna — 1 AU . . . . . . . . . . . . . . . . . . . . . . . . . 94
9.3.1 Wpływ paralaksy geocentrycznej na wspólrz˛edne równikowe . . . . . . . . . . . 96
9.4 Aberacja dobowa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
9.5 Dodatek A. Zadania . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
Streszczenie. Kierunki do ciał niebieskich można określać równnież przy założeniu, że
ciała te znajduja˛ si˛e na sferze o promieniu równym jedności. Geometria sfery jest geo-
metria˛ na zakrzywionej powierzchni dwuwymiarowej i w wielu przypadkach różni si˛e od
euklidesowej geometrii na płaszczyźnie. Elementami sfery wykorzystywanymi w astro-
nomii sa˛ koła wielkie i skonstruowane z ich pomoca˛ dwukaty ˛ i trójkaty
˛ sferyczne. W
geometrii na sferze koła wielkie pełnia˛ rol˛e analogiczna˛ do prostych w planimetrii. Suma
katów
˛ ˛ sferycznym jest zawsza wi˛eksza od π. Na sferze moga˛
wewn˛etrznych w trójkacie
istnieć trójkaty,
˛ w których wszystkie katy ˛ sa˛ katami
˛ prostymi. Elementy trójkata˛ sferycz-
nego spełniaja˛ kilka grup równań pozwalajacych˛ na rozwiazywanie
˛ wielu zagadnień z
zakresu astronomii sferycznej. Do najcz˛eściej wykorzystywanych wzorów należa˛ wzory
sinusów i cosinusów. Problemy rozwiazywane ˛ w ramach trygonometrii sferycznej daja˛
si˛e także ujać
˛ w formaliźmie wektorowym.
Położenie ciała na sferze ustalone jest z pomoca˛ dwóch katów,
˛ azymutalnego ψ i polar-
nego θ. Każdy układ współrz˛ednych sferycznych wymaga określenia bieguna układu,
wzgl˛edem którego mierzony jest kat ˛ θ, oraz koła wielkiego pełniacego
˛ rol˛e płaszczyzny
odniesienia, służacej
˛ jako poczatek˛ rachuby dwuściennego kata ˛ ψ.
W astronomii wykorzystywane sa˛ najróżniejsze układy współrz˛ednych, stad ˛ konieczna
jest umiej˛etność przeliczania współrz˛ednych pomi˛edzy dwoma układami. Transformacje
współrz˛ednych dokonuje si˛e droga˛ rozwiazania
˛ odpowiedniego trójkata
˛ sferycznego albo
za pomoca˛ macierzy obrotu, w szczególności z wykorzystaniem katów ˛ Eulera.
Słowa kluczowe: sfera niebieska, koło małe, koło wielkie, bieguny koła wielkiego, kat ˛
sferyczny, dwukat˛ sferyczny, trójkat
˛ paralaktyczny, nadmiar sferyczny, współrz˛edne sfe-
ryczne, współrz˛edne prostokatne,
˛ triada ortogonalna, transformacje współrz˛ednych, katy
˛
Eulera. a
a
[Modyfikowano AD 2011, luty, 28]
6 Położenia ciał na sferze niebieskiej
Stwierdzenie, że wektor r opisuje położenie ciała niebieskiego oznacza, że do naszej
dyspozycji sa˛ trzy składowe (x, y, z), czyli trzy liczby wyznaczone wzgl˛edem triady R.
Na piśmie wyrażamy to stwierdzenie za pomoca˛ zapisu
x
r = R y (1.1)
z
Umawiamy si˛e, że otaczajaca ˛ nas przestrzeń jest trójwymiarowa˛ przestrzenia˛ euklide-
sowa,˛ że w tej przestrzeni kartezjański układ współrz˛ednych (x, y, z) rozpi˛ety na trójce
wersorów i, j, k. 1 T˛e trójk˛e możemy połaczyć
˛ w formie wektora blokowego R nazywa-
nego triada˛ o strukturze identycznej z macierza˛ 3 × 3
R = [i, j, k] (1.2)
i=j×k
j=k×i
k=i×j (1.3)
iT i = jT j = kT k = 1
iloczyn
iT i iT j iT k 1 0 0
RT R = j T i j T j j T k = 0 1 0 = I (1.5)
kT i k T j k T k 0 0 1
1
A w jaki sposób ustalono orientacj˛e trójki wersorów i, j, k? Pytanie to wykracza poza ramy tego
wykładu, bowiem pośrednio dotyczy spraw ostatecznych podobnie jak pytanie — czy najpierw było jajo,
czy kura? Ponieważ astronomia nie zajmuje si˛e takimi problemami, dlatego poprzestajemy na tym, że
elementy tej trójki sa,˛ jakimś cudem znanymi, cosinusami kierunkowymi prostych, wzdłuż których leża˛
wersory i, j, k.
1.2 Wektorowe transformacje współrz˛ednych sferycznych 7
a)
a
cos θ
k k
k1
k 1
j1 j1
j
i
i1 θ
θ
j
sin θ
Rysunek 1.1: a) Wektor a może mieć składowe określone wzgledem wielu trójek werso-
rów i, j, k. b) Zwiazki
˛ pomi˛edzy elementami triad różniacych
˛ si˛e obrotem wokół jednej
˛ si˛e z wersorem i. Wersor j1 można otrzymać za pośrednictwem werso-
osi pokrywajacej
rów j, k odkładajac
˛ w kierunku j odcinek o długości cosθ po czym w kierunku k odcinek
sin θ.
i1 = i
j1 = (cos θ)j + (sin θ)k
k1 = (− sin θ)j + (cos θ)k
˛ θ wokół osi i. A w
co oznacza, że oba układy różnia˛ si˛e jedynie o dodatni obrót o kat
postaci dogodniejszej do umacierzowienia zapisu b˛edzie
i1 = 1i +0j +0k
j1 = 0i +(cos θ)j +(sin θ)k
k1 = 0i +(− sin θ)j +(cos θ)k
Zatem, zgodnie z definicja˛ triady, triad˛e R1 możemy określić jako iloczyn macierzowy
1 0 0
R1 = [i1 , j1 , k1 ] = R 0 cos θ − sin θ
0 sin θ cos θ
x1 x
y1 = p(θ) y (1.9)
z1 y
gdzie
1 0 0
p(θ) = 0 cos θ sin θ (1.10)
0 − sin θ cos θ
A
r z
θ
X x O
y
ψ
B Y
Z’
cos θ sin θ 0
r(θ) = − sin θ cos θ 0 (1.12)
0 0 1
W przypadku gdy układy współrz˛ednych różnia˛ si˛e dwoma lub trzema obrotami, macierz
transformacji współrz˛ednych wyznaczona jest jako iloczyn dwóch lub trzech macierzy
odpowiadajaych
˛ pojedynczym obrotom.
W praktyce możemy napotkać przypadek takich układów, dla których korzystnym b˛edzie
zastosowanie macierzy Mx , pozwalajacej˛ na lustrzane odbicie wzgl˛edem płaszczyzny
y − z czyli na zmian˛e znaku współrz˛ednej x-owej.
−1 0 0
Mx = 0 1 0 (1.14)
0 0 1
z
z1
θ x1
O
ψ
φ y
N
x
x = cos ψ cos κ
y = sin ψ cos κ (1.15)
z = sin κ
κ = arcsin z
(1.16)
ψ = arctan xy
Ale możliwe jest inne podejście, w którym wzajemna orientacja dwóch układów
współrz˛ednych określona jest za pomoca˛ tzw. katów
˛ Eulera, natomiast transformacja
składa si˛e z trzech obrotów wzgl˛edem osi w ustalonej, zawsze takiej samej kolejności.
Katy
˛ Eulera. Na rysunku 1.3 zaznaczono trzy katy
˛ Eulera wykorzystywane do trans-
formacji współrz˛ednych wyznaczonych wzgl˛edem obu układów, sa˛ to:
• kat
˛ θ zawarty pomi˛edzy osiami Z i Z1 obu układów,
• kat
˛ φ zawarty pomi˛edzy osia˛ X i linia˛ ON przeci˛ecia płaszczyzn X − Y obu
układów,
• kat
˛ ψ pomi˛edzy osia˛ X1 i linia˛ ON , liczony jako dodatni od linii ON do osi
X1 .
Sfera.
Sfera jest to powierzchnia, której punkty sa˛ jednakowo odległe od punktu zwanego środ-
kiem sfery. Jeśli w środku sfery umieścimy poczatek
˛ układu współrz˛ednych, to dla punk-
tów położonych na sferze jednostkowej b˛edzie
rT r = 1
Koło wielkie.
Każde przeci˛ecie sfery płaszczyzna˛ jest okr˛egiem. Przeci˛ecie sfery płaszczyzna˛ przecho-
˛ a˛ przez środek sfery jest kołem wielkim, np. koło AX z rysunku 2.1. Promień koła
dzac
wielkiego jest równy promieniowi sfery, a wektory r opisujace ˛ położenia punktów koła
wielkiego spełniaja˛ równanie
nT r = 0 (2.1)
Dwa punkty sfery, które nie sa˛ punktami diametralnymi jak np. A i X, wyznaczaja˛
koło wielkie jednoznacznie, bowiem łacznie
˛ ze środkiem sfery (punktem O) jednoznacz-
nie określaja˛ płaszczyzn˛e, której przeci˛ecie ze sfera˛ jest kołem wielkim. Na punktach
A, X rozpi˛ete sa˛ dwa łuki, mniejszy z nich nazywany linia˛ geodezyjna,˛ jest najkrótsza˛
2.1 Obiekty geometryczne na sferze 13
A O B
Ponieważ promień sfery r = 1, długość łuku koła wielkiego równa jest katowi
˛ w
radianach1 jaki ten łuk rozpina wzgl˛edem środka sfery.
Dwukat˛ sferyczny.
W wyniku przeci˛ecia tej samej sfery dwoma kołami wielkimi wyznaczone zostana˛ cztery
obszary (powierzchnie) zwane dwukatami ˛ sferycznymi (rysunek 2.2-a). Dwukaty˛ prze-
ciwległe sa˛ parami przystajace.
˛ Dwukat ˛ sferyczny określony jest katem
˛ sferycznym, np.
katem
˛ A na rysunku 2.2-a. Kat ˛ ten jest równy katowi
˛ płaskiemu (liniowemu) określonemu
przez płaszczyzny, na których leża˛ koła wielkie tworzace ˛ dany dwukat. ˛ Półokr˛egi tych
kół wielkich nazywamy bokami dwukata. ˛ Na danej sferze, boki wszystkich dwukatów ˛
2
sa˛ równe i maja˛ długość πr, gdzie r jest promieniem sfery. Pole pwierzchni dwukata ˛
sferycznego można obliczyć ze wzoru
S = 2r2 A
Kat
˛ pomi˛edzy płaszczyznami kół wielkich cz˛esto jest nazywany katem
˛ sferycznym.
O kacie
˛ sferycznym można też powiedzieć, że jest to kat
˛ pomi˛edzy stycznymi wystawio-
nymi w punkcie wzajemnego przeci˛ecia si˛e kół wielkich (patrz rysunek 2.2-b).
Trójkat˛ sferyczny.
Trzy koła wielkie nie przecinajace
˛ si˛e w jednej parze punktów diametralnych tworza˛ na
sferze osiem obszarów zwanych trójkatami˛ sferycznymi (rysunek 2.3-a). Znajac
˛ elementy
˛ sferyczne, np. łuki a, b, c kół wielkich i katy
jednego z nich (czyli trzy boki i trzy katy ˛
wewn˛etrzne A, B, C), łatwo wyznaczyć elementy wszystkich pozostałych trójkatów.˛ Dla-
tego zwykle rozpartuje si˛e zależności pomi˛edzy elementami tylko jednego trójkata,
˛ któ-
1
Oprócz radianów można oczywiście używać innych miar kata.
˛
2
Ze wzgl˛edu na dość cz˛este nieporozumienia, warto zapami˛etać, że dwukat
˛ sferyczny (podobie trójkat
˛
sferyczny) jest fragmentem powierzchni sfery, a nie katem.
˛
14 Elementy geometrii sferycznej
rego wszystkie boki sa˛ krótsze od połowy obwodu koła wielkiego. Taki trójkat
˛ sferyczny
nosi nazw˛e trójkata
˛ paralaktycznego lub trójkata
˛ eulerowskiego.
Boki trójkata
˛ sferycznego tradycyjnie oznaczane sa˛ małymi literami, a ich długości
mierzone sa˛ za pomoca˛ płaskich katów
˛ ˛ trójściennego OABC, np. na rysunku 2.3
kata
c = ]AOB. Wewn˛etrzne katy ˛ przy wierzchołkach trójkata˛ oznaczane dużymi literami
A, B, C mierzone sa˛ katami
˛ dwuściennymi (katami
˛ sferycznymi) tego samego kata
˛ trój-
ściennego.
A
A
O
A
A’
A’
A
c b
c b
O
a B
B C C
a
Trójkaty
˛ sferyczne eulerowskie maja˛ pewne własności wspólne z trójkatami
˛ płaskimi,
np. dowolny bok trójkata
˛ jest mniejszy od sumy, ale wi˛ekszy od różnicy dwóch boków
pozostałych. Jednak mamy też mi˛edzy nimi istotne różnice, np. w trójkacie
˛ sferycz-
nym suma katów
˛ wewn˛etrznych nie jest stała, bowiem przyjmuje wartości z przedziału
2.1 Obiekty geometryczne na sferze 15
P
θ ψ
A rM B
S
ale
E kolo m
C O D
ψ
Rysunek 2.4: Elementy sfery: koło małe AEB. Jest ono odległe o kat
˛ θ od swego bieguna
P . Wycinek AE koła małego opiera si˛e ramionach AS i ES rozwartych o kat
˛ ψ. Długość
wycinka wynosi AE = ψ sin θ.
przedziału
π < s = A + B + C < 3π
Trójkat
˛ płaski posiada tylko jeden kat
˛ prosty, trójkat
˛ sferyczny niekoniecznie, może mieć
ich dwa a nawet trzy. Różnica
s−π =ε
nazywanajest nadmiarem sferycznym.
Koło małe.
Ślad przeci˛ecia sfery płaszczyzna˛ nie przechodzac
˛ a˛ przez środek sfery jest okr˛egiem, tra-
dycyjnie zwanym kołem małym. Jego biegunami sa˛ punkty skrajne średnicy sfery wysta-
wionej prostopadle do płaszczyzny koła małego.
Promień koła małego jest zawsze mniejszy od promienia sfery. Na rysunku 2.4 wi-
dzimy koło małe AEB, równoległe doń koło wielkie CF D oraz ich wspólne w tym
przypadku bieguny P i Q.
Wyprowadzimy formuł˛e na długość łuku koła małego, cz˛esto stosowana˛ w dalszej cz˛eści
wykładu. Niech rM b˛edzie promieniem koła małego, łuk AP = θ jest miara˛ odległości
koła małego od bieguna P . Z trójkata
˛ płaskiego AOS wynika
AS = AO · sin AOS
(2.2)
rM = sin θ
Punkt E koła małego połaczmy
˛ z biegunem P łukiem koła wielkiego i powstały łuk prze-
dłużmy do przeci˛ecia z kołem CD w punkcie F . Jeśli kat
˛ sferyczny AP E oznaczymy
16 Elementy geometrii sferycznej
Z ψ Z
S
R
A θ
A
r z r
θ C
X x O
O
y
ψ
ψ
B Y B
Z’ Z’
W celu ustalenia położenia punktu na sferze, można wykorzystać różne układy współrz˛e-
dnych. Przykładowo, weźmy prawoskr˛etny prostokatny ˛ zbiór osi kartezjańskich Oxyz,
określonych trójka˛ wersorów i, j, k o poczatkach
˛ w środku O sfery jednostkowej. Do-
datnie kierunki tych osi przecinaja˛ sfer˛e w punktach X, Y, Z, natomiast koła wielkie
XY, ZX położone sa˛ w płaszczyznach xOy i zOx, odpowiednio, patrz (rysunek 2.5).
Obierzmy na sferze punkt A(xyz), wówczas. dla sfery jednostkowej prawdziwy jest
zwiazek
˛
x2 + y 2 + z 2 = 1 (2.4)
Czyli jedna z prostokatnych
˛ współrz˛ednych x, y, z punktu A jest zbyteczna, co oznacza,
że w celu ustalenia położenia ciała niebieskiego, tzn. kierunku do tego ciała, wystarczy
posłużyć si˛e dwoma liczbami.
Ponieważ dla sfery jednostkowej promień sfery r = 1, zatem dwie współrz˛edne katowe ˛
(ψ, θ) w pełni określaja˛ położenie ciała na sferze: kat
˛ θ jest długościa˛ łuku ZA, natomiast
ψ jest katem
˛ sferycznym XZA.
0≤θ≤π
(2.5)
0 ≤ ψ ≤ 2π
Pomimo nieodłacznego
˛ nadmiaru (patrz równanie (2.4)), w astronomii sferycznej warto
stosować współrz˛edne kartezjańskie. Uj˛ete w postaci uporzadkowanych
˛ trójek współrz˛e-
dne te, pozwalaja˛ nadać równaniom elegancka˛ i uniwersalna˛ wektorowa˛ form˛e.
rA = x · i + y · j + z · k (2.6)
gdzie składowe (x, y, z) sa˛ cosinusami kierunkowymi odcinka OA, obliczonymi odpo-
wiednio, wzgl˛edem osi X, Y, Z, patrz rysunek 2.5
x = cos XA
y = cos Y A (2.7)
z = cos ZA
x = sin θ cos ψ
y = sin θ sin ψ (2.8)
z = cos θ
rB = (sin c, 0, cos c)
(3.1)
rC = (sin b cos A, sin b sin A, cos b)
rB · rC = cos a
A
c
b
B a
rA
rB C
rC
90 O
rD
sin a rD = sin a(sin AD cos BAD, sin AD sin BAD, cos AD) (3.5)
Moglibyśmy teraz porównać odpowiednie składowe w równaniach (3.4) i (3.5), jednak
warto przedtem pozbyć si˛e sinusów i cosinusów katów
˛ AD i BAD.
Biorac
˛ iloczyn wektorowy dla innej pary wektorów, równanie to można uzupełnić o do-
datkowy człon, i w tej pełnej postaci nosi ono nazw˛e wzoru sinusów
sin A sin B sin C
= = (3.6)
sin a sin b sin c
Równości (3.6) wykorzystamy do wyprowadzenia kolejnych wzorów trygonometrii sfe-
˛ sferycznym BAD z rysunku 3.1, dla boku AD i kata
rycznej. W trójkacie ˛ wierzchołko-
wego ADB na mocy wzoru sinusów b˛edzie
Pozostałe pi˛eć wzorów typu (3.7) otrzymamy dzi˛eki odpowiednim permutacjom symboli
˛ sferycznym ABC. Np. zmieniajac
w trójkacie ˛ symboli aBbcbcA rol˛e B z C oraz
˛ w ciagu
b z c dostaniemy nast˛epny wzór pi˛ecioelementowy
Z wektorowego równania (3.3) nie otrzymamy już żadnych nowych formuł. Ostatni
ważny wzór trygonometrii sterycznej tzw. wzór cotangensowy (czterocz˛eściowy) można
wydedukować ze wzorów cosinusów i sinusów. W tym celu stosujmy wzór cosinusów do
boków b i c trójkata
˛ ABC (rysunek 3.1), mamy
X X’
ds
P
s
α−αο O θ
so
90−δ
0
90−δ
C s x so
U
X’
χ
L
X
γ O θ C
Rysunek 3.2: Małe przesuni˛ecie na sferze niebieskiej. a) Ciało niebieskie uległo drob-
nemu przesuni˛eciu z punktu X do X 0 , wzdłuż koła wielkiego OX. b) wersja wektorowa
małego przesuni˛ecia.
Wielkość przesuni˛ecia może być różna, m.in. jest zależna od położenia obiektu, ale
zawsze przesuni˛ecie odbywa si˛e po kole wielkim łacz ˛ acym
˛ dany obiekt z jakimś ustalo-
nym punktem sfery, wspólnym dla wszystkich obiektów. Np. małe przesuni˛ecie zwane
paralaksa˛ roczna,˛ ma miejsce zawsze wzdłuż koła wielkiego zawierajacego ˛ kierunki ku
Słońca i do danego ciała niebieskiego, przesuni˛ecie zwane aberacja˛ dobowa,˛ przebiega
po kole wielkim rozpi˛etym na kierunku do danego ciała i na kierunku do punktu wschodu
horyzontu miejsca obserwacji. Wszystkie tego typu przesuni˛ecia można traktować jako
szczególne przypadki ogólniejszego małego przesuni˛ecia opisanego poniżej.
Przyjmijmy, że np. kierunek do gwiazdy X(α, δ) z rysunku 3.2a, uległ niewielkiemu
przesuni˛eciu do punktu X 0 , oraz że odbyło si˛e to wzdłuż koła wielkiego łacz
˛ acego
˛ Xz
punktem O(α0 , δ0 ). Oznaczmy łuk OX przez θ a łuk XX 0 przez dθ. Załóżmy, że dθ jest
małym katem
˛ dodatnim. Niech dalej b˛edzie, że przesuni˛ecie dθ da si˛e wyrazić za pomoca˛
formułki
XX 0 = dθ = k sin θ (3.10)
gdzie k jest stała˛ dodatnia˛ lub ujemna˛ niezależna˛ od obranej gwiazdy, co zreszta˛ nie jest
aż tak ważne dla dalszego wywodu. W ten sposób umówiliśmy si˛e, że interesujemy si˛e
przesuni˛eciami, których wielkość dla danego k, zależy jedynie od odległości obiektu od
pewnego punktu wspólnego O.
˛ przez X 0 , przecinajace
Poprowadźmy koło małe o biegunie w P , przechodzace ˛ PX
0
w punkcie U . Możemy zatem utworzyć trójkat
˛ U XX , dla którego poszukamy wyrażeń
na niektóre jego elementy.
Poszukamy teraz wyrażeń na długości tych boków, przy czym chcemy by wyst˛epo-
wały w nich wielkości możliwie bezpośrednio zwiazane
˛ z parametrami opisujacymi
˛ małe
0 ◦
przesuni˛ecie. Oznaczmy kat ˛ sferyczny OXP przez χ, wówczas U XX = 180 − χ.
Ze wzgl˛edu na bardzo małe rozmiary w stosunku do promienia sfery, trójkat ˛ U XX 02 w
przybliżeniu można traktować jako trójkat ˛ prostym w wierzchołku U . Dla
˛ płaski, o kacie
takiego przybliżenia, mamy
Podstawiajac
˛ prawe strony tych wyrażeń do równań (3.11) ostatecznie mamy
Opis małego przesuni˛ecia dany równaniami (3.12) można również otrzymać w formie
wektorowej. Niech s b˛edzie wersorem położenia punktu X a s0 b˛edzie wersorem poło-
żenia punktu O (rysunek 3.2b). Ponieważ sa˛ to wektory jednostkowe stad ˛ wektor s × s0
ma długość sin θ, i skierowany jest do punktu L na sferze, odległego o 90◦ zarówno od O
jak i od X. Punkt L jest wi˛ec biegunem koła wielkiego OX.
s · ds = 0 (3.13)
ds = s0 − s = k s × (s × s0 ) (3.14)
Równanie to jest bardziej ogólne niż równanie (3.12), bowiem dotyczy ono dowolnego
układu współrz˛ednych, nie tylko układu równikowego.
Analizujac
˛ teoretyczne problemy na sferze, wygodnie jest wyrażać katy ˛ w radianach.
Kiedy jednak wykorzystujemy rezultaty analiz, opłaca si˛e stosować inne miary katów,
˛
np. małe katy
˛ takie jak kat
˛ paralaksy najcz˛eściej podane sa˛ w sekundach łuku.
sin θ ≈ θ
cos θ ≈ 1 (3.16)
tan θ ≈ θ
Wyrażenie to znakomicie nadaje si˛e do zamiamy radianów na sekundy łuku. Np. jeśli
zapis θ00 oznacza liczb˛e sekund w małym kacie
˛ θ, to
θ00
sin θ ≈ θ ≈ = θ00 sin 100
206265
A0 = 180◦ − a a0 = 180◦ − A
B 0 = 180◦ − b b0 = 180◦ − B
C 0 = 180◦ − c c0 = 180◦ − C
28 Astronomiczne układy współrz˛ednych
Rozdział 4
Umówiliśmy si˛e, że sfera niebieska to sfera o promieniu jednostkowym. O jej środku po-
wiedzieliśmy jedynie tyle, że jest tam gdzie znajduje si˛e obserwator, np. na powierzchni
Ziemi. Skoro tak to w dowolnym miejscu gdzie znajduje si˛e obserwator da si˛e rozpiać ˛
sfer˛e niebieska˛ i skojarzyć z nia˛ jakiś układ współrz˛ednych. Dlatego mówimy o:
• sferze niebieskiej topocentrycznej czyli o środku na powierzchni Ziemi,
• sferze niebieskiej geocentrycznej o środku w centrum Ziemi,
• sferze niebieskiej heliocentrycznej o środku w centrum Słońca,
• sferze niebieskiej lunocentrycznej o środku w centrum Ksi˛eżyca,
• etc.
Rysunek 4.2: Układ współrz˛ednych sferycznych równikowych (α, δ). Deklinacja δ jest
miara˛ wysokości obiektu nad równikiem niebieskim, rektascensja α określa katow
˛ a˛ od-
ległość astronomicznego południka obiektu, od południka przechodzacego˛ przez punkt
równonocy wiosennej Υ. Dla obserwatora znajdujacego ˛ si˛e na północnym biegunie P
świata, kierunek rachuby rektascensji jest przeciwny do ruchu wskazówek zegara. Jest to
kierunek zgodny z pozornym, rocznym ruchem Słońca na sferze.
gdzie Υ oznacza punkt na równiku niebieskim, pełni on rol˛e punktu zerowego miary rek-
tascensji. Wybrano go droga˛ konwencji tak, by znajdował si˛e możliwie blisko położenia
Słońca w momencie równonocy wiosennej (około 21 marca), kiedy to Słońce przechodzi
przez równik niebieski z półsfery południowej na półsfer˛e północna.˛ Dlatego punkt Υ
nazywany jest punktem równonocy. Ponieważ w czasach kiedy zaproponowano t˛e kon-
wencj˛e punkt Υ położony był w gwiazdozbiorze Barana, nazwano go również punktem
barana.
Współrz˛edna rektascensja mierzona jest wzdłuż równika, w kierunku zgodnym z kie-
runkiem pozornego ruchu Słońca na sferze, czyli w kierunku antyzegarowym dla ob-
serwatora znajdujacego
˛ si˛e na północnym biegunie świata P . Deklinacj˛e mierzymy w
płaszyźnie południka niebieskiego. Wartości jakie moga˛ przyjmować obie współrz˛edne
należa˛ do przedziałów
−90◦ ≤ δ ≤ 90◦
(4.5)
0◦ ≤ α ≤ 360◦
Jednak tradycyjna˛ miara˛ rektascensji nie sa˛ stopnie ale jednostki czasu. Pomi˛edzy jed-
nostkami czasowymi i katowymi
˛ mamy proste zależności, mianowicie, jeżeli przyjmiemy,
h ◦
że 24 = 360 , to łatwo sprawdzić, że
1h = 15◦ 1◦ = 4m
1m = 150 10 = 4s (4.6)
1s = 1500 100 = 1/15s
a) Z Z
b) Q
P
Rownik Rownik
R X R
W W
z z
X
O
S N S N
O A A
Horyzont Horyzont
E E
Poludnik Poludnik
Q P
Na Na
x = cos δ cos α
y = cos δ sin α (4.7)
z = sin δ
Z U
P H T
Y
R
V X
E
R
N δ S
O
L
D
W
Rownik
Poludnik
Istnieje także obszar sfery, który nigdy nie jest widoczny dla danego obserwatora. Na
mocy symetrii odpowiedni warunek ma postać
P ZX = 360◦ − A ZP X = H
ZX = z P X = 90◦ − δ
Mamy także, że P Z = 90◦ − φ, gdzie φ jest szerokościa˛ miejsca obserwacji. Stosujac
˛
˛ P ZX wzór cosinusów dostaniemy
dwukrotnie do trójkata
sin δ = cos z sin φ + sin z cos φ cos A (4.14)
cos z = sin δ sin φ + cos δ cos φ cos H (4.15)
Równania te wystarczaja˛ do przeliczenia współrz˛ednych horyzontalnych na godzinne i
odwrotnie. Problemy normalizacyjne towarzyszace˛ obliczeniom katów
˛ A i H, można
3
Czy aby na pweno?
4
Rozważania te dotycza˛ wyłacznie
˛ obserwatorów z półkuli północnej.
4.6 Współrz˛edne ekliptyczne 37
90−φ
P Z
H
360−A
90
z
−δ
X
P
ε α
K
λ 90
−β
ε
γ Rownik
tyka
lip
Ek
β = 90◦ − KX
(4.17)
λ = ΥKX
przy czym
−90◦ ≤ β ≤ 90◦
0 ≤ λ ≤ 360◦
Długości ekliptyczne rosna˛ w kierunku ruchu prostego, dla Słońca współrz˛edna ta wzrasta
monotonicznie. W przypadku planet w efekcie złożenia ich ruchów prostych z ruchem
orbitalnym Ziemi, ruch wypadkowy dla obserwatora na powierzchni Ziemi może okazać
si˛e ruchem wstecznym.
Zwiazki
˛ pomi˛edzy współrz˛ednymi ekliptycznymi i równikowymi można łatwo wy-
prowadzić rozwiazuj
˛ ac ˛ P KX z rysunku 4.6. Niech obiekt X, obok współrz˛e-
˛ trójkat
dnych ekliptycznych (λ, β) ma współrz˛edne równikowe (α, δ). Boki trójkata ˛ P KX wy-
nosza˛
KP = ε P X = 90◦ − δ KX = 90◦ − β
P KX = 90◦ − λ KP X = 90◦ + α
˛ KP X, posługujac
Z pomoca˛ tych pi˛eciu elementów trójkata ˛ si˛e standardowymi wzorami
trygonometrii sferycznej można przeliczać współrz˛edne z jednego układu do drugiego.
Bardziej ogólne i bezpośrednie podejście do tego zagadnienia wymaga zastosowa-
nia współrz˛ednych prostokatnych
˛ (patrz rysunek 6.1). Prostokatny
˛ układ współrz˛ednych
równikowych (x, y, z) określony jest przez wybór osi x w kierunku punktu Υ, osi z w
kierunku bieguna świata P oraz osi y tak by układ był prawoskr˛etny.
˛ układ współrz˛ednych ekliptycznych (ξ, η ζ) ma oś ξ skierowana˛ w kie-
Prostokatny
runku punktu Υ, oś ζ w kierunku bieguna ekliptyki K natomiast oś η skierowana jest ku
punktowi o współrz˛ednych (λ = 90◦ , β = 0).
4.7 Współrz˛edne galaktyczne 39
ξ=x x=ξ
η = y cos ε + z sin ε y = η cos ε − ζ sin ε (4.19)
ζ = −y sin ε + z cos ε z = η sin ε + ζ cos ε
x = cos δ cos α
y = cos δ sin α
z = sin δ
Podstawiajac
˛ te wyrażenia oraz (4.18) do równań (4.19) otrzymamy zwiazki
˛
θ P
G
X
l V
G G
b
Centrum Galaktyki Y
S
Równik Swia l
l ta γ
C S N
X
Płaszczyzna Galaktyki b U C Równik galaktyczny
G’ G’
b = 90◦ − GX
(4.22)
l = CGX
−90◦ ≤ b ≤ 90◦
0 ≤ l ≤ 360◦
6
Chodzi oczywiście o kierunek ku centrum Galaktyki wzgl˛edem obserwatora znajdujacego
˛ si˛e w środku
sfery.
4.7 Współrz˛edne galaktyczne 41
B1950 J2000
αG = 12h 49m αG = 12h 51m 26.282s
(4.23)
δG = 27◦ 240 δG = 27◦ 070 42.0100
θ = 123◦ θ = 122.932◦
Epoki B1950, J2000 podano tu ze wzgl˛edu na precesyjny ruch bieguna świata. Współ-
rz˛edne galaktyczne obliczone za pomoca˛ wzorów podanych niżej b˛eda˛ wi˛ec również od-
niesione do układów współrz˛ednych z odpowiedniej epoki, co zawsze należy wyraźnie
zaznaczyć.
Wyprowadzimy teraz formuły transformacyjne pomi˛edzy współrz˛ednymi (α, δ) i (l, b).
Rozważmy trójkat˛ sferyczny GP X (rysunek 4.7b), w którym punkt X ma współrz˛edne
równikwe (α, δ) a którego współrz˛edne galaktyczne wynosza˛ (b, l). Możemy łatwo usta-
lić, że
HM X = HGX + λ (5.1)
Układy godzinny i równikowy maja˛ wiele ze soba˛ wspólnego. Oba zdefiniowane sa˛
w oparciu o biegun świata P , różnia˛ si˛e natomiast wyborem koła odniesienia, poczatku
˛
rachuby współrz˛ednej azymutalnej w tych układach. W obu układach koła te przecho-
dza˛ przez bieguny świata P Q. Dla kata ˛ godzinnego płaszczyzna˛ odniesienia jest po-
łudnik obserwatora, dla rektascensji jest nia˛ koło wielkie P Υ. Ponieważ punkt Υ jest
punktem należacym
˛ do sfery niebieskiej, stad˛ nie zmienia on swego położenia wzgl˛edem
gwiazd. Oznaczajac ˛ przez α rektascensj˛e gwiazdy X z rysunku 5.1, mamy, że wynosi
ona α = ΥP X = RA X . Z powodu ruchu wirowego ziemi a wraz z nia˛ z powodu ruchu
obserwatora, kat ˛ godzinny gwiazdy X zmienia si˛e w czasie. Rektascensja tej gwiazdy
pozostaje stała.
Punkt równonocy Υ stanowi ważne odniesienie wykorzystywanej w astronomii w
koncepcji czasu. W myśl niej, czas mierzony jest za pomoca˛ obserwacji ruchu dobo-
wego gwiazd a nie Słońca jak to ma miejsce w przypadku skali czasu słonecznego to-
warzyszacej
˛ nam w życiu codziennym. Taka koncepcja czasu, (albo jak mówimy skala
5.1 Czas gwiazdowy i rektascensja 43
λ
G M
X p
g m
rotacja
γ Ziemi C
rotacja sfery q
Rysunek 5.1: Sfera niebieska i umieszczona w niej sfera ziemska. Ilustracja układów
sferycznych równikowego i godzinnego oraz układu sferycznego na powierzchni Ziemi.
CG M = CG G + λ (5.4)
Z tego co powiedziano wynika, że czas gwiazdowy musi wzrastać o 24 godziny gwiaz-
dowe, czyli po interwale czasu, w którym kat ˛ godzinny punktu równonocy wzrośnie o
2
24h . Interwał ten nazywamy doba˛ gwiazdowa,˛ nie jest on równy dobie słonecznej, bo-
wiem doba gwiazdowa trwa 23h 56m w skali czasu słonecznego. Przyczyna tej różnicy
leży w tym, że punkt odniesienia gwiazdowej skali czasu, czyli punkt Υ jest nieruchomy
wzgl˛edem tła gwiazdowego. Tymczasem punkt odniesienia słonecznej skali czasu, czyli
Słońce, nieustannie przemieszcza si˛e na sferze wzgl˛edem tła gwiazdowego.
Czas gwiazdowy jest znakomitym łacznikiem
˛ pomi˛edzy katem
˛ godzinnym i rekta-
scensja.˛ Na rysunku 5.1, dla obserwatora w punkcie m, zgodnie z definicja˛ (5.2) czas
gwiazdowy równy jest katowi˛ sferycznemu M P Υ. Ponieważ kat ˛ godzinny punktu X
dla tego obserwatora wynosi M P X, a rektascensja punktu X równa jest α = RA X =
ΥP X, stad˛
CG M = HM X + RA X (5.5)
1
Punkt ten pełni tu rol˛e gwiazdy o zerowej rektascensji i deklinacji.
2
Jest to interwał bardzo bliski okresowi jednego obrotu Ziemi wokół jej osi rotacji.
44 Czas gwiazdowy i czas słoneczny
gdzie HM to kat ˛ kodzinny słońca. Stały składnik 12h wprowadzono by poczatek˛ doby
słoneczenj przypadał w nocy. Natomiast określenie "miejscowy"ma na celu podkreślenie
lokalnego charakteru kata ˛ godzinnego, co inaczej oznacza, że obserwatorzy znajdujacy
˛
si˛e na różnych długościach geograficznych obserwuja˛ inna˛ wartość prawdziwego czasu
słonecznego.
Aby powiazać
˛ czas słoneczny z gwiazdowym wystarczy zastosować równanie (5.5),
˛ Słońce jako punkt X, wówczas dla obserwatora w m b˛edzie
traktujac
P
b) K α
a) γ
R
V
λ S D
λ δ
ε
A B α T F
S γ
U
B
C
E
Pierwsze dwa prawa wykorzystane zostana˛ natychmiast. Rysunek 5.2a, ilustruje elips˛e
orbity Ziemii, S oznacza Słońce, odcinek AB oś wielka˛ elipsy. Punkt A, perihelium, to
punkt w którym w styczniu Ziemia znajduje si˛e najbliżej Słońca, punkt B, aphelium to
punkt najwi˛ekszego oddalenia Ziemi od Słońca. Długość półosi orbity ziemskiej wybrano
na jednostk˛e długości w astronomii, tzw. jednostka astronomiczna (AU ). Odległość ta
wynosi w jednostkach układu SI 1.49597870691 · 1011 m.
Niech punkt C (rysunek 5.2a) oznacza położenie Ziemi w momencie równonocy wio-
sennej, tzn. w chwili gdy obraz Słońca na tle gwiazd znajduje si˛e w punkcie Υ. Niech
punkt E przedstawia położenie Ziemi w jakiś czas później, obraz Słońca przemieści si˛e
wówczas w położenie R. Kat ˛ ΥSR jest zatem długościa˛ ekliptyczna˛ λ Słońca. Pr˛ed-
kość katowa
˛ Ziemi w ruchu orbitalnym nie jest stała, co łatwo wydedukować z 1 i 2
prawa Keplera, a to oznacza, że długość ekliptyczna Słońca nie zmienia si˛e jednostajnie
na przestrzeni roku. Najszybciej zmienia si˛e w styczniu, najwolniej w czerwcu w czasie
przejścia Ziemi przez aphelium. Niejednorodności w tempie zmiany długości Słońca z
oczywistych wzgl˛edów sa˛ przyczyna˛ zmian w tempie przyrostu jego rektascensji, a to
z kolei pociaga
˛ nierównomierność skali prawdziwego czasu słonecznego, co łatwo wy-
wnioskować z równania (5.7).
Inna˛ przyczyna nieregularności w przyrostach rektascensji Słońca jest nachylenie eklip-
tyki do równika. Niech sfera z rysunku 5.2b, przedstawia geocentryczna˛ sfer˛e niebieska˛
z ekliptyka˛ i równikiem — koło U AΥSV i koło ΥT F odpowiednio. Punkty V i U ozna-
czja˛ położenia na ekliptyce o najwi˛ekszej i najmniejszej deklinacji (δ = ±ε)). Sa˛ to tzw.
punkty przesilenia letniego i przesilenia zimowego.
Niech Słońce znajduje si˛e na ekliptyce w punkcie S(α , δ ). Koło wielkie P S prze-
46 Czas gwiazdowy i czas słoneczny
ΥS = λ ΥT = α T S = δ ΥT S = 90◦ SΥT = ε
Stosujac
˛ tego trójkata
˛ odpowiedni wzór cotangensowy otrzymamy
Wynika stad, ˛ że rektascensja Słońca nie zmienia si˛e jednostajnie z jego długościa:
˛ przy-
rosty rektascensji sa˛ najmniejsze w okresie równonocy, najwi˛eksze w czasie przesileń.
I właśnie z powodu wyraźnie nierównomiernych przyrostów rektascensji, prawdziwe
Słońce nie nadaje si˛e jako wzorzec skali czasu cywilnego, która powinna być skala˛ o
dużej regularności.
Do tego celu wykorzystuje si˛e obiekt zwany słońcem dynamicznym definiowany po-
gladowo
˛ ˛ sposób. Niech τ oznacza moment przejścia Ziemi przez peryhe-
w nast˛epujacy
lium, na rysunku 5.2a prawdziwe Słońce znajduje si˛e wówczas na sferze w punkcie B.
˛ a˛ pr˛edkość Ziemi na orbicie, czyli n = 360◦ /rok.
Dalej niech n oznacza średnia˛ katow
Wyobraźmy sobie fikcyjny obiekt poruszjacy ˛ a˛ n, w
˛ si˛e po ekliptyce z pr˛edkościa˛ katow
taki sposób, że przez punkt B przechodzi w tej samej chwili co słońce prawdziwe. 5 Ten
fikcyjny obiekt nazwano dynamicznym słońcem średnim.
Przypuśćmy, że w pewnej chwili t prawdziwe Słońce znajduje si˛e w S a słońce dyna-
miczne jest w D (rysunek 5.2b). Jeśli czas wyrażony jest w latach to położenie słońca
dynamicznego da si˛e wyznaczyć za pomoca˛ formuły
BD = n(t − τ )
Skale czasu słonecznego miejscowego używane sa˛ bardzo rzadko. W życiu codziennym
wymagana jest synchronizacja czasu na dostatecznie dużym obszarze. Z tego wzgl˛edu
powierzchnia kuli ziemskiej podzielona została na strefy czasowe oddzielone od siebie
tzw. południkami standardowymi. Wewnatrz˛ każdej strefy obowiazuje
˛ ten sam czas sło-
neczny średni, zwany czasem strefowym,
Przemiana współrz˛ednych
z
z1
P
K
y1
a
yk y
ipt
Elk
Row
nik
γ
x, x1
N
S O x
E
Horyzont
y
˛ 180◦ ,
[x1 , y1 , z1 ]T jest złożeniem obrotów układu horyzontalnego wokół osi osi z o kat
◦
a nast˛epnie wokół nowej y o kat ˛ −(90 − φ). Transformacja ma zatem postać
x1 x
y1 = q(φ − 90◦ )r(180◦ ) y (6.3)
z1 H,δ z A,h
z z1
Rownik
S x
γ
x1
y y1
G podane jest za pomoca˛ (α, δ) albo (H, δ). Deklinacja δ nie wymaga transformacji, na-
tomiast do współrz˛ednych azymutalnych α i H możemy stosować zwiazek ˛ (5.5). W ten
sposób problem transformacji współrz˛ednych z układu równikowego do układu godzin-
nego wyczerpuje si˛e.
Ponieważ porzadek ˛ musi być, stad
˛ nic dziwnego, że dla porzadku,
˛ podajemy macie-
T
rzowa˛ wersj˛e tej transformacji. Zatem, jeśli [x, y, z] b˛eda˛ składowymi wersora kierunku
do obiektu wzgl˛edem układu godzinnego, to składowe tego wersora wzgl˛edem układu
równikowego daja˛ si˛e obliczyć jako
x1 x
y1 = r(−S)My y (6.5)
z1 α,δ z H,δ
z
δG
z1 9 0−
P Rownik galaktyczny
G
θ υ
γ x
φ
O ψ
N
90
θ−
C
x1
2. Samolot startuje w Limie kierujac ˛ si˛e wprost na Rzym. Oblicz przebyta˛ odległość
w kilometrach, a także podaj długość geograficzna˛ samolotu, w momencie gdy
przelatywał nad równikiem. Współrz˛edne geograficzne Limy i Rzymu wynosza,˛
odpowiednio:
(12◦ 100 S, 77◦ 050 W ), (41◦ 530 N, 12◦ 330 E).
3. Oblicz długość najkrótszej drogi powietrznej z San Francisco (37o 400 N, 122◦ 250 W )
do Tokio (35◦ N, 139◦ 450 E). Wyznacz kierunek w jakim samolot powinien wystar-
tować w San Francisco oraz oblicz współrz˛edne geograficzne najbardziej północ-
nego punktu tej drogi.
11. Skale czasu gwiazdowego i słonecznego maja˛ różne jednostki. Czy zatem w rów-
naniu wia˛żacym
˛ obie te skale
57
Rozdział 7
1. układ odniesienia nie może obracać si˛e wzgl˛edem Wszechświata jako całości,
2. poczatek
˛ układu odniesienia nie może poruszać si˛e ruchem przyspieszonym.
W dalszej cz˛eści wykładu rozważymy w jaki sposób można tym warunkom zadośćuczy-
nić.
x − vt t − v/c2 x
x0 = p , y 0 = y, z 0 = z, t0 = p
1 − β2 1 − β2
ψ.t P
K
P1
ekliptyka
ψ.t X
.
ε rowniki
γ
ε γ1
Rysunek 7.1: Precesja luni-solarna powoduje zmian˛e położenia bieguna świata z miejsca
P do niejsca P1 , biegun ekliptyki K nieruchomy.
słonecznego. Natomiast najwi˛eksze pary sił od planet, od Jowisza i Wenus sa˛ o czynnik
105 słabsze i w przypadku osi obrotu Ziemi najcz˛eściej bywaja˛ pomijane. Wypadkowa
para sił skr˛ecajacych
˛ od Ksi˛eżyca i Słońca nie jest stała, zmienia si˛e wraz ze zmianami w
konfiguracji i wzajemnej odległości tych ciał. I właśnie dlatego ruch bieguna P na sferze
jest tak bardzo skomplikowany. Dla wygody rozdzielono go na dwie cz˛eści: cz˛eść uśred-
niona˛ na długim interwale czasu, inaczej cz˛eść wiekowa˛ zwana˛ precesja˛ luni-solarna,˛ oraz
na okresowe oscylacje wokół pozycji średniej zwane nutacja.˛ Na rysunku 7.1, ruch pre-
cesyjny bieguna wykreślono linia˛ przerywana˛ P P1 , natomiast linia falista reprezentuje
faktyczny ruch bieguna uwzgl˛edniajacy ˛ nutacj˛e.
Przemieszczenie nutacyjne bieguna jest rz˛edu 1500 i zostało odkryte w ubiegłym stu-
leciu przez Anglika Bradley’a, który poprawnie zinterpretował drobne okresowe zmiany
deklinacji gwiazd jakie zauważył podczas obserwacji południkowych.
Efekt precesji luni-solarnej jest wi˛ekszy od nutacyjnego, a co ważniejsze kumuluje si˛e
w miar˛e upływu czasu. Precesj˛e znali już starożytni Grecy. Dwa wieki przed narodzinami
Chrystusa Hiparchus z Rodos porównywał swoje obserwacje gwiazd z wykonanymi 150
lat wcześniej. Zauważył, że szerokości ekliptyczne gwiazd nie zmieniły si˛e podczas gdy
w ich długościach była wyrażna różnica, odpowiadajaca ˛ przyrostowi około 5000 rocznie.
Niech ψ oznacza roczne tempo precesji luni-solarnej. Zatem jeśli punkt P na rysunku
7.1 odpowiada położeniu bieguna świata w epoce poczatkowej, ˛ a P1 położeniu bieguna t
˛ sferyczny P KP1 = ψ · t. Dalej, skoro KP1 = KP = ε, to nachylenie
lat później, to kat
ekliptyki do równika nie nie uległo w tym czasie żadnym zmianom. Jeśli teraz (λ1 , β1 )
b˛eda˛ współrz˛ednymi gwiazdy X w epoce późniejszej, to z równań 7.1 i 7.2 mamy
P1 KX = 90◦ − λ1 , KX = 90◦ − β1
A ponieważ P1 KX = P KX − P KP1 , możemy napisać
λ1 = λ + ψ · t (7.3)
Skoro w omawianym zjawisku punkty K, X sa˛ nieruchome to odległość KX nie zmie-
niła si˛e, a wi˛ec nie zmieniła si˛e szerokość ekliptyczna gwiazdy. Odpowiednie zmiany we
współrz˛ednych równikowych powodowane precesja˛ luni-solarna˛ sa˛ bardziej skompliko-
wane i należałoby je wyprowadzić rozważajac ˛ sferyczne KP X oraz KP1 X.
˛ trójkaty
Północny biegun świata wskutek precesji luni-solarnej zakreśla wokół bieguna eklip-
tyki koło małe w czasie około 26000 lat. Oś rotacji bryły ziemskiej zmieniajac ˛ kierunek
w przestrzeni jest jednak ciagle
˛ jednakowo nachylona do płaszczyzny ekliptyki. Wskutek
tego zjawiska w miar˛e upływu lat współrz˛edne gwiazd moga˛ ulec drastycznej zmianie.
Powyższy opis precesji jest dość grubym przybliżeniem gdyż opiera si˛e na dwóch
nieścisłych założeniach. Mianowicie, że nachylenie ekliptyki do równika oraz tempo
precesji luni-solarnej sa˛ stałe. Ponadto dotad ˛ nie wzi˛eliśmy w rachub˛e ruchu punktu K
czyli bieguna ekliptyki. Zgodnie z dynamika˛ Newtonowska˛ Ziemia porusza si˛e wokół
Słońca po orbicie keplerowskiej w stałej płaszczyźnie. Tak definiowaliśmy płaszczyzn˛e
ekliptyki. Nie jest to jednak całkiem ścisły wniosek, gdyż wyprowadzony został z od-
działywań jedynie dwóch ciał, Słońca i Ziemi. Zupełnie pomini˛eto wpływ pozostałych
planet. Wpływ ten powoduje drobne perturbacje ziemskiej keplerowskiej orbity, w rezul-
tacie czego obserwujemy małe przesuni˛ecia punktu K na sferze (rysunek 7.2).
Przypuśćmy, że biegun ekliptyki uległ przesuni˛eciu z K do K1 . Jest to nieduża
zmiana około 0.100 rocznie. Zbadamy wpływ przesuni˛ecia KK1 na współrz˛edne równi-
7.3 Układ inercjalny a precesja, nutacja iruch własny gwiazd 63
. P α
α1
K
K1
X ekliptyki
C rownik
γ
ε γ1
E1
kowe gwiazdy X zakładajac,˛ że biegun świata P jest nieruchomy. Skoro P X = 90o −δ, to
przesuni˛ecie KK1 nie ma żadnego wpływu na deklinacj˛e gwiazdy. Kat ˛ KP X = 90o + α
(rysunek 7.2) został zredukowany o kat ˛ KP K1 . W konsekwencji o taki sam kat ˛ uległa
zmniejszeniu rektascensja gwiazdy, co nie zależy od położenia gwiazdy na sferze. Dla
wszystkich gwiazd, efekt przesuni˛ecia punktu K wskutek perturbacji planetarnych jest
taki sam: każdego roku rektascensje ulegaja˛ zmniejszeniu o wartość oznaczana˛ tradycyj-
nie przez λ0 zwana˛ precesja˛ planetarna.˛ Zmianom tym towarzyszy zmniejszenie kata ˛ ε,
nachylenia równika do ekliptyki.
wówczas stosować formuły ścisłe a te sa˛ bardziej złożone. Współczynniki, które w nich
wyst˛epuja˛ daja˛ si˛e jednak wyliczyć z pomoca˛ stałych precesji ψ, λ0 , p, ich aktualne war-
tości obliczane sa˛ z pomoca˛ fotmuł:
ψ = 50.387800 + 0.004900 · T
λ0 = 0.105500 − 0.018900 · T (7.5)
p = 50.291000 + 0.022200 · T
wzgl˛edem środka sfery, a także od jej odległości. A zatem, odległe najcz˛eściej słabe
gwiazdy b˛eda˛ miały niewielki ruch własny.
Przy założeniu, że ruchy własne gwiazd maja˛ kierunki przypadkowe, z południko-
wych obserwacji daje si˛e wydzielić systematyczne efekty precesyjne. Oznacza to, że
możemy wyznaczać zarówno ruchy własne jak i stałe precesji. Jednak dokładność wy-
znaczenia jednej wielkości ogranicza dokładność określenia drugiej. A wi˛ec mimo, iż
ruchy własne moga˛ być odniesione do układu inercjalnego, to w przypadku obserwacji
południkowych b˛edzie to możliwe jedynie z dokładnościa˛ z jaka˛ znane sa˛ stałe precesji.
Ruchy własne wyznaczane sa˛ z obserwacji wykonanych w odległych od siebie epo-
kach. Potrzeba bowiem sporo czasu by przemieszczenie gwiazdy narosło do mierzalnej
wielkości, co pozwoliłoby na wyznaczenie składowych ruchu własnego w rektascensji
i deklinacji z duża˛ dokładnościa.˛ Najpowszechniej stosowane metody polegaja˛ na po-
równaniu rezultatów precyzyjnych pomiarów klisz fotograficznych wykonanych w róż-
nych epokach. Obserwacje fotograficzne gwiazd dostarczaja˛ jedynie wzgl˛ednych położeń
obiektów znajdujacych˛ si˛e na kliszy. Ich położenia wzgl˛edne moga˛ być wyznaczone z
duża˛ precyzja,˛ ale chcac˛ znać wartości absolutne współrz˛ednych (α, δ), trzeba wykorzy-
stać niektóre z gwiazd na kliszy jako gwiazdy odniesienia. Oznacza to, że ich współrz˛e-
dne a także ruchy własne sa˛ znane z góry. A zatem mimo iż wzgl˛edne położenia gwiazd
znane sa˛ bardzo dokładnie, dokładniej niż z obserwacji południkowych, przewaga ta w
dużym stopniu znika z powodu konieczności oparcia si˛e o absolutne pomiary południ-
kowe. I dlatego ruchy własne gwiazd otrzymane z klisz fotograficznych nie sa˛ wolne od
bł˛edów systematycznych zastosowanego układu odniesienia. Ta˛ ostatnia˛ trudność można
zminimalizować wybierajac ˛ jako gwiazdy odniesienia obiekty tak odległe, że ich ruchy
własne można uznać za zaniedbywalne.
Ruchy własne zwykle wyznacza si˛e wzgl˛edem heliocentrycznej sfery niebieskiej.
Zgodnie z definicja˛ ruch własny jest efektem niezerowej tangencjalnej składowej wek-
tora pr˛edkości gwiazdy wzgl˛edem Słońca. Ponieważ Galaktyka jako całość obraca si˛e,
zarówno gwiazda jak i Słońce poruszaja˛ si˛e własnymi ruchami wzgl˛edem środka Galak-
tyki. Słońce jak si˛e uważa znajduje si˛e ok. 10 kpc od centrum Galaktyki i podobnie
do gwiazd ze swego najbliższego sasiedztwa
˛ porusza si˛e po z grubsza kołowej orbicie z
szybkościa˛ liniowa˛ 200−250 km/s. Pełnego obiegu wokół centrum dokonuje wi˛ec w cza-
sie około 0.25 miliarda lat. Informacje te pochodza˛ z badań nad kinematyka˛ Galaktyki, w
których wykorzystano pomiary ruchów własnych i pr˛edkości radialnych gwiazd.
Średnia pr˛edkość wszystkich gwiazd w najbliższym sasiedztwie
˛ Słońca definiuje tzw.
lokalny standard spoczynku (LSR). Każda gwiazda ma pewna˛ pr˛edkość wzgl˛edem LSR,
a wi˛ec i Słońce. Obserwowany ruch własny gwiazdy zależy od jej pr˛edkości wzgl˛edem
Słońca i dlatego wpływ pr˛edkości własnej Słońca (tzw. lokalny ruch słoneczny) tkwi w
obserwowanych ruchach własnych gwiazd. Ale zakładajac, ˛ że poszczególne ruchy własne
gwiazd maja˛ rozkład losowy, możliwym jest na drodze statystycznych analiz wyznaczyć
lokalna˛ składowa˛ słoneczna,˛ jej wielkość i kierunek.
Również LSR ma pewna˛ pr˛edkość wzgl˛edem środka Galaktyki. Jest to główna pr˛ed-
kość rotacyjna Glaktyki w kierunku l = 90◦ , b = 0. Rotacyjna pr˛edkość galaktyczna
zmienia si˛e jednak wraz z odległościa˛ od środka, co wywiara mały systematyczny wpływ
na ruchy własne gwiazd. Wi˛eksza cz˛eść gwiazd jakie możemy obserwować, znajduje si˛e
w tym samym obszarze Galaktyki co Słońce. Sa˛ to gwiazdy odległe od Słońca nie wi˛ecej
niż kilka kiloparseków. Powinniśmy zatem opisywać ich ruch jako różnicowa˛ rotacj˛e ga-
66 Astronomiczne układy odniesienia
laktyczna˛ aniżeli rotacj˛e jako całości. Można pokazać, że różnicowa rotacja galaktyczna
wywołuje ruch własny gwiazdy w płaszczyźnie Galaktyki dany formuła˛
µ ∝ A · cos(2 · l) + B (7.6)
gdzie l, jest długościa˛ galaktyczna˛ gwiazdy. Taki ruch własny nie zależy od odległości
i powoduje wzrost długości galaktycznej. Stałe A, B nazywane sa˛ stałymi Oorta. Nie
znamy ich zbyt dokładnie, w szczególności słabo znamy stała˛ B. Obie stałe wynosza˛
około 0.01 sekund łuku na rok. Co jest tu bardzo istotne to to, że efekt opisany rów-
naniem (7.6) oznacza, że każda gwiazda na kliszy wykazuje pewien ruch własny. Stad ˛
statystyczne założenie, że gwiazdy słabe, znajdujace
˛ si˛e daleko od Słońca, maja˛ znikomy
ruch własny nie jest uzasadnione. Dlatego, bez dodatkowych poprawek, słabe gwiazdy
nie nadaja˛ si˛e do definicji dobrego układu odniesienia.
do zera z odległościa.˛ A przynajmniej dla tych obiektów, które możemy jeszcze ob-
serwować należy oczekiwać, że b˛eda˛ one bardzo małe. Oszacowanie przypadkowych
pr˛edkości transwersalnych dokonane zostało z pomoca˛ pr˛edkości własnej Słońca wzgl˛e-
dem tła Wszechświata. Z badań tła mikrofalowego wynika, że pr˛edkość ta wynosi około
400 km/s, co odpowiada zaniedbywalnemu ruchowi własnemu rz˛edu 10−4 sekundy łuku
w odległości jednego megaparseka. Dlatego układ odniesienia zdefiniowany z pomoca˛
położeń obiektów pozagalaktycznych spełnia wszystkie wymagania inercjalnego układu
odniesienia.
7.4 Poczatek
˛ układu odniesienia — środek sfery niebies-
kiej
Najcz˛eściej stosowane bywaja˛ trzy poczatki˛ układów odniesienia: miejsce obserwacji,
środek Ziemi i środek Słońca co odpowiada sferze topocentrycznej, sferze geocentrycznej
i sferze heliocentrycznej. Obserwacje z konieczności dokonywane sa˛ wzgl˛edem sfery to-
pocentrycznej, natomiast ze wzgl˛edu na kompromis dogodny dla wszystkich astronomów
zamieszkujacych
˛ kul˛e ziemska,˛ efemerydy obserwacyjne planet podawane sa˛ w odniesie-
niu do środka Ziemi. Sfera heliocentryczna wygodna jest do opisu świata gwiazd.
Rozważmy najpierw topocentryczny poczatek ˛ układu odniesienia. Jak wiadomo, uczest-
niczy on w ruchu dobowym Ziemi. Jest on zmienny co do kierunku, czyli przyspieszony,
a zatem topocentrum nie może stanowić poczatku ˛ inercjalnego układu odniesienia. Ruch
dobowy obserwatora przejawia si˛e w dwojaki sposób: po pierwsze wprowadza zmienna˛
składowa˛ do pr˛edkości radialnych, po drugie powoduje okresowe zmiany obserwowanych
z Ziemi położeń obiektów wskutek zjawisk paralaksy i aberracji. Wpływ aberracyjny oraz
zmiana pr˛edkości radialnej zależa˛ od stosunku szybkości obserwatora do szybkości świa-
tła. Maksymalna˛ pr˛edkość ma obserwator na równiku ziemskim i wynosi ona 0.465 km/s.
Paralaksa natomiast, zależy od długości wektora przemieszczenia obserwatora z jednego
poczatku
˛ do innego. W omawianym przypadku przemieszczenie jest wielkościa˛ rz˛edu
promienia Ziemi, 6378 km.
Podobne rozważania dotycza˛ geocentrycznego poczatku ˛ układu odniesienia. Z po-
wodu ruchu orbitalnego Ziemi środek geocentrycznego układu także doznaje przyspie-
szeń i dlatego nie można traktować go jako poczatku
˛ układu inercjalnego. Ruch przyspie-
szony poczatku˛ układu powoduje roczne zmiany w pr˛edkościach radialnych, paralaks˛e
roczna˛ i roczna˛ aberracj˛e. 6 Orbitalna pr˛edkość Ziemi wynosi w przybliżeniu 30 km/s
co stanowi 104 pr˛edkości światła. Przemieszczenie poczatku ˛ układu odpowiedzialne za
paralaks˛e roczna˛ jest rz˛edu jednostki astronomicznej. Zatem, układ odniesienia o po-
czatku
˛ w środku Ziemi nie realizuje układu inercjalnego, może jednak być podstawa˛ do
wyznaczania odległości do gwiazd.
Sytuacja zmienia si˛e gdy poczatek˛ układu umieścimy w środku Słońca. W trakcie
ruchu dookoła centrum Galaktyki Słońce porusza si˛e zasadniczo ze stała˛ pr˛edkościa.˛ Jego
przyspieszenie a, można wyznaczyć nast˛epujaco. ˛ Przyjmijmy promień kołowej orbity
Słońca r = 104 pc,7 , okres obiegu T = 2.5 · 108 lat. Wówczas w jednostkach SI, pr˛edkość
6
Przyspieszenie liniowe w ruchu rocznym wynosi około 6 · 10−3 ms−2 .
7
Dla wygody podajemy, że 1pc = 3.1 · 1016 m
68 Astronomiczne układy odniesienia
katowa
˛ Słońca wynosi
2π
ω= = 8.0 · 10−16 s−1
T
a przyspieszenie liniowe
Tak małe przyspieszenie jest niewykrywalne przez współczesne pomiary pozycyjne i ra-
dialne, a w czasie jednego stulecia powoduje zmian˛e w pr˛edkości mniejsza˛ niż 1 m/s. Stad
˛
w praktyce, heliocentryczny układ odniesienia można traktować jako układ inercjalny.
Ale z jednym zastrzeżeniem. To cały Układ Słoneczny znajduje si˛e w niemal jednostaj-
nym ruchu wzgl˛edem środka Galaktyki a nie Słońce z osobna. Dlatego punktem, który
należałoby adoptować jako poczatek
˛ inercjalnego układu odniesienia jest barycentrum
Układu Słonecznego. Słońce przemieszcza si˛e wokół tego punktu w zmiennej odległości
rz˛edu 106 km. Należy wi˛ec odróżniać sfer˛e heliocentryczna˛ od barycentrycznej, gdyż
tylko ta ostatnia z nieruchomym równikiem i punktem równonocy w pełni realizuje układ
inercjalny.
7.5.1 Paralaksa
Niech C oznacza jakiś poczatek
˛ standardowy, O poczatek
˛ w miejscu obserwatora. Wek-
tor położenia punktu O wzgl˛edem C oznaczmy przez R (rysunek 7.3). Położenie obiektu
wzgl˛edem C oznaczmy przaz r. Widomy, czyli obserwowany kierunek do obiektu zmiej-
sca O dany jest wektorem r0
r0 = r − R (7.7)
r=rs r0 = r0 s0 R = R so
7.5 Przesuni˛ecie paralaktyczne i aberracyjne 69
E
r’
S
p
r 0 s0 = r s − R so . (7.8)
lub
r R
s0 = 0
s − 0 so
r r
Jeśli pomnożymy je dwukrotnie, lewostronnie, wektorowo przez wersor s, to
0 r R
s × s × s = s × s × 0 s − 0 so
r r
Wykorzystujac
˛ znane zwiazki
˛ wektorowe (patrz [?]), otrzymamy
R
(s · s0 )s − (s · s)s0 = − s × (s × so )
r0
˛ R r, a wi˛ec dla małych przesuni˛eć paralaktycznych możemy podstawić
Zakładajac
r = r oraz s · s0 = 1, i w rezultacie dostajemy wzór przybliżony
0
R
ds = s0 − s = s × (s × so ) (7.9)
r
Równanie to bardzo przypomina wyprowadzona˛ w rozdziale drugim formuł˛e 3.14 na
małe przesuni˛ecie na sferze.
Na rysunku 7.3, punkt S i O reprezentuja˛ położenia obiektu i obserwatora odpowia-
dajce pewnemu momentowi czasu t. Sa˛ to tzw. położenia (miejsca) geometryczne, w
których obiekt i obserwator znajduja˛ si˛e w chwili t. Podobnie mówimy o kierunkach na
punkty S i O, sa˛ to kierunki geometryczne wzgl˛edem C.
GT (t −τ )
∆ θI
E(t)
∆ θII GA(t −τ )
G(t)
B(t)
II- ∆θII to efekt ruchu obiektu, zmiana kierunku do obiektu jest konsekwencja˛ nieze-
rowego czasu propagacji promieniowania od obiektu do obserwatora.
Gdy obserwowany obiekt jest gwiazda˛ w poprawce ∆θ = ∆θI nie uwzgl˛ednia si˛e skład-
nika od czasu propagacji promieniowania. Dlatego poprawka ta nosi miano aberracji
gwiazdowej rocznej, dobowej, . . . zależnie od tego czy w rachub˛e bierzemy ruch roczny,
dobowy, . . . . Dla obiektów położonych w Układzie Słonecznym uwzgl˛ednianie sa˛ obie
przyczyny a poprawka ∆θ = ∆θI + ∆θII nazywana jest aberracja˛ planetarna.˛
Wobec takiej konwencji, katalogowe położenia gwiazd jako uwolnione od aberra-
cji gwiazdowej, można porównać z położeniami obiektów Układu Słonecznego dodajac ˛
do efemerydy tych obiektów jedynie poprawk˛e ∆θII . Nie ma potrzeby wprowadzania
poprawki z tytułu ruchu obserwatora, bo efemeryd˛e porównuje si˛e z niewielkim polem
gwiazdowym, dla którego poprawka ∆θI dla gwiazd i planet jest niemal taka sama. Ten
dziwny kierunek, powstały porzez dodanie do miejsca geometrycznego planety jedynie
poprawki za czas propagacji nosi miano kierunku (miejsca) astrometrycznego. Jest to
kierunek do poruszajacego
˛ si˛e obiektu jaki wyznaczyłby hipotetyczny nieruchomy obser-
wator.
Wyprowadzimy teraz wzór pozwalajacy ˛ na wyznaczenie poprawki ∆θI wynikajacej ˛
z ruchu obserwatora wzgl˛edem obiektu. Rozważamy przypadek, w którym na rysunku
7.4 B b˛edzie punktem nieruchomym, natomiast punkt E ma pr˛edkość V wzgl˛edem B.
Przyjmijmy, że obserwacj˛e kierunku do pewnego obiektu wykonano w E w momencie
t. Oznaczmy przez τ czas propagacji kwantu promieniowania od obiektu do obserwa-
tora. Zatem na rysunku 7.4 punkt E reprezentuje położenie obserwatora w momencie t,
a punkt GA położenie obiektu w momencie wcześniejszym (t − τ ). 8 Linia EGA repre-
zentuje trajektori˛e fotonu ale taka˛ jaka˛ wyznaczonoby w układzie odniesienia o poczatku
˛
8
Wprowadzenie tych założeń nie narusza ścisłości równania (7.7) dla konfiguracji punktów B, E, GA .
7.6 Dodatek A. Zadania 71
7.5.3 Łaczny
˛ wpływ paralaksy i aberracji
Całkowity wpływ przemieszczenia poczatku ˛ z punktu B do ruchomego punktu E otrzy-
mamy dodajac ˛ do siebie równania (7.9) i (7.12). Gwarantuje to jedynie dokładność pierw-
szego rz˛edu w Rr i Vc , co dla wielu zastosowań w zupełności wystarcza. Przy tej dokład-
ności, sA może być zastapione
˛ po prawej stronie równania (7.12) przez s, i w rezultacie
możemy napisać
R V
sT − s = s × (s × sR ) − s × (s × n) (7.13)
r c
Podsumowujac, ˛ sT jest kierunkiem obserwowanym (topocentrycznym) w E w momencie
t, s jest astrometrycznym kierunkiem obiektu wzgl˛edem B w chwili (t − τ ). Przemiesz-
czenie poczatku˛ układu z B do E wynosi RsR , a pr˛edkość nowego poczatku ˛ wzgl˛edem
B jest równa V n.
Jeśli wymagana jest duża dokładność, rezygnujemy z formuł przybliżonych i para-
laksa musi być uwzgl˛edniona z pomoca˛ dokładnej formuły (7.7). Chcac ˛ podwyższyć
dokładność opisu przesuni˛ecia aberracyjnego należy zastosować aparat szczególnej teorii
wzgl˛edności.
Rozdział 8
Refrakcja
mikrofale
a)
00
11
350 km 11
00
00
11
00
11
00
11
00
11
00
11
00
11
fale krotkie 00
11
310 Warstwa2F
11
00
00
11
00
11
00
11
00
11
Zorze 00
11
00
11
220 00Warstwa 1F
11
Zorze
00
11
11
00
00
11
00
11
00
11 Zorze
00
11
00
11
00
11
00
11
00
11 111
000
00 125
11 000 Warstwa F
111
00
11
00
11
00
11
00 110
11
0000
1111 111
000
0000
1111
0000
1111 111 Warstwa D
000
Meteory1111
0000 90 fale srednie
55
36 0000Warstwa ozonowa
1111
28
12 Pow. Ziemi 1111 Troposfera
0000
b)
km Satelita
700
600
Jonosfera
500
400
300
20 MHz
100
10 MHz
z N =z
nN
nN−1
zi
ni
n i−1
z0 z i−1
n1
n0
Powierzchnia Ziemi
R = z − z0 (8.3)
Ponieważ zmmiany kierunku propagacji promieni świetlnych nie sa˛ duże, we wzorze po-
wyżej możemy zastosować przybliżenie małych katów,
˛ stad ˛ refrakcji R w radianach
˛ kat
dany jest jako
R = K tan z0 (8.5)
gdzie
W ramach płaskiego modelu atmosfery równanie (8.2) jest dokładne, natomiast równania
(8.4) i (8.5) sa˛ przybliżeniami pierwszego rz˛edu ze wzgl˛edu na (n0 − 1) i daja˛ dobre
rezultaty dla niedużych wartości odległości zenitalnych.2
Dla dużych odległości zenitalnych nie ma jednak sensu wprowadzanie wyższych rz˛e-
dów przybliżenia. A to dlatego, że dla dużych odległości od zenitu, dla których wy-
razy wyższych rz˛edów sa˛ znaczace,
˛ bardziej istotnym jest zmodyfikownie równania (8.5)
przez właczenie
˛ do modelu efektu krzywizny atmosfery.
W formułach (8.4), (8.5) kat ˛ refrakcji zależy tylko od obserwowanej odległości z0
oraz n0 współczynnika refrakcji przyziemnej warstwy atmosfery. Zupełnie nie intere-
suje nas stan wyższych warstw atmosfery, gdyż do obliczenia wartości współczynnika
n0 wystarczajace˛ sa˛ informacje o stanie atmosfery bezpośrednio otaczajacej
˛ obserwatora.
Właściwość ta stanowi spora˛ zalet˛e modelu płaskiej atmosfery, która znika natychmiast
po wprowadzeniu efektów krzywizny.
Współczynnik n0 zależy od lokalnych warunków atmosferycznych. Jako standardowe
przyj˛eto warunki odpowiadajace ˛ ciśnieniu 760 mm Hg i temperaturze 0o C, dla tych da-
nych współczynnik załamania wynosi
n0 = 1.0002927 (8.7)
co pociaga
˛
K = 60.4” (8.8)
Taka wartość stałej K nazywana jest stała˛ refrakcji. Dla warunków niestandardowych
K obliczana jest w oparciu o prawo Dale-Gladstone stwierdzajace, ˛ że (n0 − 1) jest pro-
porcjonalne do g˛estości powietrza. Oznacza to, że jeśli P jest ciśnieniem atmosfery w
milimetrach słupa rt˛eci a T temperatura˛ w stopniach Celsjusza, to z prawa opisujacego
˛
własności gazu wynika
P
n0 − 1 ∝ (8.9)
273 + T
Równanie (8.9) zastosowane łacznie
˛ z wartościami dla warunków standardowych daje
wzór na kat
˛ refrakcji
P/760
R = 60.4” · tan z0 (8.10)
1 + T /273
Udokładnienie równania (8.10) bez właczenia
˛ efektu zakrzywienia atmosfery mija si˛e z
celem.
Trzeba jednak jeszcze wspomnieć o dwóch ważnych kwestiach. Współczynnik zała-
mania zależy nie tylko od lokalnych wartości ciśnienia i temperatury, ale także od składu
chemicznego powietrza oraz od długości fali padajacego
˛ promieniowania elektromagne-
tycznego. Skład powietrza (patrz tabela 8.1) można tu uważać za stały poza drobnymi
zmianami ilości pary wodnej. Jest to jednak tak mały efekt, że zaniedbanie go w przybli-
żonej formule (8.10) jest zupełnie uzasadnione.
Poważniejszy problem wynika z zależności współczynnika załamiania n0 od długo-
ści fali. Standardowa wartość n0 = 1.0002927 odpowiada środkowi pasma wizualnego
2
Cz˛esto napotykana˛ w podr˛ecznikach wartościa˛ graniczna˛ stosowalności formuły (8.5) jest z0 = 45◦ .
78 Refrakcja
α0 = CGM
δ0 = φ
gdzie CGM jest miescowym czasem gwiazdowym, φ jest szerokościa˛ geograficzna˛ ob-
serwatora.
8.4 Refrakcja — model atmosfery radialnie symetrycznej 79
Z
90−φ
P t z
90−δ X’
X
.
Horyzont
dθ = −R = −K tan z0 ≈ −K tan z
k = −K sec z
K sec2 δ sin t
dα = α0 − α =
cos t + tan φ tan δ
tan φ − tan δ cos t
dδ = δ 0 − δ = K (8.13)
cos t + tan φ tan δ
Z’ Q
Z
θ
S
z0 +R ψ
ψS
A P
h0 z0
r rS
O
r0 θ
z =θ+ψ (8.14)
˛ ψ jest katem
Ponieważ kat ˛ pomi˛edzy kierunkiem radialnym punktu P i styczna˛ do krzy-
wej θ = θ(r), stad
˛
dθ
tan ψ = r (8.15)
dr
Na koniec skonstruujmy asymptot˛e trajektorii promieniowania, przecina ona półprosta˛
OZ w punkcie A. Oznaczmy wysokość punktu A nad miejscem obserwacji przez h0 .
Współczynnik refrakcji n = n(r) jest funkcja˛ odległosci od środka Ziemi, możemy
wi˛ec wyobrazić sobie atmosfer˛e jako złożona˛ ze skończonej liczby odpowiednio cienkich
3
˛ pomi˛edzy styczna˛ w P do trajektorii promieniowania i kierunkiem za zenit.
Kat
8.4 Refrakcja — model atmosfery radialnie symetrycznej 81
n i+1
ni
R ψi P ψi+1
χ
ri
C ri+1
do
Stad
˛
Oznacza to, że iloczyn r n sin ψ, dla każdej warstwy jest niezmiennikiem, w szczegól-
ności dla warstwy przyziemnej mamy
Równanie (8.17) daje nam elementarna˛ refrakcj˛e przy przejściu promienia światła przez
powierzchni˛e rozdziału dwóch sasiednich
˛ nieskończenie cieńkich kulistych warstw at-
mosfery. Całkowita˛ refrakcj˛e R przy przejściu światła przez cała˛ atmosfer˛e ziemska˛ uzy-
skamy całkujac
˛ równanie (8.17).
Zamiast odległości r, za zmienna˛ niezależna˛ wygodniej jest przyjać ˛ współczynnik
załamania n. Jeżeli tan ψ wyrazimy z pomoca˛ równania (8.16) jako
r0 n0 sin z0
tan ψ = (8.18)
(r2 n2 − r02 n20 sin2 z0 )1/2
wówczas, kładac
˛ (8.18) do równania (8.17), po scałkowaniu, otrzymamy kat
˛ pełnej re-
frakcji
Z z Z n0
dn
R= dz = r0 n0 sin z0 (8.19)
z0 1 n(r n − r0 n0 sin2 z0 )1//2
2 2 2 2
Jest to dokładna formuła zwana całka˛ refrakcji, wyst˛epuja˛ w niej dwie zmienne r i n. Po-
wiedziano wcześniej, że współczynnik załamania zmniejsza si˛e z odległościa˛ r od środka
Ziemi. A zatem gdybyśmy dysponowali funkcja˛ r = r(n), równanie (8.19) można by
˛ refrakcji R.
scałkować numerycznie, tzn. obliczylibyśmy dokładna˛ wartość kata
Obok kata˛ refrakcji można jeszcze wyznaczyć równanie trajektorii wiazki˛ promie-
niowania. W tym celu, podstawiajac ˛ prawa˛ stron˛e równania (8.18) do równania (8.15)
otrzymamy całk˛e
Z r0
dr
θ = r0 n0 sin z0 2 2 2 2 2 1//2
(8.20)
r0 r(r n − r0 n0 sin z0 )
Znajac
˛ równanie trajektorii możemy wyznaczyć jej asymptot˛e a w dalszej kolejności wy-
sokość h0 (patrz rysunek 8.4). Jednakże rozważajac ˛ ACS z tego rysunku, można
˛ trójkat
otrzymać bardzo prosty wzór na h0 . Potrzeba tylko założyć, że punkt S odpowiada ta-
kiemu miejscu na trajektorii promienia świetlnego, dla którego n = 1.
Niech rs i ψs b˛eda˛ odpowiadały takiemu właśnie punktowi. Ponieważ ZAS = z0 + R,
stad ˛ ACS oraz ze wzoru sinusów mamy
˛ z trójkata
Stosujac
˛ prawo refrakcji (8.16) do lewej strony tego równania otrzymamy
n0 sin z0
h0 = r0 −1 (8.21)
sin(z0 + R)
Dla niewielkich odległości zenitalnych wysokość h0 jest zaniedbywalnie mała. Nawet
dla umiarkowanych odległości zenitalnych wynosi zaledwia kilka metrów. Wzrasta jed-
nak szybko w pobliżu horyzontu i wynosi 1.5 km dla gwiazd znajdujacych˛ si˛e tuż przy
horyzoncie.
Sens fizyczny wysokości h0 jest nast˛epujacy:
˛ po uwzgl˛ednieniu refrakcji, topocen-
tryczne współrz˛edne ciała dotycza˛ nie tyle miejsca gdzie znajduje si˛e obserwator, ale
miejsca położonego o h0 powyżej obserwatora. Rozróżnienie to może być ważne podczas
wyznaczania poprawki z tytułu paralaksy geocentrycznej. Ze wzgl˛edu na niewielka˛ war-
tość h0 odpowiedniej korekty trzeba dokonać jedynie w przypadku obserwacji ksi˛eżyca,
8.4 Refrakcja — model atmosfery radialnie symetrycznej 83
r = r0 + h (8.22)
gdzie
Z n0
dn
R1 = n0 sin z0 (8.24)
1 n(n2 − n20 sin2 z0 )1/2
oraz
Z n0
n0 sin z0 hndn
R2 = (8.25)
r0 1 (n2 − n20 sin2 z0 )3/2
Aby rozwini˛ecie to było przydatne w praktyce, wyrażenie (n2 − n20 sin2 z0 )3/2 nie może
być wielkościa˛ mała,˛ tego samego rz˛edu co h/r0 . Pociaga
˛ to by z0 było wyraźnie mniejsze
aniżeli arcsin(1/n0 ), a co z kolei sprowadza si˛e do warunku by źródło promieniowania
nie znajdowało si˛e zbyt blisko horyzontu. Całka R1 daje si˛e scałkować dokładnie,
n0
n0 sin z0
R1 = − arcsin = arcsin(n0 sin z0 ) − z0
n 1
sin(z0 + R1 ) = n0 sin z0
84 Refrakcja
co dokładnie odpowiada równaniu (8.2) w modelu płaskiej atmosfery. Dla wygody warto
˛ R1 w szereg wzgl˛edem (n0 − 1), ale tym razem należy zachować wyrazy kwa-
rozwinać
dratowe (n0 − 1)2 , które poprzednio, w modelu płaskiej atmosfery nie zostały wzi˛ete w
rachub˛e. Rozwini˛ecie ma postać
Pierwszy wyraz po prawej stronie znika gdy ρ = ρ0 bo wówczas h = 0, da˛ży też do zera
gdy ρ → 0.
Zdefiniujmy teraz H0 jako tzw. wysokość równoważnej jednorodnej atmosfery4
1 ∞
Z
H0 = ρdh (8.28)
ρ0 0
H0
R2 = (n0 − 1) tan z0 sec2 z0 (8.29)
r0
4
Jest to wysokość atmosfery o tej samej masie masie co atmosfera rzeczywista, ale której g˛estość byłaby
stala i równa g˛estości atmosfery na poziomie morza, ρ = 1.293 · 10−3 [g · cm−3 ]. H0 ≈ 8 km.
8.5 Stałe refrakcji. Tablice refrakcji 85
X1 X
Z
φ
90− z0
P
90−δ
z’0
Y1
Łacz
˛ ac˛ równania (8.26) i (8.29) dostaniemy wyrażenie na całkowita˛ refrakcj˛e w formie
gdzie
φ − δ = z0 + A tan z0 + B tan3 z0
Przyczyny powodujace ˛ zmiany stałej refrakcji sa˛ bardzo różnorodne. Mówiliśmy o za-
leżności n0 od ciśnienia atmosferycznego, temperatury, wilgotności oraz długości fali.
Wartość skali wysokości H0 zależy od lokalnej grawitacji, a ponieważ ta głównie zależy
od odległości obserwatora od środka Ziemi, potrzebna jest wi˛ec poprawka w szerokości i
wysokości nad poziomem morza dla danego miejsca obserwacji.
Z powodu skomplikowanego charakteru zjawiska refrakcjij, od wielu lat zamiast for-
muł uwzgl˛edniajacych
˛ wszystkie te wpływy, wykorzystuje si˛e specjalne tablice refrakcji.
Kilka dużych obserwatoriów zestawiło tzw. tablice refrakcji uwzgl˛edniajace ˛ wspomnia-
ne wyżej efekty, np. Pułkowskie Tablice Refrakcji. Ale trzeba zdawać sobie spraw˛e, że
nawet najbardziej rozległe tablice pozwalaja˛ wyznaczyć refrakcj˛e jako funkcj˛e jedynie
miejscowych warunków atmosferycznych. A tymczasem ogólny wzór na refrakcj˛e (8.19)
pokazuje, że całkowita refrakcja zależy nie tylko od warunków lokalnych (przyziemnych)
ale także od zmian współczynnika załamania z wysokościa.˛ Równanie (8.30), w którym
współczynniki daja˛ si˛e określić na podstawie lokalnych warunków atmosferycznych, sta-
nowi tylko dwa pierwsze wyrazy rozwini˛ecia równania dokładnego. Daje ono zadowa-
˛ rezultaty jedynie dla odległości zenitalnych mniejszych niż 75o . Wyższe wyrazy
lajace
rozwini˛ecia zależa˛ od szczegółowej struktury atmosfery i dla z > 80o , jest to zależność
krytyczna. Dla normalnych odległości zenitalnych tablice refrakcji opieraja˛ si˛e zarówno
na teorii jak i obserwacjach. W pobliżu horyzontu tablice te opieraja˛ si˛e niemal wyłacz-˛
nie na obserwacjach. W pobliżu horyzontu kat ˛ refrakcji jest tak duży i zmienny, że w
takich przypadkkacch tablice refrakcji moga˛ dać jedynie wartości przybliżone. Wyklucza
to jakiekolwiek precyzyjne pozycyjne obserwacje obiektów znajdujacych ˛ si˛e w pobliżu
horyzontu.
Z tego co powiedziano o zjawisku refrakcji jest chyba oczywiste, że refrakcja zmienia
współrz˛edne gwiazd w sposób nie do końca dajacy ˛ si˛e modelować precyzyjnie, zwłaszcza
gdy gwiazda znajduje si˛e na dużych odległościach zenitalnych. Wynika stad, ˛ że w celu
precyzyjnego określenia położeń gwiazd, trzeba ograniczyć si˛e do niewielkich odledłości
zenitalnych. Istnieja˛ instrumenty zaprojektowane specjalnie do tego typu obserwacji.
Problem atmosferycznej refrakcji może być jednak rozwiazany ˛ w całości. W tym
celu należy przenieść optyczne teleskopy w przestrzeń kosmiczna.˛ W ostatnim wierszu
tabeli 8.2 widzimy, że nie musielibyśmy wówczas wprowadzać refrakcyjnych poprawek.
Dla teleskopów optycznych umieszczonych na okołoziemskich orbitach, pozycja obser-
wowana (pomijamy tu bł˛edy instrumentalne) b˛edzie od razu pozycja˛ topocentryczna.˛
8.6 Dodatek A. Zadania 87
K sec2 z
f=
1−K
89
Rozdział 9
Współrz˛edne geocentryczne
Linia pionu
Powierzchnie poziome
P Powierzchnia ladow
H
Poziom morza
N
Geoida
Elipsoida
Dno oceanu
x2 y2
+ =1 (9.1)
a2 a2 (1 − f )2
gdzie
b = a(1 − f ) (9.2)
Zenit geodezyjny
000000000000
00000000
111111111111
11111111 Zenit geocentryczny
000000000000
00000000
y 111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
ν
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000 00
11
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111 00O
11
000000000000
00000000
111111111111
11111111 OO’=h
11
00000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000 CO=ρ
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
11O’
ρ 00
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
b 111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
00000000
111111111111
11111111
000000000000
φ’ 00000000
111111111111
11111111
φ
000000000000
00000000
111111111111
11111111
000000001 0 x
000000000000
C 111111111111
11111111 0
1
Q a
a = 6378.140 [km]
f = 0.00335281 = 1/298.257
ρ cos φ0 = x0 = x + h cos φ
(9.3)
ρ sin φ0 = y0 = y + h sin φ
Ponieważ punkt (x, y) leży na elipsie o równaniu (9.1) a tan φ jest nachyleniem normalnej
do elipsy w tym punkcie, pociaga ˛ to
−dx
tan φ =
dy
Różniczkujac
˛ równanie (9.1) dostaniemy
Kładac
˛ jego prawa˛ stron˛e spowrotem do (9.1) b˛edziemy mieli
x2 (1 + (1 − f )2 tan2 φ) = a2
Z pomoca˛ tego równania oraz równania (9.4) można wyrazić współrz˛edne x i y poprzez
φ. Mianowicie,
x = aC cos φ
(9.5)
y = aS sin φ
gdzie
A zatem uwzgl˛edniajac
˛ (9.5), kkońcowa˛ postacia˛ równań (9.3) jest
Z’
O z’
r’ S
p
ρ
r
z
a zwłaszcza w przypadku sztucznych satelitów Ziemi, paralaksa osiaga ˛ bardzo duże war-
tości.
Na rysunku 9.3, punkt O oznacza obserwatora, C środek Ziemi a S pewne pobliskie
ciało niebieskie. Linia CO, jej przedłużenie określa kierunek na geocentryczny zenit Z 0
miejsca obserwacji. Płaszczyzna rysunku jest zdefiniowana przez trzy punkty C, O i S,
a zatem leży w płaszczyźnie koła wertykalnego przechodzacego ˛ przez gwiazd˛e. Dlatego
przekrój Ziemi pokazany na rysunku niekoniecznie musi przebiegać wzdłuż ziemskiego
południka.
Oznaczmy kat ˛ Z 0 OS przez z 0 . Jest to obserwowana odległość zenitalna odniesiona do
geocentrycznego zenitu. Kierunki zenitów geodezyjnego i astronomicznego, ogólnie nie
musza˛ leżeć w płaszczyźnie rysunku.
Niech r0 i r b˛eda˛ topocentryczna˛ i geocentryczna˛ odległościa˛ źródła promieniowania,
ρ odległościa˛ obserwatora od środka Ziemi (rysunek 9.3). Paralaksa˛ geocentryczna˛ p,
nazywamy kat ˛ OSC taki, że
z0 = z + p (9.8)
gdzie z, jest geocentryczna˛ odległościa˛ zenitalna˛ obiektu, jaka˛ obserwowanoby ze środka
Ziemi. Paralaksa geocentryczna zwi˛eksza geocentryczna˛ odległość zenitalna˛ o kat ˛ pa
skoro zmiana ta odbywa si˛e w płaszczyźnie OCS (w płaszczyźnie wertykału), azymut
geocentryczny pozostaje niezmieniony.
Stosujac ˛ OCS wzór sinusów dostaniemy
˛ do trójkata
ρ ρ
sin p = sin z 0 = 0 sin z (9.9)
r r
Wynika stad,˛ że paralaksa dla danego obiektu poza zależnościa˛ od r, zależy od jego od-
ległości zenitalnej a także od odlegości ρ obserwatora od środka Ziemi. Potrzebna jest
zatem pewna standaryzacja i jest nia˛ tzw. horyzontalna paralaksa równikowa P . Jest to
paralaksa fikcyjnego obiektu, położonego na horyzoncie obserwatora znajdujacego ˛ si˛e na
równiku ziemskim (ρ = a), czyli w warunkach gdy (z 0 = 90◦ ). Po podstawieniu do
równania (9.9), mamy
a
sin P = (9.10)
r
94 Współrz˛edne geocentryczne
Paralaksy geocentryczne planet maja˛ wyraźnie mniejsze wartości. Dla Saturna wynosi
ona około 100 , a dla najbliższej planety, dla Wenus waha si˛e w granicach 500 − 3400 . Dla
obiektu znajdujacego
˛ si˛e w odległości 1 AU , jego paralaksa nosi nazw˛e paralaksy sło-
necznej. Poza drobnymi różnicami jest ona bardzo bliska średniej paralaksie prawdzi-
wego Słońca.
11
00
00 Eros
11
M1
11
00
W1
11
00
1Z 1
0 00
11
11
00
00
11 γ
S
11
00
00
11
1
0 W2
M2 11
00
00
11
Z2
Rysunek 9.4: Rzuty orbit planet: Venus, Ziemi, Marsa oraz planetoidy Eros na płaszczyź-
nie ekliptyki.
pierwsza˛ planeta,˛ której odległość wyznaczono metoda˛ radiowa,˛ udało si˛e tego dokonać
podczas konjunkcji w roku 1959. Te i kolejne pomiary w nast˛epnych latach doprawa-
dziły do nowego określenia paralaksy słonecznej. W adoptowanym przez MUA systemie
stałych mamy, że
Jednostki astronomicznej nie definiuje si˛e już jako długości półosi wielkiej orbity Ziemi,
ponieważ półoś ta zmienia si˛e z powodu perturbacji planetarnych. Obecnie definiuje si˛e ja˛
za pomoca˛ teorii grawitacji. Stała k z równania (9.13) znana jako stała grawitacji Gaussa,
w systemie stałych zalecanych przez MUA jej wartość wynosi
k = 0.01720209895 (9.16)
Wartość ta przetrwała ostatnie zmiany systemu stałych i nic nie wskazuje by miało być
inaczej w najbliższej przyszłości. Przy takiej stałej k, jednostk˛e astronomicza˛ (AU ) defi-
niuje si˛e jako jednostk˛e długości w jakiej musi być wyrażona półoś ap z równania (9.13),
gdy okres Tp podany jest w dobach.
Wartość paralaksy słonecznej wyznaczona za pomoca˛ metod radarowych określona
jest z dokładnościa˛ do jednej mikrosekundy łuku. Jest to poza zasi˛egiem współczesnych
metod obserwacji pozycyjnych i dlatego kolejne zbliżenie Erosa do Ziemi wykorzystano
jedynie do badań jego własności topograficznych. Nawiasem mówiac, ˛ dokonano tego
technikami radarowymi.
Ze wzgl˛edu na bliska˛ odległość, paralaks˛e Ksi˛eżyca daje si˛e określić dokładniej ani-
żeli paralaksy planet. Ale i tutaj techniki radarowe wyeliminowały klasyczne metody
optyczne, te zaś po pewnym czasie, zastapiono˛ nowoczesna˛ technika˛ optyczna˛ — tech-
nika˛ laserowa.˛ Umożliwiły to misje ksi˛eżycowe Apollo, kiedy to astronauci na poczatku ˛
lat 1970-tych umieścili na Ksi˛eżycu odbłyśniki laserowe.
Jeśli chodzi o wyznaczanie odległości paralaksa geocentryczna ma obecnie mniej-
sze znaczenie, ale nadal jest istotna jako efekt pozycyjny podczas określania dokładnych
położeń ciał niebieskich.
• k = ρr ,
• α0 − α = t, kat
˛ godzinny obiektu.
Przy takich oznaczeniach, dla obiektu o (α, δ), formuły na małe przesuni˛ecie maja˛ postać
Sa˛ to formuły przybliżone (pierwszy rzad ˛ ze wzgl˛edu na (ρ/r)), dlatego nie należy ich
używać w przypadku Ksi˛eżyca praz sztucznych satelitów Ziemi. Nadaja˛ si˛e dla pozo-
stałych ciał niebieskich o ile ciała te nie znajduja˛ si˛e w zbyt bliskim sasiedztwie
˛ Ziemi,
kiedy to paralaksy geocentryczne przekraczaja˛ wartości kilku sekund łuku.
Dla obiektów bardzo dalekich jak planety zewn˛etrzne, paralaksy geocentryczne sa˛
bardzo małe i dlatego w formułach (9.18) nie ma potrzeby rozróżnienia pomi˛edzy sze-
rokościami geodezyjna˛ i geocentryczna,˛ jako ρ można do wzorów podstawić wartość a
równikowego promienia Ziemi.
Formuły przybliżone należy stosować z rozwaga,˛ a w sytuacjach watpliwych
˛ opłaca
si˛e stosowanie rozwiazań
˛ dokładnych. Mianowicie, niech r i R b˛eda˛ geocentrycznymi
wektorami położenia obserwowanego obiektu S i obserwatora O. Wówczas wektor OS ~ =
r0 od obserwatora do źródła promieniowania dany jest jako różnica
r0 = r − R (9.19)
Rysunek 9.3 odpowiada właśnie tej sytuacji. Załóżmy, że wektor R jest znany dokładnie,
jest to wektor o długości ρ skierowany na geocentryczny zenit obserwatora, a zatem we
współrz˛ednych równikowych ma on składowe
r = a csc P
r = 7205.843 [km]
α = 8h 47m 13s rownania (9.22), (9.23))
δ = 28◦ 150 3800
gdzie (α0 , δ0 ) sa˛ równikowymi współrz˛ednymi kierunku jaki wskazuje wektor jednost-
kowy n. Formuły (9.25) sa˛ ogólnymi wyrażeniami na aberracyjne zmiany rektascensji i
deklinacji.
Rozpatrzmy teraz przypadek obserwatora znajdujacego ˛ si˛e na pewnej szerokości i od-
ległości geocentrycznej φ0 , ρ. Jeśli ω oznacza katow
˛ a˛ szybkość wirowania Ziemi, liniowa
szybkość obserwatora wzgl˛edem środka Ziemi dana jest wzorem
V = ρω cos φ0 (9.26)
ρ = 6378.140 [km]
φ = φ0
2π
ω = doba gwiazdowa = 7.292 · 10−5 s−1
mamy, że w jednostkach praktycznych poprawki na aberracj˛e dobowa˛ wyrażaja˛ si˛e wzo-
rami
3. Korzystajac
˛ z rezultatów poprzedniego zadania, oblicz maksymalna˛ geodezyjna˛
wysokość jaka˛ może osiagn
˛ ać
˛ w tym miejscu obserwacji satelita, poruszajacy
˛ si˛e
po kołowej orbicie o promieniu 8798 km, nachylonej pod katem
˛ 18◦ 360 do równika.
e2 sin 2φ
tan ν =
2(1 − e2 sin2 φ)
Rozdział 10
10.1 Wst˛ep
Omówimy sposoby wyznaczania współrz˛ednych równikowych α, δ ciał niebieskich za
pomoca˛ klasycznego instrumentu astrometrycznego — koła południkowego. Nie podamy
jego pełnej teorii, ograniczymy si˛e jedynie do przedstawienia zasady obserwacji i kilku
podstawowych poprawek instrumentalnych tego narz˛edzia.
Bł˛edy instrumentalne każdego teleskopu astronomicznego dziela˛ si˛e na:
• bł˛edy pochodzace
˛ z niedoskonałości samego instrumentu,
• bł˛edy wynikajace
˛ z niedoskonałości montażu.
α0 = T + ∆T (10.1)
104 Pomiary rektascensji i deklinacji
a) Z b)
P
B
G
1
0
φ
E δ
φ h Z W
N S
W δ
Z
a) b)
P
90
a
−n
90−b
m
90−φ
P Z
90−m
90+a
N W’ S
−b
90
Horyzo
90
−n
nt
W
n ik
Row W’
Rysunek 10.2: Bł˛edy ustawienia koła południkowego: oś instrumentu przebija sfer˛e w
punkcie W’ zamiast W; położenie punktu W’ wzgl˛edem W określone jest para˛ małych
katów
˛ (a, b) lub (m, n).
a = W ZW 0 b = 90◦ − ZW 0 (10.4)
P W 0 = 90◦ − n
W 0 Z = 90◦ − b
P Z = 90◦ − φ
W 0 P Z = 90◦ − m
W 0 ZP = 90◦ + a
to za pomoca˛ cztero-elementowej formuły cotangensowej otrzymamy
P
a) b)
m P
τ
90−δ
Z
−n
90−m+τ
90
n−
90−δ
90
X
’
90−δ
W 90+c
E
W’ ik
rown X
’
90−δ 90 c
t
S horyzon W’ O
Rysunek 10.3: Z powodu niedokładnego ustawienia koła południkowego jak i bł˛edu ko-
limacji lunety, gwiazda X góruje wzgl˛edem południka instrumentalnego o interwał τ za
wcześnie.
co redukuje si˛e do
a przy założeniu, że poprawki m, n, a, b sa˛ małe, w równaniach (10.6), (10.7) można
skorzystać z przybliżenia małych katów
˛ i wówczas
m = a sin φ + b cos φ
(10.8)
n = b sin φ − a cos φ
Stałe instrumentalne m, n, c zwykle wyrażone sa˛ w mierze czasowej bowiem potrzebne sa˛
przy redukcji obserwacji współrz˛ednej rektascensji. Jedynie stała˛ poprawk˛e deklinacyjna˛
d podaje si˛e w sekundach łuku.
τ = α − α0 (10.9)
10.3 Usuwanie wpływów instrumentalnych w kole południkowym 107
P W 0 X = 90◦ − δ 0 (10.11)
Ponieważ, jak si˛e za chwil˛e przekonamy, δ różni si˛e od δ 0 jedynie wyrazami drugiego
rz˛edu, w praktyce jest wi˛ec oboj˛etne, która z tych deklinacji zostanie podstawiona do
wzoru (10.10) na poprawk˛e τ .
Ale podczas wyznaczenia wartości samego δ może być koniecznym wprowadzenie
wyrazów drugiego rz˛edu. Dlatego wyznaczymy te wyrazy, i w tym celu w równaniu
(10.12) połóżymy δ = δ 0 + ∆,
sin(δ 0 + ∆) = sin δ 0 cos ∆ + cos δ 0 sin ∆ = − sin n sin c + cos n cos c sin δ 0
rozwijajac
˛ funkcje trygonometryczne z katami
˛ n i c w szeregi pot˛egowe, biorac
˛ jedynie
po dwa pierwsze wyrazy
gdzie τ ponownie oznacza czas prawdziwego przejścia minus czas przejścia rejestrowa-
nego, ale α0 musi być teraz rozumiana jako miejscowy czas gwiazdowy plus 12 godzin.
Teoria zawarta w równaniach (10.10) do (10.14) pozwala na uwolnienie obserwowa-
nych wartości rektascensji i deklinacji z głównych bł˛edów instrumentalnych, oczywiście
pod warunkiem, że stałe instrumentalne sa˛ wcześniej znane. Wyznaczenie niektórych z
tych stałych wymaga obserwacji gwiazd. Ponieważ bardzo pożadanym ˛ jest ciagłe
˛ ka-
librowanie instrumentu astrometrycznego, dlatego wyznaczanie tych stałych, zaleca si˛e
wprowadzić jako integralna˛ cz˛eść programu obserwacji gwiazd.
Jak widać, wpływ w rektascensji jest podobny do poprawki kolimacji instrumentu w for-
mule Bessel’a — oba efekty sa˛ proporcjonalne do sec δ. Można zatem połaczyć
˛ obie
poprawki w jedna˛ poprzez podstawienie
α = α0 + m + n tan δ 0 + c∗ sec δ 0
(10.17)
δ = δ0 − R
1958.0
y −0.2’’
E 1959.5 1963.0
1960.5
180+Γ −0.1’’
λ=0 1961.5
1962.0
1959.0
P +y
1958.5 1961.0
γ 1960.0
0.1’’
1962.5
0.2’’
Γ λ0=0
N x
P0 S
0.3’’
W 0.4’’ 0.3’’ 0.2’’ 0.1’’ +x −0.1’’−0.2’’
bieguny świata na sferze niebieskiej. Bieguny świata przemieszczaja˛ si˛e na sferze wsku-
tek precesji luni-solarnej i nutacji. Obecnie tego typu ruchem biegunów si˛e nie interesu-
jemy. W dalszych rozważaniach b˛edziemy traktować oś rotacji Ziemi jako posiadajac ˛ a˛
niezmienne położenie wzgl˛edem gwiazd. Interesuje nas pozorny ruch tej osi wzgl˛edem
powierzchni Ziemi.
Wobec tego co powiedziano wyżej, skoro ruch biegunów po powierzchni Ziemi nie
wpływa na położenie biegunów świata, to nie przyczynia si˛e do zmian rektascensji i dekli-
nacji gwiazd. Wpływa jednak na proces redukcji obserwacji tych współrz˛ednych wykona-
nych z pomoca˛ instrumentów przejściowych. Jest tak gdyż montaż teleskopu ustawiony
jest wzgl˛edem powierzchni Ziemi, a to oznacza, że na skutek ruchu skorupy ziemkiej
cały zbiór gwiazd jest przesuwany wzgl˛edem instrumentu. W efekcie z powodu ruchu
biegunów zmieniaja˛ si˛e stałe instrumentalne narz˛edzia.
Na rysunku 10.4a, naszkicowano obszary polarne widziane ”z lotu ptaka”. Punkt
P oznacza chwilowy biegun rotacji, P0 jest jego średnim położeniem nazywanym nie-
zbyt trafnie biegunem figury ziemskiej. W przeciwieństwie do chwilowego, biegun figury
jest stałym punktem na powierzchni Ziemi. Przemieszczenie γ punktu P wzgl˛edem P0
wynosi około 0.00 3, co na powierzchni Ziemi odpowiada odległości bliskiej 10 m (patrz
rysunek 10.4b). Wschodnia długość Γ tego przesuni˛ecia najcz˛eściej zwi˛eksza si˛e, a za-
tem biegun chwilowy obiega P0 w kierunku przeciwnym do wskazówek zegara. Jest to
ruch skomplikowany, nie dajacy ˛ si˛e opisać dokładnie, ale jego główne składowe sa˛ znane.
Sa˛ to okres 428 dniowy i okres roczny. Przyczyny ruchów bieguna sa˛ niemal w całości
natury geofizycznej. Gdyby Ziemia była swobodnie wirujac ˛ a˛ bryła˛ sztywna,˛ odchylenie
bieguna rotacji od osi symetrii Ziemi spowodowałoby jednostajny ruch kołowy bieguna
chwilowego wokół bieguna figury z okresem 305 dniowym. Ze wzgl˛edu na plastycz-
ność Ziemi okres ten uległ silnemu wydłużeniu (okres 428 dniowy jest zmodyfikowanym
okresem swobodnego ruchu wirowego Ziemi). Składowa roczna ruchu biegunów, jest
składowa˛ wymuszona˛ narzucona˛ przez ulegajace ˛ ciagłej
˛ zmianie warunki geofizyczne.
Przemieszczenie P wzgl˛edem P0 wyrażone jest zwykle za pomoca˛ współrz˛ednych
prostokatnych
˛ (x, y), przy czym oś x skierowana jest wzdłuż południka zerowego, oś y
10.5 Ruch biegunów 111
φ0 Z λ
a) dφ 90− b)
Q φ
90− Z λ0
P0
P
0
da
−φ
90
90−
N0 W’ S0
φ
N
W0 W
Γ)
18
λ0
(λ−
0
−Γ
+Γ
λ
0−
18
P0 Γ γ Γ
λ0 P
Do Greenwich
rownolegle do siebie
x = γ cos Γ
(10.18)
y = −γ sin Γ
cos γ cos (180◦ − (λ − Γ)) = sin γ cot (90◦ − φ) − sin (λ − Γ) cot (λ0 − Γ)
a korzystajac
˛ ze znanych tożsamości trygonometrycznych mamy
Jeżeli w równaniach (10.19) i (10.20) zastosujemy wzory na sum˛e, różnic˛e sinusa i cosi-
nusa po czym wykorzystamy równanie (10.18), wówczas dostaniemy wyrażenia na chwi-
lowa˛ długość i szerokość geograficzna˛ w postaci
φ = φ0 + x cos λ0 − y sin λ0
(10.21)
λ = λ0 + (x sin λ0 + y cos λ0 ) tan φ
Możemy teraz pokazać jak ruchy bieguna wpływaja˛ na poprawki instrumentalne koła po-
łudnikowego. Punkt W 0 (rysunek 10.5a), zachodni koniec osi poziomej instrumentu z
powodu ruchu bieguna nie zmienia swego położenia na sferze. Podobnie nie zmieni si˛e
stała kolimacji narz˛edzia. A ponieważ założyliśmy, że zenit obserwatora jest nieruchomy,
wartość stałej wysokości instrumentu również pozostanie taka sama. Przemieszczenie
bieguna zmniejszy jednak wszystkie zachodnie azymuty, w tym stała˛ azymutalna˛ narz˛e-
dzia o kat˛ N0 N , a wi˛ec
da = −P0 ZP
a po zróżniczkowaniu, z wystarczajac
˛ a˛ dokładnościa˛ b˛edzie
Jest to zmiana momentu czasu obserwowacji przejścia gwiazdy przez południk spowo-
dowana ruchami biegunów ziemskich. Zmiana szerokości obserwatora, powoduje odpo-
wiednia˛ zmian˛e w odczycie deklinacji z koła podziałowego.
φ Z
90−
P H
z
90
−δ
φ
δ= X
N S
t
W horyzon
ik
rown
Podczas obserwacji w pierwszym wertykale bł˛edy instrumentalne sa˛ naturalnie takie same
jak dla koła południkowego, ale poprawki jakie trzeba wprowadzić do obserwacji sa˛
znacznie bardziej skomplikowane. Jedyna˛ przewaga˛ koła wertykalnego jest możliwość
niezależnych pomiarów szerokości obserwatora.
Zenit
S’
Ekran z>30
S
Promienie od gwiazdy S
30
S’
S’ z=30
120 S
S
z<30
S’
Horyzont rteciowy
φ Z
90− 30
90−δ
P X’ X
H
N S
t
horyzon
ficzne (umieszczona tuż pod soczewka), ˛ oraz z ciekłego lustra rt˛eciowego. Obraz gwiazdy
po odbiciu od lustra rt˛eciowego rejestrowany jest na kliszy fotograficznej.
Ważna˛ cecha˛ tego narz˛edzia jest to, że na emulsji fotograficznej leży punkt w˛ezłowy
soczewki objektywu, wówczas dobierajac ˛ odpowiednio odległość powierzchni rt˛eci otrzy-
mujemu na kliszy zogniskowane obrazy gwiazd. Takie rozwiazanie ˛ doskonale eliminuje
bł˛edy wysokości i kolimacji narz˛edzia.
Podczas ekspozycji oprawa z klisza˛ jest z odpowiednia˛ szybkościa˛ przesuwana pro-
stopadle do osi optycznej teleskopu. Przy czasach naświetlania kliszy od 10-20 sekund,
technika ta pozwala na uzyskanie punktowych obrazów gwiazd a wi˛ec na fotografowa-
nie gwiazd nawet stosunkowo slabych. W czasie obserwacji rejestrowany jest automa-
tycznie moment czasu odpowiadajacy ˛ środkowi interwału eksponowania kliszy. Jed-
nakże kompletna obserwacja za pomoca˛ FTZ wymaga czterech ekspozycji. Po każdej
ekspozycji, soczewka wraz z klisza˛ sa˛ obracane o 180◦ i w rezultacie cztery obrazy da-
nej gwiazdy tworza˛ na kliszy równoległobok, przykładowo pokazany na rysunku 10.9b.
Niech t1 , t2 , t3 , t4 b˛eda˛ momentami czterech obserwacji danej gwiazdy. Gdyby kamery
FTZ w czasie obserwacji nie obracano o 180◦ , pomijajac ˛ niewielkie zakrzywienie, cztery
obrazy leżałyby na liniach prostych. Ze wzgl˛edu na obroty, jjedynie odcinki X1 X3 i
X2 X4 odzwierciedlaja˛ ruch dobowy sfery. Przy obrocie kliszy dokładnie o kat ˛ 180◦ , od-
cinki sa˛ do siebie równoległe. W czasie obracania FTZ, punkt zenitu na kliszy pozostaje
nieruchomy. Jest on jednakowo odległy od odcinków X1 X3 i X2 X4 .
Określmy w dowolnym mmiejscu kliszy układ współrz˛ednych (x, y) , o osi OX rów-
noległej do tych odcinków. Niech (xi , yi ) sa˛ współrz˛ednymi punktów Xi , i=1,2,3,4.
Można je oczywiście zmierzyć za pomoca˛ precyzyjnych płytomierzy. Znajac ˛ skal˛e kliszy
(co ma miejsce, bowiem przykładowo, znamy odległość X1 X3 na kliszy oraz interwał
(t3 − t1 ) odpowiadajacy ˛ ruchowi dobowemu od X1 do X3 ), możemy współrz˛edne w jed-
nostkach liniowych łatwo zamienić w katowe.˛ Zauważmy też, że y1 = y3 , y2 = y4 .
Niech (x0 , y0 ) b˛eda˛ współrz˛ednymi obrazu zenitu miejsca obserwacji. Współrz˛edna
x0 jest nieznana ale y0 możemy wyznaczyć z formuły
a)
y
b)
Klisza
x3 x1
x2 x4
x
x3 x1
c)
φ−δ
Lustro rteciowe
Z
H0 = 0.25(x4 − x1 + x3 − x2 ) (10.33)
α = t0 − H0 (10.34)
Rozdział 11
X
1
0
r
r’
C00
11 00Z
11
R
r = r0 + R (11.1)
Jest to dokładne równanie, i jeśli X oznacza ciało z Układu Słonecznego taka˛ postać
równania należy zachować w zastosowaniach praktycznych. Gdy mamy do czynienia z
gwiazdami można dokonać pewnych uproszczeń.
˛ pomi˛edzy wektorami r i r0 nazywany jest paralaksa˛ roczna˛ p. W interwale jednego
Kat
roku kat ˛ CZX,
˛ ten zmienia si˛e wraz ze zmianami położenia Ziemi na orbicie. Dla trójkata
z twierdzenia sinusów wynika
R
sin p = sin Z (11.2)
r
gdzie Z elongacja gwiazdy (kat
˛ CZX). Ze wzgl˛edu na konieczność standaryzacji para-
˛ π1
laksa roczna gwiazdy definiowana jest jako kat
1
sin π = (11.3)
r
1
Oznaczenie to wprowadza nieco zamieszania i dlatego należy zachować ostrożność aby nie pomylić
paralaksy rocznej z oznaczeniem liczby niewymiernej 3.1415... .
11.1 Paralaksa roczna 121
Pr˛edkości gwiazd ujawniaja˛ si˛e przecież w ich ruchach własnych, a te sa˛ czymś bardzo
powszechnym.
A zatem poprawienie obserwowanego położenia gwiazdy na paralaks˛e i aberracj˛e me-
toda˛ dopiero co opisana,˛ nie daje geometrycznej barycentrycznej pozycji na moment ob-
serwacji powiedzmy t0 . Opisana tutaj redukcja daje położenie jakie gwiazda zajmowała
o interwał τ wcześniej. Czas τ jest czasem propagacji światła pomi˛edzy gwiazda˛ i ob-
serwatorem. W celu otrzymania geometrycznego położenia na moment t0 , musimy do
obliczonej podana˛ wyżej metoda˛ pozycji dodać rezultat iloczynu
τ · ruch wasny
Poprawka ta nazywana jest aberracja˛ wiekowa.˛ Ze wzgl˛edu na duże niepewności w po-
miarach odległości gwiazd, tym samym i duże niepewności w τ w praktyce poprawka
ta nie jest brana pod uwag˛e. Stad, ˛ przyczynek od aberracji wiekowej tkwi w tym co
rozumiemy pod poj˛eciem barycentryczne położenie gwiazdy.
Poprawka czasu propagacji światła nie może być ignorowana w przypadku ciał z
Układu Słonecznego. Dla tych ciał mówimy o tzw. aberracji planetarnej, rozumiejac ˛
przez to łaczny
˛ efekt zmiany geocentrycznego położenia, powodowany zarówno, nieze-
rowa˛ pr˛edkościa˛ obserwatora jak i źródła. Mamy wówczas pełna˛ poprawk˛e od miejsca
widomego do geometrycznego. Poprawka za aberracj˛e roczna˛ uwgl˛ednia jedynie wpływ
ruchu środka Ziemi wzgl˛edem barycentrum Układu Słonecznego.
β
λ
G’
y 0
1
11
00 x
G
Ṙ = V0 + V1 (11.15)
gdzie κ jest bezwymiarowym stosunkiem V0 /c, jest to tzw. stała aberracji. Zgodnie z
teoria˛ ruchu orbitalnego Ziemi wynosi ona
1/2
k 1+m
κ= (11.18)
c a(1 − e)2
11.4 Aberracja planetarna 125
gdzie k — stała Gaussa, m masa Ziemi, a oraz e — półoś wielka i mimośród orbity Ziemi.
Stała κ nie jest stała˛ absolutna˛ bowiem mimośród e wykazuje drobne wiekowe zmiany.
W systemie stałych MUA z roku 1976, wartościa˛ κ na epok˛e J2000.0 jest
κ = 20.00 49552 (11.19)
Analiza formuł (11.17) ukazuje, że gwiazda wokół swej pozycji heliocentrycznej, opisuje
na sferze elips˛e aberracyjna˛ o półosi wielkiej równoległej do płaszczyzny ekliptyki. Półoś
wielka elipsy równa jest κ, mimośród elipsy aberacyjnej wynosi κ sin β.
Przesuni˛ecie aberracyjne od pr˛edkości V1 otrzymać można podobnie za pomoca˛ rów-
nań (11.10), w rezultacie zmiany współrz˛ednych ekliptycznych maja˛ postać
dλ = κe sec β cos(ω̃ − λ)
(11.20)
dβ = κe sin β sin(ω̃ − λ)
Wartości te znane sa˛ jako tzw. człony E aberracji rocznej. Sa˛ one niezależne od długo-
ści Słońca i stad
˛ nie wykazuja˛ zmian rocznych. Składowe przesuni˛eć danych wzorami
(11.20) maja˛ amplitud˛e κe = 0.00 343, powoduja˛ one przesuni˛ecie całej elipsy aberracyj-
nej o stała˛ wielkość. Człony E, nie sa˛ jednak stałymi absolutnymi dla danej pary (λ, β),
wykazuja˛ bowiem pewne zmiany wiekowe.
1P’(t− τ)
0 .
0
1
111
000 τr
111
000
111111
000000 111
000
111111
000000
P(t) 00
111111
000000 11 111
000
111111
000000
111111
000000 111
000
111111
000000
111111
000000 111
000
111111
000000
111
111111
000000 000
111111
000000
111
111111
000000 000
111111
000000
111
111111
000000 000
111111
000000
111
111111
000000 000
111111
000000
111
111111
000000 000
111111
000000
111
000
111111
000000 111111
000000
111
000
111111
000000 111111
000000
111
000
111111
000000 111111
000000
111
000
111111
000000 111111
000000
111
000
111111
111111
000000 000000
111
000
111111 ρ’s’
111111
000000 000000
111
000
111111
111111
000000 000000
111
000
111111
111111
000000 000000
111
000
111111
000000
111111
000000 ρs 111
000
111111
000000
111111
000000
111111 111
000
111111
000000
000000
111111 111
000
111111
000000
000000
111111
000000 111
000
111111
000000
111111
000000 111
000
111111
000000
111111
000000
r 111
000
111111
000000
111111
000000 111
000
111111
000000
111
s*
111111
000000 000
111111
000000
111
111111
000000 000
111111
000000
111
000
111111
000000 111111
000000
111
000
111111
111111
000000 000000
111
000
111111
111111
000000 000000
111
000
111111
111111
000000 000000
111
000
111111
000000
111111
000000 111
000
111111
000000
111111
000000
111111 111
000
111111
000000
000000
111111 111
000
111111
000000
000000
111111 111
000
111111
000000
000000
111111 111
000
111111
000000
000000
111111
000000 111
000
111111
000000
111111
000000 111
000
111111
000000 1Z(t) .
0
111111
000000 R 111
000
111111
000000
111111
000000
00
11 R
111111
000000
00
11
G(t)
muje on aberracj˛e roczna˛ powodowana˛ przez pr˛edkość orbitalna˛ Ziemi, oraz poprawk˛e
powstała˛ w wyniku ruchu np. planety w czasie interwału τ , od emisji z powierzchni pla-
nety kwantu promieniowania do jego rejestracji na powierzchni Ziemi. Wpływ aberracji
rocznej dany jest równaniem (11.9), wyrażenie na poprawk˛e czasu propagacji podane
jest poniżej. Zakładamy, że obserwacji planety dokonano w momencie t. Niech na ry-
sunku 11.3 punkty G, Z, P b˛eda˛ położeniami barycentrum, Ziemi i planety odpowiednio,
wszystkie położenia w momencie t.
Niech r i R b˛eda˛ barycentrycznymi wektorami położeń planety i Ziemi, oznaczmy
jeszcze wektor ZP~ przez ρs, |s| = 1. Wektor ρs, daje geometryczny kierunek do planety
w momencie t. Oczywiście
r = ρs + R (11.21)
gdzie s0 zastapiono
˛ po prawej stronie przez s.
2 0
s otrzymano z s∗ po dodaniu doń poprawki jedynie na aberracj˛e roczna.˛
11.4 Aberracja planetarna 127
ρ0 s0 = ρs − τ ṙ (11.23)
Jak widać zmiana położenia planety w efekcie aberracji planetarnej zależy tylko od wzgl˛ed-
nej pr˛edkości Ziemi i planety.
Różniczkujac˛ równanie (11.21) mamy
ṙ − Ṙ = ρ̇s + ρṡ
s∗ = s − τ ṡ (11.26)
dα
α∗ = α − τ
dt
∗ dδ
δ =δ−τ (11.27)
dt
128 Współrz˛edne helio- i barycentryczne
s = s∗ + τ ṡ (11.28)
Główna trudność polega na tym, że czas propagacji nie jest znany a priori. Stad,
˛ zanim
b˛edzie można obliczyć poprawki z tytułu aberracji planetarnej, trzeba dokonać pewnych
oszacowań. W tym celu dokonujemy jeszcze jednego uproszczenia, mianowicie bierzemy
współrz˛edne widome takimi jakimi sa˛ i antydatujemy moment obserwacji o pewien in-
terwał τ , i dalej, z co najmniej trzech obserwacji wyznaczamy orbit˛e planety, po czym
obliczamy nowa˛ lepsza˛ wartość τ . Porces powtarzamy aż do uzyskania zbieżności roz-
wiazania
˛ na τ .
2. Udowodnij, że istnieja˛ dwa punkty na sferze niebieskiej, dla których efekt rocz-
nej aberracji znika. Pokaż, że ich współrz˛edne rownikowe dane sa˛ w przybliżeniu
wzorami
α = − arctan(cos ε cot λ )
δ = ± arcsin(sin ε cos λ
Rozdział 12
Ruch Ziemi
Streszczenie. Ziemia porusza si˛e, mi˛edzy innymi ruchem orbitalnym i obrotowym. Ruch
orbitalny niekiedy traktuje si˛e jako złożenie ruchu pary Ziemia- Ksi˛eżyc wokół ich ba-
rycentrum oraz ruchu barycentrum Ziemi i Ksi˛eżyca wokół Słońca. Płaszczyzn˛e prze-
chodzac ˛ a˛ przez środek Słońca i obracajac
˛ a˛ si˛e z w tempie wiekowej składowej szybkości
rotacji orbity barycentrum Ziemi i Ksi˛eżyca nazywamy ekliptyka.˛
Kształt Ziemi przybliżony jest trójosiowa˛ elipsoida˛ zwana˛ elipsoida˛ figury Ziemi. Wła-
sności dynamiczne bryły charakteryzowane sa˛ jej elipsoida˛ bezwładności, która ujmuje
zarówno kształt jak i mas˛e bryły. Odpowiednikiem równań Newtona w ruchu post˛epo-
wym, sa˛ w ruchu wirowym równania Eulera. Rozwiazaniem ˛ równań Eulera jest wektor
pr˛edkości katowej
˛ ω
~ , jego kierunek nazywamy chwilowa˛ osia˛ obrotu Ziemi, punkty prze-
bicia tej osi z powierzchnia˛ elipsoidy figury Ziemi to chwilowe bieguny Ziemi, punkty
przebicia ze sfera˛ niebieska˛ to prawdziwe bieguny świata.
Ruch obrotowy Ziemi rozkładany jest na ruch regularny (precesja) i okresowy (nutacja).
Oś ruchu regularnego służy do definicji średnich biegunów świata.
Tzw. ruch biegunów ziemskich dotyczy zmian w położeniu osi obrotu w ciele Ziemi.
Obejmuje on kilka drobnych niezależnych ruchów, jednym z nich jest tzw. Chandlerowski
ruch biegunów o okresie około 14-tu miesi˛ecy. Poszczególnych przyczynków w ruchu
biegunów Ziemi nie da si˛e opisać teoretycznie tak dokładnie jak precesj˛e i nutacj˛e osi
obrotu bryły ziemskiej.
Słowa kluczowe: ekliptyka (ruchoma, nieruchoma), elipsoida figury, elipsoida bezwład-
ności, równania Eulera, precesja, nutacja, ruch biegunów, biegun średni (prawdziwy, rów-
nik średni (prawdziwy), punkt równonocy średni (prawdziwy). a
a
[Modyfikowano AD 2008, kwiecień, 25]
130 Ruch Ziemi
K K
B
a) Z
Z
b) S
B
K
Z
c)
κ
N’
U’ S
U
κ
t
Plaszczyzna
orbity barycentrum
ukladu Ziemia Ksiezyc.
cie płaszyzny orbity barycentrum B Wielkość tego obrotu można opisać za pomoca˛ kata ˛
κ, liczonego od jakiegoś ustalonego nieruchomego poczatku. ˛ Wówczas, przykładowo,
zmiany kata ˛ κ w czasie wygladałyby
˛ tak jak to pokazano na rysunku 12.1c. Linia przery-
wana opisuje ciagły˛ obrót płaszczyzny orbity barycentrum B. Linia ta może być niewiel-
kim fragmentem sinusoidy o bardzo dużym okresie. Linia ciagła ˛ wskazuje na okresowe
wahania w tempie zmian kata ˛ κ.
Jeżeli wykluczymy z rozważań zmiany okresowe, wówczas opisywany ruch płasz-
czyzny barycentrum B (tylko zmiany liniowe) b˛edzie dotyczył pewnej płaszczyzny ro-
tujacej
˛ z szybkościa˛ taka˛ jak średnia pr˛edkość rotacji płaszczyzny orbity barycentrum
B. Taka konceptualna płaszczyzna nosi nazw˛e płaszczyzny ekliptyki, a ekliptyka˛ nazy-
wamy koło wielkie powstałe jako rezultat przeci˛ecia tej płaszczyzny ze sfera˛ niebieska.˛
Zdefiniowana w ten sposób ekliptyka porusza si˛e, dlatego można natknać˛ si˛e na okre-
ślenia ekliptyka ruchoma. Inne określenie to ekpliptyka chwilowa, którym określa
si˛e ekliptyk˛e ruchoma,˛ której położenie odpowiada jakiemuś szczególnemu momentowi
czasu (epoce). Ciag ˛ ekliptyk chwilowych odpowiadajacych ˛ pewnym epokom, niekiedy
nazywany jest mianem ekliptyk nieruchomych na te epoki. Np. mówimy — ekliptyka
nieruchoma 1900.0, ekliptyka nieruchoma epoki 2000.0, ekliptyka daty, czyli ekliptyka
na bieżacy˛ moment czasu.
Zaburzenia planetarne w ruchu barycentrum Ziemi i Ksi˛eżyca, nie tylko zmieniaja˛
orientacj˛e orbity barycentrum. Zmianom ulegaja˛ także parametry określajace ˛ rozmiar i
kształt orbity. Najcz˛eściej, okresowe zaburzenia w ruchu Ziemu, przedstawiane sa˛ w for-
mie sinusoidalnych poprawek do niezaburzonego ruchu średniego. Poprawek może być
bardzo dużo, każda z indywidualnym okresem i amplituda.˛ W końcu XIX wieku, ame-
rykański astronom Simon Newcomb zestawił tablice takich poprawek, uwzgl˛edniajace ˛
nieregularności w ruchu Ziemi wokół Słońca. Newcomb nadał im nazw˛e "Tablice ruchu
Ziemi wokół Słońca". W żargonie astronomów tablice te cz˛esto nazywane sa˛ "Tablicami
Słońca".
k k
r1
m1
mi
dm m2
c
c
r
m3
gdzie ρ jest g˛estościa˛ ciała, r oznacza odległość od osi k danego fragmentu masy, V jest
obj˛etościa˛ ciała.
Przez środek masy bryły da si˛e przeprowadzić nieskończenie wiele osi, wzgl˛edem
każdej z nich można obliczyć moment bezwładności oraz jego odwrotność. Odwrotno-
ści momentów bezwładności można w formie wektorów, w identycznej skali, odłożyć
wzdłuż osi wzgl˛edem których zostały obliczone. Przeprowadzajac ˛ powierzchni˛e b˛edac ˛ a˛
obwiednia˛ końców wektorów uzyskamy brył˛e b˛edac ˛ a˛ trójosiowa˛ elipsoida,˛ tzw. elipso-
ida˛ bezwładności. Jak każda elipsoida trójosiowa, elipsoida bezwładności posiada trzy
osie główne a, b, c, wzgl˛edem których mementy bezwładności nazywane sa˛ głównymi
momentami bezwładności A, B, C. Moment bezwładności C o najwi˛ekszej wartości jest
to moment obliczony wzgl˛edem najkrótszej osi elipsoidy bezwładności.
Powierzchnia fizyczna bryły ziemskiej jest bardzo skomplikowana. Pomimo to istnieje
zupełnie dobre jej przybliżenie geometryczne — trójosiowa elipsoida — zwana elipsoida˛
figury Ziemi. Jej oś biegunowa nosi miano osi figury, a prostopadła do niej płaszczyzna,
przechodzaca ˛ przez środek masy Ziemi, określa równik figury. 1
W teorii ruchu wirowego Ziemi, jako ciała doskonale sztywnego zakłada si˛e, że osie
elipsoidy figury Ziemi pokrywaja˛ si˛e z osiami elipsoidy bezwładności bryły Ziemi, wi˛e-
cej, że pokrywaja˛ si˛e najkrótsze osie obu elipsoid. W rzeczywistości, ruch mas wewnatrz ˛
bryły ziemskiej, powoduje niewielkie skr˛ecenie osi elipsoidy bezwładności Ziemi wzgl˛e-
dem elipsoidy figury, dlatego można jedynie twierdzić, że osie elipsoidy figury pokrywaja˛
si˛e z pewnym średnim położeniem osi elipsoidy bezwładności.
1
Wykorzystywane dość szeroko poj˛ecie równika geograficznego nie ma ścisłego określenia. Wydaje
si˛e, że tak naprawd˛e jego użytkownicy maja˛ na myśli równik figury.
12.2 Ruch wirowy Ziemi 133
Warto jeszcze wspomnieć, że równikowe osie ziemskiej elipsoidy bezwładności sa˛
niemal identyczne, podobnie jest w przypadku elipsoidy figury Ziemi. Stad
˛ bryła ziemska
zwykle opisywana jest za pomoca˛ dwuosiowych elipsoid obrotowych.
Biegun
ekliptyki
.
ϕ ω
υ
.
. B
ψ
rownik A N
figury
Rysunek 12.3: Wektor ω chwilowej pr˛edkości wirowania bryły Ziemi rozłożony na skła-
dowe: ψ̇ wzdłuż kierunku na bieguny ekliptyki, ϕ̇ w kierunku bieguna figury C oraz ϑ̇
wzdłuż lini w˛ezłów ekliptyki i równika figury.
biegun ekliptyki
os figury Os regularnego
obrotu
Rysunek 12.4: Na lewo ilustracja ruchu precesyjnego osi figury i chwilowej osi regular-
nego wirowania Ziemi, po prawej ilustracja ruchu rzeczywistego osi wiropwania bryły
ziemskiej.
jej ruch wirowy odbywa si˛e wokół osi bliskiej osi najwi˛ekszego momentu bezwładności,
czyli najkrótszej osi elipsoidy bezwładności bryły.
W wyrażeniach na każda˛ ze składowych ψ̇, ϑ̇, ϕ̇ tkwia:
˛ składnik stały (prawie stały)
oraz suma dużej liczby niewielkich wyrazów okresowych. Te ostatnie wyrazy nosza˛ na-
zw˛e nutacji, chwilowo zostawimy je na boku a zajmiemy si˛e regularnym ruchem wiro-
wym Ziemi.
12.2.5 Nutacja
Wyrazy nutacyjne w wyrażeniach na składowe ψ̇, ϑ̇, ϕ̇ zniekształcaja˛ nakreślony przed
chwila˛ obraz wirowania Ziemi. Nutacje w ψ̇ i ϑ̇ wywołuja˛ kołysanie osi figury i re-
gularnej chwilowej osi obrotu Ziemi wzgl˛edem ekliptyki. W rezultacie obie osie w swej
w˛edrówce wokół normalnej do ekliptyki nieustannie odchylaja˛ si˛e od pobocznic stożków,
a ślady jakie rysuja˛ w przestrzeni końce tych osi nie sa˛ okr˛egami jak na rysunku 12.4, ale
złożonymi liniami falistymi.
Nutacja w ϕ̇ sprawia niewielkie wahania wokół wartości średniej pr˛edkości dobowego
wirowania Ziemi.
Szybkości ψ̇, ϑ̇ to szybkości przemieszczania si˛e osi figury Ziemi w przestrzeni wzgl˛e-
dem tła gwiazd. W ciele Ziemi oś figury jest nieruchoma.
• Swobodny ruch biegunów o amplitudzie 0.5”, definiowany jest jako ogólne roz-
wiazanie
˛ różniczkowych równań ruchu Eulera, czyli przy MA = MB = MC = 0,
stad
˛ ruch ten nazywany bywa eulerowskim ruchem biegunów. Ze wzgl˛edu na od-
st˛epstwa Ziemi od doskonałej sztywności, realny ruch nie zawsze jest podobny do
teoretycznych przewidywań. Np. zgodnie z rozwiazaniem
˛ ogólnym, chwilowe bie-
guny powinny przemieszczać si˛e po powierzchni Ziemi z okresem bliskim 305 dni,
po okr˛egach o środkach w tzw. biegunach średnich. Tymczasem ich ruch rzeczywi-
sty trwa około 14 miesi˛ecy (okres Chandlerowski) i wykazuje zmienna˛ amplitud˛e.
• Istnieja˛ jeszcze składowe ruchu biegunów o okresie doby i o bardzo nikłej am-
plitudzie (ułamki setnych cz˛eści sekundy). Przyczyna˛ tych drobnych ruchów jest
konieczność ciagłych
˛ zmian wzajemnej orientacji osi figury i osi obrotu w ciele wi-
rujacej
˛ Ziemi (wzgl˛edem otaczajacej ˛ przestrzeni wzajemna orientacja osi figury i
obrotu praktycznie jest stała). Ruch ten wynika głównie z niezerowej wartości skła-
dowej ψ̇ pr˛edkości obrotowej Ziemi, rezultatu przyciagania
˛ Słońca i Ksi˛eżyca. Dla-
tego nazwano te ruchy wymuszonym ruchem biegunów (ruchem słoneczno-ksi˛eż-
ycowym). W wielu zagadnieniach praktycznych (jak np. rachunki redukcyjne) ze
wzgl˛edu na niewielkie amplitudy ruchu wymuszonego, mówiac ˛ o ruchu biegunów
przyczynek ruchu wymuszonego pomija si˛e.
Ruch biegunów o jakim mowa w tym rozdziale, odzielany jest od precesji i nutacji w
pewnym stopniu na zasadzie konwencji. Kierowano si˛e tu tym, że precesja i nutacja
moga˛ być określone precyzyjnie nawet na lata wprzód, czego nie da si˛e w żaden sposób
12.4 Punkt równonocy wiosennej 137
Precesja i nutacja
13.1 Wst˛ep
W rozdziale tym zajmiemy si˛e zjawiskami precesji i nutacji nieco bardziej szczegółowo
niż to miało miejsce wcześnniej. Jak wiemy w wyniku tych zjawisk ulegaja˛ zmmianie
współrz˛edne ciał niebieskich, ale sa˛ to zmiany spowodowane ruchem biegunów ekliptyki
i biegunów równika świata oraz sprz˛eżonych z nimi płaszczyzn ekliptyki i równika. Za-
tem mamy tu do czynienia ze zmianami układu odniesienia jako całości, a nie z ruchem
poszczególnych gwiazd, planet etc.
Podział na precesj˛e i nutacj˛e jest arbitralny, dokonano go kierujac
˛ si˛e wygoda˛ w
pracach pozycyjnych. Regularne długoterminowe przemieszczenie biegunów nazwano
precesja,˛ nutacja˛ nazwano przemieszczenia krótkookresowe zachodzace ˛ wokół średniego
położenia bieguna świata.
W dalszych rozważaniach zakładamy, że gwiazdy na sferze nie wykonuja˛ żadnych
ruchów a ich współrz˛edne ulegaja˛ zmianom jedynie w wyniku ruchu układu współrz˛e-
dnych.
• łuk KP = ε,
◦
• gwiazda w punkcie X ma współrz˛edne (α, δ), łuk P X = 90 − δ, natomiast kat
˛
◦
KP X = 90 + α,
◦
• w układzie ekliptycznym gwiazda ma współrz˛edne (λ, β), KX = 90 −β i P KX =
◦
90 − λ.
Rysunek 13.1: Precesja luni-solarna to ruch średniego bieguna świata P wokół nierucho-
mego bieguna ekliptyki K.
Rysunek 13.2: Zjawiskko precesji bryły ziemskiej można opisać za za pomoca˛ pary sił
F1 , F2 , b˛edacej
˛ rezultatem grawitacyjnego oddziaływania Słońca i Ksi˛eżyca na równi-
kowe wybrzuszenia Ziemi. Wypadkowy, średni moment składowych tych sił prostopa-
dłych do płaszczyzny równika, da˛ży do ustawienia równika Ziemi w płaszczyźnie eklip-
tyki.
πτ . Położenie osi obrotu N N 0 określone jest przez jej długość ekliptyczna˛ (wzgl˛edem
ekliptyki poczatkowej)
˛ i oznaczane jest przez Π. Wobec tego ΥN = Π. Zauważmy też,
że N i N 0 sa˛ biegunami łuku KK 0 jaki biegun ekliptyki zakreśla na sferze niebieskiej.
Precesja planetarna wpływa na (α, δ) gwiazd w sposób bardzo prosty. Ponieważ punkt
P jest teraz nieruchomy to dδ = 0, równonoc zaś doznaje przesuni˛ecia po łuku ΥΥ0 . Dłu-
gość tego łuku wyrażona jest w formie iloczynu λ0 τ , gdzie λ0 zwana jest roczna˛ precesja˛
planetarna.˛ Zatem w efekcie precesji planetarnej
dα = −λ0 τ
(13.6)
dδ = 0
• ΥΥ0 = λ0 τ ,
• ΥN Υ0 = πτ ,
Dla τ dostatecznie małego, gdy sin(πτ ) ≈ πτ , sin(λ0 τ ) ≈ λ0 τ oraz sin (ε + dε) ≈ sin ε
b˛edzie
sin(ε + dε) cos ε cos(λ0 τ ) − cos(ε + dε) sin ε = sin(λ0 τ ) cot Π sin(ε + dε)
144 Precesja i nutacja
sin(ε + dε)
sin dε = πτ cos Π
sin ε
a przy założeniach: sin dε ≈ dε oraz sin(ε + dε) ≈ sin ε b˛edziemy mogli napisać
dε = πτ cos Π (13.8)
Różniczkujac
˛ to równanie dostaniemy
cos β dβ = −(sin δ sin ε + cos ε cos δ sin α)dε − sin ε cos δ cos α dα
czyli
a korzystajac
˛ z równań (13.7) i (13.8), po paru przekształceniach przekonamy si˛e, że
dβ = πτ sin(Π − λ) (13.9)
Czynnik cos δ sin α już wiemy jak wyeliminować, mamy też, że dα = −λ0 τ . Z kolei
(−λ0 ) można z gracja˛ zastapić
˛ prawa˛ strona˛ równania (13.7), natomiast zamiast dβ mo-
żemy wziać
˛ prawa˛ stron˛e równania (13.9) — zatem, po podstawieniach b˛edzie
1
cos β sin λdλ = (− sin β sin ε + cos β cos ε sin λ) · (−πτ sin Π · )−
sin ε
sin β cos λ · πτ sin(Π − λ)
1
dλ = πτ tan β sin Π · − sin Π cot λ cos λ + cos Π cot λ sin λ −
sin λ
cot ε sin Π]
po wyłaczeniu
˛ sin Π z dwóch pierwszych wyrazów w nawiasach okragłych
˛
1 − cos2 λ
dλ = πτ tan β sin Π · + cos Π cos λ − sin Π cot ε
sin λ
czyli
a dalej mamy
dλ = pτ + πτ tan β cos (Π − λ)
(13.15)
dβ = πτ sin (Π − λ)
p = ψ − λ0 cos ε (13.16)
Stałe m, n, p nie sa˛ stałymi absolutnymi bowiem doznaja˛ zmian wiekowych. Na podsta-
wie podanych wcześniej formuł można napisać, że
• P0 P = θA ,
• P0 X = 90◦ − δ0 ,
• P P0 X = α + ζA ,
• P X = 90◦ − δ,
• P0 P X = 180◦ − (α − zA ).
˛ ze soba˛ współrz˛edne (α0 , δ0 ) z epoki t0 ze współrz˛ednymi (α, δ)
Możemy teraz powiazać
z epoki t. W wyprowadzonych formułach b˛eda˛ tkwiły parametry katowe ˛ θA , ζA , zA . I
tak za pomoca˛ trygonometrii sferycznej, słynnego już wzoru 5-cio elementowego, wzoru
sinusów i wzoru cosinusów, odpowiednio, mamy
Wzory (13.19) i (13.20) sa˛ ścisłe, nie dokonaliśmy w trakcie ich wyprowadzania żadnych
założeń upraszczajacych.
˛ Oczywiście by je zastosować, konieczna jest znajomość katów
˛
13.6 Precesyjna macierz obrotu 149
ζA , zA , θA , te zaś można uzyskać z teorii precesji ziemskiej osi obrotu. W praktyce sa˛
one obliczone za pomoca˛ szeregów pot˛egowych interwału czasu (t − t0 ), o wyrazach do
trzeciego rz˛edu włacznie.
˛ Współczynniki szeregów różnia˛ si˛e nieco od epoki do epoki.
Dla epoki J2000 mamy formuły
ε = 23◦ 260 21.00 448 − 46.00 815T − 0.00 001T 2 + 0.00 002T 3 (13.22)
Podobnie b˛edzie dla innej epoki t, kierunek do tej samej gwiazdy podany jest wówczas
wersorem s
x cos δ cos α
s = y = cos δ sin α (13.24)
z sin δ
150 Precesja i nutacja
s = Ps0 (13.29)
s0 = P−1 s (13.30)
a wi˛ec
s 0 = PT s (13.31)
13.7 Nutacja
Nie b˛edziemy zajmowali si˛e dynamiczna˛ teoria˛ precesji i nutacji. Sa˛ to bardzo złożone
teorie wykraczajace ˛ poza ramy naszego wykładu. Podamy jednak pewne użyteczne ko-
mentarze ilustrujace˛ podstawowe aspekty zagadnienia.
Wiemy już, że rzeczywisty ruch bieguna świata po sferze niebieskiej jest bardzo zło-
żony. Z tego też wzgl˛edu rozdzielono go na precesj˛e luni-solarna˛ i nutacj˛e. Nutacja obej-
muje wszystkie okresowe składowe zmian w położeniu prawdziwego bieguna wzgl˛edem
jego położenia średniego.
Wiemy, że przyczyna˛ precesji i nutacji jest moment skr˛ecajacy
˛ pary sił (patrz rysunek
13.2) usiłujacy
˛ ustawić płaszczyzn˛e równika ziemskiego w płaszczyźnie ekliptyki, wiemy
też, że główna˛ przyczyna˛ tego zjawiska sa˛ grawitacyjne oddziaływania pomi˛edzy Ziemia,˛
Ksi˛eżycem i Słońcem. Precsja L-S to regularny ruch średniego bieguna świata wokół
nieruchomego bieguna ekliptyki. Skad ˛ biora˛ si˛e wyrazy nutacyjne? Rozważmy najpierw
moment skr˛ecajacy ˛ pary sił pochodzacy˛ od Słońca działajacy
˛ na oś obrotu Ziemi w chwili
gdy Słońce znajduje si˛e w położeniu (αs , δs ). Z teorii wirowania doskonale sztywnej
Ziemi wiadomo nam, że wielkość momentu skr˛ecajacego ˛ okazuje si˛e być proporcjonalna
do sin 2δs . Skoro tak, to moment ten ma charakter okresowy i np. w warunkach równo-
nocy znika.
152 Precesja i nutacja
˛ moment sił można wyobrazić sobie jako wektor k, który z powodu syme-
Skr˛ecajacy
trii rozważanego zagadnienia jest prostopadły do linii Ziemia-Słońce oraz do chwilowej
osi obrotu Ziemi. A wi˛ec wektor k leży w płaszczyźnie równika i jest skierowany ku
punktowi o rektascensji αs − 90◦ . Jego składowe w układzie współrz˛ednych równiko-
wych określone sa˛ z pomoca˛ zwiazku
˛
Kładac
˛ te zwiazki
˛ do (13.32) dostaniemy
a dalej
Wektor momentu p˛edu ruchu wirowego Ziemi jest skierowany ku punktowi P , czyli ku
biegunowi świata (rysunek 13.6). Jak powiedziano wcześniej, moment skr˛ecajacy ˛ k jest
prostopadły do kierunku chwilowej osi obrotu Ziemi i dlatego nie może zmienić długości
wektora momentu p˛edu Ziemi. Może jednak zmienić jego kierunek, a wi˛ec położenie
bieguna P .
Aby móc dyskutować zmiany tego położenia, czyli śledzić ruch bieguna P , dobrym
pociagni˛
˛ eciem jest posłużenie si˛e układem współrz˛ednych określonym w oparciu o jakieś
wybrane, ustalone położenie bieguna. Niech zatem wersor s(x, y, z) b˛edzie określał po-
łożenie bieguna na sferze, wzgl˛edem osi równikowych, określonych za pomoca˛ średniego
bieguna i średniej równonocy z epoki poczatkowej,
˛ takiej kiedy to długość Słońca L = 0.
Wówczas, jako że mamy do czynienia z drobnymi ruchami, składowe wektora k b˛eda˛
proporcjonalne do (dx/dt, dy/dt, 0), czyli
dx = k0 cos ε sin ε[1 − cos 2L] · dt
dy = k0 cos ε tan ε[− sin 2L] · dt
dz = 0 · dt
albo
dx sin1 ε = ψ1 dt − ψ1 cos 2L · dt
−1
x csc ε = ψ1 t − 0.5 · ψ1 dLdt
sin 2L
gdzie ψ1 = k0 cos ε jest stała˛ zależna˛ od wielkości wirowego momentu p˛edu średniego
momentu skr˛ecajacego
˛ i nachylenia ekliptyki do równika.
Po scałkowaniu wszystkich składowych mamy, że po upływie czasu t od momentu od-
powiadajacego
˛ położeniu poczatkowemu,
˛ współrz˛edne bieguna w przybliżeniu wynosza˛
−1
x csc ε = ψ1 t − 0.5 · ψ1 dL
dt
sin 2L
dL −1
y = 0.5ψ1 tan ε dt cos 2L (13.35)
z=1
Wyrażenie x csc ε jest przemieszczeniem bieguna w długości, y natomiast opisuje przy-
rost w nachyleniu ekliptyki do równika. W równaniu (13.35) można wyróżnić różne
człony, liniowy ze wzgl˛edu na czas wyraz ψ1 stanowiacy˛ przyczynek od precesji słonecz-
nej (stanowi on około 1/3 wpływu) oraz dwa wyrazy nutacyjne, jeden w długości, drugi
w nachyleniu. Oba człony nutacyjne maja˛ okres półroczny. Jeśli czas t wyrazimy w la-
tach to pochodna dL/dt = 2π i wówczas równanie (13.35) wia˛że amplitudy wyrazów
nutacyjnych z tempem precesji słonecznej.
W równaniu (13.35) mogłyby si˛e pojawić dalsze wyrazy pochodzenia czysto słonecz-
nego. Pojawia˛ si˛e one jeśli do odpowiednich formuł wprowadzimy roczne zmiany odleg-
łości Ziemi od Słońca, czyli po uwzgl˛ednieniu mimośrodu orbity Ziemi np. do wyrazów
rz˛edu pierwszego. Wówczas w wyniku sprz˛eżeń z wyrazem precesyjnym powstana˛ nowe
wyrazy nutacyjne o okresie jednego roku. A wskutek sprz˛eżeń z istniejacymi
˛ już czło-
nami nutacyjnymi, powstana˛ dodatkowe dwa człony nutacyjne o okresach roku i czterech
miesi˛ecy, itd.
154 Precesja i nutacja
I = ε + i cos Ω (13.40)
Ponieważ i cos Ω jest wielkościa˛ mała˛ pierwszego rz˛edu, możemy położyć cos (i cos Ω) ≈
1 oraz sin (i cos Ω) ≈ i cos Ω. Odrzucajac ˛ jeszcze wyrazy drugiego rz˛edu ze wzgl˛edu na
(i cos Ω), uzyskamy
Całkujac
˛ to równanie otrzymamy wyrażennia na składowe przemieszczenia bieguna wzgl˛e-
dem jego położenia poczatkowego,
˛ kiedy to Ω = 0
dΩ −1
x csc ε = ψ10 · t + 2 i ψ10 cot 2ε
dt
sin Ω
−1
y = −iψ10 dΩdt
cos Ω (13.42)
z=1
gdzie ψ10 jest stała˛ tak dobrana˛ by reprezentowała wyraz precesyjny w długości. Skła-
dowe te wyrażone sa˛ w układzie współrz˛ednych o osiach zorientowanych zgodnie z
osiami układu równikowego.
Podobnie jak to było dla Słońca, w równaniu (13.42) można zidentyfikować ksi˛eży-
cowy człon precesyjny jak i ksi˛eżycowe wyrazy nutacjne w długości i nachyleniu, ale
tutaj maja˛ one okresy 18.6 lat. Sa˛ to najwi˛eksze człony nutacyjne, około 10 razy wi˛eksze
od 6-cio miesi˛ecznych członów słonecznych, które jeśli chodzi o wielkość sa˛ zaraz na
drugim miejscu.
Kończymy dyskusj˛e przebiegu zjawiska nutacji, to co powiedziano wyżej miało na
celu ukazanie w jaki sposób powstaja˛ najważniejsze człony precesyjne i nutacyjne. Pełna
teoria precesji i nutacji jest bardzo skomplikowana i wykracza poza ramy podstawowego
kursu astronomii sferycznej.
Jeszcze nie tak dawno teoria nutacji oparta była na modelu sztywnej Ziemi. Równania
(13.35) i (13.42) odpowiadaja˛ takiemu właśnie podejściu. W roku 1980 opublikowano
nowa˛ teori˛e nutacji, która˛ w dwa lata później zaaprobowała MUA. Teoria ta oparta jest
na bardziej realistycznym modelu Ziemi, modelu elastycznym osiowo niesymetrycznym.
Jest to tzw. pełna teoria nutacji zawierajaca ˛ po 106 wyrazów zarówno w długości jak i
w nachyleniu. W tej teorii przemieszczenie bieguna świata w długości oznaczone zostało
jako ∆ψ, przemieszczenie prostopadłe do niego przez ∆ε. Przemieszczenia te nazwano
156 Precesja i nutacja
∆ψ = 106 0
P
1 Si sin (ai L + bi L + ci F + di D + ei Ω)
P106 (13.43)
∆ε = 1 Ci cos (ai L + bi L0 + ci F + di D + ei Ω)
• D jest średnia˛ długościa˛ Ksi˛eżyca minus średnia długość Słońca, czyli średnia˛ elon-
gacja˛ Ksi˛eżyca od Słońca,
∆ψ = −17.00 1996 sin Ω − 1.00 3187 sin (2F − 2D + 2Ω) − 0.00 2274 sin (2F + 2Ω) + 0.00 2062 sin (2Ω)
∆ε = 9.00 2025 cos Ω + 0.00 5736 cos (2F − 2D + 2Ω) + 0.00 0977 cos (2F + 2Ω) − 0.00 0895 cos (2Ω)
(13.44)
Współczynniki 17.00 1996 oraz 9.00 2025, niekiedy nazywane sa˛ stałymi nutacji. War-
tości ich jak i pozostałych współczynników w równaniu (13.44) odpowiadaja˛ epoce J2000.0
.
Niekiedy wygodnym jest podział na długo i krótkookresowe człony nutacyjne. Wy-
razy długo-okresowe sa˛ to wyrazy niezależne od średniej długości Ksi˛eżyca, wszystkie te
wyrazy maja˛ okresy wi˛eksze od 90 dni. Pośród wyrazów krótkookresowych nie ma ani
jednego o okresie przekraczajacym˛ 35 dni. Zsumowane wyrazy krótkookresowe ozna-
czane sa˛ przez dψ i dε.
i stad
˛
dδ = ∆ε sin α
˛ P KX daje zwiazek
Wzór sinusów zastosowany do trójkata ˛ na rektascensj˛e
Niech s0 b˛edzie wersorem tego samego kierunku, ale określonym w oparciu o prawdziwy
równik i prwdziwy punkt równonocy na dana˛ dat˛e JD. Na rysunku 13.7b łatwo zauwa-
żyć, że transformacja składowych wersora s określonych wzgl˛edem średniego układu
współrz˛ednych do układu prawdziwego, wymaga złożenia trzech transformacji obrotu.
Mianowicie:
• obrotu wokół osi x układu równikowego średniego o kat
˛ ε0 , przechodzimy wów-
czas do średniego układu ekliptycznego daty JD,
• nast˛epnie obrotu wokół osi z średniego układu ekliptycznego o kat
˛ (−∆ψ), jeste-
śmy wówczas w prawdziwym układzie ekliptycznym daty JD,
• obrotu wokół osi x prawdziwego układu ekliptycznego o kat
˛ (−ε), co daje nam
prawdziwe współrz˛edne równikowe — czyli składowe wersora s0 na dat˛e JD.
Zatem powiazania
˛ wersorów s i s0 można dokonać za pomoca˛ rachunku macierzowego
s0 = N s (13.50)
gdzie
N = p(−ε)r(−∆ψ)p(ε0 ) (13.51)
13.9 Łaczny
˛ wpływ precesji i nutacji w formaliźmie macierzowym 159
lub
" #
cos ∆ψ, − sin ∆ψ cos ε0 , − sin ∆ψ sin ε0
N= − sin ∆ψ cos ε, cos ∆ψ cos ε cos ε0 + sin ε sin ε0 , cos ∆ψ cos ε sin ε0 − sin ε cos ε0
− sin ∆ψ sin ε, cos ∆ψ sin ε cos ε0 − cos ε sin ε0 , cos ∆ψ sin ε sin ε0 + cos ε cos ε0
13.9 Łaczny
˛ wpływ precesji i nutacji w formaliźmie ma-
cierzowym
W notacji macierzowej możemy napisać wyrażenie uwzgl˛edniajace˛ łaczny
˛ wpływ prece-
sji i nutacji. Niech s0 b˛edzie wersorem położenia gwiazdy wzgl˛edem średniego równika
i równonocy z pewnej epoki standardowej, np. J1950.0. Wówczas wersor położenia
gwiazdy wzgl˛edem średniego układu daty JD, uzyskamy za pomoca˛ formuły (13.29) jako
s = Ps0
gdzie P jest precesyjna˛ macierza˛ obrotu. Współrz˛edne prawdziwe tej gwiazdy na dat˛e JD
dostaniemy za pomoca˛ równania (13.50), a wi˛ec
s0 = N s = (NP) s0 = R s0 (13.52)
Macierz R ≡ NP jest także macierza˛ obrotu pozwalajac ˛ a˛ na jednoczesne uwzgl˛ednie-
nie precesji i nutacji. Jej elementy publikowane sa˛ w corocznych wydaniach niektórych
roczników astronomicznych.
Transformacja odwrotna, od współrz˛ednych prawdziwy daty JD do współrz˛ednych
średnich daty JD, otrzymamy bez trudu z równania (13.52) działajac ˛ na obie jego strony,
lewostronnie, macierza˛ odwrotna˛ R−1
R−1 s0 = R−1 R s0 = s0 (13.53)
A szcz˛eśliwie, wobec ortogonalności macierzy obrotu mamy
R−1 = RT = (NP)T = PT NT
PT = r(ζA )q(−θA )r(zA ) (13.54)
NT = p(−ε0 )r(∆ψ)p(ε)
4. Pokaż, że wszystkie pozadiagonalne elementy precesyjnej macierzy obrotu sa˛ małe,
diagonalne natomiast sa˛ bliskie jedności. Udowodnij, że
Rozdział 14
14.1 Streszczenie
W astronomii czas odgrywa szczególna˛ rol˛e. W odróżnieniu od np. fizyka, astronomowi
nie wystarcza tylko zegar, czy jak kto woli chronometr odmierzajacy ˛ precyzyjnie krót-
kie interwały czasu. Astronomowie podobnie jak historycy, potrzebuja˛ skali czasu czyli
sposobu przyporzadkowania
˛ zdarzeniom na osi czasu punktu, czyli daty. Do współcze-
śnie wykorzystywanych skal czasu należa: ˛ średni czas gwiazdowy w Greenwich, średni
czas słoneczny UT1, czas koordynowany UTC. W definicjach pierwszych dwóch skal
wykorzystywane jest zjawisko ruchu wirowego Ziemi i z tego powodu sa˛ to skale czasu
niejednostajnego.
W celu utworzenia skal czasu biegnacego ˛ równomiernie posłużono si˛e koncepcja˛ czasu
dynamicznego. Istnieje szereg realizacji tej koncepcji jak: czas efemerydalny ET, czas
układowy geocentryczny TCG, czas układowy barycentryczny TCB. Obecnie najważniej-
sza˛ rol˛e gra czas atomowy TAI. Jest to skala oparta o zjawiska atomowe a jej jednostka˛
jest sekunda, fundamentalna jednostka układu SI.
Pomi˛edzy skalami czasu można określić bardziej lub mniej precyzyjne zwiazki ˛ pozwa-
lajace
˛ na transformacj˛e daty z jednej skali do drugiej. Przy współczesnej dokładności
pomiaru czasu nie stwierdzono rozbieżności pomi˛edzy skalami ET, TT i TAI. Dlatego
skale te różnia˛ si˛e jedynie o stała˛ zwana˛ offsetem.
Skala TCB została tak zdefiniowana,˛ że różni si˛e od TCG wyrazami wiekowym i okreso-
wymi oraz poprawkami relatywistycznymi zależnymi od ruchu obserwatora i geocentrum
wzgl˛edem barycentrum Układu Słonecznego.
Różnica skal czasu ET i UT1 nie może być podana precyzyjnie z wyprzedzeniem. Przy-
czyna˛ tego sa˛ nieregularności w rotacji ziemskiej, których obecnie nie możemy precyzyj-
nie przewidywać.
W praktyce astronomicznej dat˛e (zwana˛ niekiedy epoka) ˛ piszemy w formie roku z ułam-
kiem. Możemy podać ja˛ jako rok Bessel’a lub rok Juliański. Inny sposób oparty na ciagłej
˛
rachubie dni nazywany jest data˛ Juliańska.˛
Słowa kluczowe: zegar, skala czasu, data, czas gwiazdowy, czas słoneczny, czas efe-
merydalny, UT1, UT2, UTC, południk efemerydalny, czas własny TDT, TT, skale czasu
układowego TCG, TDB, TCB, czas atomowy TAI, rok zwrotnikowy, rok Bessel’a, epoka
B1950, epoka J2000, data juliańska. a
a
[Modyfikowano AD 2006, Grudzień 19.]
162 Koncepcja czasu w astronomii
• jednostka˛ czasu.
Istnieja˛ zegary a wi˛ec i skale czasu o zmiennej jednostce czasu. Sa˛ nimi np. praw-
dziwy czas słoneczny, prawdziwy czas gwiazdowy, czas UT0, czas UT1. Trzeba tu roz-
różnić zmienność jednostki czasu, która daje si˛e precyzyjnie modelować od zmienności
nie poddajacej
˛ si˛e modelowaniu. W pierwszym wypadku wystarczy znaleźć odpowiednie
zależności by utworzyć papierowa˛ skal˛e czasu jednostajnego, w drugim przypadku jest to
niemożliwe.
Określenie stopnia zmienności jednostki czasu zegara wia˛że si˛e z dokładnościa˛ jej
wyznaczenia i z interwałem uśredniania. Przy bardzo małych dokładnościach, wi˛ekszość
zegarów ma stała˛ jednostk˛e.
Data. Pod poj˛eciem daty jakiegoś wydarzenia należy rozumieć podanie momentu
tego zdarzenia w ramach którejś ze skal czasu. Data jest współrz˛edna˛ na osi czasu. Jest
ona zwiazana
˛ z całka˛ funkcji jednostki czasu określajacej
˛ zegar. Stad
˛ by dokonać datowa-
nia niezb˛edne jest wi˛ec ustalenie stałej całkowania, punktu zerowego skali, tzw. epoki. 2
Interwał jednoznaczności. Zegary maja˛ na ogół urzadzenia ˛ wskaźnikowe o ogra-
niczonej pojemności cyfrowej. Może to spowodować, że odczytana wartość wskazania
zegara jest wieloznaczna. W wielu produkowanych mechanizmach zegarowych inter-
wałem jednoznaczności jest okres 12 godzin i gdy taki mechanizm wykorzystuje si˛e do
wskazywania czasu gwiazdowego bywa to dość kłopotliwe. Cz˛esto do porównania dwóch
odległych zegarów wykorzystuje si˛e elektromagnetyczne sygnały cz˛estotliwości wzorco-
wej. W tym przypadku interwałem jednoznaczności jest okres drgań sygnału.
Czas gwiazdowy
Czas gwiazdowy jest katem
˛ godzinnym punktu równonocy. Pomijajac ˛ drobne wpływy
precesyjne, skala tego czasu opiera si˛e wyłacznie
˛ o zjawisko ruchu wirowego Ziemi. Po-
nieważ tempo ruchu wirowego Ziemi nie jest stałe, skala czasu gwiazdowego nie jest
2
Z poj˛eciem daty, epoki mamy drobne zamieszanie bowiem można napotkać sytuacje, w których słowo
epoka bywa stosowane zamiast słowa data. W toku naszego wykładu epokami b˛edziemi nazywać pewne
wyróżnione punkty na osi czasu, czyli wyróżnione daty.
164 Koncepcja czasu w astronomii
Czas słoneczny
Czas słoneczny zdefiniowany jest jako kat˛ godzinny środka tarczy Słońca. Wybór
Słońca stanowi przyczyn˛e dużych nieregularności tej skali czasu. W celu ich elimina-
cji wprowadzono koncepcj˛e średniego czasu słonecznego. Pomijajac ˛ duże ale usuwalne
nieregularności, prawdziwy czas słoneczny jest hybryda˛ (krzyżówka)
˛ dwóch nie dajacych
˛
si˛e połaczyć
˛ zjawisk okresowych: dobowego ruchu wirowego Ziemi i rocznego ruchu
orbitalnego Ziemi wokół Słońca. Dlatego w praktyce astronomicznej stosowane sa˛ dwie
odmienne skale czasu słonecznego, mianowicie: czas uniwersalny UT wyłacznie ˛ zależny
od rotacji Ziemi oraz czas efemerydalny ET b˛edacy
˛ dynamiczna˛ skala˛ czasu, zależna˛ je-
dynie od ruchu orbitalnego Ziemi.
Czas UT określony jest jako, obserwowany w Greenwich kat ˛ godziny (HG ) fikcyj-
nego obiektu tzw. uniwersalnego słońca średniego (U s ). Aby poczatek
˛ doby słonecznej
przypadał o północy, do kata
˛ godzinnego musimy dodać dwanaście godzin:
U T = 12h + HG U S (14.1)
γ G
U E H
Formalnie czas ET definiowany jest poprzez równanie 14.22. Ale można go rów-
nież zdefiniować za pomoca˛ kata
˛ godzinnego (HE ) efemerydalnego słońca średniego,
mierzonego wzgl˛edem południka efemerydalnego. Zatem czas efemerydalny ET można
zdefiniować jako
ET = 12h + HE Es (14.2)
DT = ET − U T (14.3)
Wielkość ta nie może być podana ściśle z wyprzedzeniem ponieważ tkwia˛ w niej niere-
gularności ziemskiej rotacji, których jak dotad˛ nie potrafimy precyzyjnie przewidywać.
Teorie dynamiczne umożliwiaja˛ obliczenie heliocentrycznych położeń planet na do-
wolny moment ET. Jest tak również w przypadku położenia środka Ziemi, a wi˛ec można
zestawiać geocentryczne efemerydy ciał niebieskich z czasem ET jako zmienna˛ nieza-
leżna.˛
Czas atomowy
Współczesne zegary atomowe umożliwiaja˛ najbardziej dokładny i zgodny wewn˛etrz-
nie pomiar czasu. Czas atomowy, tzw. mi˛edzynarodowy czas atomowy TAI jest praktycz-
nym czasowym standardem, najbardziej zgodnym z definicja˛ jednostki czasu w układzie
SI. Czyni on zadość również potrzebom astronomów jeśli chodzi o precyzj˛e, długotrwala˛
stabilność oraz niezawodność.
˛ uśredniona˛ ze skal atomowych ATi kilku wybranych
Skala TAI jest skala˛ “papierowa”
laboratoriów czasu. Jest ona redukowana na wpływ pola grawitacyjnego do poziomu
morza i odnosi si˛e do grupy zegarów nieruchomych wzgl˛edem powierzchni Ziemi.
166 Koncepcja czasu w astronomii
Instytucja˛ uprawniona˛ do tworzenia TAI jest Bureau International des Poids et Mesu-
res (BIPM), czyli w ludzkim j˛ezyku Mi˛edzynarodowe Biuro Miar i Wag w Paryżu. 3 Czas
TAI wprowadzono w roku 1972, ale zegary atomowe pracowały już w latach 1950-tych,
dlatego możliwa˛ jest ekstrapolacja TAI aż do lipca roku 1955.
Koncepcja skali czasu atomowego polega na zliczaniu cykli sygnału o wysokiej cz˛e-
stości utrzymywanego w rezonansie z przejściami atomowymi cezu 133. Fundamentalna
jednostka czasu w układzie SI określona jest jako interwał w trakcie którego, zliczonych
zostanie 9192631770 cykli o cz˛estotliwości magnetycznego rezonansu atomowego ato-
mów cezu 133, odpowiadajacych ˛ przejściom pomi˛edzy poziomem F = 4, mf = 0 i
F = 3, mf = 0 przy zerowym zewn˛etrznym polu magnetycznym.
Istnieje zatem koncepcyjna różnica pomi˛edzy skala˛ czasu atomowego a czasem dyna-
micznym ET. Jednak przy założeniu absolutnego charakteru stałych fizycznych, skale te
powinny dać si˛e powiazać
˛ ściśle. Dlatego mimo że definicje sekund TAI i ET sa˛ formal-
nie niezleżne, numerycznie sa˛ ze soba˛ zgodne, a w takiej sytuacji skale te moga˛ różnić
si˛e jedynie stałym przesuni˛eciem (offsetem). Różnica mi˛edzy TAI i ET byłaby istotna
gdyby okazało si˛e, że np. stała grawitacji zmienia si˛e w kosmologicznej skali czasu. Ta-
kie zmiany sugeruja˛ niektóre nieortodoksyjne teorie kosmologiczne, a to oznaczałoby, że
czas dynamiczny spowalnia w porównaniu z czasem atomowym.
W fizyce klasycznej czas jest wielkościa˛ absolutna,˛ a rozróżnienie pomi˛edzy skalami
atomowa˛ i dynamiczna˛ dotyczy jedynie kwestii ich odczytu i definicji. W teorii wzgl˛ed-
ności różnica mi˛edzy tymi skalami ma wi˛eksze znaczenie i odpowiada różnicy pomi˛edzy
czasem własnym i czasem laboratoryjnym (układowym).
stu ogólnej teorii wzgl˛edności skala TDT jest skala˛ czasu własnego na powierzchni Ziemi
wskazywanego przez zegary atomowe. Oznacza to, że TDT jest skala˛ zwiazana˛ z TAI,
a tymczasem w celu obliczenia efemerydy geocentrycznej wymagana jest skala czasu
odczytywana w geocentrum. Dlatego po roku 1991 czas TDT zastapiono ˛ dwiema ska-
lami czasu: skala˛ TT zwiazan˛ a˛ z TAI, oraz ze skala˛ TCG (Geocentric Coordinate Time)
“odczytywana” ˛ w geocentrum i zwiazan˛ a˛ z efemerydami geocentrycznymi.
Druga˛ skala˛ czasu rekomendowana˛ w 1976 roku przez MUA jest barycentryczny czas
dynamiczny TDB (Barycentric Dynamical Time). Czas ten nadaje si˛e do badania ruchów
ciał wzgl˛edem barycentrum Układu Słonecznego. Nie jest to skala czasu zdefiniowana
jednoznacznie, zależy bowiem od wyboru konkretnej post-newtonowskiej teorii grawita-
cji.
Na TDB nałożono warunek by pomi˛edzy ta˛ skala˛ a TDT miały miejsce jedynie różnice
okresowe. Warunek ten może być spełniony w dowolnej teorii grawitacji. W ogólnej
teorii wzgl˛edności TDB odpowiada czasowi układowemu.
Od 1991 roku do obliczeń efemeryd barycentrycznych MUA zaleca skal˛e czasu TCB
(Barycentric Coordinate Time), zast˛epuje ona skal˛e TDB, od której różni si˛e trendem
wiekowym.
TG = HGΥ (14.4)
gdzie HGΥ to kat ˛ godzinny punktu równonocy wiosennej podany wzgl˛edem południka
Greenwich.
W równaniu tym punkt równonocy zawsze odnosi si˛e do równonocy daty, ale jeśli chodzi
o nutacj˛e to może być ona brana pod uwag˛e lub nie. I dlatego, jeżeli w równaniu (14.4)
mamy na myśli prawdziwa˛ równonoc, definiowany czas gwiazdowy określony jest mia-
nem prawdziwego czasu gwiazdowego. Jeśli Υ jest średnia˛ równonoca,˛ równanie (14.4)
definiuje średni czas gwiazdowy. Różnica pomi˛edzy tymi czasami nazywana jest równa-
niem równonocy, oznaczymy ja˛ przez DE. Korzystajac ˛ z wyrażeń na poprawki nutacyjne
w rektascensji i deklinacji można pokazać, że
DE = ∆ψ · cos ε0 (14.5)
Czas gwiazdowy średni różni si˛e od okresu rotacji Ziemi, przyczyna˛ jest zjawisko
precesji rotacyjnej osi Ziemi. Z powodu precesji rektascensja gwiazdy znajdujacej ˛ si˛e
w tym samym miejscu co punkt równonocy, w ciagu ˛ roku zwi˛ekszyłaby si˛e o ψ cos ε0 .
Odpowiada to przyrostowi dziennemu 0.s 0084 i o taka˛ właśnie wartość okres wirowania
Ziemi przewyższa długość średniej doby gwiazdowej.
168 Koncepcja czasu w astronomii
Precesja wpływa również na długość roku, a ściślej długość roku zwrotnikowego. Po-
dobnie jak i miesiac,
˛ rok można zdefiniować na wiele sposobów. Posiadajacy ˛ najwi˛eksze
znaczenie praktyczne rok zwrotnikowy jest średnim interwałem potrzebnym na przyrost
(liczony od średniego punktu równonocy) średniej długości Słońca o 360◦ . Odpowiada to
interwałowi od jednej równonocy wiosennej do drugiej, czyli jest to interwał pomi˛edzy
kolejnymi identycznymi porami roku.
Kalendarz cywilny został dopasowany właśnie do tego roku. W kalendarzu gregoriań-
skim uwzgl˛edniono wyst˛epowanie lat przest˛epnych w liczbie 97-miu na każde 400 lat.
Żadanie
˛ to sprawia, że czas trwania średniego roku kalendarzowego jest bliski długości
średniego roku zwrotnikowego. W kalendarzu juliańskim posiadajacym ˛ lata przest˛epne
co 4 lata, długość roku wynosi 365.25 — jest to rok julianski. Stulecie juliańskie trwa
zatem 36525 dni.
Okres obiegu Ziemi wokół Słońca, w odniesieniu do gwiazd stałych nazywamy ro-
kiem gwiazdowym. Perturbacje planetarne sprawiaja,˛ że jest on nieco różny od interwału
potrzebnego Ziemi na przejście od jednego perihelium do nast˛epnĺgo, zwanego rokiem
anomalistycznym.
Mamy także rok zaćmieniowy (smoczy), definiowany jako interwał mi˛edzy dwoma
kolejnymi przejściami średniego słońca przez w˛ezeł wst˛epujacy˛ średniej orbity Ksi˛eżyca.
Rok ten różni si˛e wyraźnie od pozostałych, określa średnie cz˛estości wyst˛epowania za-
ćmień Słońca i Ksi˛eżyca. Wreszcie możemy napotkać interwał nazywany rokiem gaus-
sowskim wyznaczonym droga˛ obliczeń z II-go prawa Keplera, w którym podstawiono
stała˛ Gaussa k = 0.01720209895 oraz wielka˛ półoś równa˛ a = 1. Poniżej podano długo-
ści (w średnich dobach słonecznych) odpowiadajace ˛ zdefiniowanym wcześniej interwa-
łom rocznym:
CGG = HG X + RA X (14.7)
gdzie CGG oznacza czas gwiazdowy w Greenwwich, HG X kat ˛ godziny obiektu X mie-
rzony w Greenwich natomiast RA X jest rektascensja˛ obiektu.
˛ rachuby rektascensji obierzemy średni punkt równonocy Υ, wów-
Jeśli jako poczatek
czas równanie (14.7) może posłużyć do znalazienia zwiazku
˛ mi˛edzy średnim czasem
gwiazdowym i czasem uniwerslanym. Równanie to dla obiektu X identycznego z uni-
wersalnym słońcem średnim U S , zgodnie z (14.1), ma postać
U T = CGG − RA U S + 12 (14.8)
14.4 Czas gwiazdowy i słoneczny 169
przy czym
Kładac
˛ (14.9) do (14.8) otrzymamy formuł˛e
wykorzystywana w opisie ruchu ciał w Układzie Planetarnym, uzyskała status skali jed-
nostajnej. I to z mocy samej koncepcji, w myśl której skala ta jest wykorzystywana w
równaniach newtonowskiej dynamiki — T tkwi w tych równaniach jako zmienna nieza-
leżna.
Zatem w równaniu (14.9) zamiast T należy napisać TE , a zamiast RA U S (T ) wsta-
wić RA U S (TE ) = RA ES , gdzie czas efemerydalny TE liczony jest od epoki
AD1900, sty.0.12h ET , a RA ES oznacza efemerydalne słońce średnie. Wówczas
równanie (14.9) opisuje ruch fikcyjnego punktu poruszajacego ˛ si˛e jednostajnie po rów-
niku z pr˛edkościa˛ bliska˛ średniej pr˛edkości słońca prawdziwego. Ruch tego obiektu nie
ma nic istotnie wspólnego z ruchem wirowym Ziemi, badanym za pomoca˛ obserwowa-
nych katów
˛ godzinnych gwiazd.
Na to by równanie (14.10) było spójne z obserwowanymi katami ˛ godzinnymi gwiazd
potrzeba innej kancepcji czasu służacego˛ jako miara ruchu wirowego Ziemi. Utworzono
ja˛ w bardzo prosty sposób. Otóż, ogłoszono, że równanie (14.10) definiuje nowa˛ skal˛e
czasu UT, nieregularna,˛ zwiazan ˛ a˛ z wirowaniem Ziemi, przy czym zmienna T w równa-
niu (14.10) to interwał TU w stuleciach juliańskich liczony od epoki AD1900, sty.0, 12h U T .
Różnica
DT
RA ES − RA U S = 1 + = 1.002737909 DT
365.2422
DT = ET − U T
mi˛edzy skalami ET i UT, może być wyznaczona jedynie za pomoca obserwacji, a to
oznacza, że znana b˛edzie z pewnym opóźnieniem. Ze wzgl˛edu na trudne w modelowa-
niu nieregularności skali UT, DT może myć ekstrapolowana w przyszłość z ograniczona˛
precyzja.˛
W latach 1960-1983 równanie (14.11) z TU było wykorzystywane w rocznikach as-
tronomicznych do obliczenia CGG na 0h U T
Od roku 1984, w celu powiazania
˛ skal czasu gwiazdowego i U T stosowano formuł˛e
(CGG1)0h U T 1 = 24110.s 54841 + 8640184.s 812866 TU + 0.s 093104 TU2 − 6.s 2 · 10−6 TU3
JD − J2000U T 1
TU = (14.12)
36525
J2000U T 1 = AD2000, sty.1, 12H U T 1
1
TU ∈ [±0.5, ±1.5, ±2.5 + . . .] ·
36525
Jest to formuła zgodna z położeniem i ruchem punktu równonocy zdefiniowanym w IAU
1976 ”System stałych astronomicznych”, IAU 1980 ”Teoria nutacji” oraz zgodna z poło-
żeniem i ruchem gwiazd podanych w katalogu FK5.
Analogicznie do równania (14.11) definiujacego
˛ skal˛e UT, równanie (14.12) definiujje
skal˛e czasu UT1. Definicja ta zacz˛eła obowiazywać
˛ o epoki AD1984, sty.1.0, a wyraz
stały w (14.12) dobrano w taki sposób, by zachować ciagłość
˛ pomi˛edzy stara˛ i nowa˛
koncepcja.˛ By zagwarantować ciagłość
˛ w tempie zmian skali UT1 w stosunku do UT
dobrano współczynnik przed TU . Jednak nie usuni˛eto pewnych nieciagłości
˛ majacych
˛
źródło w zmianach w teorii nutacji wykorzystywanej w definicji skali UT1.
14.4 Czas gwiazdowy i słoneczny 171
d(CGG1)
s0 = = 8640184.s 812866 + 2 · 0.s 093104 TU − 1.s 86 · 10−5 TU2 (14.15)
d(U T 1)
Ze wzgl˛edu na zmiany tempa wirowania Ziemi, długość doby gwiazdowej i doby UT1
ulegaja˛ zmianom, ale ich stosunki dane równaniami (14.17), (14.18) zawsze pozwola˛ na
przeliczenie interwału czasu wyrażonego w jednej z tych skal na interwał wyrażony w
drugiej skali. Pomijajac
˛ drobne wyrazy wiekowe, z równań tych mamy
M CGS = M CGW − DE
CGG1 = M CGS − λ
θ = CGG1 − (CGG1)0h U T 1
UT 1 = θ · r
obserwacji, trzeba dokonać redukcji stosujac˛ chwilowa˛ wartość λ. Tak uzyskana˛ skal˛e
czasu UT oznaczamy skrótem UT1. I właśnie t˛e form˛e czasu UT stosujemy w wyraże-
niach wia˛żacych
˛ UT z innymi skalami czasu. Odtad ˛ gdziekolwiek napiszemy UT, należy
rozumieć to jako skal˛e UT1. Różnica pomi˛edzy tymi skalami czasu wynosi
przyj˛eto jako definicj˛e skali ET, w której T jest interwałem czasu wyrażonym w stuleciach
juliańskich, biegnacego jednostajnie od epoki poczatkowej
˛ AD 1900 styczeń 0, 12h ET.
Równanie (14.22) poprzez współczynnik przy T , definiuje także jednostk˛e czasu efe-
merydalnego. Współczynnik ten daje przyrost średniej długości Słońca w czasie 36525
174 Koncepcja czasu w astronomii
jest natomiast sekunda SI, określona jako 9192631770 okresów radiacji odpowiadaja- ˛
cych przejściom pomi˛edzy dwoma nadsubtelnymi poziomami znajdujacego
˛ si˛e w stanie
podstawowym atomu cezu 133. Mi˛edzy tymi jednostkami nie ustalono żadnych systema-
tycznych różnic i dlatego pomi˛edzy ET i TAI mamy zwiazek
˛
ET = T AI + 32.s 184 (14.24)
W roku 1984, w celu naprawy sytuacji skal˛e czasu ET oficjalnie zastapiono
˛ nowymi
koncepcjami, skala˛ TDT oraz TDB. Na skutek tych pociagni˛ ˛ eć, w tymże roku czas ET
znika ze wszystkich roczników astronomicznych, gdzie zastapiono
˛ go nowa˛ podstawowa˛
astronomiczna˛ skala˛ czasu tzw. dynamicznym czasem ziemskim TDT. Skal˛e TDT wpro-
wadzona˛ już w roku 1977, zaprojektowano tak, by zachować ciagłość
˛ z czasem efemery-
dalnym ET. Jest ona bliższa skali TAI aniżeli skala ET gdyż sekunda˛ TDT jest sekunda
SI. A zatem by powiazać
˛ TDT z TAI, wystarczy znać wartość offsetu o jaki skale te si˛e
różnia.˛ Zwiazek
˛ ten ma postać
1977 styczen 1.0 TAI = 1977 styczen 1.0003725 TDT (14.25)
W celu zagwarantowania ciagłości
˛ z ET, jest to ten sam offset co podany w równaniu
(14.24).
W roku 1991 na konferencji MUA uporzadkowano
˛ z kolei skal˛e czasu TDT, która˛ do
tego momentu stosowano do obliczeń geocentrycznych efemeryd, a co oznacza, że miała
ona status skali dynamicznej. A tymczasem jak wynika z (14.25) jej definicja oparta jest
o skal˛e czasu atomowego odczytywanego na geoidzie a nie w geocentrum. Stad ˛ w celu
uporzadkowania
˛ poj˛eć wprowadzono dwie nowe skale czasu. Skal˛e czasu TT (Terrestial
Time) odczytywana˛ na geoidzie b˛edac˛ a˛ kontynuacja˛ skali TDT i ET. Sekunda˛ skali TT
jest sekunda TAI, a od TAI wszystkie te trzy skale różnia˛ si˛e jedynie offsetem
T T = T DT = ET = T AI + 32.s 184 (14.26)
Dla celów obliczeniowych geocentrycznych efemeryd zdefiniowano nowa˛ skal˛e, tzw.
czas układowy geocentryczny TCG (Geocentric Coordinate Time). Różni si˛e on od czasu
TT o niewielki wyraz wiekowy
T CG = T T + LG (JD − 2443144.5) · 86400 (14.27)
gdzie LG = 6.969291 · 10−10 ± 3 · 10−16 .
Do obliczeń efemerydalnych zwiazanych
˛ z barycentrum Układu Słonecznego, w roku
1976 MUA zaleciła stosowanie skali TDB (Barycentric Dynamical Time). Jest to czas
układowy odczytywany w barycentrum Układu Słonecznego, a w równaniach ruchu wzgl˛e-
dem barycentrum czas ten ma ststus zmiennej niezależnej. W praktyce skal˛e TDB otrzy-
muje si˛e ze skali TDT za pośrednictwem formuły
T DB = T DT + P
przy czym P obejmuje 500 wyrazów trygonometrycznych (patrz dygresja czas własny i
czas układowy).
W roku 1991 skal˛e TDB zastapiono
˛ czasem TCB (Barycentric Coordinate Time). Jest
to skala czasu z sekunda˛ SI (sekunda atomowa), odczytywana w barycentrum Układu
Słonecznego. Obie skale barycentryczne różnia˛ si˛e jedynie trendem wiekowym
T CB = T DB + LB (JD − 2443144.5) · 86400 (14.28)
176 Koncepcja czasu w astronomii
Możemy też obliczyć juliańska˛ epok˛e na dowolny inny moment. Przykładowo inna stan-
dardowa epoka J1950.0 oddalona jest o dokładnie 18262.5 od epoki fundamentalnej,
mamy wi˛ec
J1950.0 ≡ 1950 sty. 1.d 0
Podany wyżej nowy system epoki juliańskiej wprowadzono w 1976 roku łacznie˛ z rewizja˛
stałych astronomicznych. W tym samym czasie przedefiniowano stara˛ epok˛e bessl’owska˛
poprzez uproszczenie definicji roku Bessel’a. Odtad ˛ rok Bessel’a ma być równy co do
długości rokowi zwrotnikowemu 1900.0.
W astronomii bardzo cz˛esto operuje si˛e data˛ wyrażona˛ jedynie za pomoca˛ dni. Do
tego celu wykorzystuje si˛e poj˛ecie daty juliańskiej JD, wyrażonej z pomoca˛ liczby dni
oraz ułamka dnia jakie min˛eły od epoki 4713 pne styczeń 1.d 5. Juliańskie daty podanych
wcześniej epok fundamentalnych wynosza: ˛
B1900.0 = JD 2415020.313
B1950.0 = JD 2433282.423 (14.33)
J2000.0 = JD 2451545.0
Poczatkowo
˛ JD zdefiniowano z pomoca˛ skali UT, a kolejne dni zliczano od średniego
południa Greenwich, 1 styczeń 4713 pne. Podobna definicja oparta o skal˛e ET nosi nazw˛e
juliańskiej daty efemerydalnej. Można wykorzystać w tym celu i inne skale. W praktyce
jeżeli zachodzi obawa nieporozumienia, trzeba podać jakiego typu JD ma si˛e na myśli. W
równaniach (14.33) jest chyba naturalnym, że pierwsze dwie epoki wyrażone sa˛ w skali
ET, a ostatnia w TT.
Epoki poczatkowe
˛ daty juliańskiej dla dowolnych skal czasu formalnie sa˛ zawsze takie
same (4713 . . .) ), ale nie odpowiadaja˛ one temu samemu momentowi czasu. Niestety,
niewiele wiadomo o relacjach pomi˛edzy różnymi skalami czasu w odległej przeszłości.
Zwiazki
˛ pomi˛edzy epoka˛ Juliańska˛ i Bessel’a oraz bieżacymi
˛ JD sa˛ nast˛epujace:
˛
epoka Julianska = J2000.0 + (JD − 2451545)/365.25
(14.34)
epoka Bessel0 a = B1900.0 + (JD − 2415020.31352)/365.242198781
Cz˛eść całkowita daty juliańskiej nosi nazw˛e dnia juliańskiego. Umówiono si˛e, że dzień
juliański rozpoczyna si˛e w momencie południa, nie o północy. W celu skrócenia liczby
cyfr koniecznych przy zapisie JD, wprowadzono zmodyfikowany dzień juliański — MJD,
M JD = JD − 2400000.5 (14.35)
W tym sposobie rachuby dni, nowy dzień rozpoczyna si˛e o północy, a epoka˛ poczatkow ˛ a˛
jest 1858 listopad 17.0.
Zmodyfikowany dzień juliański nie jest jedyna˛ odmiana˛ dnia juliańskiego istnieje ich
wi˛ecej ale sa˛ to modyfikacje nie majace
˛ charakteru standardu. Przykładowo, w zagadnie-
niach satelitarnych wygodnie jest zliczać dni poczawszy
˛ od roku 1956.
E = RA U s − RA (14.37)
E ∗ = RA Es − RA (14.38)
E ∗ = HE − HE Es = HE − ET + 12
ET = 12 − E ∗ (14.39)
gdzie:
6
Duża litera oznacza, że mamy na myśli prawdziwe słońce.
14.9 Dygresja: czas własny i czas układowy 179
• m = GM/c2 ≈ 1.5 km, G to stała grawitacyjna, M jest masa˛ Słońca, c jest pr˛ed-
kościa˛ światła,
• a jest półosia˛ wielka˛ orbity ziemskiej,
• e jest mimośrodem orbity ziemskiej,
• n jest ruchem średnim orbity ziemskiej,
• E jest anomalia˛ mimośrodowa˛ Ziemi,
• t0 , s0 sa˛ odpowiednio, czasem laboratoryjnym i własnym momentu przejścia Ziemi
przez perihelium.
Ograniczajac˛ si˛e do wyrazów zawierajacych
˛ co najwyżej pierwsze pot˛egi m nieistotne jest
rozróżnienie pomi˛edzy klasycznymi i relatywistycznymi wartościami elementów orbity.
B˛edziemy identyfikowali czas dynamiczny ziemski TDT z czasem własnym. W opar-
ciu o czas laboratoryjny spróbujemy skonstruować globalny czas dynamiczny tzn. czas
TDB, wspomiany już wcześniej. Zakładamy, że mi˛edzy TDT oraz TDB moga˛ istnieć je-
dynie różnice okresowe. W równaniu (14.40) sa˛ one reprezentowane przez drugi wyraz.
A zatem jako TDB adoptujemy taki czas laboratoryjny T , dla którego mamy
−1
3m
T DB = T − T0 = 1 + (t − t0 ) (14.41)
2a
˛ równanie (14.40) przez odwrotność (1 + 3/2 · m/a) b˛edziemy mieli
Mnożac
−1 −1
3m 2m e sin E 3m
1+ (t − t0 ) = (s − s0 ) + · 1+
2a a n 2a
po rozwini˛eciu w szereg po prawej stronie wyrażenia w nawiasie kwadratowym, ograni-
˛ si˛e do wyrazów z m w pierwszej pot˛edze,
czajac
3m −1 3m 3m
(1 + ) ≈1− + ..., bowiem mamy <1
2a 2a 2a
w wyrażeniu na T DB dostaniemy
2m e sin E
T DB = T DT + (14.42)
a n
Kładac
˛ do tego równania
2m = 2.956 [km],
a = 1.496 · 108 [km],
e = 0.01671,
n = 1.991 · 10−7 [rad//sek]
dostaniemy
T DB = T DT + 0.s 001658 sin E (14.43)
albo po wyeliminowaniu anomalii mimośrodowej E, b˛edzie
T DB = T DT + 0.s 001658 sin M + 0.s 000014 sin (2M ) (14.44)
gdzie M jest anomalia˛ średnia˛ Ziemi.
180 Koncepcja czasu w astronomii
3. Udowodnij, że czas gwiazdowy średni w Greenwich dany jest dla dowolnej epoki
T wzorem
4. Oblicz czas gwiazdowy w Greenwich odpowiadajacy ˛ dacie 1996 styczeń 1.0 oraz
1997 styczeń 1.0. Rachunek wykonaj z pomoca˛ formuły (14.12) z wykładu oraz po-
danej wyżej formuły uproszczonej. Wyniki porównaj z odpowiednimi wartościami
z rocznika astronomicznego.
Bibliografia