06 Współrzędne geocentryczneid 6215

background image

Rozdział 6

Współrz˛edne geocentryczne

Streszczenie

Rzeczywisty kształt bryły ziemskiej przybli˙zany jest ró˙znymi powierzchniami np. geoid ˛

a, elipsoid ˛

a obro-

tow ˛

a, sfer ˛

a. Poło˙zenie obserwatora na rzeczywistej powierzchni Ziemi okre´slane jest za po´srednictwem tych

powierzchni. Geoida to powierzchnia zamkni˛eta, wsz˛edzie pozioma, pokrywaj ˛

aca si˛e ze ´srednim poziomem

zrównowa˙zonego grawitacynnie oceanu. Nie ppotrafimy jej opisa´c prostym wyra˙zeniem analitycznym, dlat-
ego zamiast geiod ˛

a, w praktyce posługujemy si˛e przybli˙zonym opisem kształtu Ziemi — elipsoid˛e obro-

tow ˛

a aa nawet sfer ˛

a. Wzgl˛edem elipsoidy poło˙zenie jest ustalone za pomoc ˛

a współrz˛ednych geodezyjnych

(;

;

h)

obserwatora, czyli szeroko´sci i długo´sci geodezyjnej, oraz wysoko´sci nad elipsoid ˛

a, odpowiednio.

Współrz˛edne te słu˙z ˛

a do obliczenia składowe geocentrycznego wektora poło˙zenia obserwatora. Instrumenty

astronomiczne ustawiane s ˛

a na powierzchni Ziemi tak by główna o´s instrumentu pokrywała si˛e z kierunkiem

lokalnej siły ci˛e˙zko´sci (kierunek pionu). Kierunek ten okre´sla astronomiczny zenit obserwatora, który nie
jest identyczny z zenitem geodezyjnym obserwatora (z kierunkiem jaki tworzy normalna do elipsoidy), czy
te˙z z jego zenitem geocentrycznym (z kierunkiem geocentrycznego wektora miejsca obserwacji). Wzgl˛edem
ka˙zdego z tych zenitów mo˙zna okre´sli´c układ współrz˛ednych geograficznych. Znajomo´s´c składowych geo-
centrycznego wektora poło˙zenia obserwatora jest konieczna ww celu transformacji rezultatów obserwacji z
układu odniesienia topocentrycznego do układu geocentrycznego lub odwrotnie. Pełna transformacja obej-
muje dwa zjawiska: paralaks˛e geocentryczn ˛

a i aberrecj˛e dobow ˛

a. Oba zjawiska powoduj ˛

a zmiany współrz˛e-

dnych okre´slaj ˛

acych kierunek propagacji promieniowania elektromagnetycznego. Pierwsze wynika z faktu

przeniesienia pocz ˛

atku układu współrz˛ednych z powierzchni Ziemi do ´srodka masy Ziemi, drugie to efekt

wpływu niezerowej (wzgl˛edem ´srodka Ziemi) pr˛edko´sci obserwatora. Poniewa˙z warto´s´c paralaksy zale˙zy od
odległo´sci pomi˛edzy obiektem i obserwatorem, pomiary paralaks słu˙z ˛

a do okre´slania odległo´sci ciał niebies-

kich. Wyznaczenie w pierwszej połowie XX stulecia paralaksy geocentrycznej planetoidy Eros pozwiło na
okre´slenie warto´sci jednostki astronomicznej

(1

A

U

)

w metrach. Zmiany współrz˛ednych poło˙zenia obiektu

spowodowane aberracj ˛

a dobow ˛

a osi ˛

agaj ˛

a warto´sci do

0:

00

34

, nie zale˙z ˛

a od odległo´sci obiektu od obserwatora,

zale˙z ˛

a natomiast od warto´sci ułamka

V

=

, stosunku szybko´sci obserwatora do szybko´sci ´swiatła. Efekty

relatywistyczne aberracji dobowej, ze wzgl˛edu na mał ˛

a warto´s´c stosunku

V

=

, najcz˛e´sciej s ˛

a pomijane w

transformacji topo-geo centrum. Podobnie nieistotne s ˛

a relatywistyczne efekty odchylenia kierunku propa-

gacji ´swiatła w polu grawitacyjnym Ziemi.
Słowa kluczowe: geoida, elipsoida obrotowa, współrz˛edne geodezyjnne i geocentryczne na powwierzchni
Ziemi, zenit geocentryczny, geodezyjny i zenit astronomiczny, paralaksa dobowa, aberracja dobowa.

background image

76

Współrz˛edne geocentryczne

N

H

P

Powierzchnie poziome

Linia pionu

Elipsoida

Geoida

Poziom morza

Dno oceanu

Powierzchnia ladow

Rysunek 6.1: Pogl ˛

adowa ilustracja przekroju powierzchni Ziemi. Schematycznie naniesiono frag-

menty powierzchni geoidy i elipsoidy modelujacych kształt Ziemi. Odst˛epstwo geoidy od elip-
soidy jest na tym rysunku mocno przesadzone.

6.1

Współrz˛edne geocentryczne obserwatora

Omówimy transformacj˛e współrz˛ednych z układu odniesienia topocentrycznego do układu geocentrycznego.
Ruch obserwatora wzgl˛edem ´srodka Ziemi jest przyczyn ˛

a dwóch zjawisk: aberracji i paralaksy, nazywanych

w tym przypadku aberracj ˛

a dobow ˛

a, paralaks ˛

a dobow ˛

a b ˛

ad´z paralaks ˛

a geocentryczn ˛

a.

Wyznaczenie poprawki paralaktycznej wymaga znajomo´sci poło˙zenia obserwatora wzgl˛edem ´srodka

Ziemi. Prosty sferyczny model Ziemi nie jest ju˙z dla tego celu wystarczaj ˛

acy. Rzeczywisty kształt bryły

ziemskiej jest bardzo skomplikowany i nie daje si˛e opisa´c prost ˛

a zale˙zno´sci ˛

a funkcyjn ˛

a. Kształt Ziemi defin-

iowany jest w oparciu o ´sredni poziom oceanu, znajduj ˛

acego si˛e w grawitacyjnej równowadze i pokrywa-

j ˛

acego si˛e z powierzchni ˛

a ekwipotencjaln ˛

a obserwowanej siły ci ˛

a˙zenia (patrz rysunek 6.1). W potencjale

pola sił, w którym “zanurzone” s ˛

a punkty bryły ziemskiej, poza grawitacyjnymi uwzgl˛ednia si˛e wyrazy

reprezentuj ˛

ace siły od´srodkowe b˛ed ˛

ace efektem ziemskiej rotacji wokół osi. Tak ˛

a powierzchni˛e ekwipotenc-

jaln ˛

a pokrywaj ˛

ac ˛

a powierzchni˛e oceanu, rozci ˛

agni˛et ˛

a pod masami l ˛

adowymi nazywamy geoid ˛

a.

1

Na mocy

definicji kierunek lokalnej siły ci ˛

a˙zenia jest wsz˛edzie normalny do geoidy.

Powierzchnia geoidy posiada liczne, w porównaniu do jej rozmiarów niewielkie nieregularno´sci, i dlatego
mo˙zna j ˛

a stosunkowo dokładnie przybli˙zy´c dobieraj ˛

ac odpowiedni ˛

a elipsoid˛e obrotow ˛

a o osi obrotu pokry-

waj ˛

acej si˛e z osi ˛

a rotacji Ziemi. Parametry elipsoidy, jej równikowy promie´n

a

oraz spłaszczenie biegunowe

f

, okre´slaj ˛

a tzw. standardowy sferoid wykorzystywany do celów astronomicznych i geodezyjnych. Rysunek

6.2 przedstawia południkowy przekrój standardowego sferoidu. Przekrój jest elips ˛

a o półosiach wielkiej i

małej

a

i

b

odpowiednio, a równanie elipsy ma posta´c

x

2

a

2

+

y

2

a

2

(1

f

)

2

=

1

(6.1)

gdzie

b

=

a(1

f

)

(6.2)

Decyzj ˛

a Mi˛edzynarodowej Unii Astronomicznej (MUA) z 1976 roku, jako standardowy przyj˛eto sferoid o

parametrach:

a

=

6378:140

[km℄

f

=

0:00335281

=

1=298:257

gdzie

a

jest promieniem równikowym,

f

spłaszczeniem ziemskiej elipsoidy.

Niech na rysunku 6.2 obserwator

O

znajduje si˛e na wysoko´sci

h

wzgl˛edem standardowego sferoidu.

2

Nor-

malna do sferiodu w punkcie

O

przebija sferoid w

O

0

, a jej przedłu˙zenie przecina poło˙zon ˛

a w płaszczy´znie

równika o´s

X

w punkcie

Q

. Oczywi´scie odcinek

O

O

0

=

h

.

Dla obserwatora

O

mo˙zna teraz poda´c co najmniej trzy definicje zenitu miejsca obserwacji. Rzeczywisty

kierunek pionu definiuje zenit astronomiczny, nie zaznaczono go na rysunku. Kierunek

QO

definiuje zenit

1

Wzgl˛edem geoidy podaje si˛e tzw. wysoko´s´c nad poziomem morza.

2

Pomijaj ˛

ac niewielkie (do kilku metrów) poprawki geodezyjne, odpowiada to wysoko´sci nad poziomem morza.

background image

6.1 Współrz˛edne geocentryczne obserwatora

77

00

00

11

11

00

00

11

11

00

00

11

11

00

11

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

000000000

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

111111111

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

000000000000

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

111111111111

b

a

Q

O’

O

y

x

ρ

Zenit geodezyjny

Zenit geocentryczny

ν

φ

φ

OO’=h

C

CO=

ρ

Rysunek 6.2: Południkowy przekrój ziemskiej elipsoidy. Spłaszczenie bryły ziemskiej na tym
rysunku jest silnie przesadzone. Nie narysowano kierunku pionu (zenit astronomiczny), bowiem
niekoniecznie musi on le˙ze´c w tej samej płaszczy´znie co zenit geocentryczny i geodezyjny.

geodezyjny, który byłby identyczny z zenitem astronomicznym gdyby geoida ´sci´sle pokrywała si˛e ze sferoi-
dem standardowym. K ˛

at pomi˛edzy tymi dwoma zenitami (efekt istnienia anomalii grawitacyjnej) zwany jest

odchyleniem pionu. Trzeci punkt zenitu to zenit geocentryczny powstały w wyniku przeci˛ecia sfery niebies-
kiej półprost ˛

a

C

O

. K ˛

at



pomi˛edzy kierunkami na zenit geocentryczny i geodezyjny nazywany jest k ˛

atem

wertykału.

W oparciu o ka˙zdy z tych punktów mo˙zna zdefiniowa´c odmienn ˛

a szeroko´s´c obserwatora. Szeroko´s´c jest

k ˛

atem jaki kierunek na zenit tworzy z płaszczyzn ˛

a równika ziemskiego. Zatem, mamy szeroko´s´c geodezyjn ˛

a



i szeroko´s´c geocentryczn ˛

a



0

. Szeroko´sci astronomicznej nie pokazano na rysunku 6.2, bowiem zenit

astronomiczny zwykle nie le˙zy w płaszczy´znie, któr ˛

a ilustruje ten rysunek. Druga współrz˛edna — długo´s´c

geocentryczna i geodezyjna, jak łatwo pokaza´c, s ˛

a sobie równe i tradycyjnie oznaczane przez



.

Niech



b˛edzie geocentryczn ˛

a odległo´sci ˛

a obserwatora,



=

C

O

. Wówczas współrz˛edne geocen-

tryczne

(;



0

;

)

w pełni okre´slaj ˛

a poło˙zenie obserwatora wzgl˛edem ´srodka Ziemi. Dla rozwi ˛

azania wielu

zagadnie´n praktycznych musimy dysponowa´c formułami umo˙zliwiaj ˛

acymi wzajemn ˛

a transformacj˛e pomi˛e-

dzy współrz˛ednymi geocentrycznymi i współrz˛ednymi geodezyjnymi

(;

;

h)

. Niech

(x

0

;

y

0

)

b˛ed ˛

a kartez-

ja´nskimi współrz˛ednymi obserwatora, natomiast

(x;

y

)

współrz˛ednymi punktu

O

0

. Na rysunku ?? mo˙zemy

zauwa˙zy´c, ˙ze



os



0

=

x

0

=

x

+

h

os





sin



0

=

y

0

=

y

+

h

sin



(6.3)

Poniewa˙z punkt

(x;

y

)

le˙zy na elipsie o równaniu (??) a

tan



jest nachyleniem normalnej do elipsy w tym

punkcie, poci ˛

aga to

tan



=

dx

dy

Ró˙zniczkuj ˛

ac równanie (6.1) dostaniemy

y

=

x(1

f

)

2

tan



(6.4)

Kład ˛

ac jego praw ˛

a stron˛e spowrotem do (6.1) b˛edziemy mieli

x

2

(1

+

(1

f

)

2

tan

2

)

=

a

2

Z pomoc ˛

a tego równania oraz równania (6.4) mo˙zna wyrazi´c współrz˛edne

x

i

y

poprzez



. Mianowicie,

x

=

aC

os



y

=

aS

sin



(6.5)

background image

78

Współrz˛edne geocentryczne

Z’

O

C

z

z’

p

r’

r

S

ρ

Rysunek 6.3: Wzgl˛edem miejsca

O

na powierzchni Ziemi obiekt

S

oddalony od geocentrycznego

zenitu

Z

0

o k ˛

at

z

0

. Wzgl˛edem ´srodka

C

Ziemi k ˛

at ten wynosi

z

. Ró˙znica

z

0

z

=

p

nazywana

jest paralaks ˛

a geocentryczn ˛

a (parralaks ˛

a dobow ˛

a).

gdzie

C

=

[ os

2



+

(1

f

)

2

sin

2

℄

1=2

S

=

(1

f

)

2

C

(6.6)

A zatem uwzgl˛edniaj ˛

ac (6.5), kko´ncow ˛

a postaci ˛

a równa´n (6.3) jest



os



0

=

a

os





(C

+

h=a)



sin



0

=

a

sin





(S

+

h=a)

(6.7)

Jest to zale˙zno´s´c dokładna, pozwalaj ˛

aca na obliczenie współrz˛ednych geocentrycznych je´sli tylko dysponu-

jemy współrz˛ednymi geodezyjnymi danego miejsca na powierzchni Ziemi.

6.2

Paralaksa geocentryczna

Warto´s´c paralaksy geocentrycznej zale˙zy od odległo´sci obiektu i jest całkowicie zaniedbywalna dla ciał spoza
Układu Słonecznego. Dla obiektów w pobli˙zu Ziemi jak Ksi˛e˙zyc a zwłaszcza w przypadku sztucznych
satelitów Ziemi, paralaksa osi ˛

aga bardzo du˙ze warto´sci.

Na rysunku 6.3, punkt

O

oznacza obserwatora,

C

´srodek Ziemi a

S

pewne pobliskie ciało niebieskie.

Linia

C

O

, jej przedłu˙zenie okre´sla kierunek na geocentryczny zenit

Z

0

miejsca obserwacji. Płaszczyzna

rysunku jest zdefiniowana przez trzy punkty

C ;

O

i

S

, a zatem le˙zy w płaszczy´znie koła wertykalnego prze-

chodz ˛

acego przez gwiazd˛e. Dlatego przekrój Ziemi pokazany na rysunku niekoniecznie musi przebiega´c

wzdłu˙z ziemskiego południka.

Oznaczmy k ˛

at

Z

0

O

S

przez

z

0

. Jest to obserwowana odległo´s´c zenitalna odniesiona do geocentrycznego

zenitu. Kierunki zenitów geodezyjnego i astronomicznego, ogólnie nie musz ˛

a le˙ze´c w płaszczy´znie rysunku.

Niech

r

0

i

r

b˛ed ˛

a topocentryczn ˛

a i geocentryczn ˛

a odległo´sci ˛

a ´zródła promieniowania,



odległo´sci ˛

a obser-

watora od ´srodka Ziemi (rysunek ??). Paralaks ˛

a geocentryczn ˛

a

p

, nazywamy k ˛

at

O

S

C

taki, ˙ze

z

0

=

z

+

p

(6.8)

gdzie

z

, jest geocentryczn ˛

a odległo´sci ˛

a zenitaln ˛

a obiektu, jak ˛

a obserwowano by ze ´srodka Ziemi. Paralaksa

geocentryczna zwi˛eksza geocentryczn ˛

a odległo´s´c zenitaln ˛

a o k ˛

at

p

a skoro zmiana ta odbywa si˛e w płaszczy´znie

O

C

S

(w płaszczy´znie wertykału), azymut geocentryczny pozostaje niezmieniony.

Stosuj ˛

ac do trójk ˛

ata

O

C

S

wzór sinusów dostaniemy

sin

p

=



r

sin

z

0

=



r

0

sin

z

(6.9)

Wynika st ˛

ad, ˙ze paralaksa dla danego obiektu poza zale˙zno´sci ˛

a od

r

, zale˙zy od jego odległo´sci zenitalnej

a tak˙ze od odlego´sci



obserwatora od ´srodka Ziemi. Potrzebna jest zatem pewna standaryzacja i jest ni ˛

a

background image

6.3 Dygresja. Jednostka astronomiczna — 1 AU

79

tzw. horyzontalna paralaksa równikowa

P

. Jest to paralaksa fikcyjnego obiektu, poło˙zonego na horyzoncie

obserwatora znajduj ˛

acego si˛e na równiku ziemskim

(

=

a)

, czyli w warunkach gdy

(z

0

=

90

Æ

)

. Po

podstawieniu do równania (6.9), mamy

sin

P

=

a

r

(6.10)

Wielko´s´c ta w skrócie zwana paralaks ˛

a horyzontaln ˛

a jest identyczna z odwrotno´sci ˛

a odległo´sci ´zródła promieniowa-

nia. Dlatego w niektórych rocznikach astronomicznych w taki wła´snie sposób stabelaryzowano odległo´sci
Ksi˛e˙zyca. Mamy zatem, ˙ze dla dowolnego obserwatora paralaksa dobowa kierunku do obiektu obser-
wowanego na odległo´sci zenitalnej

z

0

wyra˙za si˛e wzorem

sin

p

=



a

sin

P

sin

z

0

(6.11)

Paralaksy geocentryczne posiadaj ˛

a znacz ˛

ace warto´sci jedynie dla ciał Układu Słonecznego. Jednak warto´sci

paralaks horyzontalnych tych obiektów, ze wzgl˛edu na ruch orbitalny, nie s ˛

a stałe. Np. dla eliptycznej orbity

Ksi˛e˙zyca, jego paralaksa horyzontalna oscyluje pomi˛edzy

54

0

61

0

. W takich przypadkach Mi˛edzynarodowa

Unia Astronomiczna rekomenduje ´srednie warto´sci paralaksy i dla Ksi˛e˙zyca zaleca warto´s´c

P

0

podan ˛

a w

1983

roku przez Murray’a

sin

P

0

=

3422:485



sin

1

00

albo co jest równowa˙zne

P

0

=

57

0

02:

00

6050

(6.12)

Paralaksy geocentryczne planet maj ˛

a wyra´znie mniejsze warto´sci. Dla Saturna wynosi ona około

1

00

, a dla

najbli˙zszej planety, dla Wenus waha si˛e w granicach

5

00

34

00

. Dla obiektu znajduj ˛

acego si˛e w odległo´sci

1

AU

, jego paralaksa nosi nazw˛e paralaksy słonecznej. Poza drobnymi ró˙znicami jest ona bardzo bliska

´sredniej paralaksie prawdziwego Sło´nca.

6.3

Dygresja. Jednostka astronomiczna — 1 AU

Pomiary paralaksy geocentrycznej umo˙zliwiaj ˛

a wyznaczenie odległo´sci pomiedzy ciałami Układu Słonecz-

nego. Jednak obecnie nie jest to ju˙z podstawowy sposób pomiaru odległo´sci bowiem w wielu przypadkacch
zast ˛

apiła go technika radarowa. Ze wzgl˛edów historycznych warto po´swi˛eci´c mu nieco uwagi.

Pomiary pozycyjne planet s ˛

a interpretowane przez mechanik˛e nieba w oparciu o prawa dynamiki graw-

itacyjnej. Je˙zeli w stosunku do masy Sło´nca zaniedbamy masy planet, to mo˙zemy pomin ˛

a´c skomplikowany

opis ruchu planet zast˛epuj ˛

ac go prostym ruchem keplerowskim. Trzecie prawo Keplera powiada wówczas,

˙ze sze´sciany półosi orbit

a

p

planet s ˛

a proporcjonalne do kwadratów ich okresów obiegu

T

p

,

k

2

a

3

p

=

T

2

p

(6.13)

gdzie

k

2

jest stał ˛

a.

Pozycyjne obserwacje w długich interwałach czasu pozwalaj ˛

a dokładnie wyznaczy´c okres orbitalny plan-

ety. Z równania (6.13) daje si˛e wówczas obliczy´c wzgl˛edne rozmiary orbit planet. Aby wyznaczy´c rozmiary
absolutne, potrzeba jeszcze warto´sci stałej

k

2

, ta za´s wyra˙zona jest m.in. poprzez nieznan ˛

a mas˛e Sło´nca. A

zatem, mo˙zna skonstruowa´c model całego Układu Słonecznego w pewnej skali, na to jednak by był to model
np. w kilometrach konieczny jest pomiar obległo´sci cho´cby do jednej planety.

Odległo´sci do planet wyznaczane z obserwacji ich geocentrycznych paralaks, które z natury rzeczy

s ˛

a nieprecyzyjne. Główn ˛

a przyczyn ˛

a s ˛

a tu bardzo małe warto´sci samej paralaksy. Je˙zeli zale˙zy nam na

mo˙zliwie dokładnym pomiarze trzeba dokona´c wyboru odpowiedniego do tego celu obiektu, a wi˛ec obiektu
najbli˙zszego Ziemi czyli o najwi˛ekszej paralaksie. Zagwarantuje to najwy˙zsz ˛

a procentow ˛

a precyzj˛e wyz-

naczonej odległo´sci, co pozwoli dobrze wyskalowa´c rozmiary Układu Słonecznego. Na rysunku ?? na
płaszczy´znie orbity Ziemi narysowano rzuty orbit Wenus Marsa i planetoidy Eros. Zaznaczono równie˙z
poło˙zenia tych obiektów dla dwóch wybranych epok. Punkty

Z

1

;

W

1

;

M

1

oznaczaj ˛

a konfiguracj˛e, w której

Mars znajduje si˛e w opozycji ze Sło´ncem. W par˛e miesi˛ecy pó´zniej mamy konfiguracj˛e

Z

2

;

W

2

;

M

2

, w której

background image

80

Współrz˛edne geocentryczne

00

00

11

11

00

00

11

11

00

00

11

11

0

1

00

11

0

1

00

00

11

11

00

00

11

11

W

Z

Eros

γ

S

W

Z

M

M

1

1

1

2

2

2

Rysunek 6.4: Rzuty orbit planet: Venus, Ziemi, Marsa oraz planetoidy Eros na płaszczy´znie ek-
liptyki.

Wenus jest w konjunkcji dolnej wzgl˛edem Ziemi. Jasne jest, ˙ze dla tych wła´snie konfiguracji mamy najlep-
sz ˛

a okazj˛e do pomiarów paralaksy. Warto jeszcze zauwa˙zy´c, ˙ze ze wzgl˛edu na spory mimo´sród orbity Marsa

(e



0:1)

jego odległo´s´c od Ziemi w momencie kolejnych opozycji zmienia si˛e i to do´s´c wyra´znie. Wenus

posiada bardziej kołow ˛

a orbit˛e dlatego w jej wypadku takich zmian nie obserwujemy. Jednak w momencie

konjunkcji dolnej Wenus jest trudna do obserwacji bowiem przeszkadza tu ´swiatło słoneczne a obserwacja
jest mo˙zliwa jedynie podczas przej´scia planety przez tarcz˛e Sło´nca. Warto´s´c paralaksy daje si˛e wówczas
wyznaczy´c w oparciu o pomiary czasowe zjawiska przej´scia Wenus przez tarcz˛e Sło´nca, obserwowanego z
kilku punktów na powierzchni Ziemi. Niestety g˛esta atmosfera Wenus bardzo utrudnia dokładne pomiary i
była przyczyn ˛

a niepowodzenia kampanii obserwacyjnej podczas przej´scia Wenus przez dysk Sło´nca w ko´ncu

XIX stulecia.

Dlatego wzrosło zainteresowanie obserwacjami Marsa a tak˙ze małej planety Eros, któr ˛

a odkryto pok

koniec XIX wieku. Jak wida´c na rysunku ?? orbita Erosa jest silnie ekscentryczna, sama za´s planetka mo˙ze
znacznie zbli˙zy´c si˛e do Ziemi (na odległo´s´c

0:16

AU

). Dlatego główne programy pomiaru paralaksy geocen-

trycznej w

1901

i

1931

roku po´swi˛econo tej planetce. Obserwacje zako´nczyły si˛e sukcesem co w nast˛epst-

wie doprowadziło do rewizji skali Układu Słonecznego. Wyznaczena w roku

1931

paralaksa Erosa stanowiła

podstaw˛e do wyznaczania warto´sci paralaksy słonecznej a˙z do lat

1960

-tych.

Dopiero technika radarowa umo˙zliwiła bardziej bezpo´srednie pomiary odległo´sci do planet. Stosunek

sygnału do szumu w radarowym echu jest bardzo czuły na odległo´s´c, jest bowiem odwrotnie proporcjonalny
do czwartej pot˛egi odległo´sci obiektu. Wenus była pierwsz ˛

a planet ˛

a, której odległo´s´c wyznaczono metod ˛

a

radiow ˛

a, udało si˛e tego dokona´c podczas konjunkcji w roku

1959

. Te i kolejne pomiary w nast˛epnych latach

doprawadziły do nowego okre´slenia paralaksy słonecznej. W adoptowanym przez MUA systemie stałych [?]
mamy, ˙ze

1

[A

U℄

=

A

=

1:49597870



10

11

[m℄

P

0

=

8:

00

794148

(paralaksa

sone zna)

(6.14)

Wielko´sci te s ˛

a ze sob ˛

a zwi ˛

azane, z równania (6.10) mamy, ˙ze

sin

P

0

=

a=

A

(6.15)

Jednostki astronomicznej nie definiuje si˛e ju˙z jako długo´sci półosi wielkiej orbity Ziemi, poniewa˙z póło´s ta
zmienia si˛e z powodu perturbacji planetarnych. Obecnie definiuje si˛e j ˛

a z pomoc ˛

a teorii grawitacji. Stała

k

z równania (6.13) znana jako stała grawitacji Gaussa, w systemie stałych zalecanych przez MUA jej warto´s´c
wynosi

k

=

0:01720209895

(6.16)

Warto´s´c ta przetrwała ostatnie zmiany systemu stałych i nic nie wskazuje by miało by´c inaczej w najbli˙zszej
przyszło´sci. Przy takiej stałej

k

, jednostk˛e astronomicz ˛

a

(AU

)

definiuje si˛e jako jednostk˛e długo´sci w jakiej

musi by´c wyra˙zona póło´s

a

p

z równania (6.13), gdy okres

T

p

podany jest w dobach.

background image

6.3 Dygresja. Jednostka astronomiczna — 1 AU

81

Warto´s´c paralaksy słonecznej wyznaczona za pomoc ˛

a metod radarowych okre´slona jest z dokładno´sci ˛

a

do jednej mikrosekundy łuku. Jest to poza zasi˛egiem współczesnych metod obserwacji pozycyjnych i dlatego
kolejne zbli˙zenie Erosa do Ziemi wykorzystano jedynie do bada´n jego własno´sci topograficznych. Nawiasem
mówi ˛

ac, dokonano tego technikami radarowymi.

Ze wzgl˛edu na blisk ˛

a odległo´s´c, paralaks˛e Ksi˛e˙zyca daje si˛e okre´sli´c dokładniej ani˙zeli paralaksy planet.

Ale i tutaj techniki radarowe wyeliminowały klasyczne metody optyczne, te za´s po pewnym czasie, zast ˛

api-

ono nowoczesn ˛

a technik ˛

a optyczn ˛

a — technik ˛

a laserow ˛

a. Umo˙zliwiły to misje ksi˛e˙zycowe Apollo, kiedy to

astronauci na pocz ˛

atku lat

1970

-tych umie´scili na Ksi˛e˙zycu odbły´sniki laserowe.

Je´sli chodzi o wyznaczanie odległo´sci paralaksa geocentryczna ma obecnie mniejsze znaczenie, ale nadal

jest istotna jako efekt pozycyjny podczas okre´slania dokładnych poło˙ze´n ciał niebieskich.

Wpływ paralaksy geocentrycznej na wspólrz˛edne równikowe

Powiedziano wcze´sniej, ˙ze paralaksa geocentryczna powi˛eksza geocentryczn ˛

a odległo´s´c zenitaln ˛

a obiektu

nie zmieniaj ˛

ac jego azymutu. Je˙zeli paralaksa jest dostatecznie mała, zmian˛e odległo´sci zenitalnej mo˙zna

wyrazi´c jako

dz

=

z

0

z

=



r

sin

z

(6.17)

A zatem nic nie stoi na przeszkodzie by i w tym wypadku zastosowa´c formuły na małe przesuni˛ecie na sferze
(patrz paragraf 2.6 rozdział 2), podstawiamy zatem:



k

=



r

,



(

0

;

Æ

0

)

identyfikujemy z geocentrycznym zenitem, czyli

Æ

0

=



0

, natomiast

0

=

T

, miejscowy

czas gwiazdowy,



0

=

t

, k ˛

at godzinny obiektu.

Przy takich oznaczeniach, dla obiektu o

( ;

Æ

)

, formuły na małe przesuni˛ecie maj ˛

a posta´c

d

=

0

=



r

os



0

sin

t

se

Æ

=

Æ

0

Æ

=



r

( os



0

os

t

sin

Æ

sin



0

os

Æ

)

(6.18)

S ˛

a to formuły przybli˙zone (pierwszy rz ˛

ad ze wzgl˛edu na

(=r

)

), dlatego nie nale˙zy ich u˙zywa´c w przypadku

Ksi˛e˙zyca praz sztucznych satelitów Ziemi. Nadaj ˛

a si˛e dla pozostałych ciał niebieskich o ile ciała te nie

znajduj ˛

a si˛e w zbyt bliskim s ˛

asiedztwie Ziemi, kiedy to paralaksy geocentryczne przekraczaj ˛

a warto´sci kilku

sekund łuku.

Dla obiektów bardzo dalekich jak planety zewn˛etrzne, paralaksy geocentryczne s ˛

a bardzo małe i dlatego

w formułach (6.18) nie ma potrzeby rozró˙znienia pomi˛edzy szeroko´sciami geodezyjn ˛

a i geocentryczn ˛

a, jako



mo˙zna do wzorów podstawi´c warto´s´c

a

równikowego promienia Ziemi.

Formuły przybli˙zone nale˙zy stosowa´c z rozwag ˛

a, a w sytuacjach w ˛

atpliwych opłaca si˛e stosowanie

rozwi ˛

aza´n dokładnych, które mo˙zna bez trudu wyprowadzi´c. Mianowicie, umie´s´cmy pocz ˛

atek układu odnie-

sienia w ´srodku Ziemi, niech

r

i

R

b˛ed ˛

a wektorami poło˙zenia obserwowanego obiektu

S

i obserwatora

O

.

Wektor

~

O

S

od obserwatora do ´zródła promieniowania dany jest jako ró˙znica

r

0

=

r

R

(6.19)

Rysunek ?? odpowiada wła´snie tej sytuacji. Załó˙zmy, ˙ze wektor

R

jest znany dokładnie, jest to wektor o

długo´sci



skierowany na geocentryczny zenit obserwatora, a zatem we współrz˛ednych równikowych ma on

składowe

R

=

( os



0

os

T

;

os



0

sin

T

;

sin



0

)

(6.20)

gdzie

T

jest miejscowym czasem gwiazdowym w momencie obserwacji. Równanie (6.19) mo˙ze by´c wyko-

rzystane w obie strony, zale˙znie od tego, który wektor jest znany.

Przyjmijmy, ˙ze np. dla Ksi˛e˙zyca, z rocznika astronomicznego na pewien moment czasu zaczerpn˛eli´smy

jego geocentryczne współrz˛edne równikowe

( ;

Æ

)

oraz paralaks˛e horyzontaln ˛

a

P

. Wówczas korzystaj ˛

ac z

równania (6.10) mo˙zemy obliczy´c geocentryczn ˛

a odległo´s´c ksi˛e˙zyca

r

=

a

s

P

background image

82

Współrz˛edne geocentryczne

Geocentryczny wektor poło˙zenia Ksi˛e˙zyca ma składowe

r

=

a

s

P

( os

Æ

os

;

os

Æ

sin

;

sin

Æ

)

(6.21)

Oznaczmy przez

(x

0

;

y

0

;

z

0

)

składowe wektora

r

0

, opisuj ˛

acego obserwowane poło˙zenie Ksi˛e˙zyca. Odpowia-

daj ˛

ace im współrz˛edne sferyczne oznaczymy przez

(

0

;

Æ

0

)

. Wówczas równanie (6.19) w postaci skalarnej

b˛edzie dane jako układ

x

0

=

r

0

os

Æ

0

os

0

=

a

s

P

os

Æ

os



os



0

os

T

y

0

=

r

0

os

Æ

0

sin

0

=

a

s

P

os

Æ

sin



os



0

sin

T

z

0

=

r

0

sin

Æ

0

=

a

s

P

sin

Æ



sin



0

(6.22)

Dysponuj ˛

ac składowymi

(x

0

;

y

0

;

z

0

)

łatwo obliczymy

(

0

;

Æ

0

)

0

=

ar tan (y

0

=x

0

)

Æ

0

=

ar tan(z

0

=

p

(x

02

+

y

02

)

(6.23)

Równanie (6.23) wymaga pewnej ostro˙zno´sci podczas normowania rektascencji do odpowiedniej ´cwiartki.

Przykład.

W stacji o szeroko´sci geodezyjnej

39

Æ

42

0

48

00

dokonano obserwacji sztucznego satelity Ziemi zarówno

radiowo jak i optycznie. Wysoko´s´c stacji wynosi

456

metrów nad poziomem morza. Z obserwacji otrzymano

nast˛epuj ˛

ace równikowe współrz˛edne satelity:

r

0

=

1735:87

k

m

0

=

7

h

12

m

19

s

Æ

0

=

21

Æ

42

0

21

00

T

=

9

h

17

m

34

s

(miejs o

wy

zas

gwiazdo

wy

momen

tu

obserw

a ji)

Oblicz geocentryczne miejsce i odległo´s´c satelity.

Rozwi ˛

azanie wymaga kilku kroków.

1. Nale˙zy obliczy´c składowe wektora

R

poło˙zenia obserwatora na moment obserwacji. W tym celu

przyjmujemy

a

=

6378:14

km,

f

=

3:35281



10

3

i zamieniamy jednostki



=

39:

o

7133

,

T

=

139:

o

3917

.

Mo˙zemy teraz kolejno obliczy´c:

h=a

=

7:15



10

5

C

()

=

1:0013693

(ro

wnanie

(6:6))

S

(f

)

=

0:9946658

(ro

wnanie

(6:6))



os



0

=

4913:459

[km℄

(ro

wnanie

(6:7))



sin



0

=

4053:845

[km℄

(ro

wnanie

(6:7))

R

=

(

3730:183;

3198:095;

4053:845)

[km℄

(ro

wnanie

(6:20))

2. Wykorzystuj ˛

ac wyniki obserwacji obliczymy składowe wektora

r

0

r

0

=

1735:87

[km℄

0

=

108:0792

Æ

=

21:7058

r

0

=

(

500:498;

1533:162;

641:997)

[km℄

(ro

wnania

(6:22))

3. Geocentryczny wektor

r

ma zatem składowe

r

=

r

0

+

R

=

(

4230:681;

4731:257;

3411:849)

[km℄

4. A konwersja do współrz˛ednych sferycznych daje geocentryczne współrz˛edne sferyczne

(r

;

;

Æ

)

:

r

=

7205:843

[km℄

=

8

h

47

m

13

s

ro

wnania

(6:22);

(6:23))

Æ

=

28

Æ

15

0

38

00

background image

6.4 Aberacja dobowa

83

6.4

Aberacja dobowa

Podczas transformacji współrz˛ednych topocentrycznych w geocentryczne, paralaksa jest wa˙zn ˛

a poprawk ˛

a

jedynie dla obiektów z Układu Słonecznego. W przeciwie´nstwie do niej poprawka aberracyjna — aberracja
dobowa — jest niezale˙zna od odległo´sci ´zródła promieniowania i musi by´c uwzgl˛edniona w pozycji ka˙zdego
ciała niebieskiego.

Nie jest to jednak bardzo du˙za poprawka. Poniewa˙z liniowa szybko´s´c obserwatora na równiku ziemskim

stanowi jedynie

1:6



10

6

szybko´sci ´swiatła, przemieszczenie poło˙zenia ciała z powodu aberracji dobowej

nie przekracza

0:

00

33

łuku. Dlatego mo˙zna tu stosowa´c bez zastrze˙ze´n przybli˙zenie małych przesuni˛e´c jak i

pomin ˛

a´c aberracyjne efekty relatywistyczne.

Powtórzmy tu (patrz paragraf 4.5 rozdział 3) nieco zmnienione równania na aberacyjn ˛

a zmian˛e poło˙zenia

´zródła dla obserwatora poruszaj ˛

acego si˛e z szybko´sci ˛

a

V

w kierunku wektora jednostkowego n

ds

=

(V

= )

s



(s



n)

(6.24)

Je˙zeli interesuj ˛

a nas zmiany z tytułu aberracji we współrz˛ednych równikowych w oparciu o formuły na małe

przesuni˛ecie na sferze (2.25) mo˙zemy wyprowadzi´c wzory

d

=

(V

= )

se

Æ

os

Æ

0

sin(

0

)

=

(V

= )( os

Æ

sin

Æ

0

sin

Æ

os

Æ

0

os (

0

))

(6.25)

gdzie

(

0

;

Æ

0

)

s ˛

a równikowymi współrz˛ednymi kierunku jaki wskazuje wektor jednostkowy

n

. Formuły

(6.25) s ˛

a ogólnymi wyra˙zeniami na aberracyjne zmiany rektascencji i deklinacji.

Rozpatrzmy teraz przypadek obserwatora znajduj ˛

acego si˛e na pewnej szeroko´sci i odległo´sci geocen-

trycznej



0

;



. Je´sli

!

oznacza k ˛

atow ˛

a pr˛edko´s´c wirowania Ziemi, liniowa szybko´s´c obserwatora wzgl˛edem

´srodka Ziemi dana jest wzorem

V

=

!

os



0

(6.26)

Ruch dobowy obserwatora oczywi´scie odbywa si˛e w kierunku na wschód co oznacza, ˙ze wektor

n

skierowany

jest na punkt wschodu horyzontu obserwatora, b˛edziemy zatem mieli

0

=

T

+

6

h

Æ

0

=

0

(6.27)

Co po podstawieniu do równania (6.25), oraz zauwa˙zeniu, ˙ze

(

0

)

=

t

daje formuły

d

=

0

=

(!

os



0

= )

se

Æ

os

t

=

Æ

0

Æ

=

(!

os



0

= )

sin

Æ

sin

t

(6.28)

S ˛

a to wystarczaj ˛

aco dokładne formuły dla wszystkich przypadków. Ze wzgl˛edu na niewielki wpływ tej aber-

racji czynione s ˛

a dalsze uproszczenia, mianowicie, pomijana jest niesferyczno´s´c kształtu Ziemi i wówczas

po podstawieniach



=

6378:140

[km℄



=



0

!

=

2

doba

gwiazdo

w

a

=

7:292



10

5

s

1

mamy, ˙ze w jednostkach praktycznych poprawki na aberracj˛e dobow ˛

a wyra˙zaj ˛

a si˛e wzorami

d

=

0:

s

0213

os



se

Æ

os

t

=

0:

00

320

os



sin

Æ

sin

t

(6.29)

Formuły te daj ˛

a ró˙znice pomi˛edzy współrz˛ednymi topocentrycznymi a tymi, które wyznnaczyłby ten sam

obserwator wzgl˛edem ´srodka niewiruj ˛

acej Ziemi. Aby otrzyma´c współrz˛edne odniesione do ´srodka Ziemi

trzeba jeszcze dokona´c transformacji uwzgl˛edniaj ˛

acej wpływ paralaksy geocentrycznej.

Je˙zeli paralaksa i aberracja s ˛

a małe, porz ˛

adek uwzgl˛ednienia poprawek nie jest istotny. W wypadku

du˙zej paralaksy, a natura rei, trzeba aberracj˛e usuwa´c przed zastosowaniem ´scisłych formuł na paralaks˛e
geocentryczn ˛

a.

background image

84

Współrz˛edne geocentryczne

Powiedziano, ˙ze wpływy relatywistyczne mog ˛

a by´c pomini˛ete ze wzgl˛edu na niewielki stosunek

V

=

.

Uwaga ta dotyczyła relatywistycznych efektów w aberracji wynikaj ˛

acych ze szczególnej teorii wzgl˛edno´sci.

Pomini˛ecie relatywistycznego efektu odchylenia ´swiatła w polu grawitacyjnym Ziemi wymaga dalszego

uzasadnienia. Odchylenie to, powiedzmy

Æ

, nie przekracza

2m=

radianów, gdzie

m

jest połow ˛

a Schwarzschildowskiego

promienia Ziemi. Poniewa˙z

m

=

GM



=

2

=

4:4

[mm℄

(6.30)

gdzie

G

jest stał ˛

a grawitacyjn ˛

a, daje to

Æ

=

2m=

<

0:

00

0003

(6.31)

co usprawiedliwia pomini˛ecie wpływu grawitacji na kierunak propagacji ´swiatła w pobli˙zu Ziemi.

6.5

Zadanka na ´cwiczenia

1. Je´sli

a

i

b

s ˛

a równikowym i biegunowym promieniem sferoidy ziemskiej, poka˙z, ˙ze najwi˛eksza

warto´s´c k ˛

ata wertykalnego ma miejsce dla szeroko´sci geodezyjnej

ar tan (a=b)

.

2. Oblicz geocentryczne: odległo´s´c, szeroko´s´c oraz k ˛

at wertykalny dla obserwatora znajduj ˛

acego si˛e na

poziomie morza i szeroko´sci geodezyjnej

52

Æ

.

3. Korzystaj ˛

ac z rezultatów poprzedniego zadania, oblicz maksymaln ˛

a geodezyjn ˛

a wysoko´s´c jak ˛

a mo˙ze

osi ˛

agn ˛

a´c w tym miejscu obserwacji satelita, poruszaj ˛

acy si˛e po kołowej orbicie o promieniu

8798

km,

nachylonej pod k ˛

atem

18

Æ

36

0

do równika.

4. Podaj definicje zenitu astronomicznego, geodezyjmego i geocentrycznego. Wyja´snij co oznacza poj˛e-

cie k ˛

at vertykału ? Udowodnij, ˙ze k ˛

at ten okre´slony jest formuł ˛

a

tan



=

e

2

sin

2

2(1

e

2

sin

2

)

gdzie

e

jest mimo´srodem standardowego ziemskiego sferoidu,



jest szeroko´sci ˛

a geodezyjn ˛

a.

5. Atmosferyczna refrakcja oraz paralaksa, obie zmieniaj ˛

a odległo´sci zenitalne ciał. Wzgl˛edem jakich

punktów zenitu zmiany te s ˛

a okre´slone?


Wyszukiwarka

Podobne podstrony:
311[10] Z1 06 Stosowanie rachunku współrzędnych w obliczeniach geodezyjnych
MT st w 06
Kosci, kregoslup 28[1][1][1] 10 06 dla studentow
06 Kwestia potencjalności Aid 6191 ppt
06 Podstawy syntezy polimerówid 6357 ppt
06
06 Psych zaburz z somatoformiczne i dysocjacyjne
GbpUsd analysis for July 06 Part 1
Probl inter i kard 06'03
06 K6Z4
06 pamięć proceduralna schematy, skrypty, ramyid 6150 ppt
Sys Inf 03 Manning w 06

więcej podobnych podstron