Combining asteroid models derived by lightcurve inversion with asteroidal occultation silhouettes

Asteroid sizes can be directly measured by observing occultations of stars by asteroids. When there are enough observations across the path of the shadow, the asteroid’s projected silhouette can be reconstructed. Asteroid shape models derived from photometry by the lightcurve inversion method enable us to predict the orientation of an asteroid for the time of occultation. By scaling the shape model to fit the occultation chords, we can determine the asteroid size with a relative accuracy of typically ∼10%. We combine shape and spin state models of 44 asteroids (14 of them are new or updated models) with the available occultation data to derive asteroid effective diameters. In many cases, occultations allow us to reject one of two possible pole solutions that were derived from photometry. We show that by combining results obtained from lightcurve inversion with occultation timings, we can obtain unique physical models of asteroids.     

R-band light curve of Comet 9P/Tempel 1 during the Deep Impact event

We present high-speed CCD photometry of Comet 9P/Tempel 1 during the Deep Impact event on 2005 July 4 UT. Approximately 2 h and 50 min of R-band data were acquired at Mount Laguna Observatory with a temporal resolution of 5.5 s. The flux increased by 9% in the first minute after impact. This was followed by a more gradual two-part linear rise, with a change in slope at 9.2 min post-impact, at which time the rate of brightening increased from ∼ to ∼. An analysis of the light curve obtained with the guide camera on the United Kingdom Infrared Telescope and yields very similar results. These findings are mildly in disagreement with the 3-part linear rise found by Fernández et al. (2007) in that we do not find any evidence for a change at 4 min post-impact. We interpret the linear rise phase as due to solar illumination of the edge of an expanding optically thick dust ejecta plume. After approximately 20 min, the light curves begin to flatten out, perhaps coincident with the start of the transition to becoming optically thin. In the large apertures (>10^″) the light curve continues to gradually rise until the end of the observations. In smaller apertures, the light curves reach a peak at approximately 50 min, then decrease back towards the pre-impact flux level. The drop in flux in the smaller apertures may be caused by the ejecta expanding beyond the edge of the photometric aperture, and if so, we can use this timescale to infer an expansion velocity of ∼, consistent with previous published estimates.     

Probable swirls detected as photometric anomalies in Oceanus Procellarum

Images of the lunar nearside obtained by telescopes of Maidanak Observatory (Uzbekistan) and Simeiz Observatory (Crimea, Ukraine) equipped with Canon CMOS cameras and Sony CCD LineScan camera were used to study photometric properties of the lunar nearside in several spectral bands. A wide range of lunar phase angles was covered, and the method of phase ratios to assess the steepness of the phase function at different phase angles is applied. We found several areas with photometric anomalies in the south-west portion of the lunar disk that we refer to as Oceanus Procellarum anomalies. The areas being unique on the lunar nearside do not obey the inverse correlation between albedo and phase-curve slope, demonstrating high phase-curve slopes at intermediate albedo. Low-Sun images acquired with Lunar Orbiter IV and Apollo-16 cameras do not reveal anomalous topography of the regions, at least for scales larger than several tens of meters. The areas also do not have any thermal inertia, radar (70 and 3.8 cm), magnetic, or chemical/mineral peculiarities. On the other hand they exhibit a polarimetric signature that we interpret to be due to the presence of a porous regolith upper layer consisting of dust particles. The anomalies may be interpreted as regions of very fresh shallow regolith disturbances caused by impacts of meteoroid swarms consisting of rather small impactors. This origin is similar to one of the hypotheses for the origin of lunar swirls like the Reiner-γ formation. The photometric difference between the shallow and pervasive (Reiner-γ class) swirls is that the latter appear to have a significant amount of immature soils in the upper surface layers.     

(47171) 1999 TC36, A transneptunian triple

We present new analysis of HST images of (47171) 1999 TC₃₆ that confirm it as a triple system. Fits to the point-spread function (PSF) consistently show that the apparent primary is itself composed of two similar-sized components. The two central components, A1 and A2, can be consistently identified in each of nine epochs spread over 7 years of time. In each instance, the component separation, ranging from 0.023 ± 0.002 to 0.031 ± 0.003 arcsec, is roughly one half of the Hubble Space Telescope’s diffraction limit at 606 nm. The orbit of the central pair has a semi-major axis of a ∼ 867 km with a period of P ∼ 1.9 days. These orbital parameters yield a system mass that is consistent with M_sys = 12.75 ± 0.06 × 10¹⁸ kg derived from the orbit of the more distant secondary, component B. The diameters of the three components are . The relative sizes of these components are more similar than in any other known multiple in the Solar System. Taken together, the diameters and system mass yield a bulk density of . HST photometry shows that component B is variable with an amplitude of ?0.17 ± 0.05 magnitudes. Components A1 and A2 do not show variability larger than 0.08 ± 0.03 magnitudes approximately consistent with the orientation of the mutual orbit plane and tidally distorted equilibrium shapes. The system has high specific angular momentum of J/J′ = 0.93, comparable to most of the known transneptunian binaries.     

Determination of a precise rotation period for the Deep Space 1 target, Comet 19P/Borrelly

The rotation period derived by Mueller and Samarasinha (Mueller, B.E.A., Samarasinha, N.H. [2002]. Earth Moon Planets 90, 463-471) of the Deep Space 1 (DS1) mission target, Comet 19P/Borrelly, using ground-based data from July 28 to August 1, 2000, is improved by two orders of magnitude. This precision is reached in a multistep process.Combining all available ground-based data in 2000 decreases the error by an order of magnitude. Next, assuming that the rotation period did not change between 2000 and 2001, constraints from the HST 2001 data (Weaver, H.A., Stern, S.A., Parker, J.Wm. [2003]. Astron. J. 126, 444-451) yield three possible rotation periods: P = 1.088 ± 0.003 days, P = 1.108 ± 0.002 days, and P = 1.135 ± 0.003 days, which are consistent with our initial derivation of P = 1.08 ± 0.04 days (Mueller, B.E.A., Samarasinha, N.H. [2002]. Earth Moon Planets 90, 463-471).These three periods are further refined and the error bars further improved by another order of magnitude by linking the combined ground-based data from 2000 to the nuclear orientation of Borrelly at the DS1 encounter in 2001 (see Table 2). Due to aliasing, there are seven possible rotation periods around P = 1.088 days, five possible periods around P = 1.108 days, and six possible periods around P = 1.135 days, with precisions of the order of 0.0002 days (≈17 s).     

A three-parameter magnitude phase function for asteroids

We develop a three-parameter H, G₁, G₂ magnitude phase function for asteroids starting from the current two-parameter H, G phase function. We describe stochastic optimization of the basis functions of the magnitude phase function based on a carefully chosen set of asteroid photometric observations covering the principal types of phase dependencies. We then illustrate the magnitude phase function with a chosen set of observations. It is shown that the H, G₁, G₂ phase function systematically improves fits to the existing data and considerably so, warranting the utilization of three parameters instead of two. With the help of the linear three-parameter phase function, we derive a nonlinear two-parameter H, G₁₂ phase function, and demonstrate its applicability in predicting phase dependencies based on small numbers of observations.     

A new look at photometry of the Moon

We use ROLO photometry (Kieffer, H.H., Stone, T.C. [2005]. Astron. J. 129, 2887-2901) to characterize the before and after full Moon radiance variation for a typical highlands site and a typical mare site. Focusing on the phase angle range 45° < α < 50°, we test two different physical models, macroscopic roughness and multiple scattering between regolith particles, for their ability to quantitatively reproduce the measured radiance difference. Our method for estimating the rms slope angle is unique and model-independent in the sense that the measured radiance factor I/F at small incidence angles (high Sun) is used as an estimate of I/F for zero roughness regolith. The roughness is determined from the change in I/F at larger incidence angles. We determine the roughness for 23 wavelengths from 350 to 939 nm. There is no significant wavelength dependence. The average rms slope angle is 22.2° ± 1.3° for the mare site and 34.1° ± 2.6° for the highland site. These large slopes, which are similar to previous “photometric roughness” estimates, require that sub-mm scale “micro-topography” dominates roughness measurements based on photometry, consistent with the conclusions of Helfenstein and Shepard (Helfenstein, P., Shepard, M.K. [1999]. Icarus 141, 107-131). We then tested an alternative and very different model for the before and after full Moon I/F variation: multiple scattering within a flat layer of realistic regolith particles. This model consists of a log normal size distribution of spheres that match the measured distribution of particles in a typical mature lunar soil 72141,1 (McKay, D.S., Fruland, R.M., Heiken, G.H. [1974]. Proc. Lunar Sci. Conf. 5, Geochim. Cosmochim. Acta 1 (5), 887-906). The model particles have a complex index of refraction 1.65-0.003i, where 1.65 is typical of impact-generated lunar glasses. Of the four model parameters, three were fixed at values determined from Apollo lunar soils: the mean radius and width of the log normal size distribution and the real part of the refraction index. We used FORTRAN programs from Mishchenko et al. (Mishchenko, M.I., Dlugach, J.M., Yanovitskij, E.G., Zakharova, N.T. [1999]. J. Quant. Spectrosc. Radiat. Trans. 63, 409-432; Mishchenko, M.I., Travis, L.D., Lacis, A.A. [2002]. Scattering, Absorption and Emission of Light by Small Particles. Cambridge Univ. Press, New York. <http://www.giss.nasa.gov/staff/mmishchenko/books.html>) to calculate the scattering matrix and solve the radiative transfer equation for I/F. The mean single scattering albedo is ω = 0.808, the asymmetry parameter is 〈cos Θ〉 = 0.77 and the phase function is very strongly peaked in both the forward and backward scattering directions. The fit to the observations for the highland site is excellent and multiply scattered photons contribute ?80% of I/F. We conclude that either model, roughness or multiple scattering, can match the observations, but that the strongly anisotropic phase functions of realistic particles require rigorous calculation of many orders of scattering or spurious photometric roughness estimates are guaranteed. Our multiple scattering calculation is the first to combine: (1) a regolith model matched to the measured particle size distribution and index of refraction of the lunar soil, (2) a rigorous calculation of the particle phase function and solution of the radiative transfer equation, and (3) application to lunar photometry with absolute radiance calibration.     

Titan haze distribution and optical properties retrieved from recent observations

The VIMS instrument onboard Cassini observed the north polar region of Titan at 113° phase angle, 28 December 2006. On this spectral image, a vast polar cloud can be seen northward to 62°N, and elsewhere, the haze appears as the dominant source of scattering. Because the surface does not appear in the wavelength range between 0.3 and , this spectro-image is ideal to study airborn scatterers both in methane bands and windows. In this work, we study this image, along with another image taken at 13° phase angle. This image probe both the atmosphere and the surface from pole to pole. First, we characterise the spatial distribution of the haze layer above 100 km between 80°S and 70°N. We find a north south asymmetry with a haze opacity increasing by a factor 3 from the south pole to the equator, then a constant value up to about 30°N and a decrease of a factor 2 between 30°N and about 60°N. Beyond 60°N, we can see the influence of the north polar cloud, even in the band, but no polar haze accumulation. The fact that the north polar region is still in the polar night is a possible explanation. No haze accumulation is observed in the southern polar region. Secondly, we partly identify the origin of spectral features in the 2.8-μm methane window, which are found to be due to deuterated methane (CH₃D). This allows the analyse of this window and to retrieve the opacity of scatterers layer below 80 km (presumably made of aerosols and condensate droplets) between 35°N and 60°N. Finally, we constrained the values and the spectral behaviour of the imaginary part of the aerosol refractive index in the range between 0.3 and . To do so, we studied the 2.8-μm window with the image taken at 113° phase angle. To complete the analysis, we studied the transmission through the haze layer in the 3.4-μm band observed in solar occultation mode with VIMS, and we analysed the single scattering albedo retrieved with DISR instrument between 0.4 and . The imaginary part of the refractive index that we find for Titan aerosols follows Khare et al. (Khare, B.N. et al. [1984]. Icarus 60, 127-137) optical constant up to and becomes constant beyond this wavelength at least up to . It also has a prominent peak at and a secondary peak at , which indicates material rich in C-H bonds, with much less N-H bonds than in Khare et al. (1984) tholins.     

Sunskirting comets discovered with the LASCO coronagraphs over the decade 1996–2008

In addition to an unprecedented number of Kreutz sungrazing comets, the LASCO coronagraphs have discovered some 238 unrelated “sunskirting” comets over the 12 years from 1996 to 2008. This new class is organized in several groups, and at least two comets have further been found periodic. This article presents the photometry and the heliocentric light curves of these 238 sunskirting comets. The bulk of them exhibit a continuous increase of the brightness as the comet approaches the Sun, reach a peak before perihelion and then progressively fade with a large variety of brightness gradients. However some of them have peak brightness either at or post-perihelion, whereas a quite large number are approximately flat. Likewise for the sungrazers, we find a color effect prominent between 8 and 40R_⊙ (solar radii) which we interpret as resulting from the emission lines of the Na I doublet (D lines). We finally characterize the different groups of sunskirters on the basis of their cumulative distribution function of the peak brightness and of their fragmentation history.     

A tale of two very different comets: ISO and MSX measurements of dust emission from 126P/IRAS (1996) and 2P/Encke (1997)

We present the characteristics of the dust comae of two comets, 126P/IRAS, a member of the Halley family (a near-isotropic comet), and 2P/Encke, an ecliptic comet. We have primarily used mid- and far-infrared data obtained by the ISOPHOT instrument aboard the Infrared Space Observatory (ISO) in 1996 and 1997, and mid-infrared data obtained by the SPIRIT III instrument aboard the Midcourse Space Experiment (MSX) in 1996. We find that the dust grains emitted by the two comets have markedly different thermal and physical properties. P/IRAS's dust grain size distribution appears to be similar to that of fellow family member 1P/Halley, with grains smaller than 5 microns dominating by surface area, whereas P/Encke emits a much higher fraction of big (20 μm and higher) grains, with the grain mass distribution being similar to that which is inferred for the interplanetary dust population. P/Encke's dearth of micron-scale grains accounts for its visible-wavelength classification as a “gassy” comet. These conclusions are based on analyses of both imaging and spectrophotometry of the two comets; this combination provides a powerful way to constrain cometary dust properties. Specifically, P/IRAS was observed preperihelion while 1.71 AU from the Sun, and seen to have a 15-arcmin long mid-infrared dust tail pointing in the antisolar direction. No sunward spike was seen despite the vantage point being nearly in the comet's orbital plane. The tail's total mass at the time was about 8×10⁹ kg. The spectral energy distribution (SED) is best fit by a modified greybody with temperature T=265±15 K and emissivity ε proportional to a steep power law in wavelength λ: ε∝λ^−α, where α=0.50±0.20(2σ). This temperature is elevated with respect to the expected equilibrium temperature for this heliocentric distance. The dust mass loss rate was between 150-600 kg/s (95% confidence), the dust-to-gas mass loss ratio was about 3.3, and the albedo of the dust was 0.15±0.03. Carbonaceous material is depleted in the comet's dust by a factor of 2-3, paralleling the C₂ depletion in P/IRAS's gas coma. P/Encke, on the other hand, observed while 1.17 AU from the Sun, had an SED that is best fit by a Planck function with T=270±15 K and no emissivity falloff. The dust mass loss rate was 70-280 kg/s (95% confidence), the dust-to-gas mass loss ratio was about 2.3, and the albedo of the dust was about 0.06±0.02. These conclusions are consistent with the strongly curved dust tail and bright dust trail seen by Reach et al. (2000; Icarus 148, 80) in their ISO 12-μm imaging of P/Encke. The observed differences in the P/IRAS and P/Encke dust are most likely due to the less evolved and insolated state of the P/IRAS nuclear surface. If the dust emission behavior of P/Encke is typical of other ecliptic comets, then comets are the major supplier of the interplanetary dust cloud.