Advanced Thermophysical Model Research Articles

High‐Resolution Thermophysical Analysis of the OSIRIS‐REx Sample Site and Three Other Regions of Interest on Bennu

AbstractThe OSIRIS‐REx (Origins, Spectral Interpretation, Resource Identification, and Security–Regolith Explorer) spacecraft sampled asteroid (101955) Bennu on 20 October 2020 and will return the collected regolith to Earth in 2023. Before sample collection, spectral observations of four regions of interest on Bennu's surface were acquired at high spatial resolution (2–9 m per spectrometer spot) to identify the most suitable site for sampling and provide contextual information for the returned sample. In this study, we investigate thermal‐infrared (6–50 μm) observations of these four regions, including the site that OSIRIS‐REx ultimately sampled, using the Advanced Thermophysical Model with input digital terrain models derived from laser altimetry. From model‐to‐measurement comparisons, we find that the observed brightness temperatures depend strongly on small‐scale topography, local variations in thermal inertia, and the observation phase angle. Thermal inertia mapping reveals spatial variations that distinguish the different boulder types found on Bennu. A boulder bearing carbonate veins has higher thermal inertia than average, suggesting that cementation processes reduced its porosity. The thermal inertia of the site sampled is 190 ± 30 J m−2 K−1 s−1/2, which is consistent with observations of a fine‐grained regolith mixed with porous rocks. Thermophysical modeling of the site sampled predicts that the maximum temperatures experienced by the collected sample while on Bennu were 357 ± 3 and 261 ± 3 K for the surface and 50 cm depth, respectively. We predict that OSIRIS‐REx will return a sample with thermophysical properties unique from those of meteorites.

Thermophysical Modeling of 20 Themis Family Asteroids with WISE/NEOWISE Observations

Abstract Themis family is one of the largest and oldest asteroid populations in the main belt. Water ice may widely exist on the parent body (24) Themis. In this work, we employ the Advanced Thermophysical Model as well as midinfrared measurements from NASA’s Wide-field Infrared Survey Explorer to explore thermal parameters of 20 Themis family members. Here we show that the average thermal inertia and geometric albedo are 39.5 ± 26.0 J m−2 s−1/2 K−1 and 0.067 ± 0.018, respectively. The family members have a relatively moderate roughness fraction on their surfaces. We find that the relatively low albedos of Themis members are consistent with the typical values of B-type and C-type asteroids. As mentioned, the Themis family bears a very low thermal inertia, which indicates a fine and mature regolith on their surfaces. The resemblance of thermal inertia and geometric albedo of Themis members may reveal their close connection in origin and evolution. In addition, we present the compared results of thermal parameters for several prominent families.

Implications for Ice Stability and Particle Ejection From High‐Resolution Temperature Modeling of Asteroid (101955) Bennu

AbstractThe finding by the OSIRIS‐REx (Origins, Spectral Interpretation, Resource Identification, and Security‐Regolith Explorer) mission that its target (101955) Bennu is an active asteroid has raised questions as to whether the observed particle ejection events are driven by temperature. To investigate sublimation of water ice and rock thermal fracture as possible temperature‐driven causes, we modeled the global temperatures of Bennu and searched for correlations with the identified ejection points on the asteroid surface. We computed temperatures with the Advanced Thermophysical Model and the 75‐cm‐resolution global shape model of Bennu derived by the OSIRIS‐REx mission. We find that ~1,856 m2 of Bennu's polar regions have orbit‐averaged temperatures that are sufficiently cold to enable water ice, if buried within the top few meters of the surface, to remain stable over geological timescales. Millimeter thick layers of surface water ice are also stable over ~103‐year timescales within polar centimeter‐scale cold traps. However, we do not find evidence of conditions enabling ice stability in the warmer equatorial regions, where ejection events have been observed, implying that sublimation of water ice is not the cause of particle ejection. Conversely, rock thermal fracture remains a possible mechanism of particle ejection. We find high amplitudes of diurnal temperature variation, a proxy for the efficacy of thermal fracturing, at all latitudes on Bennu due to its extreme ruggedness. Therefore, if rock thermal fracture is the mechanism, particles could be ejected from any latitude, which is consistent with the continued observations of particle ejection by OSIRIS‐REx.

Determination of Size, Albedo, and Thermal Inertia of 10 Vesta Family Asteroids with WISE/NEOWISE Observations

Abstract In this work, we investigate the size, thermal inertia, surface roughness, and geometric albedo of 10 Vesta family asteroids using the Advanced Thermophysical Model, based on the thermal-infrared data acquired by mainly NASA’s Wide-field Infrared Survey Explorer. Here, we show that the average thermal inertia and geometric albedo of the investigated Vesta family members are 42 J m−2 s−1/2 K−1 and 0.314, respectively, where the derived effective diameters are less than 10 km. Moreover, the family members have a relatively low roughness fraction on their surfaces. The similarity in thermal inertia and geometric albedo among the V-type Vesta family members may reveal their close connection in origin and evolution. As the fragments of the cratering event of Vesta, the family members may have undergone a similar evolutionary process, thereby leading to very close thermal properties. Finally, we estimate their regolith grain sizes with different volume filling factors.

Observing the variation of asteroid thermal inertia with heliocentric distance

Thermal inertia is a useful property to characterise a planetary surface since it can be used as a qualitative measure of the regolith grain size. It is expected to vary with heliocentric distance because of its dependence on temperature. However, no previous investigation has conclusively observed a change in thermal inertia for any given planetary body. We have addressed this by using NEOWISE data and the Advanced Thermophysical Model to study the thermophysical properties of the near-Earth asteroids (1036) Ganymed, (1580) Betulia, and (276049) 2002 CE26 as they moved around their highly eccentric orbits. We confirm that the thermal inertia values of Ganymed and 2002 CE26 do vary with heliocentric distance, although the degree of variation observed depends on the spectral emissivity assumed in the thermophysical modelling. We also confirm that the thermal inertia of Betulia did not change for three different observations obtained at the same heliocentric distance. Depending on the spectral emissivity, the variations for Ganymed and 2002 CE26 are potentially more extreme than that implied by theoretical models of heat transfer within asteroidal regoliths, which might be explained by asteroids having thermal properties that also vary with depth. Accounting for this variation reduces a previously observed trend of decreasing asteroid thermal inertia with increasing size, and suggests that the surfaces of small and large asteroids could be much more similar than previously thought. Furthermore, this variation can affect Yarkovsky orbital drift predictions by a few tens of per cent.

Thermophysical characteristics of the large main-belt asteroid (349) Dembowska

(349) Dembowska, a large, bright main-belt asteroid, has a fast rotation and oblique spin axis. It may have experienced partial melting and differentiation. We constrain Dembowska's thermophysical properties, e.g., thermal inertia, roughness fraction, geometric albedo and effective diameter within 3$\sigma$ uncertainty of $\Gamma=20^{+12}_{-7}\rm~Jm^{-2}s^{-0.5}K^{-1}$, $f_{\rm r}=0.25^{+0.60}_{-0.25}$, $p_{\rm v}=0.309^{+0.026}_{-0.038}$, and $D_{\rm eff}=155.8^{+7.5}_{-6.2}\rm~km$, by utilizing the Advanced Thermophysical Model (ATPM) to analyse four sets of thermal infrared data obtained by IRAS, AKARI, WISE and Subaru/COMICS at different epochs. In addition, by modeling the thermal lightcurve observed by WISE, we obtain the rotational phases of each dataset. These rotationally resolved data do not reveal significant variations of thermal inertia and roughness across the surface, indicating the surface of Dembowska should be covered by a dusty regolith layer with few rocks or boulders. Besides, the low thermal inertia of Dembowska show no significant difference with other asteroids larger than 100 km, indicating the dynamical lives of these large asteroids are long enough to make the surface to have sufficiently low thermal inertia. Furthermore, based on the derived surface thermophysical properties, as well as the known orbital and rotational parameters, we can simulate Dembowska's surface and subsurface temperature throughout its orbital period. The surface temperature varies from $\sim40$ K to $\sim220$ K, showing significant seasonal variation, whereas the subsurface temperature achieves equilibrium temperature about $120\sim160$ K below $30\sim50$ cm depth.

Surface thermophysical properties on the potentially hazardous asteroid (99942) Apophis

We investigate the surface thermophysical properties (thermal emissivity, thermal inertia, roughness fraction and geometric albedo) of asteroid (99942) Apophis, using the currently available mid-infrared observations from CanariCam on Gran Telescopio CANARIAS and far-infrared data from PACS on Herschel, based on the Advanced Thermophysical Model. We show that the thermal emissivity of Apophis should be wavelength dependent from 8.70 μm to 160 μm, and the maximum emissivity may appear around 20 μm, similar to that of Vesta. Moreover, we further derive the thermal inertia, roughness fraction, geometric albedo and effective diameter of Apophis within a possible 1 scale of , , and , and 3 scale of , , and . The derived low thermal inertia but high roughness fraction may imply that Apophis could have regolith on its surface, where stronger space weathering but weaker regolith migration has happened in comparison with asteroid Itokawa. Our results show that small-size asteroids could also have fine regolith on the surface, and further infer that Apophis may have been delivered from the Main Belt by the Yarkovsky effect.

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 NASA's 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 characterised by an RMS slope of 38 $\pm$ 8{\deg} at the 0.5-cm spatial scale associated with its thermal-infrared beaming effect. This is slightly greater than the RMS slope of 25 $\pm$ 5{\deg} 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{\deg}, and if the asteroids are observed at phase angles of less than 40{\deg}. They will likely exhibit NEATM beaming parameters that are lower than expected for a typical asteroid at all phase angles up to 100{\deg}.

Physical characterisation of near-Earth asteroid (1620) Geographos

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 Geographos's 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 Geographos's 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.

Shape, thermal and surface properties determination of a candidate spacecraft target asteroid (175706) 1996 FG3

In this paper, a 3D convex shape model of (175706) 1996 FG3, which consists of 2040 triangle facets and 1022 vertices, is derived from the known lightcurves. The best-fit orientation of the asteroid's spin axis is determined to be $\lambda =237.7^\circ$ and $\beta=-83.8^{\circ}$ considering the observation uncertainties, and its rotation period is $\sim$ 3.5935 h . Using the derived shape model, we adopt the so-called advanced thermophysical model (ATPM) to fit three published sets of mid-infrared observations of 1996 FG3 \citep{Wolters2011,Walsh2012}, so as to evaluate its surface properties. Assuming the primary and the secondary bear identical shape, albedo, thermal inertia and surface roughness, the best-fit parameters are obtained from the observations. The geometric albedo and effective diameter of the asteroid are reckoned to be $p_{\rm v}=0.045\pm0.002$, $D_{\rm eff}=1.69^{+0.05}_{-0.02}$ km. The diameters of the primary and secondary are determined to be $D_{1}=1.63^{+0.04}_{-0.03}$ km and $D_{2}=0.45^{+0.04}_{-0.03}$ km, respectively. The surface thermal inertia $\Gamma$ is derived to be a low value of $80\pm40\rm Jm^{-2}s^{-0.5}K^{-1}$ with a roughness fraction $f_{\rm R}$ of $0.8^{+0.2}_{-0.4}$. This indicates that the primary possibly has a regolith layer on its surface, which is likely to be covered by a mixture of dust, fragmentary rocky debris and sand. The minimum regolith depth is estimated to be $5\sim20\rm mm$ from the simulations of subsurface temperature distribution, indicating that 1996 FG3 could be a very suitable target for a sample return mission.

The nucleus of Comet 67P/Churyumov-Gerasimenko

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 P sid = 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 H R (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 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.

The influence of rough surface thermal-infrared beaming on the Yarkovsky and YORP effects

It is now becoming widely accepted that photon recoil forces from the asymmetric reflection and thermal re-radiation of absorbed sunlight are, together with collisions and gravitational forces, primary mechanisms governing the dynamical and physical evolution of asteroids. The Yarkovsky effect causes orbital semi-major axis drift and the YORP effect causes changes in the rotation rate and pole orientation. We present an adaptation of the Advanced Thermophysical Model (ATPM) to simultaneously predict the Yarkovsky and YORP effects in the presence of thermal-infrared beaming caused by surface roughness, which has been neglected or dismissed in all previous models. Tests on Gaussian random sphere shaped asteroids, and on the real shapes of asteroids (1620) Geographos and (6489) Golevka, show that rough surface thermal-infrared beaming enhances the Yarkovsky orbital drift by typically tens of percent but it can be as much as a factor of two. The YORP rotational acceleration is on average dampened by up to a third typically but can be as much as one half. We find that the Yarkovsky orbital drift is only sensitive to the average degree, and not to the spatial distribution, of roughness across an asteroid surface. However, the YORP rotational acceleration is sensitive to the surface roughness spatial distribution, and can add significant uncertainties to the predictions for asteroids with relatively weak YORP effects. To accurately predict either effect the degree and spatial distribution of roughness across an asteroid surface must be known.

Physical characterization of low delta-V asteroid (175706) 1996 FG3★

Asteroid (175706) 1996 FG3 is a binary asteroid and the baseline target for the proposed MarcoPolo-R sample return mission. We present thermal-infrared photometry obtained with the European Southern Observatory (ESO) Very Large Telescope using the VISIR instrument, together with optical photometry obtained with the ESO New Technology Telescope using the EFOSC2 instrument. An absolute visual magnitude HV= 17.833 ± 0.024 and phase parameter G=-0.041 ± 0.005 are derived. The near-Earth asteroid thermal model has been fitted to the measured fluxes to derive a geometric visual albedo ?v= 0.046 p± 0.014, effective diameter at the observed aspect Deff= 1.68 ± 0.25 km and beaming parameter ? = 1.15 for phase angle ?±= 117. The advanced thermophysical model (ATPM) has been fitted to the measured fluxes to derive a more accurate effective diameter Deff= 1.71 ± 0.07km and albedo ?v= 0.044 ± 0.004. Based on the ATPM results, assuming the same albedo for primary and secondary, we derive a primary mean spherical diameter D?= 1.69+0.18 - 0.12km, secondary diameter Ds= 0.51 ± 0.03 km and a secondary orbital semimajor axis a= 2.8+1.7 -0.7km. A low surface thermal inertia ?= 120 ± 50Jm-2s-1/2K-1 was also derived, suggesting a dusty surface and raising questions as to the binary formation mechanism of this asteroid. These physical properties are used to predict a Yarkovsky drift in semimajor axis of -60+31 -45myr-1. © 2011 The Authors Monthly Notices of the Royal Astronomical Society © 2011 RAS.

Directional characteristics of thermal-infrared beaming from atmosphereless planetary surfaces - a new thermophysical model

We present a new rough-surface thermophysical model (Advanced Thermophysical Model or ATPM) that describes the observed directional thermal emission from any atmosphereless planetary surface. It explicitly incorporates partial shadowing, scattering of sunlight, self-heating and thermal–infrared beaming (re-radiation of absorbed sunlight back towards the Sun as a result of surface roughness). The model is verified by accurately reproducing ground-based directional thermal emission measurements of the lunar surface using surface properties that are consistent with the findings of the Apollo missions and roughness characterized by an rms slope of ∼32°. By considering the wide range of potential asteroid surface properties, the model implies a beaming effect that cannot be described by a simple parameter or function. It is highly dependent on the illumination and viewing angles as well as surface thermal properties and is predominantly caused by macroscopic rather than microscopic roughness. Roughness alters the effective Bond albedo and thermal inertia of the surface as well as moving the mean emission away from the surface normal. For accurate determination of surface properties from thermal–infrared observations of unresolved bodies or resolved surface elements, roughness must be explicitly modelled, preferably aided with thermal measurements at different emission angles and wavelengths.

10 Hygiea: ISO Infrared Observations

Observations of emissivity features of 10 Hygiea have been made for the first time in the relatively unexplored thermal-infrared wavelength region with the ISO (Infrared Space Observatory) satellite. Spectrophotomer (PHT-S) and short wavelength spectrometer (SWS) spectra of 10 Hygiea, obtained at 5.8–11.6 and 7–45 μm, respectively, are presented. In order to remove the thermal emission continuum, an advanced thermo–physical model has been applied to the observational data. To better interpret the spectral features above the thermal emission continuum, we compared the ISO observations with laboratory spectra available in the literature. Several laboratory experiments on minerals and meteorites have been performed to complete the analysis and to study the spectral behavior at various grain sizes. A possible spectral similarity with CO carbonaceous chondrites at small grain size is demonstrated.