Constraining Venus's interior with gravity and topography predictions from geodynamic models

One of the most informative ways of studying the interior structure and geodynamics of terrestrial planets is the joint investigation of gravity and topography data. In the case of Venus, this is in fact one of the only sources of information about the planet's interior, along with estimations of the tidal Love number [1], moment of inertia factor [2], and the absence of an internally generated magnetic field [3].

Early gravity and topography studies that made use of Pioneer Venus and later Magellan data revealed unique properties of Venus's interior. They showed that Venus has a notably higher gravity-topography correlation for long wavelengths compared to Earth and a globally large apparent depth of compensation [4]. Considering these characteristics and analyzing the wavelength-dependent ratio between gravity and topography, the so-called spectral admittance, [5] concluded that the long-wavelength topography of the planet was supported dynamically, i.e., through convection in the mantle.

Since then, several studies have investigated the gravity and topography signature of Venus with the goal of understanding the planet’s interior by constraining geophysical properties of the mantle, such as the mantle viscosity. Some estimated the dynamic geoid contribution from three-dimensional geodynamic simulations [6,7]. Others adopted the analytical viscous flow model by [11] to study the viscosity structure of Venus's mantle [8,9,10]. The main advantage of the latter method is that it is computationally inexpensive, allowing for the performance of inversions. However, it is a simplified model which neglects lateral viscosity variations.

In this study, we estimate the dynamic topography and geoid signatures from the geodynamical models by [12], which include a strongly temperature-dependent viscosity, hence lateral viscosity variations, to evaluate the influence of different parameters such as the increase of viscosity with depth, the presence of viscosity jumps, and the ratio between intrusive and extrusive magmatism. In a second step, we plan to systematically evaluate the influence of lateral viscosity variations on the geoid and topography and thus to quantify the importance of the simplifications adopted in the analytical model.

Our first results show that the increase of viscosity with depth should be no more than 2 orders of magnitude, since larger values strongly decrease the spectral correlation and admittance at long wavelengths which is inconsistent with the observations. In addition, scenarios where extrusive magmatism dominates tend to overestimate the admittance due to the generation of thick thermal lithospheres in excess of 300 km thickness. These results underscore the importance of gravity and topography analyses for deciphering the geodynamical evolution and tectonic style of Venus.

[1] Konopliv and Yoder (1996) GRL, 23; [2] Margot et al. (2021) Nat Astron; [3] Phillips and Russel (1987) JGR, 92; [4] Sjogren et al. (1980) JGR, 85;  [5] Kiefer et al. (1986) GRL, 13; [6] Huang et al. (2013) EPSL, 362; [7] Rolf et al. (2018) Icarus, 313; [8] Pauer et al. (2006) JGR, 111; [9] Steinberger et al. (2010) Icarus, 207; [10] Maia et al. (2023) GRL, 50; [11] Hager & Clayton (1989) Mantle Convection; [12] Plesa et al. (2023) EGU.