Reproducing LMRv2.1 with PAGES2k and LinkedEarth Tools

Reproducing LMRv2.1 with PAGES2k and LinkedEarth Tools#

By Tanaya Gondhalekar, Deborah Khider & Julien Emile-Geay.

Climate field reconstruction is the task of estimating variations in one or more climate fields (e.g. surface temperature or precipitation) from a collection of paleoclimate observations (aka “proxies”). Many statistical methods are available for doing so ; a relatively new and impactful one has been offline data assimilation, as implemented in the Last Millennium Reanalysis Hakim et al. (2016). Part of PReSto’s mission is to democratize these tools and enable a wider variety of actors, from seasoned researchers to students or citizen scientists, to generate their own reconstructions based on available code and data.

The purpose of this repository is to showcase how to use tools from the LinkedEarth Python research ecosystem (and broader scientific Python stack) to emulate and expand on the Last Millennium Reanalysis, version 2.1 (Tardif et al, 2019), which used the offline data assimilation method of Hakim et al. (2016) together with the PAGES 2k database, version 2.0.0.

The reconstruction workflow is broken down into 3 major steps, some of which have variants:

  1. Data assembly: gathering, selection and cleaning

  2. Data assimilation, which blends proxy observations with calibration data and the model prior

  3. Validation and comparison to other relevant reconstructions

Here we offer two different ways to carry out Step 1, both of which result in a netCDF file that can be used in Step 2:

  • Step 1a illustrates how to get a proxy database derived from the two pickle files (metadata and data) provided by (Tardif et al, 2019).

  • Step 1b illustrates how to get the proxy database directly from what is built into the cfr codebase.

Step 2 is common to all workflows, and only one instance is illustrated here (results obviously vary as a function of how Step 1 is implemented). Step 3 focuses here on emulating LMR v2.1, with validating both our results from Step 2, as well as comparing different proxy databases from Step 1 for forensics purposes. Newer reconstructions based on updates to the PAGES 2k database will be the subject of upcoming notebooks/publications.