Shape and spin distributions of asteroid populations from brightness variation estimates and large databases

Context. Many databases on asteroid brightnesses (e.g. ALCDEF, WISE) are potential sources for extensive asteroid shape and spin modelling. Individual lightcurve inversion models require several apparitions and hundreds of data points per target. However, we can analyse the coarse shape and spin distributions over populations of at least thousands of targets even if there are only a few points and one apparition per asteroid. This is done by examining the distribution of the brightness variations observed within the chosen population. Aims. Brightness variation has been proposed as a population-scale rather than individual-target observable in two studies so far. We aim to examine this approach rigorously to establish its theoretical validity, degree of ill-posedness, and practical applicability. Methods. We model the observed brightness variation of a target population by considering its cumulative distribution function (CDF) caused by the joint distribution function of two fundamental shape and spin indicators. These are the shape elongation and the spin latitude of a simple ellipsoidal model. The main advantage of the model is that we can derive analytical basis functions that yield the observed CDF as a function of the shape and spin distribution. The inverse problem can be treated linearly. Even though the inaccuracy of the model is considerable, databases of thousands of targets should yield some information on the distribution. Results. We establish the theoretical soundness and the typical accuracy limits of the approach both analytically and numerically. Using simulations, we derive a practical estimate of the model distribution in the (shape, spin)-plane. We show that databases such as Wide-field Infrared Survey Explorer (WISE) yield coarse but robust estimates of this distribution, and as an example compare various asteroid families with each other.Comment: 16 pages, 21 figures, manuscript accepted in Astronomy & Astrophysics, to be published in section 10. Planets and planetary system

Additional spectra of asteroid 1996 FG3, backup target of the ESA MarcoPolo-R mission

Near-Earth binary asteroid (175706) 1996 FG3 is the current backup target of the ESA MarcoPolo-R mission, selected for the study phase of ESA M3 missions. It is a primitive (C-type) asteroid that shows significant variation in its visible and near-infrared spectra. Here we present new spectra of 1996 FG3 and we compare our new data with other published spectra, analysing the variation in the spectral slope. The asteroid will not be observable again over the next three years at least. We obtained the spectra using DOLORES and NICS instruments at the Telescopio Nazionale Galileo (TNG), a 3.6m telescope located at El Roque de los Muchachos Observatory in La Palma, Spain. To compare with other published spectra of the asteroid, we computed the spectral slope S', and studied any plausible correlation of this quantity with the phase angle (alpha). In the case of visible spectra, we find a variation in spectral slope of Delta S' = 0.15 +- 0.10 %/10^3 A/degree for 3 < alpha < 18 degrees, in good agreement with the values found in the literature for the phase reddening effect. In the case of the near-infrared, we find a variation in the slope of Delta S' = 0.04 +- 0.08 %/10^3 A/degree for 6 < alpha < 51 degrees. Our computed variation in S' agrees with the only two values found in the literature for the phase reddening in the near-infrared. The variation in the spectral slope of asteroid 1996 FG3 shows a trend with the phase angle at the time of the observations, both in the visible and the near-infrared. It is worth noting that, to fully explain this spectral variability we should take into account other factors, like the position of the secondary component of the binary asteroid 1999 FG3 with respect to the primary, or the spin axis orientation at the time of the observations. More data are necessary for an analysis of this kind.Comment: 4 pages, 3 figures, Accepted in A&A 25 June 201

The Moon at thermal infrared wavelengths: A benchmark for asteroid thermal models

Thermal-infrared measurements of asteroids are crucial for deriving the objects' sizes, albedos, and also the thermophysical properties of the surface material. Depending on the available data, a range of simple to complex thermal models are applied to achieve specific science goals. However, testing these models is often a difficult process and the uncertainties of the derived parameters are not easy to estimate. Here, we make an attempt to verify a widely accepted thermophysical model (TPM) against unique thermal infrared (IR), full-disk, and well-calibrated measurements of the Moon. The data were obtained by the High-resolution InfraRed Sounder (HIRS) instruments on board a fleet of Earth weather satellites that serendipitously scan over the Moon. We found 22 Moon intrusions, taken in 19 channels between 3.75 micron and 15.0 micron, and over a wide phase angle range from -73.1 deg to +73.8 deg. The similarity between these Moon data and typical asteroid spectral-IR energy distributions allows us to benchmark the TPM concepts and to point out problematic aspects. The TPM predictions match the HIRS measurements within 5% (10% at the shortest wavelengths below 5 micron when using the Moon's known properties (size, shape, spin, albedo, thermal inertia, roughness) in combination with a newly established wavelength-dependent hemispherical emissivity. In the 5-7.5 micron and in the 9.5 to 11 micron ranges, the global emissivity model deviates considerably from the known lunar sample spectra. Our findings will influence radiometric studies of near-Earth and main-belt asteroids in cases where only short-wavelength data (from e.g., NEOWISE, the warm Spitzer mission, or ground-based M-band measurements) are available. The new, full-disk IR Moon model will also be used for the calibration of IR instrumentation on interplanetary missions (e.g., for Hayabusa-2) and weather satellites.Comment: 21 pages, 9 figures, 7 tables, accepted for publication in Astronomy & Astrophysics in March 202

The visible and near-infrared spectra of asteroids in cometary orbits

We study the visible and near-infrared (NIR) spectral properties of different ACO populations and compare them to the independently determined properties of comets. We select our ACOs sample based on published dynamical criteria and present our own observational results obtained using the 10.4m Gran Telescopio Canarias (GTC), the 4.2m William Herschel Telescope (WHT), the 3.56m Telescopio Nazionale Galileo (TNG), and the 2.5m Isaac Newton Telescope (INT), all located at the El Roque de los Muchachos Observatory (La Palma, Spain), and the 3.0m NASA Infrared Telescope Facility (IRTF), located at the Mauna Kea Observatory, in Hawaii. We include in the analysis the spectra of ACOs obtained from the literature. We derive the spectral class and the visible and NIR spectral slopes. We also study the presence of hydrated minerals by studying the 0.7 $\mu$m band and the UV-drop below 0.5 $\mu$m associated with phyllosilicates. We present new observations of 17 ACOs, 11 of them observed in the visible, 2 in the NIR and 4 in the visible and NIR. We also discuss the spectra of 12 ACOs obtained from the literature. All but two ACOs have a primitive-like class spectrum (X or D-type). Almost 100\% of the ACOs in long-period cometary orbits (Damocloids) are D-types. Those in Jupiter family comet orbits (JFC-ACOs) are $\sim$ 60\% D-types and $\sim$ 40\% X-types. The mean spectral slope $S'$ of JFC-ACOs is 9.7 $\pm$ 4.6 \%/1000 \AA \ and for the Damocloids this is 12.2 $\pm$ 2.0 \%/1000 \AA . No evidence of hydration on the surface of ACOs is found from their visible spectra. The slope and spectral class distribution of ACOs is similar to that of comets. The spectral classification, the spectral slope distribution of ACOs, and the lack of spectral features indicative of the presence of hydrated minerals on their surface, strongly suggest that ACOs are likely dormant or extinct comets.Comment: 11 pages, 10 Figures, published in A&

Near-Earth asteroid (3200) Phaethon. Characterization of its orbit, spin state, and thermophysical parameters

The near-Earth asteroid (3200) Phaethon is an intriguing object: its perihelion is at only 0.14 au and is associated with the Geminid meteor stream. We aim to use all available disk-integrated optical data to derive a reliable convex shape model of Phaethon. By interpreting the available space- and ground-based thermal infrared data and Spitzer spectra using a thermophysical model, we also aim to further constrain its size, thermal inertia, and visible geometric albedo. We applied the convex inversion method to the new optical data obtained by six instruments and to previous observations. The convex shape model was then used as input for the thermophysical modeling. We also studied the long-term stability of Phaethon's orbit and spin axis with a numerical orbital and rotation-state integrator. We present a new convex shape model and rotational state of Phaethon: a sidereal rotation period of 3.603958(2) h and ecliptic coordinates of the preferred pole orientation of (319$^{\circ}$, $-$39$^{\circ}$) with a 5$^{\circ}$ uncertainty. Moreover, we derive its size ($D$=5.1$\pm$0.2 km), thermal inertia ($\Gamma$=600$\pm$200 J m$^{-2}$ s$^{-1/2}$ K$^{-1}$), geometric visible albedo ($p_{\mathrm{V}}$=0.122$\pm$0.008), and estimate the macroscopic surface roughness. We also find that the Sun illumination at the perihelion passage during the past several thousand years is not connected to a specific area on the surface, which implies non-preferential heating.Comment: Astronomy and Astrophysics. In pres

Haumea's thermal emission revisited in the light of the occultation results

A recent multi -chord occultation measurement of the dwarf planet (136108) Haumea (Ortiz et al., 2017) revealed an elongated shape with the longest axis comparable to Pluto's mean diameter. The chords also indicate a ring around Haumea's equatorial plane, where its largest moon, Hi'iaka, is also located. The Haumea occultation size estimate (size of an equal-volume sphere(1) D-equ = 1595 km) is larger than previous radiometric solutions (equivalent sizes in the range between 1150 and 1350 km), which lowers the object's density to about 1.8 g/cm3, a value closer to the densities of other large TNOs. We present unpublished and also reprocessed Herschel and Spitzer mid- and far-infrared measurements. We compare 100 and 160 pm thermal lightcurve amplitudes - originating from Haumea itself - with models of the total measured system fluxes (ring, satellite, Haumea) from 24-350 pm. The combination with results derived from the occultation measurements allows us to reinterpret the object's thermal emission. Our radiometric studies show that Haumea's crystalline water ice surface must have a thermal inertia of about 5 J K-3 m-s(-1) (combined with a root mean square of the surface slopes of 0.2). We also have indications that the satellites (at least Hi'iaka) must have high geometric albedos > 0.5, otherwise the derived thermal amplitude would be inconsistent with the total measured system fluxes at 24, 70, 100, 160, 250, and 350 pm. The high albedos imply sizes of about 300 and 150 km for Hi'iaka and Namaka, respectively, indicating unexpectedly high densities > 1.0 g cm(-3) for TNOs this small, and the assumed collisional formation from Haumea's icy crust. We also estimated the thermal emission of the ring for the time period 1980-2030, showing that the contribution during the Spitzer and Herschel epochs was small, but not negligible. Due to the progressive opening of the ring plane, the ring emission will be increasing in the next decade when JWST is operational. In the MIRI 25.5 pm band it will also be possible to obtain a very high-quality thermal lightcurve to test the derived Haumea properties

The non-convex shape of (234) Barbara, the first Barbarian

Asteroid (234) Barbara is the prototype of a category of asteroids that has been shown to be extremely rich in refractory inclusions, the oldest material ever found in the Solar System. It exhibits several peculiar features, most notably its polarimetric behavior. In recent years other objects sharing the same property (collectively known as "Barbarians") have been discovered. Interferometric observations in the mid-infrared with the ESO VLTI suggested that (234) Barbara might have a bi-lobated shape or even a large companion satellite. We use a large set of 57 optical lightcurves acquired between 1979 and 2014, together with the timings of two stellar occultations in 2009, to determine the rotation period, spin-vector coordinates, and 3-D shape of (234) Barbara, using two different shape reconstruction algorithms. By using the lightcurves combined to the results obtained from stellar occultations, we are able to show that the shape of (234) Barbara exhibits large concave areas. Possible links of the shape to the polarimetric properties and the object evolution are discussed. We also show that VLTI data can be modeled without the presence of a satellite.Comment: 10 pages, 6 figure