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.     

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.     

Dark red debris from three short-period comets: 2P/Encke, 22P/Kopff, and 65P/Gunn

We present observations of the extended dust structures near the orbits of three short-period comets: 2P/Encke, 22P/Kopff, and 65P/Gunn. The dust trails were originally discovered by the Infrared Astronomical Satellite (IRAS). Our observations were made using wide-field optical CCD cameras on the University of Hawaii 2.24-m telescope, the Canada-France-Hawaii 3.6-m telescope, and the Kiso 1.05-m Schmidt telescope. We compared the observed images with models and found that the extended structures seen around 2P/Encke and 22P/Kopff before perihelion passage were most likely “dust trails,” whereas images taken after perihelion passage show a high contamination by recently released particles (i.e., particles in Neck-Line structures are visible). We could not confirm the existence of a dust trail from 65P/Gunn within the field of view of the camera used. The effective sizes of the particles responsible for the scattered light were estimated at 1-100 mm (2P/Encke), 1-10 mm (22P/Kopff), and 100 μm-1 mm (65P/Gunn), respectively, which is consistent with previous studies of dust trails made with infrared space telescopes and optical telescopes. We evaluated the mass loss rates of these comets, averaged over their orbits, as reaching (2P/Encke), (22P/Kopff), and (65P/Gunn). These values are consistent with previous work. Therefore, the total amount of material ejected from these three comets is , which would contribute a considerable fraction of the lost within 1 AU that needs to be replaced if the zodiacal cloud is to be maintained in a steady state. We also found that the particles in the dust structures are significantly redder than the Sun and the zodiacal light, and might be redder than the average short-period comet nuclei. Specifically, the reflectivity gradients of 2P/Encke, 22P/Kopff, and 65P/Gunn are 13±7 (% ¹⁰³ Å⁻¹), 20±5 (% ¹⁰³ Å⁻¹), and 15±4 (% ¹⁰³ Å⁻¹), respectively. We examined the change in color with distance from the nucleus. No clear correlation was detected for 2P/Encke or 22P/Kopff to an accuracy of 3-11%, while the 65P/Gunn tail did show color variation, becoming redder with increasing distance from the nucleus. This dark red material, consisting of particles of sand-cobble size, has marginally escaped from the nuclei and will evolve into finer-grained interplanetary dust particles after subsequent collisions.     

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 .     

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.     

Environmental predictions for the Beagle 2 lander, based on GCM climate simulations

The Mars climate database (MCD) is a database of statistics based on output from physically consistent numerical model simulations which describe the climate and surface environment of Mars. It is used here to predict the meteorological environment of the Beagle 2 lander site. The database was constructed directly on the basis of output from multiannual integrations of two general circulation models, developed jointly at Laboratoire de Météorologie Dynamique du Centre National de la Recherche Scientifique, France, and the University of Oxford, UK. In an atmosphere with dust opacities similar to that observed by Mars Global Surveyor, the predicted surface temperature at the time of landing (L_s=322°, 13:00 local time), is , and varying between ∼186 and over the Martian day. The predicted air temperature at above the surface, the height of the fully extended Beagle 2 robot arm, is at the time of landing. The expected mean wind near the surface on landing is north-westerly in direction, becoming more southerly over the mission. An increase in mean surface pressure is expected during the mission. Heavy global dust storm predictions are discussed; conditions which may only occur in the extreme as the expected time of landing is around the end of the main dust storm period. Past observations show approximately a one in five chance of a large-scale dust storm in a whole Mars year over the landing region, Isidis Planitia. This statistic results from observations of global, encircling, regional and local dust storms but does not include any small-scale dust “events” such as dust devils.     

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.     

Fragment-collision model for compound chondrule formation: Estimation of collision probability

We propose a new scenario for compound chondrule formation named as “fragment-collision model,” in the framework of the shock-wave heating model. A molten cm-sized dust particle (parent) is disrupted in the high-velocity gas flow. The extracted fragments (ejectors) are scattered behind the parent and the mutual collisions between them will occur. We modeled the disruption event by analytic considerations in order to estimate the probability of the mutual collisions assuming that all ejectors have the same radius. In the typical case, the molten thin () layer of the parent surface will be stripped by the gas flow. The stripped layer is divided into about 200 molten ejectors (assuming that the radius of ejectors is 300 μm) and then they are blown away by the gas flow in a short period of time (). The stripped layer is leaving from the parent with the velocity of depending on the viscosity, and we assumed that the extracted ejectors have a random velocity Δv of the same order of magnitude. Using above values, we can estimate the number density of ejectors behind the parent as . These ejectors occupy ∼9% of the space behind the parent in volume. Considering that the collision rate (number of collisions per unit time experienced by an ejector) is given by Rcoll=σcollneΔv, where σcoll is the cross-section of collision [e.g., Gooding, J.K., Keil, K., 1981. Meteoritics 16, 17-43], we obtain by substituting above values. Since most collisions occur within the short duration () before the ejectors are blown away, we obtain the collision probability of Pcoll∼0.36, which is the probability of collisions experienced by an ejector in one disruption event. The estimated collision probability is about one order of magnitude larger than the observed fraction of compound chondrules. In addition, the model predictions are qualitatively consistent with other observational data (oxygen isotopic composition, textural types, and size ratios of constituents). Based on these results, we concluded that this new model can be one of the strongest candidates for the compound chondrule formation. It should be noted that all collisions do not necessarily lead to the compound chondrule formation. The formation efficiency and the future works which should be investigated in the forthcoming paper are also discussed.     

Gravitational tidal waves in Titan's upper atmosphere

Tidal waves driven by Titan's orbital eccentricity through the time-dependent component of Saturn's gravitational potential attain nonlinear, saturation amplitudes (|T^′|>10 K, , and ) in the upper atmosphere (?500 km) due to the approximate exponential growth as the inverse square root of pressure. The gravitational tides, with vertical wavelengths of ∼100-150 km above 500 km altitude, carry energy fluxes sufficient in magnitude to affect the energy balance of the upper atmosphere with heating rates in the altitude range of 500-900 km.     

The effect of using different scale lengths on the production rates of Comet 46P/Wirtanen

The variations in production rates for Comet 46P/Wirtanen for the species H₂O and the parents of C₂ and CN are examined from the point of view of a variety of commonly used scale lengths. The calculations are carried out as a function of heliocentric distance. It is shown that, by using a common set of scale lengths, the results of various investigators can be brought into acceptable accord. The resulting production rates of H₂O and the parents of C₂ and CN versus heliocentric distance are recalculated and plotted versus the heliocentric distance r_H. The curves show reasonable agreement with a slope of ∼r_H⁻⁴. The water production rate near perihelion of 46P/Wirtanen is close to .     

Near infrared imaging spectroscopy of Venus with the Anglo-Australian Telescope

We have obtained full-disk spatially resolved spectra of the Venus nightside at near-infrared wavelengths during July 2007 using the Anglo-Australian Telescope and Infrared Imager and Spectrograph 2 (IRIS2). The data have been used to map the intensity and rotational temperature of the O₂(a¹Δ_g) airglow band at . The temperatures agree with those obtained in earlier IRIS2 observations and are significantly higher than expected from the Venus International Reference Atmosphere (VIRA) profile. We also report the detection of the corresponding ν=0-1O₂ airglow band at with a similar spatial distribution to the ν=0-0 band. Observations in the thermal window have been used to image surface topography using two different methods of cloud correction. We have also obtained images that can be used to study cloud motion.     

CCD photometry of 23 minor planets

We present CCD photometric observations of 23 main-belt asteroids, of which 8 have never been observed before; thus, the data of these objects are the first in the literature. The majority showed well-detectable light variations, exceeding 0^m1. We have determined synodic periods for 756 Lilliana (936), 1270 Datura (34), 1400 Tirela (1336), 1503 Kuopio (998), 3682 Welther (359), 7505 Furushu (414) and 11436 1969 QR (123), while uncertain period estimates were possible for 469 Argentina (123), 546 Herodias (104) and 1026 Ingrid (53). The shape of the lightcurves of 3682 Welther changed on a short time-scale and showed dimmings that might be attributed to eclipses in a binary system. For the remaining objects, only lower limits of the periods and amplitudes were concluded.     

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.     

Interaction between eddies and mean flow in Jupiter's atmosphere: Analysis of Cassini imaging data

Beebe et al. [Beebe, R.F., et al., 1980. Geophys. Res. Lett. 17, 1-4] and Ingersoll et al. [Ingersoll, A.P., et al., 1981. J. Geophys. Res. 86, 8733-8743] used images from Voyagers 1 and 2 to analyze the interaction between zonal winds and eddies in Jupiter's atmosphere. They reported a high positive correlation between Jupiter's eddy momentum flux, , and the variation of zonal velocity with latitude, . This correlation implied a surprisingly high rate of conversion of energy from eddies to zonal flow: , a value more than 10% of Jupiter's thermal flux emission. However, Sromovsky et al. [Sromovsky, L.A., et al., 1982. J. Atmos. Sci. 39, 1413-1432] argued that possible biases in the analysis could have caused an artificially high correlation. In addition, significant differences in the derived eddy flux between datasets put into question the robustness of any one result. We return to this long-standing puzzle using images of Jupiter from the Cassini flyby of December 2000. Our method is similar to previous analyses, but utilizes an automatic feature tracker instead of the human eye. The number of velocity vectors used in this analysis is over 200,000, compared to the 14,000 vectors used by Ingersoll et al. We also find a positive correlation between and and derive a global average power per unit mass, , ranging from . Utilizing Ingersoll et al.'s estimate of the mass per unit area involved in the transport, this would imply a rate of energy conversion of . We discuss the implications of this result and employ several tests to demonstrate its robustness.