Here we provide time series of three isotope-enabled model simulations of the iHadCM3 Model in supplement to Bühler et al. (2021). Simulation output was extracted at SISAL v.2. sites for the last millennium (810-1850CE). The spatial resolution in the atmosphere is 2.5 by 3.75 degrees in the horizontal with 19 vertical levels, and 1.25 by 1.25 degrees in the ocean, with 20 irregularly spaced vertical levels. Details of the isotope-enabled model extension can be found in Tindall et al. (2009). The model simulations were performed with constant orbital forcing, and time-varying prescribed land cover changes, solar irradiance, volcanic aerosol optical depth and greenhouse gases (CO2, CH4, NO2). Changes in greenhouse gas forcings occur 50 years earlier than other forcings, which should not cause issues in interpretation due to the absence of large trends prior to 1800. We provide output for surface temperature, total precipitation and the precipitation-weighted oxygen isotope ratio extracted by bilinear interpolation at the cave site. To facilitate comparisons we provide both the temporally aggregated output (i.e. model output aggregated to the same temporal resolution as the speleothem data) as well as annual means.

These files contain the output of three ensemble members ('xnapa', 'xnapb' and 'xnapc') extracted for cave locations within the SISAL v.2. database (https://researchdata.reading.ac.uk/256/, Comas-Bru et al. (2020)). We include output for sites that pass the resolution and dating screening, meaning that have at least 2 radiometric dates (or are lamina-counted) and provide ten oxygen isotope ratio measurements within the last millennium. See the README.txt files for more details.