The paleoclimate simulations used here come from the HadCM3 version of the UK Met Office Unified Model General Circulation Model (GCM). This is a well-established coupled ocean atmosphere climate model, having contributed to the last three Intergovernmental Panel on Climate Change (IPCC) Assessment Reports (AR3, AR4 and AR5), and used to simulate climate for nearly 20 years. HadCM3 has a horizontal resolution of 2.5° in latitude and 3.75° in longitude, and higher resolution ocean of 1.25° x 1.25° regular long-lat grid, with 19 vertical levels in atmosphere and 20 in the ocean. The atmospheric component has a time-step of 30 minutes, and is coupled to the ocean every day. Typically, the climatology is output every month, and the mean annual and monthly climate are calculated from these data. As the name GCM suggests, this class of climate model is able to reproduce the major circulations in the both the atmosphere and ocean, as well as major drivers of inter-annual variability. The resolution also allows for synoptic weather patterns to be simulated, along with key climate oscillations, but may not simulate well local extremes or regions with high gradients (e.g. extreme convective events). HadCM3 is in the middle of the range of overall climate sensitivities exhibited by the IPCC-class climate models.
If simulations represent a composite of a continuous period (i.e. the data from PaleoView), mean/min/max climate values were calculated in 50 year window incriments, which were subsequently processed for the whole focal time period (e.g. for the Heinrich Stadial 1, which occured from 17.0-14.7 ka, all 50 yr composites within this period were aveaged).
The paleoclimate simulations are mathematical models of the general circulation of a planetary atmosphere or ocean. They are massively complex simulations of weather on Earth during as the result of global conditions at that time period, including: location and shape of terrestrial landmasses, CO2/N2O/CH4 levels, and orbital parameters
We employed the Change-Factor method to downscale the paleoclimatic climatologies. This approach creates high-resolution layers by quantifying the differences between the paleo and current (control) climatologies for each raw variable, at the native model-specific spatial resolution. This functions as a calibration step to measure the raw climate anomalies at the coarser spatial scale climate model. Once this step is completed, the difference layers (commonly called delta layers, change-factor differences, or climate change anomalies) are downscaled to high-resolutions (here 5-20km) and summed to a matching high-resolution current climate variable. This method is relatively quick, requiring less than a day of computational time per raster layer, and can be efficiently applied to global datasets. A major benefit of the Change-Factor method relative to other methods of downscaling is its ability to incorporate small-scale topographic nuances in regional climatologies that are often not captured in climate models, but present in the high-resolution current datasets. Examples include climatic differences in mountainous regions such as differences between valleys, mid-elevation ranges, and their peaks.
Here, we created global delta layers by subtracting the raw temperature and precipitation values of each snapshot paleoclimatic simulation from corresponding HadCM3 control simulations that represent the pre-industrial era. The delta layer represents the pixel-by-pixel changes from pre-industrial conditions, within the constraints of each snapshot climate simulation. The delta layers were downscaled 60 fold from 2.5 arc-degrees to 2.5 arc-minutes (ca. 5km) using a tensioned spline in ArcGIS 10.5 (sampling=12 nearest observations to a focal point, weight of 0.1, ESRI 2018). A spline is a deterministic interpolation method that ihas been commonly considered as appropriate for interpolation environmental variables. We used a tensioned spline (instead of a regularized spline) to avoid extraneous inflection points, and more generally to preserve shape properties, such as monotonicity and convexity, of a set of data points – and to do so without sacrificing smoothness . Spline approaches are based on requirement that the interpolation function passes through the data points, but also yield the smoothest transition as possible.
The high-resolution delta layers were then summed to a corresponding current monthly temperature or precipitation climate layers from the Climatologies at High-Resolution for the Earth’s Land Surface Areas (CHELSA) database, at the same resolution (see download links here). Though rare in our analyses, negative precipitation values were converted to zero. To reduce pixel-depth and file sizes of final products, all monthly temperature raster layers were multiplied by 10 and converted to integers. Prior to the creation of bioclimate layers, final monthly layers were adjusted to the mean sea-level of paleoclimatic period, based on adjustments to a contemporary bathymetry dataset.
Downscaling creates high resolution datasets useful for biological modeling from the coarse data output from paleoclimatic simulations
From the high-resolution monthly temperature and precipitation values, we calculated a set of derived parameters broadly used in ecological applications. These bioclimatic variables are derived from the monthly mean temperature (or minimum and maximum temperature, depending on their availability) and precipitation values. They are specifically developed for species distribution modelling and related ecological applications . The procedure for generating bioclimatic variables followed WorldClim and used the ‘biovars’ function of the R package dismo. Output bioclimate layers were saved as individual GeoTiffs (*tif) and projected in the WGS 1984 projection.
Bio_1=Annual Mean Temperature [°C*10]
Bio_2=Mean Diurnal Range [°C]
Bio_4=Temperature Seasonality [standard deviation*100]
Bio_5=Max Temperature of Warmest Month [°C*10]
Bio_6=Min Temperature of Coldest Month [°C*10]
Bio_7=Temperature Annual Range [°C*10]
Bio_8=Mean Temperature of Wettest Quarter [°C*10]
Bio_9=Mean Temperature of Driest Quarter [°C*10]
Bio_10=Mean Temperature of Warmest Quarter [°C*10]
Bio_11=Mean Temperature of Coldest Quarter [°C*10]
Bio_12=Annual Precipitation [mm/year]
Bio_13=Precipitation of Wettest Month [mm/month]
Bio_14=Precipitation of Driest Month [mm/month]
Bio_15=Precipitation Seasonality [coefficient of variation]
Bio_16=Precipitation of Wettest Quarter [mm/quarter]
Bio_17=Precipitation of Driest Quarter [mm/quarter]
Bio_18=Precipitation of Warmest Quarter [mm/quarter]
Bio_19=Precipitation of Coldest Quarter [mm/quarter]
Bioclimatic variables, or bioclims, are derived from the monthly temperature and rainfall values in order to generate more biologically meaningful variables
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25, 1965-1978 (2005).
Gordon, C. et al. The simulation of SST, sea ice extents and ocean heat transports in a version of the Hadley Centre coupled model without flux adjustments. Climate dynamics 16, 147-168 (2000).
Pope, V., Gallani, M., Rowntree, P. & Stratton, R. The impact of new physical parametrizations in the Hadley Centre climate model: HadAM3. Climate Dynamics 16, 123-146 (2000).
Randall, D. A. et al. Climate Models and Their Evaluation. (2007).
Kendon, E. J. et al. Do convection-permitting regional climate models improve projections of future precipitation change? Bulletin of the American Meteorological Society 98, 79-93 (2017).
Lüthi, D. et al. High-resolution carbon dioxide concentration record 650,000–800,000 years before present. Nature 453, 379 (2008).
Loulergue, L. et al. Orbital and millennial-scale features of atmospheric CH 4 over the past 800,000 years. Nature 453, 383 (2008).
Spahni, R. et al. Atmospheric methane and nitrous oxide of the late Pleistocene from Antarctic ice cores. Science 310, 1317-1321 (2005).
Laskar, J. et al. A long-term numerical solution for the insolation quantities of the Earth. Astronomy & Astrophysics 428, 261-285 (2004).
Bragg, F., Lunt, D. & Haywood, A. Mid-Pliocene climate modelled using the UK hadley centre model: PlioMIP experiments 1 and 2. Geoscientific Model Development 5, 1109 (2012).
Hill, D. J. The non-analogue nature of Pliocene temperature gradients. Earth and Planetary Science Letters 425, 232-241 (2015).
Dowsett, H. et al. The PRISM3D paleoenvironmental reconstruction. Stratigraphy 7, 123-139 (2010).
Salzmann, U., Haywood, A., Lunt, D., Valdes, P. & Hill, D. A new global biome reconstruction and data‐model comparison for the middle Pliocene. Global Ecology and Biogeography 17, 432-447 (2008).
Sohl, L. et al. PRISM3/GISS topographic reconstruction: US Geological Survey Data Series 419. US Geological Survey, Reston VA (2009).
Haywood, A. M. et al. On the identification of a Pliocene time slice for data–model comparison. Phil. Trans. R. Soc. A 371, 20120515 (2013).
Lisiecki, L. E. & Raymo, M. E. A Pliocene‐Pleistocene stack of 57 globally distributed benthic δ18O records. Paleoceanography 20 (2005).
Dolan, A. M. et al. Modelling the enigmatic Late Pliocene Glacial Event—Marine Isotope Stage M2. Global and Planetary Change 128, 47-60 (2015).
Tan, N. et al. Exploring the MIS M2 glaciation occurring during a warm and high atmospheric CO2 Pliocene background climate. Earth and Planetary Science Letters 472, 266-276 (2017).
Bartoli, G., Hönisch, B. & Zeebe, R. E. Atmospheric CO2 decline during the Pliocene intensification of Northern Hemisphere glaciations. Paleoceanography 26 (2011).
Haywood, A. M. et al. Large-scale features of Pliocene climate: results from the Pliocene Model Intercomparison Project. Clim. Past 9, 191-209, doi:10.5194/cp-9-191-2013 (2013).
Wilby, R. et al. Guidelines For Use of Climate Scenarios Developed From Statistical Downscaling Methods. (2004).
Lima-Ribeiro, M. et al. EcoClimate: a database of climate data from multiple models for past, present, and future for macroecologists and biogeographers. Vol. 10 (2015).
Laslett, G. M. Kriging and Splines: An Empirical Comparison of Their Predictive Performance in Some Applications. Journal of the American Statistical Association 89, 391-400, doi:10.2307/2290837 (1994).
Dubrule, O. Comparing splines and kriging. Computers & Geosciences 10, 327-338, doi:https://doi.org/10.1016/0098-3004(84)90030-X (1984).
Hutchinson, M. F. Interpolating mean rainfall using thin plate smoothing splines. International Journal of Geographical Information Systems 9, 385-403, doi:10.1080/02693799508902045 (1995).
Laslett, G. M., McBratney, A. B., Pahl, P. J. & Hutchinson, M. F. Comparison of several spatial prediction methods for soil pH. Journal of Soil Science 38, 325-341, doi:10.1111/j.1365-2389.1987.tb02148.x (1987).
Karger, D. N. et al. Climatologies at high resolution for the earth’s land surface areas. Scientific Data 4, 170122, doi:10.1038/sdata.2017.122 (2017).
Weatherall, P. et al. A new digital bathymetric model of the world’s oceans. Earth and Space Science 2, 331-345, doi:10.1002/2015EA000107 (2015).
Thompson, R. S. & Fleming, R. F. Middle Pliocene vegetation: reconstructions, paleoclimatic inferences, and boundary conditions for climate modeling. Marine Micropaleontology 27, 27-49, doi:https://doi.org/10.1016/0377-8398(95)00051-8 (1996).
Markwick, P. J. The palaeogeographic and palaeoclimatic significance of climate proxies for data-model comparisons. (2007).
De Boer, B., Van de Wal, R., Bintanja, R., Lourens, L. & Tuenter, E. Cenozoic global ice-volume and temperature simulations with 1-D ice-sheet models forced by benthic δ 18 O records. Annals of Glaciology 51, 23-33 (2010).
Fordham, D. A., Saltré, F. , Haythorne, S. , Wigley, T. M., Otto‐Bliesner, B. L., Chan, K. C. and Brook, B. W. (2017), PaleoView: a tool for generating continuous climate projections spanning the last 21 000 years at regional and global scales. Ecography, 40: 1348-1358. doi:10.1111/ecog.03031