Skip to content

Elastic thickness of the lithosphere, yield strength envelope, and heat flow calculations.

License

Notifications You must be signed in to change notification settings

AB-Ares/Te_HF_Conversion

Repository files navigation

Code style: black DOI

Te_HF_Conversion

Elastic thickness of the lithosphere, yield strength envelope, and heat flow calculations.

Description

Te_HF_Conversion (TeHF) is a simple code that allows to:

(1) Convert the elastic thickness of the lithosphere to planetary heat flow (and a yield strength envelope) given several input parameters including crustal thickness, strain rate, or radiogenic heating.

(2) Retrieve the elastic thickness (or mechanical thickness, and a yield strength envelope) of the lithosphere based on a temperature profile.

The model makes use of the equating bending moment approach of McNutt (1984) and has been used in Broquet et al (2020,2021), or Maia & Wieczorek (2022).

Benchmarks

Heat flow calculations have been benchmarked to various studies (e.g., McNutt, 1984 or Solomon & Head, 1990). I thank Julia Maia for having performed some further benchmarks to the literature.

Methods

Curv_Moment Determine the bending moment given the input yield strength envelope and curvature.

Conversion_Te_HF Determine the surface, crustal, and mantle heat flows, mechanical thickness, and thermal gradients from input rheology and elastic parameters.

Conversion_Tprofile_Te Determine yield strength envelope, mechanical thickness given the input temperature profile. Elastic thickness will also be output from the assumed plate curvature.

Brittle_Strength Compute the brittle part of the yield strength envelope based on planetary constants.

Example scripts

Mars_YSE Determine the surface, crust, and mantle heat flows for a given elastic thickness on Mars assuming a wet rheology for the diabase crust and olivine mantle. Plot the yield strength envelope.

Mars_T_profile_Te Determine the yield strength envelope, mechanical and elastic thickness for an assumed temperature profile that is linear but has a cosine tapered temperature anomaly between 20 and 50 km depth.

Venus_YSE Determine the surface, crust, and mantle heat flows for a given elastic thickness on Venus assuming a dry rheology for the diabase crust and olivine mantle. Plot the yield strength envelope.

How to install and run Te_HF_Conversion

If you would like to modify the source code, download the Displacement_strain_planet repository and install using pip (or pip3 depending on your installation).

    git clone https://github.com/AB-Ares/Te_HF_Conversion.git
    cd Te_HF_Conversion/
    pip install .

To run the example scripts

    cd examples
    python Mars_YSE.py 
    python Venus_YSE.py 

Author

Adrien Broquet ([email protected])

Cite

If you use the package, please cite the latest release as: Adrien Broquet. AB-Ares/Te_HF_Conversion: 0.2.3 (Version 0.2.3). Zenodo. http://doi.org/10.5281/zenodo.4973893

This package was created for Broquet et al., 2020, which you can also cite as

    @article{Broquet2020,
    author = {Broquet, A. and Wieczorek, M. A. and Fa, W.},
    title = {Flexure of the Lithosphere Beneath the North Polar Cap of Mars: Implications for Ice Composition and Heat Flow},
    journal = {Geophysical Research Letters},
    volume = {47},
    number = {5},
    doi = {https://doi.org/10.1029/2019GL086746},
    year = {2020}}