Last Millennium Reanalysis (LMR) Project Data Page

LMR version 2 is complete and available for download. We are working with our partners at NOAA to distribute data from their servers, but for now you can access the data from this page.

LMR version 2 reference:

Tardif, R., G. J. Hakim, W. A. Perkins, K. A. Horlick, M. P. Erb, J. Emile-Geay, D. M. Anderson, E. J. Steig, and D. Noone, 2018: Last Millennium Reanalysis with an Expanded Proxy Database and Seasonal Proxy Modeling. Clim. Past, 14, submitted.

LMR version 1 reference:

Hakim, G. J., J. Emile-Geay, E. J. Steig, D. Noone, D. M. Anderson, R. Tardif, N. J. Steiger, and W. A. Perkins, 2016: The Last Millennium Climate Reanalysis Project: Framework and first results. J. Geophys. Res. Atmos. 121,  10.1002/2016JD024751

The source code for the LMR is open on Github. Web-based documentation describing the code also provides a page with access to data to perform basic reconstructions.

The data provided in the following links are for the LMR version 2 reconstruction in netcdf format. File and variable naming conventions follow those for the 20th Century Reanalysis.

Gridded fields have the format: (time, MCrun, lat, lon) where time is the year, lat is the latitude index, lon is the longitude index, and MCrun indicate the Monte Carlo iteration index. There are in fact 20 LMR reconstructions contained in these arrays. They differ in the climate model ensemble prior to assimilation (random draws from the CCSM4 Last Millennium simulation) and the proxies that were drawn randomly for the reconstruction (75% of all available proxies). You probably want to start with the "grand mean," which is an average over the MCrun dimension. NOTE: the variance on the MCrun index does NOT represent the uncertainty (error) in the reconstructions! It only represents the piece from the prior and proxies. The major portion of the error is in the full ensemble, which is summarized by the spread field for each iteration.

All fields are anomalies from the 1951--1981 time-mean.

Example plotting code coming soon (September 2018)

Two-meter air temperature:

Sea-surface temperature:

500 hPa geopotential height:


Mean-sea-level pressure:

Palmer drought severity index (PDSI):

Precipitable water:

Global-mean two-meter air temperature (format: (time, MCrun, members); this has all ensemble members):