js1019 / PlanetaryModels

Build 3D planetary models on deformable meshes
GNU Lesser General Public License v3.0
19 stars 9 forks source link
earth-science planetary-science unstructured-meshes

Planetary Model Builder

This repository provides scripts to build a planetary model on a deformable tetrahedal mesh as well as its reference gravity. It supports and provides input files for the repository "normal mode computation at planetary scales".

The above are two animations for compressional wave speed of PREM and shear wave speed of PREM.

Contents

How to make it work for you?

This repository uses several packages, including

We've put most of the original codes under different folders with some minor changes for this application.

You may follow the README.md under packages to compile what you need. If you only need to build a tetrahedral mesh, please compile TetGen; if you are a Linux user, just go to packages/tetgen1.5.0 and type make. If you need to compute the reference gravity, please check the readme under packages/ and compile packages/fmmlib3d-1.2.

Build your planetary models

Please check the scripts under Demos/CONST for a constant ball model, Demos/PREM for a standard Earth model, Demos/Moon for a Moon model and Demos/RT_MarsDWAK for a Mars model. Under these folders, you may run

run {CONST/PREM/M6Ltopo/MarsDWAK}_mesh; run Gravity;

to obtain the model and its reference gravity. Note that the 3D crusts used in this work are based on the results of several NASA missions, including the Lunar Orbiter Laser Altimeter (LOLA), Gravity Recovery and Interior Laboratory (GRAIL) and the Mars Orbiter Laser Altimeter (MOLA). Please see discontinuities/ for more details.

Illustration of other terrestrial planets

Here, we show a Moon shear wave speed model in the top right figure, the Mars topography in the top middle figure and a Mars shear wave speed model in the top right figure. These models can be reproduced by the scripts provided in the Demos.

The corresponding Moon and Mars models can be found via this IEEE data port.

Include more heterogeneity

You can build more realistic models using similiar ideas and following the above tips! You may need to understand what distmesh can do. It is quite simple, please see its Demos.

The top right figure illustrates an Earth compressional wave speed model based on MIT tomographic results (Burdick et al. 2017) and crust 1.0 (Laske et al. 2013). The top left one shows the topography of the Moho discontinuity under Tibetan Plateau. If you would like to reproduce these models, you may need to download the model data and make some minor changes in the scripts to obtain and visualize your results.

A few remarks

Postprocess for visualizing the modes

You can use scripts in visual/ to visualize your computed normal modes from the repository "normal mode computation at planetary scales". Please check visual/README.md for more details.

Reference

Other related references can be found in modelbuilder/, packages/, radialmodels/ and discontinuities/.

Report

Please let us know any issues of this repository. Contributions are welcome.