GlobalTe

Simple tools to extract and plot a global model of effective elastic thickness (<i>T_e</i>) of the lithosphere. This <i>T_e</i> model was obtained from the inversion of the wavelet coherence between the topography and Bouguer gravity anomalies. See Audet and Burgmann (2011) for details.

See also the software PlateFlex, which describes how to produce regional grids of <i>T_e</i> using the wavelet transform.

Installation

Dependencies

The current version was developed using Python3.7 Also, the following packages are required:

• numpy
• matplotlib
• cartopy

Conda environment

We recommend creating a custom conda environment where GlobalTe can be installed along with its dependencies.

• Create a environment called te and install all dependencies:

conda create -n te python=3.7 numpy matplotlib cartopy -c conda-forge

• Activate the environment:

conda activate te

• Clone the GlobalTe repository and install using pip:

git clone https://github.com/paudetseis/GlobalTe.git
cd GlobalTe
pip install .

Usage

Once installed, you can produce a <i>T_e</i> map with three lines in Python:

>>> from globalte import TeModel as TM
>>> model = TM()
>>> model.plot_global()

You can also specify to load a data set that does not include data points biased by gravitational noise. You can also check out a different global projection and a different colormap:

>>> model_nobias = TM(nobias=True)
>>> model_nobias.plot_global(proj='IGH', cmap='magma')

The corresponding <i>T_e</i> grids are available as attributes of the model object:

>>> model.te_global
array([[nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan],
       ...,
       [nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan],
       [nan, nan, nan, ..., nan, nan, nan]])

---

NOTE

Most data points are NaN since <i>T_e</i> is only calculated over dry areas (including shallow shelf areas).

---

It is also possible to extract a <i>T_e</i> value at any arbitrary point on the globe:

>>> model.get_te_point(62.3, -125.7)
44.9015197754

Reference

• Audet, P. and Burgmann, R. Dominant role of tectonic inheritance in supercontinent cycles, Nat. Geosci., 4, 184-187, 2011