Comet Shoemaker-Levy 9 Collision with Jupiter
16-22 July 1994

An animation of one impact fragment

Terry Robb
<terryr@wolfram.com>

Between 16 and 22 July 1994 comet P/Shoemaker-Levy 9 will collide with Jupiter. On the comet's previous pass around Jupiter in July 1992 it came so close that it moved within the Roche limit and subsequently shattered into many fragments. These fragments, now numbering 21 major bodies, range in size between 1 to 5 kilometres and all will eventually impact on Jupiter. This is the first major collision between any two solar system bodies that has ever been predicted in advance. At last count there will be seven spacecraft and a large proportion of Earth-based observatories watching these events, both before and after they occur.

The 21 fragments have been labeled A through W, and they will all impact on Jupiter at a point just out of direct view from Earth, but near the edge on the far side of Jupiter. However within 10 or 15 minutes after impact the impact region will rotate into view of Earth. Also one impact – the bright fragment labelled 'K=12' – may cause a flash to be reflected off the surface of Europa which will be in eclipse during the impact. It is possible that this flash might be detected from Earth.

In this Notebook I have produced a simple animation to illustrate the K=12 impact. The animation was produced using the Scientific Astronomer package for Mathematica, and it will give you an idea of the timescale of the impact, and also things to look out for during the event. With good binoculars you should be able to see the four Galilean moons orbiting Jupiter, and with a good telescope you should be able to see the cloud features on Jupiter including the Great Red Spot. It is unlikely that you will be able to see any impact, but you might be able to see an after effect if you have a very good telescope and if the impact causes a disturbance in the cloud patterns.

For those of you who don't know where to look to see Jupiter, it should be very easy to find during July 1994. In the mid-evening sky Jupiter will be the brightest object high in the sky. Earlier in the evening the only object brighter will be Venus which will be following the Sun and setting in the western horizon.

Impact times

Impact times have been calculated by Paul Chodas et al. There is still some uncertainly in the precise times, so a colunm of uncertainty is also given. In general the smaller the uncertainty, the bigger the fragment.

In the table the angle E-J-F is the Earth-Jupiter-Fragment angle at impact. Values greater than 90˚ indicate a farside impact. All impacts will be just on the farside as viewed from Earth, but later impacts will be closer to the limb. Each impact will however rotate into sight of Earth after about 15 minutes though.

The labelling system looks a bit chaotic because the letters I and O were never used, and also because the J and M fragments have faded from view. The initially named P and Q fragments have also subsequently split into two smaller fragments each.

Earth viewpoint

The K=12 impact occurs at about 19 July 1994 10:21:00 (Universal Time), which in Mathematica Date notation could be written as {1994,7,19,10+TimeZone[],21,0}. This is equivlavent to 8:21pm Melbourne time (Melbourne, Australia is 10 hours ahead of GMT or Universal Time).

Using the Scientific Astronomer function PlanetPlot3D, we can see a view of Earth as seen from Jupiter during the K=12 impact:

PlanetPlot3D[Earth,
        {1994,7,19,10+TimeZone[],30,0},
        ViewPoint->Jupiter];

The image shows that Australia is best placed to view the event because not only is it facing Jupiter during the K=12 impact, but also because it is in darkness at that time. Local sunset occurs at:

SunSet[{1994,7,19}]

which is about 5:21pm. It's winter in Melbourne during July, hence the Sun sets early.

And from Melbourne, at the time of the K=12 impact, Jupiter will be well placed in the sky:

HorizonCoordinates[Jupiter,
    {1994,7,19,10+TimeZone[],30,0}]

The above says that Jupiter will be about 55˚ above the horizon at 8:30pm. What's more the position of the Moon is not too close:

Separation[Jupiter, Moon,
    {1994,7,19,10+TimeZone[],30,0}]

Although the next full moon will be four days away:

FullMoon[{1994,7,19}]

so the Moon will be quite bright.

The K=12 impact has one special property. It will accur when the moon Europa is visible from Earth, but eclipsed from the Sun. A reflected flash might therefore be visible. Of the 21 impacts it is the only one that is certain to happen during a Jovian lunar eclipse.

Comet impact

To produce an animation of the K=12 impact I have used the standard Scientific Astronomer PlanetPlot3D function to show a graphic of Jupiter along with it's moons and cloud features (the Great Red Spot and the three lesser known White Spots). Labelling and comet trajectory were added as an after thought by using the Epilog option. Additional options such as PlotRange and PlotRegion were used to get a slightly bigger field of view centred just to the left of Jupiter where all the impacts occur.

So here is the animation of Jupiter, with a frame every 15 minutes, from 6pm on 19 July 1994 until 4am the next morning (Melbourne time):

Do[PlanetPlot3D[
     Jupiter, {1994,7,19, 10+TimeZone[],30+m,0},
     ViewPoint  -> Earth,
     PlotRange  -> {{-1,1},
                    {-1,1},
                    {-1,1}}*72000*2.5,
     PlotRegion -> {{-0.5,2}, {-1,2}},
     Epilog     -> annotate[m]],
  {m, -150, 435, 15}]

The two moons visible during this animation are Io and Europa. Both are shown as green dots, but Europa is the dot appearing from behind Jupiter. Due to the Sun being at an angle – slightly towards the right but behind the viewpoint above – the moon Europa will be in eclipse (and shaded as a darker green in the graphic) for about two hours or so. During that eclipse the K=12 impact will occur. Later the shadow of Io can be seen passing across the Jovian disk. In the animation the blue comet trajectory line, as well as the blue impact dots were added with the Epilog option. Local clock time is shown with white text.

At the time of impact, I have also used yellow lines to indicate the flash that will occur. In the first version of this animation I used a Rectangle[mins, maxs, soundgraphic] function to generate a crashing sound when the impact frame appeared – but this got annoying after a while so I took it out.

Finally the reason I chose the duration of the animation to be 10 hours, was because the Great Red Spot approximately repeats one full orbit around Jupiter in that time, and so if you run the animation in a loop it will appear to repeat without a glitch.

Here's the definition of annotate[m_] as passed to Epilog

Based on some orbital elements I found for the Q=7 fragment of the comet, I was able to add a new object to the Scientific Astronomer:

SetOrbitalElements[Q, ViewPoint->Jupiter, Date->{1994,7,15,TimeZone[],0,0},
  OrbitalMeanMotion    -> 1/930.*Revs/Day,
  OrbitalSemiMajorAxis -> 34776.7/(1-0.9987338)/149579000*AU,
  OrbitalEccentricity  -> 0.9987338,
  OrbitalInclination   -> 94.2333*Degree,
  MeanLongitude        -> (43.47999-360/930.*(20.7846-15))*Degree,
  PerigeeLongitude     -> (43.47999)*Degree,
  AscendingLongitude   -> (290.87450+55.19080)*Degree];

From the above elements I worked out the final impact trajectory and hence found the comet's position relative to Jupiter in 5 minute intervals to be:

comet={{0.451195, 0.269150}, {0.450219, 0.283792}, 
       {0.448266, 0.298433}, {0.445338, 0.314051}, 
       {0.443386, 0.329669}, {0.44241 , 0.346263}, 
       {0.441434, 0.362856}, {0.440457, 0.381402}, 
       {0.440457, 0.400925}, {0.440457, 0.420447}};

I then used the above to add the blue annotations to the PlanetPlot3D animation.

annotate[m_] := Block[{d,h,m0=(10+TimeZone[])*60+30+m},
  d = Floor[m0/60/24];
  h = Floor[(m0-24*60*d)/60];
  m0= Floor[(m0-24*60*d)-60*h];
  {
(* If[m===0, Rectangle[{0,0},{1,1}, {}(*crash sound*)], {}], *)
  {RGBColor[0,0,1], Thickness[0], Line[comet],
   PointSize[0.01],
      Which[m==0, {RGBColor[1,1,0],
                  Line[{{.425,.440},{.4, .44}}],
                  Line[{{.435,.432},{.42,.412}}],
                  Line[{{.446,.420},{.45,.40}}]},
       -50<=m< 0, Point[comet[[m/5]]],
        300>m>10, Point[{0.5-0.07*Cos[0.01*(m+10)],
                         0.437}],
            True, {}]},
  {GrayLevel[1], Text["K=12 Impact", {.4,.70}]},
  {GrayLevel[1], Text[{1994,7,19+d,
    PaddedForm[h, 2],PaddedForm[m0, 2],0},
                                     {.4,.67}]}}]

Galilean moons

During the period over which the comet fragments will impact, the inner moons of Jupiter will make several orbits. The following chart shows the position of the four Galilean moons during the period 17 July until 23 July. Time is shown down the right side, and is measured in hours after local midnight Melbourne time.

JupiterMoonChart[{1994,7,17}, {1994,7,23}];

The black, red, green and blue lines represent the positions, as seen from Earth, of the Galilean moons Io, Europa, Ganymede and Callisto respectively. And the white line wrapping round the purple coloured Jovian disk represents the position of the Great Red Spot.

The package

Scientific Astronomer

Contact astronomer@wolfram.com to learn more about the package.

See also the web address:

  http://store.wolfram.com/view/astronomer/animation.html
  http://www.wolfram.com/~terryr/Astronomer/index.html
  http://www.wolfram.com/~terryr/apps/AsteroidTrack.html

Star Charts: Five types of charts are defined in Scientific Astronomer, including two wide field star charts. With the star charts you can zoom into any portion of the sky. All the charts have options to show star spectral colors, mesh lines, a skyline, the horizon line, and the Milky Way; and to label constellations, stars, planets, deep sky objects, and so on.

Planet Plots: Planet plotting is done in two- and three-dimensional forms. Surface features for the Earth, the Moon, Mars, and Jupiter are shown on the plots. Moons and their shadows are displayed for the Earth and Jupiter. Related functions allow you to produce planet position finder charts and planet rise/set timing charts.

Eclipses: Several functions are provided for dealing with eclipses. These functions provide information about both solar and lunar eclipses, and are general enough to handle Galilean moon eclipses, occultation of stars by the Moon, and transits of Mercury or Venus across the solar disk. You can produce umbra and penumbra track plots and perform eclipse prediction.

Satellite Tracking: Satellite tracking is another feature of Scientific Astronomer. You can create track plots, make visibility predictions, and project satellite tracks onto star charts.

Miscellaneous: Miscellaneous other features are available, such as producing planisphere plates, planet charts, and solar system plots. In addition, sunrise, moonrise, and full moon functions are provided, as well as functions for adding new objects, such as comets and satellites.

Scientific Astronomer is Mathematica 3.0 compatible. It has palettes and buttons and is fully integrated into the Help Browser system.