Hybrid Open-Access
Refine
Has Fulltext
- no (3)
Document Type
- Article (3)
Language
- English (3)
Is part of the Bibliography
- yes (3)
Keywords
- methods: data analysis (3) (remove)
Institute
Thermophysical modelling and parameter estimation of small solar system bodies via data assimilation
(2020)
Deriving thermophysical properties such as thermal inertia from thermal infrared observations provides useful insights into the structure of the surface material on planetary bodies. The estimation of these properties is usually done by fitting temperature variations calculated by thermophysical models to infrared observations. For multiple free model parameters, traditional methods such as least-squares fitting or Markov chain Monte Carlo methods become computationally too expensive. Consequently, the simultaneous estimation of several thermophysical parameters, together with their corresponding uncertainties and correlations, is often not computationally feasible and the analysis is usually reduced to fitting one or two parameters. Data assimilation (DA) methods have been shown to be robust while sufficiently accurate and computationally affordable even for a large number of parameters. This paper will introduce a standard sequential DA method, the ensemble square root filter, for thermophysical modelling of asteroid surfaces. This method is used to re-analyse infrared observations of the MARA instrument, which measured the diurnal temperature variation of a single boulder on the surface of near-Earth asteroid (162173) Ryugu. The thermal inertia is estimated to be 295 +/- 18 Jm(-2) K-1 s(-1/2), while all five free parameters of the initial analysis are varied and estimated simultaneously. Based on this thermal inertia estimate the thermal conductivity of the boulder is estimated to be between 0.07 and 0.12,Wm(-1) K-1 and the porosity to be between 0.30 and 0.52. For the first time in thermophysical parameter derivation, correlations and uncertainties of all free model parameters are incorporated in the estimation procedure that is more than 5000 times more efficient than a comparable parameter sweep.
The observation of the non-Keplerian behavior of propeller structures in Saturn's outer A ring raises the question: how does the propeller respond to the wandering of the central embedded moonlet? Here, we study numerically how the structural imprint of the propeller changes for a libration of the moonlet. It turns out that the libration induces an asymmetry in the propeller, which depends on the libration period and amplitude of the moonlet. Further, we study the dependence of the asymmetry on the libration period and amplitude for a moonlet with a 400 m Hill radius, which is located in the outer A ring. This allows us to apply our findings to the largest known propeller Blériot, which is expected to be of a similar size. For Blériot, we can conclude that, supposing the moonlet is librating with the largest observed period of 11.1 yr and an azimuthal amplitude of about 1845 km, a small asymmetry should be measurable but depends on the moonlet's libration phase at the observation time. The longitude residuals of other trans-Encke propellers (e.g., Earhart) show amplitudes similar to Blériot, which might allow us to observe larger asymmetries due to their smaller azimuthal extent, allowing us to scan the whole gap structure for asymmetries in one observation. Although the librational model of the moonlet is a simplification, our results are a first step toward the development of a consistent model for the description of the formation of asymmetric propellers caused by a freely moving moonlet.
Context. On 27 April 2015, when comet 67P/Churyumov-Gerasimenko was at 1.76 au from the Sun and moving toward perihelion, the OSIRIS and VIRTIS-M instruments on board the Rosetta spacecraft simultaneously observed the evolving dust and gas coma during a complete rotation of the comet. Aims. We aim to characterize the spatial distribution of dust, H2O, and CO2 gas in the inner coma. To do this, we performed a quantitative analysis of the release of dust and gas and compared the observed H2O production rate with the rate we calculated using a thermophysical model. Methods. For this study we selected OSIRIS WAC images at 612 nm (dust) and VIRTIS-M image cubes at 612 nm, 2700 nm (H2O emission band), and 4200 nm (CO2 emission band). We measured the average signal in a circular annulus to study the spatial variation around the comet, and in a sector of the annulus to study temporal variation in the sunward direction with comet rotation, both at a fixed distance of 3.1 km from the comet center. Results. The spatial correlation between dust and water, both coming from the sunlit side of the comet, shows that water is the main driver of dust activity in this time period. The spatial distribution of CO2 is not correlated with water and dust. There is no strong temporal correlation between the dust brightness and water production rate as the comet rotates. The dust brightness shows a peak at 0 degrees subsolar longitude, which is not pronounced in the water production. At the same epoch, there is also a maximum in CO2 production. An excess of measured water production with respect to the value calculated using a simple thermophysical model is observed when the head lobe and regions of the southern hemisphere with strong seasonal variations are illuminated (subsolar longitude 270 degrees-50 degrees). A drastic decrease in dust production when the water production (both measured and from the model) displays a maximum occurs when typical northern consolidated regions are illuminated and the southern hemisphere regions with strong seasonal variations are instead in shadow (subsolar longitude 50 degrees-90 degrees). Possible explanations of these observations are presented and discussed.