Thermophysical modelling for high-resolution digital terrain models

A method is presented for efficiently calculating surface temperatures for highly resolved celestial body shapes. A thorough investigation of the necessary conditions leading to reach model convergence shows that the speed of surface temperature convergence depends on factors such as the quality of initial boundary conditions, thermal inertia, illumination conditions, and resolution of the numerical depth grid. The optimization process to shorten the simulation time while increasing or maintaining the accuracy of model results includes the introduction of facet-specific boundary conditions such as pre-computed temperature estimates and pre-evaluated simulation times. The individual facet treatment also allows for assigning other facet-specific properties such as local thermal inertia. The approach outlined in this paper is particularly useful for very detailed digital terrain models in combination with unfavourable illumination conditions such as little-to-no sunlight at all for a period of time as experienced locally on comet 67P/Churyumov–Gerasimenko. Possible science applications include thermal analysis of highly resolved local (landing) sites experiencing seasonal, environment, and lander shadowing. In combination with an appropriate roughness model, the method is very suitable for application to disc-integrated and disc-resolved data. Further applications are seen where the complexity of the task has led to severe shape or thermophysical model simplifications such as in studying surface activity or thermal cracking.

[1]  B. Altieri,et al.  Thermal and shape properties of asteroid (21) Lutetia from Herschel observations around the Rosetta flyby , 2012 .

[2]  J. Blum,et al.  Outgassing of icy bodies in the Solar System – II: Heat transport in dry, porous surface dust layers , 2011, 1111.0535.

[3]  T. Encrenaz,et al.  Subsurface properties and early activity of comet 67P/Churyumov-Gerasimenko , 2015, Science.

[4]  S. Debei,et al.  Insolation, erosion, and morphology of comet 67P/Churyumov-Gerasimenko , 2015 .

[5]  J. Blum,et al.  Outgassing of icy bodies in the Solar System – I. The sublimation of hexagonal water ice through dust layers , 2011, 1101.2518.

[6]  Line Drube,et al.  THERMAL TOMOGRAPHY OF ASTEROID SURFACE STRUCTURE , 2016, 1608.06839.

[7]  Giampiero Naletto,et al.  The rotation state of 67P/Churyumov-Gerasimenko from approach observations with the OSIRIS cameras on Rosetta , 2014 .

[8]  S. Squyres,et al.  Methods for computing comet core temperatures. , 1986, Icarus.

[9]  H. Keller,et al.  MIRO observations of subsurface temperatures of the nucleus of 67P/Churyumov-Gerasimenko , 2015 .

[10]  E. Kührt.  Temperature profiles and thermal stresses in cometary nuclei , 1984 .

[11]  Giampiero Naletto,et al.  Shape model, reference system definition, and cartographic mapping standards for comet 67P/Churyumov-Gerasimenko Stereo-photogrammetric analysis of Rosetta/OSIRIS image data , 2015 .

[12]  Jens Biele,et al.  Thermophysical modeling of Didymos’ moon for the Asteroid Impact Mission , 2017 .

[13]  H. Senshu,et al.  Feasibility and Accuracy of Thermophysical Estimation of Asteroid 162173 Ryugu (1999 JU3) from the Hayabusa2 Thermal Infrared Imager , 2017 .

[14]  H. Rickman,et al.  Surface roughness and three-dimensional heat conduction in thermophysical models , 2014 .

[15]  Paul Hartogh,et al.  Continuum and spectroscopic observations of asteroid (21) Lutetia at millimeter and submillimeter wavelengths with the MIRO instrument on the Rosetta spacecraft , 2012 .

[16]  Akira Fujiwara,et al.  Pole and Global Shape of 25143 Itokawa , 2006, Science.

[17]  Tilman Spohn,et al.  Line heat-source measurements of the thermal conductivity of porous H2O ice, CO2 ice and mineral powders under space conditions , 1996 .

[18]  J. Bandfield,et al.  Interpretation of thermal emission. I. The effect of roughness for spatially resolved atmosphereless bodies , 2015 .

[19]  T. Spohn,et al.  Three-dimensional illumination and thermal model of the Abydos region on comet 67P/Churyumov-Gerasimenko , 2017 .

[20]  Jens Biele,et al.  Rosetta Lander - Landing and operations on comet 67P/Churyumov-Gerasimenko , 2016 .

[21]  S. Erard,et al.  How pristine is the interior of the comet 67P/Churyumov-Gerasimenko? , 2017 .

[22]  S. Green,et al.  Directional characteristics of thermal–infrared beaming from atmosphereless planetary surfaces – a new thermophysical model , 2011, 1211.1844.

[23]  Sukhan Lee,et al.  Interpretation of combined infrared, submillimeter, and millimeter thermal flux data obtained during the Rosetta fly-by of Asteroid (21) Lutetia , 2012 .

[24]  M. Banaszkiewicz,et al.  Thermal and mechanical properties of the near-surface layers of comet 67P/Churyumov-Gerasimenko , 2015, Science.

[25]  S. Hasegawa,et al.  (25143) Itokawa : The power of radiometric techniques for the interpretation of remote thermal observations in the light of the Hayabusa rendezvous results (New Insights into Near-Earth and Main-Belt Asteroids) , 2014, 1404.5842.

[26]  S. Hasegawa,et al.  Hayabusa-2 mission target asteroid 162173 Ryugu (1999 JU3): Searching for the object's spin-axis orientation , 2016, 1611.05625.

[27]  T. Fusco,et al.  3D shape of asteroid (6) Hebe from VLT/SPHERE imaging: Implications for the origin of ordinary H chondrites , 2017, 1705.10515.

[28]  Hajime Yano,et al.  Regolith Migration and Sorting on Asteroid Itokawa , 2007, Science.

[29]  S. Debei,et al.  The global meter-level shape model of comet 67P/Churyumov-Gerasimenko , 2017 .