Ben Rozitis

The internal structure of asteroid (25143) Itokawa as revealed by detection of YORP spin-up

Context. Near-Earth asteroid (25143) Itokawa was visited by the Hayabusa spacecraft in 2005, resulting in a highly detailed shape and surface topography model. This model has led to several predictions for the expected radiative torques on this asteroid, suggesting that its spin rate should be decelerating. Aims. To detect changes in rotation rate that may be due to YORP-induced radiative torques, which in turn may be used to investigate the interior structure of the asteroid. Methods. Through an observational survey spanning 2001 to 2013 we obtained rotational lightcurve data at various times over the last five close Earth-approaches of the asteroid. We applied a polyhedron-shape-modelling technique to assess the spin-state of the asteroid and its long term evolution. We also applied a detailed thermophysical analysis to the shape model determined from the Hayabusa spacecraft. Results. We have successfully measured an acceleration in Itokawa’s spin rate of dω/dt = (3.54 ± 0.38) × 10-8 rad day-2, equivalent to a decrease of its rotation period of ~45 ms year-1. From the thermophysical analysis we find that the centre-of-mass for Itokawa must be shifted by ~21 m along the long-axis of the asteroid to reconcile the observed YORP strength with theory. Conclusions. This can be explained if Itokawa is composed of two separate bodies with very different bulk densities of 1750 ± 110 kg m-3 and 2850 ± 500 kg m-3, and was formed from the merger of two separate bodies, either in the aftermath of a catastrophic disruption of a larger differentiated body, or from the collapse of a binary system. We therefore demonstrate that an observational measurement of radiative torques, when combined with a detailed shape model, can provide insight into the interior structure of an asteroid. Futhermore, this is the first measurement of density inhomogeneity within an asteroidal body, that reveals significant internal structure variation. A specialised spacecraft is normally required for this.

The nucleus of Comet 67P/Churyumov-Gerasimenko - A new shape model and thermophysical analysis

Context. Comet 67P/Churyumov-Gerasimenko is the target of the European Space Agency Rosetta spacecraft rendez-vous mission. Detailed physical characteristation of the comet before arrival is important for mission planning as well as providing a test bed for ground-based observing and data-analysis methods. Aims. To conduct a long-term observational programme to characterize the physical properties of the nucleus of the comet, via ground-based optical photometry, and to combine our new data with all available nucleus data from the literature. Methods. We applied aperture photometry techniques on our imaging data and combined the extracted rotational lightcurves with data from the literature. Optical lightcurve inversion techniques were applied to constrain the spin state of the nucleus and its broad shape. We performed a detailed surface thermal analysis with the shape model and optical photometry by incorporating both into the new Advanced Thermophysical Model (ATPM), along with all available Spitzer 8–24 μm thermal-IR flux measurements from the literature. Results. A convex triangular-facet shape model was determined with axial ratios b/a = 1.239 and c/a = 0.819. These values can vary by as much as 7% in each axis and still result in a statistically significant fit to the observational data. Our best spin state solution has Psid = 12.76137 ± 0.00006 h, and a rotational pole orientated at Ecliptic coordinates λ = 78◦(±10◦), β = +58◦(±10◦). The nucleus phase darkening behaviour was measured and best characterized using the IAU HG system. Best fit parameters are: G = 0.11 ± 0.12 and HR(1,1,0) = 15.31 ± 0.07. Our shape model combined with the ATPM can satisfactorily reconcile all optical and thermal-IR data, with the fit to the Spitzer 24 μm data taken in February 2004 being exceptionally good. We derive a range of mutually-consistent physical parameters for each thermal-IR data set, including effective radius, geometric albedo, surface thermal inertia and roughness fraction. Conclusions. The overall nucleus dimensions are well constrained and strongly imply a broad nucleus shape more akin to comet 9P/Tempel 1, rather than the highly elongated or “bi-lobed” nuclei seen for comets 103P/Hartley 2 or 8P/Tuttle. The derived low thermal inertia of <15 J m−2 K−1 s−1/2 is comparable with that measured for other comets scaled to similar heliocentric distances, and implies a surface regolith finer than lunar surface material.

Granular convection in microgravity

We investigate the role of gravity on convection in a dense granular shear flow. Using a microgravity-modified Taylor-Couette shear cell under the conditions of parabolic flight microgravity, we demonstrate experimentally that secondary, convective-like flows in a sheared granular material are close to zero in microgravity and enhanced under high-gravity conditions, though the primary flow fields are unaffected by gravity. We suggest that gravity tunes the frictional particle-particle and particle-wall interactions, which have been proposed to drive the secondary flow. In addition, the degree of plastic deformation increases with increasing gravitational forces, supporting the notion that friction is the ultimate cause.

The influence of global self-heating on the Yarkovsky and YORP effects

In addition to collisions and gravitational forces, there is a growing amount of evidence that photon recoil forces from the asymmetric reflection and thermal re-radiation of absorbed sunlight are primary mechanisms that are fundamental to the physical and dynamical evolution of small asteroids. The Yarkovsky effect causes orbital drift, and the Yarkovsky–O’Keefe–Radzievskii–Paddack (YORP) effect causes changes in the rotation rate and pole orientation. We present an adaptation of the Advanced Thermophysical Model to simultaneously predict the Yarkovsky and YORP effects in the presence of global self-heating that occurs within the large concavities of irregularly shaped asteroids, which has been neglected or dismissed in all previous models. It is also combined with rough surface thermal-infrared beaming effects, which have been previously shown to enhance the Yarkovsky orbital drift and dampen on average the YORP rotational acceleration by orders of several tens of per cent. Tests on all published concave shape models of near-Earth asteroids, and also on 100 Gaussian random spheres, show that the Yarkovsky effect is sensitive to shadowing and global self-heating effects at the few per cent level or less. For simplicity, Yarkovsky models can neglect these effects if the level of accuracy desired is of this order. Unlike the Yarkovsky effect, the YORP effect can be very sensitive to shadowing and global self-heating effects. Its sensitivity increases with decreasing relative strength of the YORP rotational acceleration, and does not appear to depend greatly on the degree of asteroid concavity. Global self-heating tends to produce a vertical offset in an asteroids YORP-rotational-acceleration versus obliquity curve which is in opposite direction to that produced by shadowing effects. It also ensures that at least one critical obliquity angle exists at which zero YORP rotational acceleration occurs. Global self-heating must be included for accurate predictions of the YORP effect if an asteroid exhibits a large shadowing effect. If global self-heating effects are not included, then it is found in ~75 per cent of cases that better predictions are produced when shadowing is also not included. Furthermore, global self-heating has implications for reducing the sensitivity of the YORP effect predictions to detailed variations in an asteroids shape model.

Physical and dynamical characterisation of the unbound asteroid pair 7343-154634

Context. Models have shown that asteroids can undergo fission if their rate of rotation is steadily increased. The forces acting to pull the asteroid apart exceed the material strength and gravitational force holding the asteroid together and material can escape from the surface of the asteroid. Initially forming a binary asteroid system, the components are capable of decoupling at low relative velocity from their mutual orbit if their mass ratio is less than 0.2. A number of asteroids with very similar orbital elements have been shown to have had very recent (<1 Myr) encounters at distances smaller than the Hill sphere radius of the larger of the asteroids. The mass ratio of the asteroids in each pair is estimated to be less than 0.2, suggesting that these unbound pairs are the result of rotational fission. Aims. We determine whether the asteroids in one such unbound pair, (7343) Ockeghem and (154 634) 2003 XX28, share a common composition, indicative of asteroids formed from a common parent and further constrain a likely formation age for this pair. Methods. We have obtained spectroscopic observations of each asteroid covering the wavelength range 0.45 to 1.0 microns. Using thermal observations we have measured the size and albedo of (7343) Ockeghem. Combined with optical lightcurve data of both asteroids, we have constrained the size and density of the asteroids and estimated the strength of the Yarkovsky force experienced by both. This improved physical information has been used in new dynamical simulations of the asteroids’ orbits to better constrain a formation time of this pair. Results. We find that the asteroids have very similar spectra consistent with an S-type taxonomy. The geometric albedo of (7343) Ockeghem, 0.20 ± 0.06 is consistent with this classification. The mass ratio range of the asteroids assuming an equal density, 0.007 to 0.065, is consistent with models of unbound asteroid pair formation. A new dynamical analysis has indicated that an absolute lower limit for the age of this pair is 400 kyr with a more likely age around 560 kyr, lower than a previous estimate of 800 kyr.

Physical characterisation of near-Earth asteroid (1620) Geographos - Reconciling radar and thermal-infrared observations

The Yarkovsky (orbital drift) and YORP (spin state change) effects play important roles in the dynamical and physical evolution of asteroids. Thermophysical modelling of these observed effects, and of thermal-infrared observations, allows a detailed physical characterisation of an individual asteroid to be performed. We perform a detailed physical characterisation of near-Earth asteroid (1620) Geographos, a potential meteor stream source and former spacecraft target, using the same techniques as previously used in Rozitis et al. (2013) for (1862) Apollo. We use the advanced thermophysical model (ATPM) on published light-curve, radar, and thermal-infrared observations to constrain the thermophysical properties of Geographos. The derived properties are used to make detailed predictions of the Yarkovsky orbital drift and YORP rotational acceleration, which are then compared against published measurements to determine Geographoss bulk density. We find that Geographos has a thermal inertia of 340 +140/-100 J m-2 K-1 s-1/2, a roughness fraction of >50%, and a bulk density of 2100 +550/-450 kg m-3 when using the light-curve-derived shape model with the radar-derived maximum equatorial diameter of 5.04 +/- 0.07 km. It is also found that the radar observations had overestimated the z-axis in Geographoss shape model because of their near-equatorial view. This results in a poor fit to the thermal-infrared observations if its effective diameter is kept fixed in the model fitting. The thermal inertia derived for Geographos is slightly higher than the typical values for a near-Earth asteroid of its size, and its derived bulk density suggests a rubble-pile interior structure. Large uncertainties in shape model z-axes are likely to explain why radar and thermal-infrared observations sometimes give inconsistent diameter determinations for other asteroids.

Granular shear flow in varying gravitational environments

Despite their very low surface gravities, asteroids exhibit a number of different geological processes involving granular matter. Understanding the response of this granular material subject to external forces in microgravity conditions is vital to the design of a successful asteroid sub-surface sampling mechanism, and in the interpretation of the fascinating geology on an asteroid. We have designed and flown a Taylor–Couette shear cell to investigate granular flow due to rotational shear forces under the conditions of parabolic flight microgravity. The experiments occur under weak compression. First, we present the technical details of the experimental design with particular emphasis on how the equipment has been specifically designed for the parabolic flight environment. Then, we investigate how a steady state granular flow induced by rotational shear forces differs in varying gravitational environments. We find that the effect of constant shearing on the granular material, in a direction perpendicular to the effective acceleration, does not seem to be strongly influenced by gravity. This means that shear bands can form in the presence of a weak gravitational field just as on Earth.

A thermophysical analysis of the (1862) Apollo Yarkovsky and YORP effects

Context. The Yarkovsky effect, which causes orbital drift, and the YORP effect, which causes changes in rotation rate and pole orientation, play important roles in the dynamical and physical evolution of asteroids. Near-Earth asteroid (1862) Apollo has strong detections of both orbital semimajor axis drift and rotational acceleration. Aims. To produce a unified model that can accurately match both observed effects using a single set of thermophysical properties derived from ground-based observations, and to determine Apollo’s long term evolution. Methods. We use light-curve shape inversion techniques and the Advanced Thermophysical Model (ATPM) on published light-curve, thermal-infrared, and radar observations to constrain Apollo’s thermophysical properties. The derived properties are used to make detailed predictions of Apollo’s Yarkovsky and YORP effects, which are then compared with published measurements of orbital drift and rotational acceleration. The ATPM explicitly incorporates 1D heat conduction, shadowing, multiple scattering of sunlight, global self-heating, and rough surface thermal-infrared beaming in the model predictions. Results. We find that ATPM can accurately reproduce the light-curve, thermal-infrared, and radar observations of Apollo, and simultaneously match the observed orbital drift and rotational acceleration using: a shape model with axis ratios of 1.94:1.65:1.00, an effective diameter of 1.55 ± 0.07 km, a geometric albedo of 0.20 ± 0.02, a thermal inertia of 140 +140-100 J m-2 K-1 s-1/2, a highly rough surface, and a bulk density of 2850 +480-680 kg m-3. Using these properties we predict that Apollo’s obliquity is increasing towards the 180 degree YORP asymptotic state at a rate of 1.5 +0.3-0.5 degrees per 105 yr. Conclusions. The derived thermal inertia suggests that Apollo has loose regolith material resting on its surface, which is consistent with Apollo undergoing a recent resurfacing event based on its observed Q-type spectrum. The inferred bulk density is consistent with those determined for other S-type asteroids, and suggests that Apollo has a fractured interior. The YORP effect is acting on a much faster timescale than the Yarkovsky effect and will dominate Apollo’s long term evolution. The ATPM can readily be applied to other asteroids with similar observational data sets.

Simulating regoliths in microgravity

Despite their very low surface gravities, the surfaces of asteroids and comets are covered by granular materials – regolith – that can range from a fine dust to a gravel-like structure of varying depths. Understanding the dynamics of granular materials is, therefore, vital for the interpretation of the surface geology of these small bodies and is also critical for the design and/or operations of any device planned to interact with their surfaces. We present the first measurements of transient weakening of granular material after shear reversal in microgravity as well as a summary of experimental results recently published in other journals, which may have important implications for small-body surfaces. Our results suggest that the force contact network within a granular material may be weaker in microgravity, although the influence of any change in the contact network is felt by the granular material over much larger distances. This could mean that small-body surfaces are even more unstable than previously imagined. However, our results also indicate that the consequences of, e.g., a meteorite impact or a spacecraft landing, may be very different depending on the impact angle and location, and depending on the prior history of the small-body surface.

The Surface Roughness of (433) Eros as Measured by Thermal-Infrared Beaming

In planetary science, surface roughness is regarded to be a measure of surface irregularity at small spatial scales, and causes the thermal-infrared beaming effect (i.e. re-radiation of absorbed sunlight back towards to the Sun). Typically, surface roughness exhibits a degeneracy with thermal inertia when thermophysical models are fitted to disc-integrated thermal-infrared observations of asteroids because of this effect. In this work, it is demonstrated how surface roughness can be constrained for near-Earth asteroid (433) Eros (i.e. the target of NASAs NEAR Shoemaker mission) when using the Advanced Thermophysical Model with thermal-infrared observations taken during an ‘almost pole-on’ illumination and viewing geometry. It is found that the surface roughness of (433) Eros is characterized by an rms slope of 38 ± 8° at the 0.5-cm spatial scale associated with its thermal-infrared beaming effect. This is slightly greater than the rms slope of 25 ± 5° implied by the NEAR Shoemaker laser ranging results when extrapolated to this spatial scale, and indicates that other surface shaping processes might operate, in addition to collisions and gravity, at spatial scales under one metre in order to make asteroid surfaces rougher. For other high-obliquity asteroids observed during ‘pole-on’ illumination conditions, the thermal-infrared beaming effect allows surface roughness to be constrained when the sub-solar latitude is greater than 60°, and if the asteroids are observed at phase angles of less than 40°. They will likely exhibit near-Earth asteroid thermal model beaming parameters that are lower than expected for a typical asteroid at all phase angles up to 100°.