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 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 data assimilation method, the Ensemble Square Root Filter, to 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 \pm 18$ $\mathrm{J\,m^{-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 $\mathrm{W\,m^{-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 which is more than 5000 times more efficient than a comparable parameter sweep.

Heat flow is an important constraint on planetary formation and evolution. It has been suggested that Martian obliquity cycles might cause periodic collapses in atmospheric pressure, leading to corresponding decreases in regolith thermal conductivity (which is controlled by gas in the pore spaces). Geothermal heat would then build up in the subsurface, potentially affecting present-day heat flow - and thus the measurements made by a heat-flow probe such as the InSight HP$^{3}$ instrument. To gauge the order of magnitude of this effect, we model the diffusion of a putative heat pulse caused by thermal conductivity changes with a simple numerical scheme and compare it to the heat-flow perturbations caused by other effects. We find that an atmospheric collapse to 300 Pa in the last 40 kyr would lead to a present-day heat flow that is up to $2-8\%$ larger than the average geothermal background. Considering the InSight mission with expected $5-15\%$ error bars on the HP$^{3}$ measurement, this perturbation would only be significant in the best-case scenario of full instrument deployment, completed measurement campaign, and a well-modelled surface configuration. The prospects for detecting long-term climate perturbations via spacecraft heat-flow experiments remain challenging.

Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on.

The evolution and escape of the martian atmosphere and the planet's water inventory can be separated into an early and late evolutionary epoch. The first epoch started from the planet's origin and lasted $\sim$500 Myr. Because of the high EUV flux of the young Sun and Mars' low gravity it was accompanied by hydrodynamic blow-off of hydrogen and strong thermal escape rates of dragged heavier species such as O and C atoms. After the main part of the protoatmosphere was lost, impact-related volatiles and mantle outgassing may have resulted in accumulation of a secondary CO$_2$ atmosphere of a few tens to a few hundred mbar around $\sim$4--4.3 Gyr ago. The evolution of the atmospheric surface pressure and water inventory of such a secondary atmosphere during the second epoch which lasted from the end of the Noachian until today was most likely determined by a complex interplay of various nonthermal atmospheric escape processes, impacts, carbonate precipitation, and serpentinization during the Hesperian and Amazonian epochs which led to the present day surface pressure.

Strange mode instabilities in the envelopes of massive stars lead to shock waves, which can oscillate on a much shorter timescale than that associated with the primary instability. The phenomenon is studied by direct numerical simulation using a, with respect to time, implicit Lagrangian scheme, which allows for the variation by several orders of magnitude of the dependent variables. The timestep for the simulation of the system is reduced appreciably by the shock oscillations and prevents its long term study. A procedure based on domain decomposition is proposed to surmount the difficulty of vastly different timescales in various regions of the stellar envelope and thus to enable the desired long term simulations. Criteria for domain decomposition are derived and the proper treatment of the resulting inner boundaries is discussed. Tests of the approach are presented and its viability is demonstrated by application to a model for the star P Cygni. In this investigation primarily the feasibility of domain decomposition for the problem considered is studied. We intend to use the results as the basis of an extension to two dimensional simulations.

The evolution of strange mode instabilities into the non linear regime has been followed by numerical simulation for an envelope model of a massive star having solar chemical composition, M=50M_sun, T_eff=10^4K and L=1.17*10^6 L_sun. Contrary to previously studied models, for these parameters shocks are captured in the H-ionisation zone and perform rapid oscillations within the latter. A linear stability analysis is performed to verify that this behaviour is physical. The origin of an instability discovered in this way is identified by construction of an analytical model. As a result, the stratification turns out to be essential for instability. The difference to common stratification instabilities, e.g., convective instabilities, is discussed.