Current technology is not able to map the topography of rocky exoplanets, simply because the objects are too faint and far away to resolve them. Nevertheless, indirect effect of topography should be soon observable thanks to photometry techniques, and the possibility of detecting specular reflections. In addition, topography may have a strong effect on Earth-like exoplanet climates because oceans and mountains affect the distribution of clouds \citep{Houze2012}. Also topography is critical for evaluating surface habitability \citep{Dohm2015}. We propose here a general statistical theory to describe and generate realistic synthetic topographies of rocky exoplanetary bodies. In the solar system, we have examined the best-known bodies: the Earth, Moon, Mars and Mercury. It turns out that despite their differences, they all can be described by multifractral statistics, although with different parameters. Assuming that this property is universal, we propose here a model to simulate 2D spherical random field that mimics a rocky planetary body in a stellar system. We also propose to apply this model to estimate the statistics of oceans and continents to help to better assess the habitability of distant worlds.

Topography is the expression of both internal and external processes of a planetary body. Thus hypsometry (the study of topography) is a way to decipher the dynamic of a planet. For that purpose, the statistics of height and slopes may be described by different tools, at local and global scale. We propose here to use the multifractal approach to describe fields of topography. This theory both encompass height and slopes and other statistical moment of the field, tacking into account the scale invariance. Contrary to the widely used fractal formalism, multifractal is able to describe the intermittency of the topography field. As we commonly observe the juxtapostion of rough and smooth at given scale, the multifractal framework seems to be appropriate for hypsometric studies. Here we analyze the data at global scale of the Earth, Mars, Mercury and the Moon and find that the statistics are in good agreement with the multifractal theory for scale larger than 10km. Surprisingly, the analysis shows that all bodies have the same fractal behavior for scale smaller than 10km. We hypothesized that dynamic topography of the mantle may be the explanation at large scale, whereas the smaller scales behavior may be related to elastic thickness.