Create a crustal thickness map of a planet
Project description
Description
ctplanet provides several functions and example scripts for generating crustal thickness maps of a planet from gravity and topography data, and the calculation of hydrostatic relief along density interfaces beneath the lithosphere.
Methods
pyMoho Calculate relief using a constant density crust and mantle.
pyMohoRho Calculate relief using a constant density mantle and a variable density crust.
HydrostaticShapeLith Calculate the relief of hydrostatic interfaces beneath the lithosphere along with the predicted gravity, taking into account rotation and/or tides using the approach of Wieczorek et al. (2019).
HydrostaticShape Calculate the relief of hydrostatic interfaces and predicted gravity of a rotating hydrostatic planet using the approach of Wieczorek et al. (2019).
InertiaTensor_from_shape Calculate the inertia tensor given a radial density profile and shape of each interface.
InertiaTensor_from_C Calculate the inertia tensor given the polar moment of inertia and the gravitational potential coefficients.
moi Calculate the mean, normalized, moment of inertia up to index n.
ReadRefModel Read the reference interior model file.
Example scripts
Crust_Moon A script that demonstrates how to calculate the thickenss of the lunar crust using either a constant or variable density crust. The latter can be used to reproduce the results presented in Wieczorek et al. (2013).
Crust_Mars A script that demonstrates how to calculate the thickenss of the Martian crust using either a constant or variable density crust. For the variable density crust, the density is assumed to change discontinuously across the dichotomy boundary.
Mars_crust_thick_test Create a crustal thickness map of Mars from gravity and topography and compare how results change if hydrostatic interfaces are not taken into account.
Crust_Mars_InSight Create a crustal thickness map of Mars from gravity and topography, using the InSight crustal thickness constraint.
Crust_Mars_InSight_dichotomy Create a crustal thickness map of Mars from gravity and topography, using the InSight crustal thickness constraint and different densities across the dichotomy boundary.
mars_fcn Compute the free core nutation period of Mars.
mars_figs Create images related to Mars in Wieczorek et al. (2019).
mars_j2 Compute the contribution to the gravitational J2 of Mars from hydrostatic interfaces beneath the lithosphere.
Core-Moon Calculate the hydrostatic relief of the lunar core accounting for the non-hydrostatic potential that comes from the lithosphere.
Earth_test Compute hydrostatic relief of Earth using PREM.
ceres Calculate the hydrostatic shape of Ceres.
How to install and run ctplanet
Download the ctplanet repository and install using pip
git clone https://github.com/MarkWieczorek/ctplanet.git
pip install .
To execute a script
cd examples
python Crust_Moon.py
Depending on how your system is set up, it might be necessary to use explicitly python3 and pip3 instead of python and pip in the above commands.
Reference
Wieczorek, M. A., G. A. Neumann, F. Nimmo, W. S. Kiefer, G. J. Taylor, H. J. Melosh, R. J. Phillips, S. C. Solomon, J. C. Andrews-Hanna, S. W. Asmar, A. S. Konopliv, F. G. Lemoine, D. E. Smith, M. M. Watkins, J. G. Williams, M. T. Zuber (2013), The crust of the Moon as seen by GRAIL, Science, 339, 671-675, doi:10.1126/science.1231530.
Wieczorek, M. A., M. Beuthe, A. Rivoldini, and T. Van Hoolst (2019), Hydrostatic interfaces in bodies with nonhydrostatic lithospheres, Journal of Geophysical Research: Planets, 124, doi:10.1029/2018JE005909.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.