Skip to content

Extract and plot a global map of effective elastic thickness values

License

Notifications You must be signed in to change notification settings

paudetseis/GlobalTe

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

12 Commits
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

GlobalTe

Simple tools to extract and plot a global model of effective elastic thickness (Te) of the lithosphere. This Te 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 Te using the wavelet transform.

Installation

Dependencies

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

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 Te 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 Te 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 Te is only calculated over dry areas (including shallow shelf areas).


It is also possible to extract a Te value at any arbitrary point on the globe:

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

Reference

About

Extract and plot a global map of effective elastic thickness values

Topics

Resources

License

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages