Apply a simple ice-volume correction to marine-isotope proxy records in Python
The correction uses the LR04 benthic stack from Lisiecki and Raymo, 2005. The stack is scaled so that the LGM-to-present change is assumed to be 1.0 ‰ in accordance with the pore-water estimate of Schrag et al. 1996. The package is adapted from the ice-volume correction in Tierney et al. 2017.
Start by importing
import erebusfall as ef import numpy as np
...and creating a hypothetical proxy time series...
age_ka = np.arange(0, 20, 1) proxy = np.random.normal(loc=-2.0, size=len(age_ka))
Now we can plug this into
ef.icevol_correction() with a few key options:
proxy_adjusted = ef.icevol_correction(age_ka, proxy, proxytype='d18o', timeunit='ka')
We first plug in the proxy age and proxy values. The
indicates that we're dealing with δ18O. We can set
proxytype='dd' for a δD
timeunit argument indicates that
age_ka is in thousands of
years before present. There are also options for for "years BP" and
"million years BP". The output from the function,
proxy_adjusted, is the
corrected isotope proxy as a numpy array.
help(ef.icevol_correction) for more documentation.
You can install the package from PyPI with
pip install erebusfall
If you want to use
conda install erebusfall -c sbmalev
Development and Support
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
|Filename, size & hash SHA256 hash help||File type||Python version||Upload date|
|erebusfall-0.0.3-py3-none-any.whl (31.8 kB) Copy SHA256 hash SHA256||Wheel||py3|
|erebusfall-0.0.3.tar.gz (20.2 kB) Copy SHA256 hash SHA256||Source||None|