Astronomy - The Solar System

We begin with a tour (in visible light) of some of the well-photographed bodies in our Solar System (all photographs are from The Planetary PhotoJournal at the Jet Propulsion Laboratory, although some have been cropped and/or color-corrected). We've seen many of these already, but it is instructive to see them in the context of their brethren. Both the planets and their moons are ordered with increasing orbital distance; asteroids and comets appear alphabetically:

MercuryGaspra SaturnIapetus Comet Wild 2
VenusIda MimasUranus
EarthMathilde EnceladusMiranda
Earth's Moon (Luna)Jupiter TethysAriel
MarsIo DioneTitania
PhobosEuropa RheaNeptune
DeimosGanymede TitanTriton
ErosCallisto HyperionComet Tempel 1

Pay attention to atmospheric features such as density, coloration and movement.

Yellow and orange coloration are largely due to sulfur compounds, although in Titan's case its coloration comes from a heavy hydrocarbon haze. Jupiter's ammonia content also affects its color, and can be seen in the greenish tints reflected from the night sides of the Galilean moons. Methane absorbs light in the red and yellow wavelengths, leading to Uranus' (with less methane) and Neptune's (with more methane) coloration.
Also look closely for surface features which might give clues to past tectonic, volcanic (source), geyser or meteor activity. Distinguish between volcanic and meteor craters by looking for impact debris and ejecta (source) (some craters could be calderas). Look for volcanic rilles (trenches) (source), lava flows (source) and domes (source), as well as ridges and scarps (cliffs) formed by cooling and shrinking. Note the presence of liquids (source), dust or ice, and evidence of major collisions. What do these tell you about the surface history?
To get a feel for how difficult this business is, read this analysis of Europa.
Here is a radar map of Venus, and one of Titan, to help you identify surface features. We also include one of Mars because our other image was chosen to highlight atmospheric and polar features.

This image shows the planets to scale:

(source), and this one shows some of the moons to scale:

(source). This image shows the orbits of the outer planets:

(source).

Have you ever seen the far side of the Moon? (source)

(View GOES 2/26/98 eclipse images (source).)

(View Cosmos DVD 4, episode 6, on Voyager at Europa.)

(View Cosmos DVD 3, episode 4, comet simulations, on the Tunguska Event.)

(View Solar System Animations from the Minor Planet Center.)

Here is a panorama of the meteor crater near Winslow, AZ, and a movie of the Peekskill fireball (source).

Some of the smallest constituents of the Solar System are dust particles.


Portfolio Exercise 4: Categorize each body in the list above as to the distinguishing feature(s) of its landscape: tectonic, volcanic, geyser or meteoric. You may have some bodies in which multiple features are present.

For any characterization other than meteoric, find NASA images which support your categories; include them (with source references) in your portfolio. Annotate the pictures, circling and identifying the features which support your categorization.


Solar System Observation

With careful telescopic observation, we can learn a number of important parameters about some of these bodies in our Solar System:
BodyPerihelionAphelionAng. Diam. @OrbitalOrbitalRotationalAxial Tilt
(AU)(AU)ca. (as)Period (Yrs)InclinationPeriod (Days)
Sun19182.2 * 10825.387.25
Mercury0.30750.466713.030.241758.650
Venus0.7180.72865.440.6153.39-243177.3
Earth0.9831.01715. * 10-50.997323.45
The Moon0.002430.0027119680.0755.14527.326.68
Mars1.3811.66625.731.8811.81.02625.19
Phobos6.175 * 10-56.363 * 10-50.10478.731 * 10-411.026
Deimos1.567 * 10-41.569 * 10-40.058240.0034561.81.026
Gaspra0.2917
Ida0.193
Jupiter4.9525.45550.111.861.3050.41353.12
Io0.002810.0028321.27610.0048440.041.769
Europa0.0044450.0045261.0970.0097230.473.551
Ganymede0.0071380.0071671.8460.019590.217.155
Callisto0.01250.012681.6840.045690.5116.69
Saturn9.02110.0520.7629.42.4840.44426.73
Mimas0.0012150.0012650.067530.002581.530.9424
Enceladus0.0015840.0015980.08510.00375101.37
Tethys0.001970.001970.18260.0051691.861.888
Dione0.0025170.0025280.19290.0074930.022.737
Rhea0.003520.0035270.26320.012370.354.518
Titan0.0079290.0084060.88710.043660.3315.95
Hyperion0.008870.010930.030620.05830.4313(chaotic)
Iaepetus0.023130.024480.24740.217214.7279.33
Uranus18.2920.14.08684.020.77-0.719697.86
Miranda8.658 * 10-48.705 * 10-40.037690.003874.21.413
Ariel0.0012720.0012810.092460.00690.32.52
Titania0.002910.0029230.1260.023840.148.706
Neptune29.8130.332.372164.81.7690.671229.58
Triton0.0023710.0023710.1295-0.01609157.3-5.877

("ca." stands for "closest approach".)

Since we telescopically measure only angular positions in the sky, we cannot measure perihelion and aphelion directly; we require an independent measurement of distance. This is done, for example, by radar ranging of Venus at closest approach. Thus the perihelion and aphelion values are actually computed and not directly observed. Planetary orbits are elliptical, and orbital parameters such as perihelion and aphelion are obtained by fitting careful measurements to elliptical orbits.

The data in the tables in this section come from NASA. There is an Excel spreadsheet available containing some of this information. It will serve as a starting point for your efforts to duplicate the conclusions below.


Orbits

An ellipse is parametrized by a pair of lengths. In the definition of an ellipse, one is the length between two reference points (called "foci"), and the other defines the size of the ellipse: for every point on the ellipse, the sum of the distances between that point and the two foci is a constant. We can trade these lengths for two of more physical interest: the perihelion and aphelion. But we often prefer a different pair of values: the length of the semimajor axis (denoted "a") and the eccentricity (denoted "e"), which is the ratio of the distance between the foci to the semimajor axis length:

In this plot, the sum of the lengths of the two red lines is a constant for every point on the ellipse.
Using the above information and the equations

and

we can compute the following parameters:

Bodya(AU)EccentricityR(m)R(Earth)
Sun1.62 * 1096.955 * 108109
Mercury0.38710.20562.44 * 1060.3825
Venus0.72330.00686.052 * 1060.9488
Earth10.016716.378 * 1061
The Moon0.002570.05491.734 * 1060.2719
Mars1.5240.09343.397 * 1060.5326
Phobos6.269 * 10-50.0151.3816 * 1040.002166
Deimos1.568 * 10-45. * 10-476880.001205
Gaspra84280.001321
Ida1.57 * 1040.002462
Mathilde2.971 * 1040.004658
Jupiter5.2030.048397.149 * 10711.21
Io0.0028210.0041.821 * 1060.2855
Europa0.0044850.0091.565 * 1060.2454
Ganymede0.0071530.0022.634 * 1060.413
Callisto0.012590.0072.403 * 1060.3768
Saturn9.5370.054156.027 * 1079.449
Mimas0.001240.02021.96 * 1050.03073
Enceladus0.0015910.004522.47 * 1050.03873
Tethys0.0019705.3 * 1050.0831
Dione0.0025230.002235.6 * 1050.0878
Rhea0.0035230.0017.64 * 1050.1198
Titan0.0081670.029192.575 * 1060.4037
Hyperion0.00990.1048.878 * 1040.01392
Iaepetus0.023810.028287.18 * 1050.1126
Uranus19.190.047172.559 * 1074.012
Miranda8.681 * 10-40.00272.36 * 1050.037
Ariel0.0012760.00345.79 * 1050.09078
Titania0.0029160.00227.889 * 1050.1237
Neptune30.070.008592.476 * 1073.883
Triton0.0023711.6 * 10-51.352 * 1060.212
Tempel 131034.865 * 10-4
Wild 227504.312 * 10-4

(Here and in the following, "(Earth)" means the values are multiples of that value for the Earth.)

Of course, a and e are sufficient to describe any one of these orbits, but in order to understand how the planets are oriented with respect to each other we need 4 additional values:

Together these 6 values are called orbital elements. From them we can predict the past and future positions of any of the planets, to reasonable accuracy, within about 20 years of the time mentioned above. Beyond those dates, gravitational interactions between the planets must be taken into account.

There are subtle variations in eccentricity and tilt which cause long term cycles in the amount of sunlight received in the northern hemisphere, and which drive cyclic climate change:

These cycles were discovered by Milutin Milankovitch, a Serbian mathematician, and account for the ice ages which occur in 100000 and 41000 year cycles, as well as the smaller variations that occur in 19000 to 23000 year cycles (source).

They do not account for the increases in average global temperature since 1880.


Portfolio Exercise 5: Compute the following characteristics for each body in the table above: volume (in units of Earth's volume), surface area (in units of Earth's surface area) and average orbital velocity (in km/hour). Include details on the formulas used to compute the radius of each body.

Mass and Motion

Using Newton's version of Kepler's Third Law:
orbital period = 2 * p * (semimajor axis3 / (G * mass))1/2
we can compute the mass of any body with an orbiting satellite using the satellite's orbital parameters. But since we are able to send probe satellites to these bodies, we are able to measure their masses more accurately by measuring their gravitational influence on the paths of our probes and using Newton's Equation:
force due to gravity = G * mass1 * mass2 / distance2
The relationship between Newton's Law of Gravity and Kepler's Laws of Motion is discussed in a Mathematica Notebook. You need Mathematica or Mathematica Player to view Mathematica Notebooks.

Gauss' Law tells us that the gravitational force due to a spherically symmetric mass is equivalent to the force due to a point particle of the same mass, located at the center. Therefore we will usually treat planets (and stars!) as if all their mass was located at a point at their center.

Using these masses, we can compute a number of interesting parameters: their average densities, their surface gravities (by setting "mass2" to 1 and "distance" to the radius), their orbital angular momenta ("L", using the equation

angular momentum = 2 * p * mass * rotation rate * orbital radius2)
and their central pressures. The last we can estimate from dimensional considerations, using the formula
pressure = G * mass2 / radius4
Note that this value of the Earth's central pressure is 14.23 million times the atmospheric pressure on the surface.
BodyM(kg)M(Earth)Densityg(Earth)L (kgm2/s) L(Earth)Central Pres.Central Pres.
(g/cm3)(Pa)(Earth)
Sun1.989 * 10303.33 * 1051.411281.424 * 1056 5.349 * 10151.128 * 1015784.1
Mercury3.302 * 10230.055285.4290.37789.149 * 1038 0.034372.054 * 10110.1427
Venus4.868 * 10240.8155.2440.90531.846 * 1040 0.69341.179 * 10120.8195
Earth5.974 * 102415.49612.662 * 1040 11.439 * 10121
The Moon7.348 * 10220.01233.3650.16642.884 * 1034 1.083 * 10-63.986 * 10100.0277
Mars6.418 * 10230.10743.9090.37883.531 * 1039 0.13262.064 * 10110.1435
Phobos1.063 * 10161.779 * 10-90.96233.792 * 10-42.132 * 1026 8.009 * 10-152.07 * 1051.438 * 10-7
Deimos2.38 * 10153.984 * 10-101.252.742 * 10-47.545 * 1025 2.835 * 10-151.082 * 1057.52 * 10-8
Jupiter1.899 * 1027317.81.242.531.932 * 1043 725.89.209 * 10126.4
Io8.932 * 10220.014953.5310.18346.538 * 1035 2.456 * 10-54.841 * 10100.03364
Europa4.8 * 10220.0080352.9890.13354.426 * 1035 1.663 * 10-52.563 * 10100.01781
Ganymede1.482 * 10230.024811.9360.14551.725 * 1036 6.479 * 10-53.044 * 10100.02116
Callisto1.076 * 10230.018011.8510.12691.662 * 1036 6.245 * 10-52.317 * 10100.0161
Saturn5.685 * 102695.170.621.0667.837 * 1042 294.41.635 * 10121.136
Mimas3.75 * 10196.278 * 10-61.1890.0066489.96 * 1031 3.742 * 10-96.359 * 1074.419 * 10-5
Enceladus7. * 10191.172 * 10-51.1090.0078142.105 * 1032 7.907 * 10-98.785 * 1076.105 * 10-5
Tethys6.27 * 10201.05 * 10-41.0050.01522.097 * 1033 7.879 * 10-83.325 * 1082.311 * 10-4
Dione1.1 * 10211.841 * 10-41.4950.023894.163 * 1033 1.564 * 10-78.21 * 1085.706 * 10-4
Rhea2.31 * 10213.867 * 10-41.2370.026951.033 * 1034 3.881 * 10-71.045 * 1097.263 * 10-4
Titan1.346 * 10230.022521.8810.13829.161 * 1035 3.442 * 10-52.748 * 10100.0191
Hyperion8 * 10171.339 * 10-70.27296.912 * 10-45.993 * 1030 2.251 * 10-106.875 * 1054.778 * 10-7
Iaepetus1.6 * 10212.678 * 10-41.0320.021141.86 * 1034 6.989 * 10-76.428 * 1084.467 * 10-4
Uranus8.685 * 102514.541.2370.90321.696 * 1042 63.731.174 * 10120.8157
Miranda6.6 * 10191.105 * 10-51.1990.008075.728 * 1031 2.152 * 10-99.371 * 1076.512 * 10-5
Ariel1.35 * 10212.26 * 10-41.660.027421.42 * 1033 5.336 * 10-81.082 * 1097.52 * 10-4
Titania3.53 * 10215.909 * 10-41.7160.038635.613 * 1033 2.109 * 10-72.147 * 1090.001492
Neptune1.024 * 102617.151.611.1382.505 * 1042 94.091.862 * 10121.294
Triton2.14 * 10220.0035822.0670.07973-3.333 * 1034 -1.252 * 10-69.146 * 1090.006356

Several densities of interest:

  • liquid water at 4 C - 1 g / cm3
  • gaseous H2 at 0 C and 1 atm - .00009 g / cm3
  • metallic H - 1.3 to 1.6 g / cm3
  • silicate minerals - 3 to 4 g / cm3
  • ferric/ferrous minerals - 7 to 8 g / cm3
Two other quantities worth noting are the orbital velocity (obtained from equating the gravitational acceleration to the acceleration v2/r required to keep a body moving in a circular orbit):
orbital velocity = (G * mass / distance)1/2
and the escape velocity (obtained from equating the kinetic energy mv2/2 to the gravitational energy Gm/r):
escape velocity = (2 * G * mass / distance)1/2
In both of these equations, the mass is the mass of the body we are orbiting or escaping from; the results are independent of the mass of our spacecraft. The distance in each equation is the distance from the center of the body (its radius in the case of the escape velocity).

Planetary magnetic fields are believed to be the product of rotation and liquid metallic cores. We have been able to measure the magnetic fields of some of these bodies using space probes:

BodyMagnetic Field (Earth)
Sun2
Mercury0.011
Venus0.001
Earth1
Mars0.001
Jupiter13.89
Saturn0.67
Uranus0.74
Neptune0.43
Here are images of auroras on Earth in optical wavelengths:

(source), Jupiter in x-ray wavelengths:

(source) and Saturn in ultraviolet wavelengths:

(source). There is also a good movie of terrestrial auroras from the POLAR mission (source).


Portfolio Exercise 6: Use the relation between the orbital period, the semimajor axis length and the total mass, to compute the masses of the Sun and those planets with moons.
Recall from our derivation of Kepler's Laws from Newton's that the mass in this relation is the combined mass of the parent and the satellite. For this reason, you will get the best results if you choose the satellite with the lowest mass (ie., use Mercury's orbital parameters to compute the mass of the Sun, Deimos' to compute Mars' mass, etc.).
Compare these values to the measured values quoted above: in terms of a percentage, how far off are your values from those in the table above?

For each body in the table above, compute the orbital velocity at a distance of 400 km from the surface, and the escape velocity from the surface. Compute each value in both km/hr and in multiples of the Earth's value.


Gravity and Tides

Using the masses above, it is interesting to compare the strengths of the gravitational forces among some of these bodies. This table gives the gravitational force between some of the planets and their moons, expressed as a fraction of the gravitational force between the Sun and that planet:
PlanetBodyFgrav (Sun)
EarthThe Moon0.005594
JupiterIo0.1528
JupiterEuropa0.03248
JupiterGanymede0.03943
JupiterCallisto0.009244
SaturnJupiter0.004623
SaturnTethys0.00739
SaturnDione0.007904
SaturnRhea0.008511
SaturnTitan0.09224
UranusMiranda0.01622
UranusAriel0.1534
UranusTitania0.07685
NeptuneTriton1.73
We can see that many of the moons have a significant perturbative influence on the orbital motion of their planets.

It is also interesting to compare the tidal accelerations caused by some of these bodies on others. We will define the tidal acceleration as the difference between the gravitational acceleration on the near side of the body and that at its center:

accelerationtidal = G * mass1 * (1 / (distance - radius2)2 - 1 / distance2)
where the distance is the distance between the bodies, and the subscripts refer to their labels below. These values have been expressed as fractions of the tidal acceleration of the Moon on the Earth (1.286 * 10-6 m/s2), and as fractions of the body's own surface gravity.
Here are some time lapse photos of tidal phenomena caused by the Moon.
In all cases, the length of the orbital semimajor axis has been used to compute the distance:
Body 1Body 2atidal (Moon on Earth)atidal / g
SunMercury2.9559.008 * 10-7
SunVenus1.1231.429 * 10-7
SunEarth0.44815.161 * 10-8
SunThe Moon0.12188.43 * 10-8
EarthThe Moon21.71.502 * 10-5
The MoonEarth11.152 * 10-7
MarsPhobos12740.387
MarsDeimos45.220.01899
JupiterIo54760.003438
JupiterEuropa11670.001007
JupiterGanymede484.53.837 * 10-4
JupiterCallisto80.967.349 * 10-5
IoJupiter13.336.069 * 10-7
IoEuropa1.0819.327 * 10-7
IoGanymede0.10298.146 * 10-8
EuropaJupiter1.5937.252 * 10-8
EuropaIo0.67684.25 * 10-7
EuropaGanymede0.23771.882 * 10-7
GanymedeJupiter1.1355.168 * 10-8
GanymedeIo0.11787.395 * 10-8
GanymedeEuropa0.43433.748 * 10-7
CallistoJupiter0.14446.574 * 10-9
SaturnMimas20670.03581
SaturnEnceladus12330.01818
SaturnTethys13960.01058
SaturnDione701.90.003384
SaturnRhea351.60.001503
SaturnTitan95.27.935 * 10-5
SaturnHyperion1.8373.062 * 10-4
SaturnIaepetus1.0695.825 * 10-6
TethysSaturn0.24782.677 * 10-8
TethysEnceladus0.10151.496 * 10-6
DioneSaturn0.192.053 * 10-8
DioneTethys0.12299.31 * 10-7
RheaSaturn0.13521.461 * 10-8
TitanSaturn0.56736.13 * 10-8
TitanHyperion0.081131.352 * 10-5
UranusMiranda11090.01583
UranusAriel8580.003604
UranusTitania97.812.917 * 10-4
MirandaUranus0.12751.626 * 10-8
ArielUranus0.739.309 * 10-8
ArielMiranda0.16632.374 * 10-6
TitaniaUranus0.14091.797 * 10-8
NeptuneTriton368.95.33 * 10-4
TritonNeptune1.5651.585 * 10-7
All of the moons in this table are locked in synchronous orbits: their rotational period is equal to their orbital period. This is due to tidal deformations (bulges toward the planet) which hold the same face toward the planet at all times. Mercury is also in a sychronous orbit, but in a more complicated fashion because of its higher eccentricity: it rotates 3 times for every 2 orbits around the Sun.

In addition, it is interesting to compare the tidal accelerations to that of Europa on Io, which is thought to account for much of Io's internal heating, and subsequent active volcanism. Note that values in the final column are all less than 1. If the tidal acceleration were equal to or greater than the surface gravity, the body would be gravitationally unstable and would break up. This happens to bodies lying withing the Roche Limit. Saturn's Rings lie within Saturn's Roche Limit:

(source)
In this image, green hues indicate particle sizes less than 5 cm, while blue indicates less than 1 cm. Purple regions had pieces larger than 5 cm (up to several meters), and white indicates an area so dense that it was opaque to the radio waves used to probe the particle sizes.


Temperature and Atmosphere

By measuring the wavelength of maximum power output and using Wien's Law
wavelength of maximum thermal outputAngstroms = 2.9 * 107 Angstroms / temperature
we can measure the effective (black body) temperature of many of these bodies. Using Stefan's Law
power emitted or absorbed per unit area = 5.67 * 10-8 * temperature4
and the fact that power per unit area (intensity) decreases with the square of the distance, we can compare the intensity of the radiation they receive to that which Earth receives (1365 W/m2), and the ratio of their output to input intensity (for moons, we include the input both from the Sun and from their host planet).

In addition, for some of these bodies we have been able to measure the surface temperature using space probes. Using these temperatures, we can learn something about possible atmospheres. If the average molecular speed of a gas

average molecular speed = 157 * (temperature / molecular weightH = 1)1/2
is less than 1/6 of the planet's escape velocity at the surface
escape velocity = (2 * G * mass / radius)1/2
there is a high probability that the gas still exists in an atmosphere near the surface. We can use these two equations to find the minimum molecular weight of any atmosphere for these bodies (the majority atmospheric component has been provided for comparison purposes, when known):
BodyTeffISun (Earth)Iout/IinMax Tsurface Min Mol. Wt.Major Atmos. Comp.Surface Atmos. Pressure (Earth)
Sun57770H
Mercury5336.6741.005700342 * 10-12
Venus2271.9110.23087356C O290
Earth25510.70253312N21
The Moon38713.727396621 * 10-12
Mars2170.43080.85532689C O2.006
Jupiter1250.036941.0980H2
Io1280.036941.19817
Europa1280.036941.2042710-7
Ganymede1280.036941.20615
Callisto1280.036941.20719
Saturn950.01011.2310H2
Titan850.01010.78839411N21.6
Uranus570.0027150.6460H2
Miranda860.0027153.3272044
Ariel840.0027153.038239
Neptune590.0011061.820H2
Triton380.0011060.312615N21.4 * 10-5
Some possible atmospheric components of interest are:
Atom or MoleculeMolecular Weight
H22
He4
C H4 (methane)16
N H3 (ammonia)17
H2 O18
C O28
N228
C2 H6 (ethane)30
O232
Ar40
C O244
C3 H8 (propane)44

(View GOES Isabel mov (source); Saturn mov (source), Cosmos DVD 4, episode 6, on Jupiter, the Red Spot and Neptune.)

After viewing these large-scale weather patterns on various planets, it is interesting to observe something much smaller on Earth, and on Mars:
(View Martian Dust Devils as photographed by Spirit (source))
Compare two dust storms, one on Earth off the northwest coast of Africa, and one on Mars in the north polar regions (source):

Dust is responsible for heat transport on Mars much as water is here on Earth.

Gravity Waves are seldom noticed in action because of their long periods; this movie is a time lapse over 40 minutes:

(View movie of gravity wave (source))

Solar System Formation

Temperatures in the protoplanetary disc are much higher in the central regions where the star is forming, due to the increased density and pressure. Further away the temperatures cool quickly, with the outer regions cooling faster. The general makeup of planets as a function of distance from the protostar can be understood by the various melting points of the materials involved:
MaterialMelting Point (K)
Iron1808
Silicon1684
Ice273
Ammonia195
Methane91
AU Microscopii is a red dwarf star approximately 12 million years old with a debris disc; the radius of the cleared central region is approximately 7 AU:

(source)

HD 107146 is a yellow star much like our Sun but between 30 and 250 million years old. The radius of the cleared central region is approximately 40 AU:

(source)

In comparison, the Kuiper Belt (home to short period comets) extends to about 50 AU from the Sun. The Oort Cloud (home to long period comets) has a radius of approximately 100000 AU.

The mass of Comet Halley is about 1014 kg. It is about 40% water and about 20% organic compounds. Assuming that it is an average comet, and that comets were responsible for the presence of water on the Earth, we can compute the frequency of cometary impact in the Earth's early history. These impacts occurred mainly during a period of about 200 million years, approximately 4 billion years ago. If the total amount of water on the Earth is 1.65 * 1021 kg, it would have taken 41.25 million cometary impacts, or an average of about 5 impacts a year. These impacts would have deposited a further 8.25 * 1020 kg of organic material.
Here is a Hubble animation revealing a protoplanetary disc in the Orion Nebula (source). Note that during formation, accretion and fragmentation are both occurring, as well as frequent collisions.
(View Animations of Solar System Formation (source) (source))


©2008, Kenneth R. Koehler. All Rights Reserved. This document may be freely reproduced provided that this copyright notice is included.

Please send comments or suggestions to the author.