.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "api/gallery/glaciology/plot_03_jamieson_glacial_isostasy.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_api_gallery_glaciology_plot_03_jamieson_glacial_isostasy.py>`
        to download the full example code or to run this example in your browser via JupyterLite.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_api_gallery_glaciology_plot_03_jamieson_glacial_isostasy.py:


Jamieson's glacial isostasy (1865)
==================================

Thomas Jamieson explained the raised beaches of Scotland, now far above the
sea, by the weight of the ice sheet that once covered the country: the land
sank under the ice and has risen since it melted. In equilibrium the mantle
gives way until the displaced rock weighs as much as the ice, so the bed sinks
by the ice thickness times the density ratio, about 0.28.

.. GENERATED FROM PYTHON SOURCE LINES 11-32



.. image-sg:: /api/gallery/glaciology/images/sphx_glr_plot_03_jamieson_glacial_isostasy_001.png
   :alt: Jamieson: the land sinks under the ice
   :srcset: /api/gallery/glaciology/images/sphx_glr_plot_03_jamieson_glacial_isostasy_001.png
   :class: sphx-glr-single-img





.. code-block:: Python


    import matplotlib.pyplot as plt
    import numpy as np

    from earthsciencekit.constants import ICE_DENSITY
    from earthsciencekit.glaciology import isostatic_depression

    mantle_density = 3300.0
    distance = np.linspace(-600e3, 600e3, 601)
    ice = 2000.0 * np.sqrt(np.clip(1 - (distance / 500e3) ** 2, 0, None))
    depression = isostatic_depression(ice, ICE_DENSITY, mantle_density)
    assert np.isclose(depression.max(), 2000 * ICE_DENSITY / mantle_density)

    fig, ax = plt.subplots()
    ax.fill_between(distance / 1e3, -depression, ice - depression, color="lightblue", label="ice")
    ax.plot(distance / 1e3, -depression, color="saddlebrown", label="depressed bed")
    ax.axhline(0, color="0.5", ls=":", label="bed before glaciation")
    ax.set(xlabel="distance [km]", ylabel="elevation [m]", title="Jamieson: the land sinks under the ice")
    ax.legend()

    plt.show()


.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (0 minutes 0.095 seconds)


.. _sphx_glr_download_api_gallery_glaciology_plot_03_jamieson_glacial_isostasy.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: lite-badge

      .. image:: images/jupyterlite_badge_logo.svg
        :target: ../../../lite/lab/index.html?path=api/gallery/glaciology/plot_03_jamieson_glacial_isostasy.ipynb
        :alt: Launch JupyterLite
        :width: 150 px

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: plot_03_jamieson_glacial_isostasy.ipynb <plot_03_jamieson_glacial_isostasy.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: plot_03_jamieson_glacial_isostasy.py <plot_03_jamieson_glacial_isostasy.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: plot_03_jamieson_glacial_isostasy.zip <plot_03_jamieson_glacial_isostasy.zip>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_