Skip to main content

A Python package for paleoclimate data analysis

Project description

PyPI PyPI license DOI NSF-1541029


Python Package for the Analysis of Paleoclimate Data

Table of contents

What is it?

Pyleoclim is a Python package primarily geared towards the analysis and visualization of paleoclimate data. Such data often come in the form of timeseries with missing values and age uncertainties, so the package includes several low-level methods to deal with these issues, as well as high-level methods that re-use those within scientific workflows.

High-level modules assume that data are stored in the Linked Paleo Data (LiPD) format and makes extensive use of the LiPD utilities. Low-level modules are primarily based on NumPy arrays or Pandas dataframes, so Pyleoclim contains a lot of timeseries analysis code (e.g. spectral analysis, singular spectrum analysis, wavelet analysis, correlation analysis) that can apply to these more common types as well. See the example folder for details.

The package is aware of age ensembles stored via LiPD and uses them for time-uncertain analyses very much like GeoChronR.

Current capabilities:

  • binning
  • interpolation
  • standardization
  • plotting maps, timeseries, and basic age model information
  • paleo-aware correlation analysis (isopersistent, isospectral and classical t-test)
  • weighted wavelet Z transform (WWZ)
  • age modeling through Bchron

Future capabilities:

  • paleo-aware singular spectrum analysis (AR(1) null eigenvalue identification, missing data)
  • spectral analysis (Multi-Taper Method, Lomb-Scargle)
  • cross-wavelet analysis
  • index reconstruction
  • climate reconstruction
  • causality
  • ensemble methods for most of the above

If you have specific requests, please contact


Python v3.6+ is required.

We recommend using Anaconda, with an environment dedicated to Pyleoclim. See the documentation for details.

To install Pyleoclim, first install numpy and Cartopy through Anaconda (conda)

conda install numpy
conda install -c conda-forge cartopy

Then install pyleoclim via pip

pip install pyleoclim

Note that the pip command line above will trigger the installation of (most of) the dependencies, as well as the local compilation of the Fortran code for WWZ with the GNU Fortran compiler gfortran. If you have the Intel's Fortran compiler ifort installed, then further accerlation for WWZ could be achieved by compiling the Fortran code with ifort, and below are the steps:

  • download the source code, either via git clone or just download the .zip file
  • modify by commenting out the line of extra_f90_compile_args for gfortran, and use the line below for ifort
  • run python build_ext --fcompiler=intelem && python install

Some functionalities require R.

Version Information

Current Version

0.4.10: Support local compilation of the Fortran code for WWZ; precompiled .so files have been removed.

Past Versions

0.4.9: Major bug fixes; mapping module based on cartopy; compatibility with latest numpy package
0.4.8: Add support of f2py WWZ for Linux
0.4.7: Update to coherence function
0.4.6: Fix an issue when copying the .so files
0.4.5: Update to to include proper .so file according to version
0.4.4: New fix for .so issue
0.4.3: New fix for .so issue
0.4.2: Fix issue concerning download of .so files
0.4.1: Fix issues with tarball
0.4.0: New functionalities: map nearest records by archive type, plot ensemble time series, age modelling through Bchron
0.3.1: New functionalities: segment a timeseries using a gap detection criteria, update to summary plot to perform spectral analysis
0.3.0: Compatibility with LiPD 1.3 and Spectral module added 0.2.5: Fix error on loading (Looking for Spectral Module)
0.2.4: Fix load error from init
0.2.3: Freeze LiPD version to 1.2 to avoid conflicts with 1.3
0.2.2: Change progressbar to tqdm and add standardization function
0.2.1: Update package requirements
0.2.0: Restructure the package so that the main functions can be called without the use of a LiPD files and associated timeseries objects.
0.1.4: Rename function using camel case and consistency with LiPD utilities version
0.1.3: Compatible with LiPD utilities version Function openLiPD() renamed openLiPDs()
0.1.2: Compatible with LiPD utilities version Uses basemap instead of cartopy
0.1.1: Freezes the package prior to version of LiPD utilities

Quickstart guide

  1. Install Pyleoclim.

  2. Wait for installation to complete, then:

    3a. Import the package into your favorite Python environment (we recommend the use of Spyder, which comes standard with the Anaconda package)

    3b. Use Jupyter Notebook to go through the tutorial contained in the PyleoclimQuickstart.ipynb Notebook, which can be downloaded here. The folder also contains a collection of LiPD files. More LiPD files available here.

  3. Help with functionalities can be found in the Documentation.


Tested with:

  • LiPD 0.2.7
  • pandas v0.25.0
  • numpy v1.16.4
  • matplotlib v3.1.0
  • Cartopy v1.17.0
  • scipy v1.3.1
  • statsmodel v0.8.0
  • seaborn 0.9.0
  • scikit-learn 0.21.3
  • tqdm 4.33.0
  • pathos 0.2.4
  • rpy2 3.0.5

The installer will automatically check for the needed updates.

Known Issues

  • Some of the packages supporting Pyleoclim do not have a build for Windows
  • Known issues with proj4 v5.0-5.1, make sure your environment is set up with v5.2

Further information

Python and Anaconda:
Jupyter Notebook:


Please report issues to


The project is licensed under the GNU Public License. Please refer to the file call license. If you use the code in publications, please credit the work using this citation.


This material is based upon work supported by the National Science Foundation under Grant Number ICER-1541029. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the investigators and do not necessarily reflect the views of the National Science Foundation.

This research is funded in part by JP Morgan Chase & Co. Any views or opinions expressed herein are solely those of the authors listed, and may differ from the views and opinions expressed by JP Morgan Chase & Co. or its affilitates. This material is not a product of the Research Department of J.P. Morgan Securities LLC. This material should not be construed as an individual recommendation of for any particular client and is not intended as a recommendation of particular securities, financial instruments or strategies for a particular client. This material does not constitute a solicitation or offer in any jurisdiction.

Project details

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Files for pyleoclim, version 0.5.0
Filename, size File type Python version Upload date Hashes
Filename, size pyleoclim-0.5.0.tar.gz (98.6 kB) File type Source Python version None Upload date Hashes View

Supported by

Pingdom Pingdom Monitoring Google Google Object Storage and Download Analytics Sentry Sentry Error logging AWS AWS Cloud computing DataDog DataDog Monitoring Fastly Fastly CDN DigiCert DigiCert EV certificate StatusPage StatusPage Status page