.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "api/gallery/glaciology/plot_08_nye_crevasse_depth.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_08_nye_crevasse_depth.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_08_nye_crevasse_depth.py:


Nye's crevasse depth (1955)
===========================

John Nye argued that a crevasse stays open only down to the depth where the
pressure of the overlying ice equals the tensile stress from stretching,
which Glen's law gives from the strain rate. Dry crevasses therefore reach
only a few tens of metres. Water in a crevasse adds its own pressure and
drives it deeper, which is how meltwater can fracture an entire ice shelf.

.. GENERATED FROM PYTHON SOURCE LINES 11-31



.. image-sg:: /api/gallery/glaciology/images/sphx_glr_plot_08_nye_crevasse_depth_001.png
   :alt: Nye: crevasses stop where the ice closes them
   :srcset: /api/gallery/glaciology/images/sphx_glr_plot_08_nye_crevasse_depth_001.png
   :class: sphx-glr-single-img





.. code-block:: Python


    import matplotlib.pyplot as plt
    import numpy as np

    from earthsciencekit.constants import GRAVITY, ICE_DENSITY, SECONDS_PER_YEAR
    from earthsciencekit.glaciology import nye_crevasse_depth, water_filled_crevasse_depth

    rate = np.logspace(-3, -1, 100) / SECONDS_PER_YEAR
    dry = nye_crevasse_depth(rate, 2.4e-24, 3, ICE_DENSITY, GRAVITY)
    assert np.all((dry > 5) & (dry < 60)) and np.isclose(dry[-1] / dry[0], 100 ** (1 / 3))
    fig, ax = plt.subplots()
    ax.semilogx(rate * SECONDS_PER_YEAR, dry, label="dry")
    for water in (10, 30):
        wet = water_filled_crevasse_depth(rate, 2.4e-24, 3, water, ICE_DENSITY, 1000, GRAVITY)
        ax.semilogx(rate * SECONDS_PER_YEAR, wet, label=f"{water} m of water")
    ax.invert_yaxis()
    ax.set(xlabel="stretching rate [1/year]", ylabel="crevasse depth [m]", title="Nye: crevasses stop where the ice closes them")
    ax.legend()

    plt.show()


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

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


.. _sphx_glr_download_api_gallery_glaciology_plot_08_nye_crevasse_depth.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_08_nye_crevasse_depth.ipynb
        :alt: Launch JupyterLite
        :width: 150 px

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

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

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

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

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

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


.. only:: html

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

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