Holocene DA Validation Report

Validation of the custom Holocene Data Assimilation reconstruction against published Holocene global temperature reconstructions. Reference datasets are discovered at runtime from reference_data/ — any CSV with columns age_BP, median, q05, q95 or age_BP, anomaly, uncertainty_1sigma will be used.

0.067
6 ka ΔT (°C, area-weighted)
0.967
GMST R vs erb2022_gmst
0.927
GMST CE vs erb2022_gmst
0.803
GMST R vs kaufman2020_temp12k
-1.169
GMST CE vs kaufman2020_temp12k
0.955
Spatial R vs erb2022
0.091
Proxy R (temp12k_v102_proxies, N=1174)

GMST Validation Metrics

Pearson correlation (R) and Nash–Sutcliffe coefficient of efficiency (CE) of the ensemble-median GMST against each reference over their common age range. CE = 1 is perfect; CE = 0 equals climatology; CE < 0 is worse than climatology.

ReferenceOverlapRCE
erb2022_gmst 99.0–11899.0 BP (60 pts) 0.9671 0.9271
kaufman2020_temp12k 99.0–11899.0 BP (60 pts) 0.8030 -1.1691

Spatial Anomaly at 6 ka

Mean temperature anomaly for 5500–6500 BP relative to 0–1000 BP baseline. Robinson projection.

6 ka spatial anomaly map

Spatial Comparison at 6 ka

Side-by-side 6 ka temperature anomaly for the custom reconstruction and each spatial reference, plus their difference. Pattern correlation (R) and area-weighted RMSE are computed after nearest-neighbor regridding of the reference onto the reconstruction grid.

ReferencePattern RRMSE (°C)Recon geo meanRef geo mean
erb20220.95500.0647+0.067 °C+0.085 °C
6 ka spatial comparison vs erb2022

Proxy-Site Comparison (6 ka)

For each proxy site, the reconstruction’s 6 ka anomaly is sampled via nearest-neighbor on the model grid, then compared to the proxy’s own 6 ka anomaly (5500–6500 BP relative to 0–1000 BP baseline, falling back to the record mean if the modern baseline is absent). R, RMSE, and bias (recon − proxy) summarize the match.

DatasetNRRMSE (°C)Bias (°C)
temp12k_v102_proxies11740.09071.8643-0.2413
Proxy comparison (temp12k_v102_proxies)
Proxy Database Comparison vs Temp12k 1.0.2 (Erb et al. 2022) (shared 956, custom-only 59, reference-only 320)

Comparison of the proxy records in this run’s lipd_legacy.pkl against Temp12k 1.0.2 (Erb et al. 2022), the version used by the published Erb et al. 2022 reconstruction. Records are matched on paleoData_TSid after the same paleoData_inCompilation == Temp12k + paleoData_units == degC filter that da_load_proxies.py applies at runtime.

1,276
Temp12k 1.0.2 (Erb et al. 2022) records
1,015
Custom records
956
Shared
59
Custom-only
320
Temp12k 1.0.2 (Erb et al. 2022)-only

Side-by-side statistics

StatisticCustomTemp12k 1.0.2 (Erb et al. 2022)
Records1,0151,276
Distinct archive types68
Earliest record start (yr BP)1,277,3461,277,346
Latest record end (yr BP)-240-275
Median record length (yr)1169011403
Median observations per record7371

Records by archive type

ArchiveSharedCustom-only Temp12k 1.0.2 (Erb et al. 2022)-onlyTotal
LakeSediment51010184704
MarineSediment2624295399
Peat163428195
GlacierIce103518
Midden80210
Speleothem3025
Wood0033
Ice-other0011
Total956593201335

Spatial distribution

Proxy spatial comparison

Temporal coverage

Records covering each 250-yr age bin, partitioned by which database they belong to.

Temporal coverage comparison

Records exclusive to the custom run (59)

Records present in this reconstruction’s proxy database but absent from Temp12k 1.0.2 (Erb et al. 2022) — typically records added in later Temp12k versions or pulled from filtered queries.

TSIDArchiveDatasetAge range (BP)n_obs
GH2d963782MarineSedimentBS79_38.Cacho.2001867–24,51771
GH6edd69deMarineSedimentGeoB33131.Lamy.2002283–6,99075
LPD52d54dc4MarineSedimentGEOFAR_KF16_MgCa.Repschlager.2016524–11,25898
LPD5d932ca5MarineSedimentPP10-07.Mary.2016422–10,150206
LPD71287822MarineSedimentTTR17_434G.Rodrigo-Gamiz.201427–13,72892
LPD73260700MarineSedimentPP10-07.Mary.2016422–10,150206
LPD74bc6584LakeSedimentLittlePondRoyalston.Marsicek.20130–11,000221
LPD779b21dcLakeSedimentBloodPond.Marsicek.20130–11,000221
LPD7abab671-dupGlacierIceLawDome.Dahl-Jensen.1999-45–8,69235
LPD92d845a8MarineSedimentPP10-07.Mary.2016422–10,150206
LPD9c79485dxxxxMarineSedimentODP658C_UK37.Zhao.1995137–22,589203
LPDa0c55c63MarineSedimentPP10-07.Mary.2016422–10,150206
LPDadaae75fMarineSedimentSN-6.Tiwari.2015154–4,77236
LPDc82d344dMarineSedimentGEOFAR_KF16_MgCa.Repschlager.2016524–11,25898
PYT917S9C8JMarineSedimentGIK17748-2.Kim.20021,080–16,07054
PYTDSG9BNO2MarineSedimentGIK23258.Sarnthein.2003698–13,987285
PYTH91MXU0ZMarineSedimentP1003.Sejrup.2011-48–7,881649
PYTRS33LNIKMarineSedimentGIK23258.Sarnthein.2003698–13,987285
PYTUQPUNFJTPeatCorser.Nichols.2014103–11,11636
R60yDfufGWZPeatStaroselskyMoch.Novenko.2018-19–9,79256
R7erF84XRXyMarineSedimentDANA12_11_2_GC01.VanNieuwenhove.2018-54–11,81080
RCP5UgOcOkRLakeSedimentGunnarsfjorden.Allen.2007-40–9,124104
RGIi8rzDMmjMarineSedimentGeoB10042_1.Setiawan.20154,200–39,420106
RI2vCrHM2fwMarineSedimentM35003_4.R_hlemann.19990–24,54091
RJ64ej30dlDMarineSedimentMD03_2607-Assemblage.LopesdosSantos.20131,100–134,800172
RKz0xLW2LFhMarineSedimentMD06_3067.Bolliet.20113,390–158,200484
RMpGVpFWptOLakeSedimentDuranunlak.EPD368–11,50136
RNAQJ20Sh7jMarineSedimentRR1313_23PC.MoffaSanchez.20191,280–16,186120
RPZj5YKrFr0LakeSedimentPupukePollen.VandenBos.20181,001–16,13862
RQ6N1Ktr5AiLakeSedimentDuranunlak.EPD368–11,50136

⬇ Download full CSV (59 rows)

Records exclusive to Temp12k 1.0.2 (Erb et al. 2022) (320)

Records used by the published reconstruction but missing from this custom run — typically records dropped by the user’s filter or removed from later Temp12k versions.

TSIDArchiveDatasetAge range (BP)n_obs
GH06e59526MarineSedimentMD79-257.Levi.20071,554–15,69547
GH194ddf83PeatHomesteadScarp.McGlone.2010-39–12,247170
GH73726ccePeatMountHoney.McGlone.2010-18–12,563109
LPD27f34021PeatYakumo.Leipe.201340–5,50638
LPD3267fe5dMarineSedimentPP10-07.Mary.2016422–10,150206
LPD575ea390LakeSedimentSihailongwanMaar.Stebich.2015150–11,992206
LPD5c8693b4MarineSedimentMD01-2416.Sarnthien.20045,727–11,92041
LPD5c89a1b5LakeSedimentHAP0212.vanderBilt.2016534–11,11760
LPD6e0eacd1LakeSedimentLagoDellAccesa.Finsinger.20104,479–11,55286
LPD7abab671GlacierIceLawDome.Dahl-Jensen.1999-45–8,69235
LPD801dea7fMarineSedimentMD01-2446.Salguero.20142,290–20,01255
LPD8d84bf64WoodGreatBasin.Salzer.2014-56–4,5244,581
LPD91355717LakeSedimentLagoVerdarolo.Samartin.2011-59–13,72050
LPDbf7cf949PeatSouthernCalifornia.Ohlwein.2012408–11,44628
LPDbfac52bbLakeSedimentQinghai.Hou.20160–15,629184
LPDd2a984feLakeSedimentAMP112.vanderBilt.20161,144–12,35366
LPDf97892fcLakeSedimentXingyunLake.Wu.2018-27–14,19267
M2L1fcqp9rpPLakeSedimentEC1.Gajewski.199358–6,61032
M2L1r9snju2PLakeSedimentAtocas.Gauthier.19810–11,97741
M2L4ptir6l2PLakeSedimentBanksIsland-12.Gajewski.2000175–10,18333
M2L68fmqq0lPLakeSedimentAtocas.Gauthier.19810–11,97741
M2L68iw0l2ePLakeSedimentNorthPond.Whitehead.1978-105–13,08355
M2L6qbln98nPLakeSedimentLoneFoxLake.MacDonald.1985-32–12,63837
M2L726ed2xjPLakeSedimentCH2.Garralla.19921–8,33289
M2L7thzxje2PLakeSedimentHams.Bennett.1987-29–11,94582
M2L8sc7gn3mPLakeSedimentHams.Bennett.1987-29–11,94582
M2LdjjxzkgtPLakeSedimentQueens.Moser.19900–8,84033
M2Lfut6kmdqPLakeSedimentQueens.Moser.19900–8,84033
M2Lg77nq75zPLakeSedimentEC1.Gajewski.199358–6,61032
M2Lh4fhttxhPLakeSedimentCH2.Garralla.19921–8,33289

⬇ Download full CSV (320 rows)

GMST Time Series

Custom reconstruction ensemble spread alongside reference medians. X-axis runs from the oldest age (left) to present (right).

GMST time series

GMST Ensemble Members (52 total)

Every ensemble member plotted individually (subsampled to 200 for readability), showing the full reconstruction spread.

GMST ensemble members

GMST Difference (Custom − erb2022_gmst)

Year-by-year difference between the custom reconstruction median and the primary reference. Red = warmer, Blue = cooler.

GMST difference plot

← Back to results