základy astronomie a astrofyziky rotace Země, precese, nutace, aberace, refrakce 4.1 Rozměry Země První historická měření rozměru Země vycházela z Eratosthénovy metody, která svazovala úhel mezi dvěma body na povrchu a vzdálenost, která je mezi nimi. Takto určili již Arabové délku 1° mezi 111,7-113,3 km. První novodobá měření podobnou metodou provedl v roce 1528 J. Fernel (1497—1558), který pomocí otáček kola na povrchu stanovil délku tzv. Pařížského poledníku mezi Paříží a Amiens. Od této metody se postupně přešlo k přesnější metodě triangulační. Na povrchu Země se vytýčily body, pro které platilo, že nejbližší mezi sebou byly viditelné. Měřením úhlů a vzdáleností mezi nimi se pak postupně proměřovaly jednotlivé trojúhelníky, které měly strany cca. mezi 30-40 km (příklad je na obrázku 56). Od roku 1669 proběhlo měření napříč celou Francii, bylo zjištěno, že 1° měří různě. Vysvětlení je prosté, Země není ideální koule, ale směrem k pólům se více projeví její zploštění. Teoretické odvození provedl Isaac Newton kolem roku 1686, příčina zploštění Země tkví v její rotaci. Další měření provedená v Laponskú a v Peru pak toto tvrzení jen potvrdila. Na základě těchto měření pak byla definována jednotka 1 m jako 1/10000000 poledníkového kvadrantu Země. Mezi roky 1733-1740 Jacques (1677—1756) a jeho syn César (1714—1784) Cas-siniové uskutečnili první triangulaci celé Francie včetně přepočtení délky jednoho stupně meridiánu, což vedlo v roce 1745 k publikaci první mapy Francie na rigorózních základech. Ve Francii byl také metr uznán jako oficiální délková míra v roce 1799, v Rakou-sko-Uhersku se tak stalo až v roce 1876. Prototyp metru je umístěn v Sévres u Paříže. V současné době je metr definován pevně stanovenou hodnotou rychlosti světla ve vakuu, kterou toto světlo urazí během časového intervalu 1/299792458 sekundy. Dnešní geodetická měření se provádějí pomocí družic Země. Měření jsou mnohonásobně přesnější než pomocí trigonometrické sítě, získáváme velmi podrobné Fig. 4. Die rheinifch-hefTifche Kelte und das nieder-rh ein i ich e Dreiecksnetz. Obrázek 56: Historická mapa Německa z 19. století se znázorněnými triangulačními body v oblasti Porýní-Hesenska [E21]. informace o tvaru, rozměrech Země a také o rozložení hmoty v ní. To vseje možné díky přímým pozorováním a měřením vzdáleností pomocí laserových odražečů ä rádiových měření, relativních rychlostí určených z Dopplerova posuvu, interfero-metrických měření a altimetrií z družic. Přesnou definicí trajektorie družic a pohybu po ní nám umožňují speciální geodetické družice, které jsou schopny potlačit negravitační poruchy, což nám pak umožňuje zjistit přesný potenciál Země. 4.2 TVar Země Snaha o popsání tvaru Země naráží na složitou aproximaci velmi členitého zemského povrchu tělesem jednodušších vlastností. Realizace ekvipotenciální plochy nejtěsněji přiléhající ke střední klidné hladině moří a oceánů, ke které se vztahuje tzv. nadmořská výška, je tvarem velice komplikovaným pro praxi. Tato plocha reprezentuje tzv. geoid (obr. 57). Pro praktické potřeby se zavádí zjednodušení, kterého se využívá v astronomii a kartografii. Prvním jednodušším modelem tvaru Země je rotační elipsoid, který lze ve válcových souřadnicích se středem v centru Země popsat vztahem kde p — sjx1 + y2 je vzdálenost od osy otáčení, a je rovníkový poloměr a & poloměr EIGEN-CG01C Geoid Obrázek 57: Geoid telné [E22]. tvar Země, výškové rozdíly jsou schválně zvětšeny, aby byly vidi- polární. Zploštění i je definováno jako poměr i = (193) a Referenční elipsoid, který slouží jako (co nejlepší) náhrada geoidu, vznikl na základě družicových měření. Současným globálním standardem mezi referenčními elipsoidy je tzv. WGS-84, jehož střed leží ve středu Země a poloosy jsou a = 6378137 m, b = 6356752,3 m. Zploštění je i=l:298,257. Tento elipsoid se používá např. při satelitní navigaci GPS. Ostatní referenční elipsoidy lze definovat pomocí posunutí jejich středu vůči středu WGS-84, případně i změnou délky poloos (změnou délky hlavní polosy a rozdílem zploštění). Zploštění Země není pouhým okem rozpoznatelné a proto má Země z vesmíru tvar dokonalé koule. Přesnějším vystižením tvaru geoidu je trojosý rotační elipsoid popsaný rovnicí (Bursa 1995) s rovníkovým zploštěním a — b 1, (194) (195) Polární zploštění je definováno obdobně jako ! p (196) Parametry elipsoidu jsou a = 6378173 m, iR = 1 : 94000, íP = 1 : 297,787, přitom hlavni osa x protíná rovník na zeměpisné délce Xa = 14,8°. Tomuto modelu se nejvíce přibližuje rotační elipsoid s parametry a = 6378139 m a i = 1 : 297,257. Skutečný tvar Země je znázorněn na obrázku 57, ze kterého je patrno, že se severní polokoule liší od jižní, tvarem připomíná hrušku. Odchylky od ideálního tvaru jsou cca. 50 metrů. Zeje Země přibližně kulatého tvaru lze dokázat jednak pozorováním jejího stínu při zatměních Měsíce, různou výškou Polárky nad obzorem v závislosti na zeměpisné šířce nebo přímým pozorováním z družic. 4.3 Dohlednost a deprese obzoru Stojíme-li na zemském povrchu, čím se nacházíme ve větší nadmořské výšce, tím se nám vzdálenost, do které vidíme objekty na Zemi zvětšuje. Je to dáno jejím zakřivením, jak je znázorněno na obrázku 58. X 1 geodetický V'V obzor yf R / \X / \ x ■ Obrázek 58: Geodetický obzor pro výsku pozorovatele h nad povrchem Země. Pro vzdálenost geodetického obzoru vyjdeme z obrázku 58, ze kterého platí, že LQ= [(R + hy-Ri]1/2 = a pro tangentu pak (a je tzv. úhel deprese) 1/2 = (2/v7?)1/2. (198) 2Rh 1 2R V2Kh (197) tan a ■■ (R+h)2-R2 R2 Pro případ, že by Země byla bez atmosféry, pak vychází pro výšku h — 10 m úhel deprese a = 6' 5" a vzdálenost geodetického obzoru L0 = 11,29 km. Vlivem refrakce je však skutečná vzdálenost L geodetického obzoru větší, L > L0. Dohlednost za normálních podmínek je větší přibližně o 6,5 %. Hodně však záleží na teplotním rozvrstvení přízemní atmosféry, která má na refrakci velký vliv. 4.4 Šířka geocentrická, geodetická a astronomická Vzhledem k tomu, že je tvar Země odlišný od koule a popisujeme jej pomocí rotačního elipsoidu, rozlišujeme různé typy zeměpisných šířek. Šířka geocentrická
+-(-) sin20. (216) Pro Zemi je poměr v/c roven přibližně 1/10000, což dává aberační úhel cca 20,5". V okolí pólu ekliptiky budeme během roku pozorovat aberační kružnici o poloměru a = 20,49552", budeme-li se blížit rovině ekliptiky, kružnice přejde do elipsy, jejíž malá poloosa bude mít úhlovou velikost b = 20,49552" sin/3. V rovině ekliptiky pak budeme pozorovat pouze aberační přímku. Vektor v se mění během roku, tím se mění i velikost a směr aberace, která leží na hlavní kručnici procházející tělesem a apexem. Složky hvězdné aberace můžeme rozdělit na a) denní - souvisí s rotací Země (v ~ 0,46km/s, b) roční - pohyb Země kolem Slunce {v ~ 30 km/s), (214) (215) c) pohyb Země vůči těžišti soustavy Země-Měsíc zanedbáváme (v ~ 0, 01 km/s), d) sekulární - pohyb těžiště sluneční soustavy prostorem (v Galaxii). V praxi sekulární aberaci nelze odlišit od opravy na šíření světla, tuto aberaci neurčujeme, je součástí tzv. střední polohy hvězd. Korekce souřadnic pro denní aberaci můžeme vypočítat podle vztahů Aa — 0, 0213s cosip cosi sec S . . AS = 0,320" sin
je geocentrická a
c = p/acosip = (C + h/a)cosíp, ^ ' kde h je nadmořská výška v metrech v místě pozorování. Parametry S a, C jsou dány vztahy (Seidelmann 2006) S (1 - if \jcos2 tp + (1 — i)2 sin2 ij? Ze vztahů (202) a (193) získáme relaci e2 = což po dosazení do rovnic (221) dává 1 _ c2 a C- cos2 ^> ■ ■ (1 — i)2 sin ■ 2z ■ (221) (222) e2 sin2
0 cos 0 v bodě 1 a A«r — v® cos0 v bodě 3. Vůči radiální heliocentrické rychlosti hvězdy vl0 tak variace radiální rychlosti způsobí rozkmit, který je roven 2-00 005/3. Rovněž rozdíl časů, kdy přijde signál z pozorovaného objektu k Zemi a Slunci, se mění, pouze v bodech 1 a 3 a v pólech ekliptiky, Slunci a hvězdě je heliocentrická korekce rovna nule. Maximální hodnoty korekce jsou pak v bodech 2 a 4 a jejich velikosti jsou Aí = —t cos 0 v bodě 2 a Aí = — t cos 0 pro bod 4. Velikost heliocentrické korekce může být maximálně rovna rcos/3 . Ar =-= t cos p c kde t = 499 s je oprava na šíření světla. Existence paralaxy je jedním z hlavních důkazů pro platnost heliocentrického systému. Pokusy o její změření selhávaly ať již z důvodu neexistence vhodné metody či z důvodů nevhodných přístrojů. Tycho Brahe nezměřil paralaxu, což jej vedlo k podpoře geocentrické soustavy, stejně tak neuspěl ani Koperník. Problémem bylo, že jejich měření byla absolutní. Galileo Galilei i Christian Huyghens navrhli, že je lepší použít měření pozic relativních vůči hvězdám v pozadí, které mohou posloužit jako opěrné body. Tím se zvýší přesnost astrometrie a navíc se bližší hvězdy projeví změnou své polohy vůči vzdálenějším hvězdám na pozadí. Ale ani to nepřineslo úspěch, i když se o heliocentrickém modelu již nepochybovalo. Příčinou byla velká vzdálenost hvězd, u kterých je paralaxa menší než 7r -c 1". V roce 1838 Friedrich Bessel (1784-1846) určil paralaxu hvězdy 61 Cygni na tt = 0, 314", dnes udávaná hodnota je 7r = 0,286". Souběžně určil paralaxu Vegy Friedrich Ge-org Wilhelm von Struve (1793-1864) a Thomas James Henderson (1798-1844) pak vzdálenost Tolimanu (alfa Centauri). Měření vzdálenějších hvězd je obtížnější, paralaktické úhly jsou velmi malé, spolehlivost byla do 20-30 pc. V této vzdálenosti se nachází ale málo důležitých a opěrných typů hvězd (ty jsou většinou hodně jasné a ve velkých vzdálenostech, v prostoru málo zastoupeny), což vede ke složitějšímu způsobu navazování vzdáleností. Velkým přínosem byla astrometrická družice HIPPARCOS (High Precision PARallax COllecting Satellite), která vystartovala v roce 1989 a po 4 letech změřila polohy a vlastní pohyby více než 100 tisíc hvězd s přesností 0,002 ". Všechny výše vzpomínané efekty (aberace, paralaxa, variace radiální rychlosti, heliocentrická korekce i varicace v počtu meteorů) nám dávají důkazy o pohybu Země kolem těžiště sluneční soustavy. 4.16 Důsledky pohybu Země kolem Slunce Země vzhledem ke hvězdám oběhne kolem Slunce za rok siderický, který je dlouhý 365,25636 dne. Rok tropický je dán dvěma po sobě následujícími průchody Slunce jarním bodem, tj. dobou od jarní do jarní rovnodennosti, trvá o něco méně než rok siderický a to 365,24219 dne. Tento rozdíl souvisí s tím, že se jarní bod pomalu posouvá po ekliptice proti pohybu Slunce. Do původní polohy se opětovně dostane až za necelých 26000 let a tato perioda se nazývá rokem Platónským. Posun jarního bodu souvisí s precesí zemské osy. Doba mezi dvěma průchody Země perihéliem se nazývá rok anomalistický, který trvá 365,25964 dne. Přímka apsid se stáčí ve směru pohybu Slunce po ekliptice, proto je anomalistický rok delší a do stejné pozice vůči hvězdám se přímka apsid vrátí po 110000 rocích. Velice důležitým rokem je rok drakonický, který je definován dobou mezi dvěma průchody Slunce výstupným uzlem měsíční dráhy. Uzlová přímka Měsíce se stáčí proti pohybu Slunce po ekliptice, proto je délka drakonického roku rovna 346,62006 dnům, do stejného bodu v prostoru se pak dostane za 19,6 let. (227) Oběh Země kolem Slunce se rovněž projeví ve střídání ročních období. Ani sklon rotační osy Země vůči rovině ekliptiky není konstantní a mění se v čase dle e(ť) = 23°26'21,4"- 0,4684(í- 2000), (228) kde t je letopočet. Sluneční paprsky dopadají během roku na obě polokoule pod měnícím se úhlem, který souvisí se změnou deklinace Slunce na hvězdné obloze. Tato změna má pak za příčinu různé teploty, které souvisejí s tzv. středními klimatickými pásy. Ty vznikají tak, že jsou v rozdílných zeměpisných šířkách a) různé délky dne a noci, b) rozdílné výšky Slunce nad obzorem během dne, c) různě zmírňovány teplotní setrvačností a d) proměnlivou vzdáleností Země a Slunce během roku. S tím souvisí množství energie, které od Slunce dopadá na zemský povrchu na jednotku plochy. Toto množství si můžeme vyjádřit jako K K W = —-cos z = — sin/i, (229) kde K je sluneční konstanta, která je pro vzdálenost 1 astronomické jednotky rovna 1360 W/m2, r je vzdálenost od Slunce v astronomických jednotkách, z je zenitová vzdálenost a h pak výška Slunce nad obzorem. Insolace (oslunění) je pak energie, která dopadá za jeden den na 1 m2 plochy na horní hranici zemské atmosféry a je dána vztahem rt=h K £ — i — sin h at Jt=-t =-E0 T sin/i = sin sin í + cos
1, Slunce nezapadá a t0 = ir. Naopak pro tan
1 Slunce nevychází a to = 0.
Astronomická roční období jsou vzhledem k různé rychlosti Země v dráze různě dlouhá. Nejdelší je astronomické léto (pro severní polokouli), nejkratší pak astronomická zima. Letní půlrok je o 7,5 dne delší než zimní. Klimaticky je tak podnebí na severní polokouli mírnější, vlivem precese se však za 13000 let polokoule prohodí a bude tomu naopak. Celkové roční úhrny insolace jsou uvedeny v tabulce 7. Rozdíly nejsou až tak extrémní, je možné osídlit i vyšší zeměpisné šířky zejména v období krátkého léta.
Parametry zemské trajektorie nejsou konstantní a mění se v čase. Sklon zemské osy kolísá mezi e — 22°04ř a e — 24°34' s periodou 41000 let. Při rostoucím sklonu
Obrázek 66: Vliv apsidálního pohybu (změna délky perihélia) na změnu ročních období (upraveno dle fE29]).
rotační osy bude i změna deklinace Slunce větší, takže léta pak budou teplejší (a zimy chladnější). Pro zjištění sklonu ekliptiky můžeme použít vztahu
e = 23°26ř21, 448" - 46, 8150fT - 0, 00059řT2 + 0,001813"T3, (232)
kde T = {JD — 2451545,0)/36525. Změna délky perihélia II se děje s periodou 21000 let a v této periodě se Země dostává nejbllže Slunci v různých částech roku. Je-li II = 0°, nastává jarní rovnodennost v přísluní, II = 90° je Země v pHsluní v době zimního slunovratu (kratší a teplejší zima, delší a chladnější léto), II = 180° je podzimní rovnodennost v přísluní a II = 270° je Země v přísluní v době letního slunovratu (kratší a teplejší léto a delší a chladnější zima). S touto periodou se vyměňují role severní a jižní polokoule. Změny ve výstřednosti dráhy Země jsou s periodou cca 100000 let, excentricita se mění v rozmezí 0, 0007 -f- 0,0658. V současnosti je rovna 0,0165. Při vyšší excentricite trajektorie je Země větší část oběhu ve větší vzdálenosti, což má za následek doby ledové. Délka velké poloosy se mění zcela nepatrně.
Všechny tyto periody se snaží popsat Milankovičova teorie (Milutin Milankovic (1879-1958)), která předpokládá, že změny klimatu Země souvisejí s periodami změn výstřednosti zemské trajektorie, délkou perihélia, sklonem zemské osy a její precesí (viz obr. 67).
MilankOvičovy Cykly
i-1-1-1-1-1-1-1-1-1-1-1-1-1-r
-íůí -700 -600 -500 -400 -300 -200 -100 0 100 200 300 400 500 KŮD 7W> $«3
v tisících letech
Obrázek 67: Minulé a budoucí Milankovičovy cykly. Nahoře je znázorněna změna sklonu zemské osy, pod ní pak změna excentricity zemské trajektorie, sinus délky perihélia a násobek excentricity a sinu délky perihélia. černou barvou je pak závislost střední hodnota insolace pro zeměpisnou šířku 65 stupňů pro horní část atmosféry. Poslední dvě křivky dole znázorňují změny zastoupení izotopu kyslíku lsO v ukládaných sedimentech a změny teploty na růstu antarktického ledu ze stanice Vostok [E30].
Dosazením obou efektů získáme aproximaci s přesností na cca. 0,3 minuty. Časová rovnice pak nabude tvaru
TPS - Tss = 0,46 cos t — 7,34 sin t - 3,36 cos 2t - 9, 33 sin 2r, (248)
kde t = 3gg°°äD a rozdíl pravého a středního slunečního času vychází v minutách. D je opětovně počet dní od začátku roku.
4.19 Precese a nutace
Precesi objevil kolem roku 130 př.n.l. Hipparchos, který porovnával svá pozorování s pozorováním řeckých astronomů. Zjistil přírůstek ekliptikálních délek, které vysvětlil posunem jarního bodu proti směru pohybu Slunce. Ptolemaios precesi změřil na 36" za rok. Arabští astronomové mezi 10. a 11. stoletím precesi určili mezi 48"-54" za rok. V roce 1260 perský hvězdář Nassir Edin pak určil velice přesnou hodnotu precese 51" za rok, která je velice blízká současné hodnotě 50,3"/rok. Od dob Hipparcha se jarní bod posunul zhruba o 30 stupňů.
Fyzikální vysvětlení precese podal až Isaac Newton, který popsal rotaci Země jako chování roztočeného setrvačníku s rotační osou, na kterou působí vnější síly a pokud mají vůči tělesu určitý moment, pak vyvolávají precesi (viz obr. 72). Moment sil D působených Měsícem a Sluncem na zploštělou Zemi chce osu Země
y
Měsíc
Obrázek 72: Výsledný moment sil D, který vzniká působením dvojice sil Měsíce a Slunce za zploštělou Zemi.
narovnat vůči rovině ekliptiky. Výsledkem tak je pohyb zemské osy okolo pólu ekliptiky. Velikost momentu sil je dán vektorovým součinem
Ď = C [ď x