Perihelion motion of small irregular dust particles due to radiation forces

Perihelion motion, i.e. a secular change of longitude of perihelion, of interplanetary dust particles is investigated under the action of solar gravity and solar electromagnetic radiation. As for spherical particle [Kla?ka, J., 2004. Electromagnetic radiation and motion of a particle. Cel. Mech. Dynam. Astron. 89, 1-61]: (i) perihelion motion is of the order ( is heliocentric velocity of the meteoroid and c is the speed of light in vacuum), if a component of electromagnetic radiation acceleration is considered as a part of central acceleration; (ii) perihelion motion is of the first order in if the total electromagnetic radiation force is considered as a disturbing force. The new facts presented in this paper concern irregular dust particles. Detailed numerical calculations were performed for the grains ejected at aphelion of comet Encke. Perihelion motion for irregular interplanetary dust particles exists already in the first order in for both cases of central accelerations. Moreover, perihelion motion of irregular particles exhibits both positive and negative directions during the particle orbital motion. Irregularity of the grains causes not only perihelion motion, but also dispersion of the dust in various directions, also normal to the orbital plane of the parent body.     

Non-gravitational force modeling of Comet 81P/Wild 2: I. A nucleus bulk density estimate

The nucleus of Comet 81P/Wild 2 is modeled by assuming various smooth triaxial ellipsoidal or irregular body shapes, having different rotational periods, spin axis orientations, and thermophysical properties. For these model nuclei, a large number of surface activity patterns (e.g., maps of active and inactive areas) are studied, and in each case the resulting water production rate and non-gravitational force vector versus time are calculated. By requiring that the model nuclei simultaneously reproduce certain properties of the empirical water production curve and non-gravitational changes of the orbit (focusing on changes of the orbital period and in the longitude of perihelion), constraints are placed on several properties of the nucleus. The simulations suggest that the mass of Comet 81P/Wild 2 is , resulting in a rather low bulk density, (depending on the assumed nucleus volume), and that the nucleus rotation is prograde rather than retrograde. The active area fraction is difficult to constrain, but at most 60% of the nucleus is likely to have near-surface ice.     

Temporal deconvolution of the hydrogen coma: II. Pre- and post-perihelion activity of Comet Hyakutake (1996 B2)

Comet 1996 B2 (Hyakutake) displayed strong evidence for break-up, with a prominent antisunward dust spike and fragments traveling antisunward for many days after an eruptive event in late March 1996. Because of its high orbital inclination and rapid southward motion after perihelion, its post-perihelion activity was not well monitored from the ground. The SWAN all-sky Lyman-alpha camera on the SOHO spacecraft was ideally placed for long-term monitoring of the hydrogen coma of Comet Hyakutake both before and after perihelion. The SWAN images were analyzed with a new time-resolved model (TRM) that provides daily averages of the water production rate and an estimate of the hydrogen atom lifetime (dominated by charge exchange with solar wind protons) during extended periods throughout the apparition. A long-term variation of water production rate of , where r is the heliocentric distance in AU was found. The daily average values of the production rate covered the March 19 outburst and two more outbursts seen in the April before perihelion, which had progressively shorter durations at respectively smaller heliocentric distances. The long-term variation of the production rate was found to be consistent with the seasonal effect predicted by the jet rotation model of Schleicher and Woodney [2003. Analyses of dust coma morphology of Comet Hyakutake (1996 B2) near perigee: Outburst behavior, jet motion, source region locations, and the nucleus pole orientation. Icarus 162, 190-213] when added to a more steady source that is about two-thirds of the maximum of the jet source. The seasonal effect in their model found the dust jet source largely not illuminated after perihelion, coinciding with somewhat reduced overall activity and the absence of outbursts and fragmentation. The locations of the dust jets appear to be responsible for the outbursts and fragmentation before perihelion. The erratic behavior of the pre-perihelion jet sources as contrasted with the smoother variation from the rest of the surface after perihelion indicates there is a strong heterogeneity in the physical make-up of active areas on the nucleus.     

On the population, physical decay and orbital distribution of Jupiter family comets: Numerical simulations

We study the Jupiter family comet (JFC) population assumed to come from the Scattered Disk and transferred to the Jupiter’s zone through gravitational interactions with the Jovian planets. We shall define as JFCs those with orbital periods and Tisserand parameters in the range 2<T?3.1, while those comets coming from the same source, but that do not fulfill the previous criteria (mainly because they have periods ) will be called ‘non-JFCs’. We performed a series of numerical simulations of fictitious comets with a purely dynamical model and also with a more complete dynamical-physical model that includes besides nongravitational forces, sublimation and splitting mechanisms. With the dynamical model, we obtain a poor match between the computed distributions of orbital elements and the observed ones. However with the inclusion of physical effects in the complete model we are able to obtain good fits to observations. The best fits are attained with four splitting models with a relative weak dependence on q, and a mass loss in every splitting event that is less when the frequency is high and vice versa. The mean lifetime of JFCs with radii and is found to be of about 150-200 revolutions (∼. The total population of JFCs with radii within Jupiter’s zone is found to be of 450±50. Yet, the population of non-JFCs with radii in Jupiter-crossing orbits may be ∼4 times greater, thus leading to a whole population of JFCs + non-JFCs of ∼2250±250. Most of these comets have perihelia close to Jupiter’s orbit. On the other hand, very few non-JFCs reach the Earth’s vicinity (perihelion distances ) which gives additional support to the idea that JFCs and Halley-type comets have different dynamical origins. Our model allows us to define the zones of the orbital element space in which we would expect to find a large number of JFCs. This is the first time, to our knowledge, that a physico-dynamical model is presented that includes sublimation and different splitting laws. Our work helps to understand the role played by these erosion effects in the distribution of the orbital elements and lifetimes of JFCs.     

Asymmetries in the distribution of H2O and CO2 in the inner coma of Comet 9P/Tempel 1 as observed by Deep Impact

Prior to the impact event, Deep Impact monitored the ambient inner coma of Comet 9P/Tempel 1 at high spatial resolution in July 2005. Gaseous H₂O and CO₂ are unambiguously detected in the infrared spectra collected with the HRI-IR spectrometer aboard Deep Impact. Detailed distribution maps of these volatiles in the inner coma, within 60 km from the nucleus, are produced from the integrated emission bands of H₂O (2.66 μm) and CO₂ (4.26 μm). Uncorrelated asymmetries are determined in the spatial distribution of both species indicating chemical heterogeneities within the nucleus. Although present at some abundance surrounding the entire nucleus, H₂O has a pronounced enhancement in abundance in the sunward direction rotational phases, evidence that the dominant process of subliming water ice from the nucleus is solar heating. In contrast, CO₂ is enhanced in the regions near the negative rotational pole of the nucleus, suggesting localized outgassing there. Both species show an increase in radiance above the limb of the nucleus toward Ecliptic North. The distribution maps also suggest that the process of dust removal from the nucleus is strongly connected to the outgassing of volatiles. Detailed study of these coma asymmetries gives insight to the relative abundances of the dominant molecular components of the inner coma, source regions of the native volatiles, anisotropic outgassing of the nucleus, and the formation and evolution of the nucleus. A quiescent water production rate for Tempel 1 on July 4, 2005, is estimated to be .     

The dust trail of Comet 67P/Churyumov-Gerasimenko between 2004 and 2006

We report on observations of the dust trail of Comet 67P/Churyumov-Gerasimenko (CG) in visible light with the Wide Field Imager at the ESO/MPG 2.2 m telescope at 4.7 AU before aphelion, and at with the MIPS instrument on board the Spitzer Space Telescope at 5.7 AU both before and after aphelion. The comet did not appear to be active during our observations. Our images probe large dust grains emitted from the comet that have a radiation pressure parameter β<0.01. We compare our observations with simulated images generated with a dynamical model of the cometary dust environment and constrain the emission speeds, size distribution, production rate and geometric albedo of the dust. We achieve the best fit to our data with a differential size distribution exponent of −4.1, and emission speeds for a β=0.01 particle of 25 m/s at perihelion and 2 m/s at 3 AU. The dust production rate in our model is on the order of 1000 kg/s at perihelion and 1 kg/s at 3 AU, and we require a dust geometric albedo between 0.022 and 0.044. The production rates of large (>) particles required to reproduce the brightness of the trail are sufficient to also account for the coma brightness observed while the comet was inside 3 AU, and we infer that the cross-section in the coma of CG may be dominated by grains of the order of .     

Observations of the inner coma of C/1995O1 (Hale-Bopp)—gas and dust production

Hale-Bopp (C/1995 O1) was the most productive recent comet observed in terms of gas and dust output. Since its discovery in 1995 at a distance of 7.14 AU from the Sun, the comet has been well observed, revealing the dynamics of a rare and large comet. Hale-Bopp showed strong emissions of the principle cometary gases CN, C₃, and C₂, as well as an abundance of dust. The production rates of these gases were found to be 1.45×10²⁸, 1.71×10²⁸, and , respectively, with dust production, in terms of Afρ, , as measured in the green continuum (5260 Å). The observations for this paper are presented in two groups spanning 10 days each, one group centered near 32 days prior to and the other 21 days after perihelion. The averages of dust and gas production rates show a slightly higher value for each prior to perihelion than after perihelion, consistent with a possible peak in production a few weeks prior to perihelion passage.     

Non-constant rotation period of Comet C/2001 K5 (LINEAR)

The article presents results of CCD photometry in R-band of a dynamically new Comet C/2001 K5 (LINEAR), obtained at a heliocentric distance of about 5.6 AU, after the perihelion passage. Being so distant from the Sun, this comet was extremely active (Afρ close to 2000 cm), exhibiting quite well developed dust coma and tail. During the observations, general photometric behavior of the comet with heliocentric distance r was well described by the 2.5nlog(r) function with coefficient n=5. The radial profiles of the coma were found to be undulated, with mean slope of the dependence between cometary magnitude and 2.5log of aperture radius (at comet distance) equal to . The light curve of Comet LINEAR exhibited short-term variability which we attributed to cyclic changes of dust emission, induced by nucleus rotation. Model computations by some authors have revealed that active comets can change their spin status quite substantially even during a single orbital revolution. Thus, attempting to search for a rotation frequency, we have modified the classical PDM approach by including the spin acceleration term. Such DynamicalPDM (DPDM) method revealed the most reliable solution for the frequency f₀=0.019048±0.000013 h⁻¹ and its first time-derivative (index “zero” denotes reference to the mid time of the whole observing run), indicating a rapid spin-down of the nucleus. These parameters are equivalent to the rotation period of 52.499±0.036 h and its relative increment of 0.02729±0.00013. We present the most probable evolution of the rotation frequency of Comet LINEAR, based on the results of periodicity analysis and a simple, almost parameter independent, dynamical model of nucleus rotation. It is also shown that the DPDM may be an effective tool for determination of a nucleus radius, which provided us with the value of 1.53±0.25 km for Comet LINEAR.     

Seasonal and inter-annual changes of volume density of martian CO2 snow from time-variable elevation and gravity

The martian atmosphere seasonally exchanges CO₂ with the surface by repeating condensation and sublimation, causing seasonal growth and decay of the polar CO₂ snowcaps. These processes leave two kinds of geodetic signatures, i.e. seasonal changes of the martian gravity field and of surface elevation of the snow-covered regions. Here we study gradual increase of the volume density of the martian snow due to compaction, by combining these two data sets during 1999-2001 covering three martian winters. We found that light fresh snow of slowly becomes denser reaching or more immediately before it thaws. The maximum snow density varies slightly from year to year, and between hemispheres. In the second southern winter, the density became as high as . This might have been caused by a dust storm activity, e.g. increased mixing of silicate particles and/or enhancement of sintering.     

Cometary dust trail associated with Rosetta mission target: 67P/Churyumov-Gerasimenko

A thin, bright dust cloud, which is associated with the Rosetta mission target object (67P/Churyumov-Gerasimenko), was observed after the 2002 perihelion passage. The neckline structure or dust trail nature of this cloud is controversial. In this paper, we definitively identify the dust trail and the neckline structure using a wide-field CCD camera attached to the Kiso 1.05-m Schmidt telescope. The dust trail of 67P/Churyumov-Gerasimenko was evident as scattered sunlight in all images taken between September 9, 2002 and February 1, 2003, whereas the neckline structure became obvious only after late 2002. We compared our images with a semi-analytical dynamic model of dust grains emitted from the nucleus. A fading of the surface brightness of the dust trail near the nucleus enabled us to determine the typical maximum size of the grains. Assuming spherical compact particles with a mass density of 10³ kg m⁻³ and an albedo of 0.04, we deduced that the maximum diameter of the dust particles was approximately 1 cm. We found that the mass-loss rate of the comet at the perihelion was on or before the 1996 apparition, while the mass-loss rate averaged over the orbit reached . The result is consistent with the studies of the dust cloud emitted in the 2002/2003 return. Therefore, we can infer that the activity of 67P/Churyumov-Gerasimenko has showed no major change over the past dozen years or so, and the largest grains are cyclically injected into the dust tube lying along the cometary orbit.     

Deep Impact: Optical spectroscopy and photometry obtained at MIRA

We present spectroscopic and photometric observations, spanning the optical UV to the far red, before, during, and after the NASA Deep Impact event of July 4, 2005. The inner 2000 km of the pre- and post-impact coma was about 0.3 magnitude redder in B-R than in the outer coma. The pre-impact spectrum was a faint reflected solar spectrum dominated by molecular emissions extending >40,000 km from the nucleus. The post-impact light curve in R and I showed a rapid rise consistent with an expanding optically thick cloud during the first 18 min after impact. During the next 8 min the cloud became optically thin. Sixty minutes after impact the impact R-band flux reached a plateau at , the comet brightening by a factor of ∼4.3 above its pre-impact value observed in a 15″ aperture. The mean expansion velocity of the grains during the first 49 min was . The spectrum became dominated by scattered sunlight during the first hour after impact. The volume scattering function (VSF) observed 32 min after impact shows strong reddening. At 49 min, however, the VSF shows an additional twofold increase in the blue but only a 20% increase at 5500 Å. Post-impact spectra and R-I photometry showed rapid reddening. The particle size distribution, dominated by 1-2.5 μm particles shortly after impact, changed dramatically during the first hour due to sublimation of water-ice particles of this size. On the night following impact the comet was still substantially brighter than before impact, but R-I had returned to its pre-impact value. B-R remained significantly redder. The ejecta 25 h after impact was fan-shaped subtending ∼180° roughly symmetrical about position angle 225°. The mean expansion velocity 90° from the direction to the Sun was .     

Non-gravitational force modeling of Comet 81P/Wild 2: II. Rotational evolution

In this paper, we have studied both the dynamical and the rotational evolution of an 81P/Wild 2-like comet under the effects of the outgassing-induced force and torque. The main aim is to study if it is possible to reproduce the non-gravitational orbital changes observed in this comet, and to establish the likely evolution of both orbital and rotational parameters. To perform this study, a simple thermophysical model has been used to estimate the torque acting on the nucleus. Once the torque is calculated, Euler equations are solved numerically considering a nucleus mass directly estimated from the changes in the orbital elements (as determined from astrometry). According to these simulations, when the water production rate and changes in orbital parameters for 1997, as well as observational rotational parameters for 2004 are imposed as constraints, the change in the orbital period of 81P/Wild 2, , will decrease so that to , which is similar to the actual tendency observed from 1988 up to 1997. This nearly constant decreasing can be explained as due to a slight drift of the spin axis orientation towards larger ecliptic longitudes. After studying the possible spin axis orientations proposed for 1997, simulations suggest that the spin obliquity and argument (I,Φ)=(56°,167°) is the most likely. As for rotational evolution, changes per orbit smaller than 10% of the actual spin velocity are probable, while the most likely value corresponds to a change between 2 and 7% of the spin velocity. Equally, net changes in the spin axis orientation of 4°-8° per orbit are highly expected.     

Radio OH observations of 9P/Tempel 1 before and after Deep Impact

We observed 18-cm OH emission in Comet 9P/Tempel 1 before and after Deep Impact. Observations using the Arecibo Observatory 305 m telescope took place between 8 April and 9 June, 2005, followed by post-impact observations using the National Radio Astronomy Observatory 100 m Green Bank Telescope 4-12 July, 2005. The resulting spectra were analyzed with a kinematic Monte Carlo model which allows estimation of the OH production rate, neutral gas outflow velocity, and distribution of the out-gassing from the nucleus. We detected typically 24% variability from the overall OH production rate trend in the two months leading up to the impact, and no dramatic increase in OH production in the days post-impact. Generally, the coma is well-described, within uncertainties, by a symmetric model with OH production rates from 1.6 to , and mean water outflow velocity of . At these low production rates, collisional quenching is expected to occur only within 20,000 km of the nucleus. However, our best-fit average quenching radius is 64,200 ± 22,000 km in April and May.     

Submillimeter Wave Astronomy Satellite observations of Comet 9P/Tempel 1 and Deep Impact

On 4 July 2005 at 5:52 UT the Deep Impact mission successfully completed its goal to hit the nucleus of 9P/Tempel 1 with an impactor, forming a crater on the nucleus and ejecting material into the coma of the comet. NASA's Submillimeter Wave Astronomy Satellite (SWAS) observed the ₁₁₀-₁₀₁ ortho-water ground-state rotational transition in Comet 9P/Tempel 1 before, during, and after the impact. No excess emission from the impact was detected by SWAS and we derive an upper limit of 1.8×¹⁰⁷ kg on the water ice evaporated by the impact. However, the water production rate of the comet showed large natural variations of more than a factor of three during the weeks before and after the impact. Episodes of increased activity with alternated with periods with low outgassing (). We estimate that 9P/Tempel 1 vaporized a total of N∼4.5×¹⁰³⁴ water molecules (∼1.3×¹⁰⁹ kg) during June-September 2005. Our observations indicate that only a small fraction of the nucleus of Tempel 1 appears to be covered with active areas. Water vapor is expected to emanate predominantly from topographic features periodically facing the Sun as the comet rotates. We calculate that appreciable asymmetries of these features could lead to a spin-down or spin-up of the nucleus at observable rates.     

Experimental study of the sublimation of ice through an unconsolidated clay layer: Implications for the stability of ice on Mars and the possible diurnal variations in atmospheric water

We have studied the sublimation of ice and water vapor transport through various thicknesses of clay (<63 μm grain size). We experimentally demonstrate that both adsorption and diffusion strongly affect the transport of water, and that the processes of diffusion and adsorption can be separately quantified once the system comes to a steady state. At shallow depths of clay, water vapor transport is determined by diffusion through both the atmosphere and the clay layer, whereas at greater depth the rate of sublimation of the ice is governed only by diffusion through the clay. Using two different models, we determine the diffusion coefficient for water vapor through unconsolidated clay layer to be 1.08±0.04×¹⁰⁻⁴ and . We also determined the adsorption isotherms for the clay layer, which follow the Langmuir theory at low water vapor pressure (<100 Pa, where a monolayer of water molecules forms on the surface of the clay) and the BET theory at higher pressure (where multiple water layers form). From our analysis of both types of isotherms we determined the adsorption constants to be and c=30±10, respectively, and specific surface areas of 1.10±0.2×¹⁰⁵ and , respectively. Finally, we report a theoretical kinetic model for the simultaneous diffusion and adsorption from which we determine adsorption kinetic constants according to the Langmuir theory of and . If the martian regolith possesses diffusive properties similar to those of the unconsolidated montmorillonite soil we investigated here, it would not represent a significant barrier to the sublimation of subsurface ice. However, at the low subsurface temperatures of high latitude (180 K on average), ice could survive from the last glaciation period (about 300 to 400,000 years ago). Higher subsurface temperatures in the equatorial regions would prevent long-timescale survival of ice in the shallow subsurface. In agreement with previous work, we show that adsorption of water by a clay regolith could provide a significant reservoir of subsurface water and it might account for the purported diurnal cycle in the water content of the atmosphere.