Category Archives: Comets

Breaking an asteroid

Hi there! Asteroids, these small bodies in the Solar System, are fascinating by the diversity of their shapes. This is a consequence of their small sizes. Another consequence is their weakness, which itself helps to split some of them into different parts, sometimes creating binary objects, asteroids families… The study I present you today, Internal gravity, self-energy, and disruption of comets and asteroids, by Anthony R. Dobrovolskis and Donald G. Korycansky, proposes an accurate computation of the required energy to provoke this break-up, at any place of the asteroid, i.e. you are more efficient when you hit at a given location. This study has recently been accepted for publication in Icarus.

Shapes of asteroids

Please allow me, in this context, to call asteroid a comet, a comet being a small body, i.e. like an asteroid, but with a cometary activity. The important thing is that the involved bodies are small enough.

Beyond a given size, i.e. a diameter of ~400 km, a planetary body is roughly spheroidal, i.e. it is an ellipsoid with it two equatorial axes almost equal and the polar one smaller, because of its rotation. For a tidally despun body, like the Moon, or a satellite of a giant planet, the shape is more triaxial, since the tidal (gravitational) action of the parent planet tends to elongate the equatorial plane. The same phenomenon affects Mercury.

However, for smaller bodies, the self-gravitation is not strong enough to make the body look more or less like a sphere. As a consequence, you can have almost any shape, some bodies are bilobate, some are contact binaries, i.e. two bodies which permanently touch together, some others are rubble piles, i.e. are weak aggregates of rocks, with many voids.

These configurations make these bodies likely to undergo or have undergone break-up. This can be quantified by the required energy to extract some material from the asteroid.

The energies involved

For that, an energy budget must be performed. The relevant energies to consider are:

  • The impact disruption energy: the minimum kinetic energy of an impactor, to shatter the asteroid and remove at least half of its mass,
  • The shattering energy: the minimum energy needed to shatter the asteroid into many small pieces. It is part of the impact disruption energy. This energy is roughly proportional to the mass of the asteroid. It represents the cohesion between the adjacent pieces.
  • The binding energy: this energy binds the pieces constituting the asteroid. In other words, once you have broken an asteroid (don’t try this at home!), you have to make sure the pieces will not re-aggregate… because of the binding energy. For that, you have to bring enough energy to disperse the fragments.
  • The self-gravitational energy: due to the mutual gravitational interaction between the blocks constituting the asteroids. Bodies smaller than 1 km are strength-dominated, i.e. they exist thanks to the cohesion between the blocks, which is the shatter energy. However, larger bodies are gravity-dominated.
  • The kinetic energy of rotation: the spin of these bodies tends to enlarge the equatorial section. In that sense, it assists the break-up process.

This study addresses bodies, which are far enough from the Sun. This is the reason why I do not mention its influences, i.e. the tides and the thermic effects, which could be relevant for Near-Earth Objects. In particular, the YORP effect is responsible for the fission of some of them. I do not mention the orbital kinetic energy of the asteroid either. Actually the orbital motion is part of the input energy brought by an impact, since the relative velocity of the impactor with respect to the target is relevant in this calculation.

I now focus on the two cases studied by the authors to illustrate their theory: the asteroid Kleopatra and the comet 67P/Churyumov-Gerasimenko.

2 peculiar cases: Kleopatra and Churyumov-Gerasimenko

216 Kleopatra is a Main-Belt asteroid. Adaptive optics observations have shown that is is constituted of two masses bound by material, giving a ham-bone shaped. As such, it can be considered as a contact binary. It is probably a rubble pile. Interestingly, observations have also shown that Kleopatra has 2 small satellites, Alexhelios and Cleoselene, which were discovered in 2008.

Reconstruction of the shape of Kleopatra. © NASA
Reconstruction of the shape of Kleopatra. © NASA

However, 67P Churyumov-Gerasimenko is a Jupiter-family comet, i.e. its aphelion is close to the orbit of Jupiter, while its perihelion is close to the one of the Earth. It has an orbital period of 6.45 years, and was the target of the Rosetta mission, which consisted of an orbiter and a lander, Philae. Rosetta orbited Churyumov-Gerasimenko between 2014 and 2016. The shape of this comet is sometimes described as rubber ducky, with two dominant masses, a torso and a head, bound together by some material, i.e. a neck.

Churyumov-Gerasimenko seen by Rosetta. © ESA
Churyumov-Gerasimenko seen by Rosetta. © ESA
216 Kleopatra 67P/Churyumov-Gerasimenko
Semimajor axis 2.794 AU 3.465 AU
Eccentricity 0.251 0.641
Inclination 13.11° 7.04°
Spin period 5.385 h 12.761 h
Mean radius 62 km 2.2 km
Magnitude 7.30 11.30
Discovery 1880 1969

The irregular shapes of these two bodies make them interesting targets for a study addressing the gravitation of any object. Let us see now how the authors addressed the problem.

Numerical modeling

Several models exist in the literature to address the gravity field of planetary bodies. The first approximation is to consider them as spheres, then you can refine in seeing them as triaxial ellipsoids. For highly irregular bodies you can try to model them as cuboids, and then as polyhedrons. Another way is to see them as duplexes, this allows to consider the inhomogeneities dues to the two masses constituting bilobate objects. The existence of previous studies allow a validation of the model proposed by the authors.

And their model is a finite-element numerical modeling. The idea is to split the surface of the asteroid into small triangular planar facets, which should be very close to the actual surface. The model is all the more accurate with many small facets, but this has the drawback of a longer computation time. The facets delimit the volume over which the equations are integrated, these equations giving the local self-gravitational and the impact disruption energies. The authors also introduce the energy rebate, which is a residual energy, due to the fact that you can remove material without removing half of it. This means that the impact disruption energy, as it is defined in the literature, is probably a too strong condition to have extrusion of material.
The useful physical quantities, which are the gravitational potential, the attraction, and the surface slope, are propagated all along the body thanks to a numerical scheme, which accuracy is characterized by an order. This order quantifies the numerical approximation which is made at each integration step. A higher order is more accurate, but is computationally more expensive.

Once the code has been run on test cases, the authors applied it on Kleopatra and Churyumov-Gerasimenko, for which the shape is pretty well known. They used meshes of 4,094 and 5,786 faces, respectively.

Results

The validation phase is successful. The authors show that with a 3rd order numerical scheme, they recover the results present in the literature for the test cases with an accuracy of ~0.1%, which is much better than the accuracy of the shape models for the real asteroids. Regarding Kleopatra and Churyumov-Gerasimenko, they get the gravity field at any location, showing in particular excesses of gravity at the two lobes.

Such a study is particularly interesting for further missions, which would determine the gravity field of asteroids, which would then be compared with the theoretical determination by this code. Other applications are envisaged, the authors mentioning asteroid mining.

The study and its authors

And that’s it for today! Please do not forget to comment. You can also subscribe to the RSS feed, and follow me on Twitter and Facebook.

And Merry Christmas!

Rotation and activity of a comet

Hi there! We, Earthians, are regularly visited by periodic comets, the most famous one being probably 1P/Halley, which will visit us in 2061. Since we cannot wait, we study others of that kind. Today I tell you about 49P / Arend-Rigaux. This is the opportunity for me to present you The rotation and other properties of Comet 49P/Arend-Rigaux, 1984 – 2012, by Nora Eisner, Matthew M. Knight and David G. Schleicher. This study has recently been published in The Astronomical Journal.

The comet 49P / Arend-Rigaux

The comet 49P / Arend-Rigaux has been discovered in February 1951 at the Royal Observatory of Belgium, by Sylvain Arend and Fernand Rigaux. It is a periodic comet of the Jupiter family, i.e. with a period smaller than 20 years. Its period is actually 6.71 years, its semimajor axis 3.55 AU (astronomical units, 1 AU being 150 millions km, i.e. the Sun-Earth distance), its eccentricity 0.6, and its orbital inclination 19°, with respect to the ecliptic. These numbers are extracted from the JPL Small-Body Database Browser, and are calculated at the date Apr 6, 2010. I have plotted below the distances Sun-comet and Earth-comet.

Distance to the Sun.
Distance to the Sun.
Distance to the Earth.
Distance to the Earth.

The distance to the Sun clearly shows the periodic variations. The orbit of the Earth is at 1 AU, the one of Mars at 1.5 AU, and the one of Jupiter at 5.2 AU. Every 6.71 years, the comet reaches its perihelion, i.e. minimizes its distance to the Sun. This proximity warms the comet and provokes an excess of cometary activity, i.e. sublimation of dirty ice. At these occasions, the distance with the Earth is minimized, but with variations due to the orbital motion of the Earth. We can see for instance that the comet gets pretty close to the Earth in 1951 (when it was discovered), in 1984, and in early 2032. These are favorable moments to observe it. The paper I present you today is mainly (but not only) based on photometric observations made between January and May 2012, at Lowell Observatory.

Observations at Lowell Observatory

Lowell Observatory is located close to Flagstaff, AZ (USA). It was founded by the famous Percival Lowell in 1894, and is the place where Clyde Tombaugh discovered Pluto, in 1930. Among its facilities is the 4.28 m Discovery Channel Telescope, but most of the data used in this study were acquired with the 1.1 m Hall telescope, which is devoted to the study of comets, asteroids, and Sun-like stars. The authors also used a 79 cm telescope. The observations were made in the R(ed) band.

The data

Besides these 33 observation nights during the first half of 2012, the authors used data acquired close to the 1984 and 2005 perihelion passages, even if the 2005 ones revealed unusable. The observations consists to measure the magnitude (somehow, the luminosity) of the comet, in correcting for atmospheric problems, so as to be able to detect the variations of this magnitude. You can find below an example of data:

Magnitude of 49P / Arend-Rigaux measured in April 2012.
Magnitude of 49P / Arend-Rigaux measured in April 2012.

Of course, the data have holes, since you cannot observe during the day. Moreover, the comet needs to be visible from Arizona, otherwise it was just impossible to observe it and make any measurements.

We can see a kind of periodicity in the magnitude, this is a signature of the rotation of the comet.

Measuring the rotation

Most of the planetary bodies are kinds of triaxial ellipsoids. Imagine we are in the equatorial plane of one of them. We see an alternation of the long and short axes of its equatorial section. If the albedo of the surface element we face depends mainly on its curvature (it depends on it, but mainly may be an overstatement), then we should see two peaks during a period. As a consequence, the period of the lightcurve we observe should be half the rotation period of the comet.

In combining all the measurements, the authors managed to derive a rotation period of 13.45 ± 0.01 hour. For that, they used two different algorithms, which gave very close results, giving the authors confidence in their conclusions. The first one, Phase Dispersion Minimization (PDM), consists to assume a given period, split the measurements into time intervals of this period, and overlap them. The resulting period gives to the best overlap. The other algorithm is named Lomb-Scargle, following its authors. It is a kind of Discrete Fourier Transform, but with the advantage of not requiring uniformly sampled data.

In addition to this rotation period, the authors detected an increasing trend in the 2012 data, as if the spin of the comet accelerated. This is in agreement with an alteration of the measured rotation from the Earth, which moves, and reveals a retrograde rotation, i.e. an obliquity close to 180°. In other words, this is an illusion due to the motion of the observer, but this illusion reveals the obliquity.

Moreover, in comparing the 2012 data with the ones of 1984, the authors managed to detect a variation in the rotation period, not larger than 54 seconds. This is possible regarding the fact that the comet is altered by each perihelion passage, since it outgasses. In this case, that would imply a change of at the most 14 seconds of the rotation period between two passages. Such variations have also been detected for at least 4 other comets (2P/Encke, 9P/Tempel 1, 10P/Tempel 2, and 103P/Hartley 2, see Samarinha and Mueller (2013)).

Comet Period (h) Variation (s)
2P/Encke 11 240
9P/Tempel 1 41 -840
10P/Tempel 2 9 16.2
103P/Hartley 2 18 7200
49P/Arend-Rigaux 13.45 -(>14)

Finally, since the lightcurve is a signature of the shape as well, the authors deduced from the amplitude of variation that the axial ratio of the nucleus, i.e. long axis / short axis, should be between 1.38 and 1.63, while an independent, previous study found 1.6.

Cometary activity

49P / Arend-Rigaux has a low activity. Anyway, the authors detected an event of impulse-type outburst, which lasted less than 2 hours. The analysis of the coma revealed an excess of cyanides with respect to the 1984 passage. Moreover, 49P / Arend-Rigaux is the first comet to show hydroxyde.

The study and its authors

That’s it for today! Please do not forget to comment. You can also subscribe to the RSS feed, and follow me on Twitter and Facebook.

The activity of the comet C/2015 ER61

Hi there! Today’s post is on the comet C/2015 ER61. Behind this weird name is a small object orbiting the Sun on a highly elongated orbit, which currently shows us a tail. The associated study is Beginning of activity in long-period comet C/2015 ER61 (PANSTARRS), by Karen J. Meech, Charles A. Schambeau, Kya Sorli, Jan T. Kleyna, Marco Micheli, James Bauer, Larry Denneau, Jacqueline V. Keane, Elizabeth Toller, Richard Wainscoat, Olivier Hainaut, Bhuwan Bhatt, Devendra Sahu, Bin Yang, Emily Kramer and Gene Magnier. It has recently been published in The Astronomical Journal.

C/2015 ER61‘s facts

This comet was discovered two years ago, in March 2015, by the telescope Pan-STARRS 1, located on the Haleakalā, Hawai’i. Its distance to the Sun was then 8.44 Astronomical Units, its absolute magnitude about 12, and no tail was visible. As such, it was supposed to be a Manx object, a Manx being a tailless cat. A Manx object would be a comet, which had no activity anymore, as if the lighter elements had already gone.

From its magnitude, it was guessed that its radius was about 10 km. Its apparent lack of activity triggered enough interest for the object to be followed, this in particularly permitted to determine its orbit, and showed that it had a huge eccentricity, i.e. some 0.998. When the eccentricity reaches 1, then the orbit is parabolic, so the orbit of C/2015 ER61 is almost parabolic. Further observations showed the beginning of a period of activity, proving that C/2015 ER61 (I would appreciate a funnier nickname…) is actually not a Manx. This period is not done yet, and the activity is actually increasing, as the comet is approaching the Sun. At its smallest distance, i.e. the perihelion, its distance to the Sun is 1.04 AU, i.e. it almost crosses the orbit of the Earth (don’t worry, I said “almost”). So, observing this comet today reveals a tail.

We are actually pretty lucky to be able to observe it, since its orbital period is some 10,000 years. This comet is considered to belong to the Oort cloud, which is a reservoir of comets at the edge of our Solar System.

Cometary outgassing

Since the comet model by Fred L. Whipple, published between 1950 and 1955, a comet is seen as a kind of dirty snowball, with a nucleus, and icy elements, which tend to sublimate when approaching the Sun, because of the elevation of the temperature. This hypothesis was confirmed in 1986 when we were visited by the well-known comet 1P/Halley (you know, Halley’s comet).
The idea is this: you have some water ice, some CO, some CO2, trapped on the comet. When it is warm enough, it sublimates.

But the intensity of the sublimation depends on several parameters:

  • the thermal inertia of the comet: how does the temperature elevate?
  • its albedo: which fraction of the incident Solar light flux is reflected?
  • its density
  • the quantity of elements, which are likely to be sublimated
  • their depth: if they are not at the surface, the heat needs to be conducted deep enough for them to sublimate
  • the distance to the Sun (of course)
  • etc.

This means that observing and measuring this outgassing gives some physical properties of the comet.

The observation facilities

To conduct this study, several observation facilities were used:

  • Pan-STARRS1 (PS1): This stands for Panoramic Survey Telescope and Rapid Response System. This is a 1.8m wide-field telescope,
  • Gemini North: this is a 8.19 m telescope, which is based in Hawai’i. It has a twin brother, Gemini South, which is based in Chile,
  • Canada-France-Hawai’i Telescope (CFHT): this 3.58m telescope is part of the Mauna Kea Observatory. For this study, the MegaPrime/Megacam wide-field imager was used, which gives of fied of view of 1°,
  • ATLAS: (for Asteroid Terrestrial-impact Last Alert System). This will be a network of two 0.5m-telescopes, both based in Hawai’i. At this time, only the ATLAS-Haleakalā has begun full operation,
  • Himalayan Chandra Telescope (HCT): this is a 2.01 m optical-infrared telescope, which is part of the Indian Astronomical Observatory, which stands on Mount Saraswati, Digpa-ratsa Ri, Hanle, India,
  • Wide-field Infrared Survey Explorer (WISE): this is an infrared space telescope, on a Sun-synchronous polar orbit. It is used in the program NEOWISE, NEO standing for Near-Earth Objects.

The diversity of observation facilities explains the numbers of authors signing this study. The observations span from February 2014 to February 2017, which means that there are pre-discovery observations. It is always easier to find an object when you know where it is, which permitted to find C/2015 ER61 on images, which were taken before its discovery.

Results

These observations (see the Figure) has shown a variation of the magnitude, which could be expected since the comet approached the Earth, but too large to be explained by its trajectory. Actually, it is enhanced by the activity of the comet, more precisely by the sublimation of CO and CO2, starting in early 2015.

The measured apparent magnitude of the comet, with respect to the date and the distance to the Sun. We can see that the comet is brighter when closer to the Sun, because of the outgassing. The measurements have some uncertainties, which are not represented here. This figure is drawn for the Tab.1 Observation Log of the paper.

The authors modeled the warming of the comet and the sublimation of the elements, in using the well-known heat equation. The observed tail suggests a radius of the nucleus of about 9 km, which is consistent with previous guesses. Moreover, they suggest that the CO2 is present at a depth of about 0.4 m. If it were present at the surface, then sublimation would have been observed even when the comet was 20 AU away from the Sun.

The closest approach of the comet with the Earth was on April 4, and with the Sun on May 10, which would result in a peak of activity… probably with some delay, please give the comet a chance to warm!

To know more

That’s all for today! Please do not forget to comment. You can also subscribe to the RSS feed, and follow me on Twitter.