Geoid and topography of Venus in various thermal convection models

Abstract

Important though indirect information about the internal structure of Venus is provided by its topography and geoid. In the last decades this information has been used to constrain the Venus mantle viscosity structure and its dynamic regime. Recently, the geodynamic inversion of the Venus’ geoid and topography resulted in a group of best fitting viscosity profiles. We use these viscosity models here as an input to our mantle convection code. We carry out simulations of the Venus’ mantle evolution in a 3D spherical shell with depth dependent viscosity and check whether the character of the dynamic topography and the geoid represented by their power spectra fits the observed quantities. We compare the results with several other models obtained for different viscosity stratifications (constant, constant with highly viscous lithosphere, linear increase of viscosity). Further, we estimate the effect of other factors such as internal heating and varying Rayleigh number. We use a 2D spherical axisymmetric convection code to study the effect of lateral viscosity variations. In these 2D models we monitor the topography and the geoid developing above the axisymmetric plume and compare them with the observed elevations of Venus’ geoid and topography in several Regia. Though none of the models fits observed data perfectly, we can generally conclude, that the best fit between the observed and predicted quantities is reached for viscosity profiles with 200 km thick lithosphere followed by a gradual increase of viscosity with depth and with the upper mantle viscosity of 2 × 10 21 Pa s. For all viscosity profiles the predicted geoid and topography spectra match the observed ones only up to the degree 40, thus indicating other than dynamic origin of these quantities for higher degrees.