Comparison of spatial downscaling methods of general circulation models to study climate variability during the Last Glacial Maximum
Abstract. The extent to which climate conditions influenced the spatial distribution of hominin populations in the past is highly debated. General Circulation Models (GCMs) and archaeological data have been used to address this issue. Most GCMs are not currently capable of simulating past surface climate conditions with sufficiently detailed spatial resolution to distinguish areas of potential hominin habitat, however. In this paper we propose a Statistical Downscaling Methods (SDM) for increasing the resolution of climate model outputs in a computationally efficient way. Our method uses a generalized additive model (GAM), calibrated over present-day data, to statistically downscale temperature and precipitation from the outputs of a GCM simulating the climate of the Last Glacial Maximum (19–23 000 BP) over Western Europe. Once the SDM is calibrated, we first interpolate the coarse-scale GCM outputs to the final resolution and then use the GAM to compute surface air temperature and precipitation levels using these interpolated GCM outputs and fine resolution geographical variables such as topography and distance from an ocean. The GAM acts as a transfer function, capturing non-linear relationships between variables at different spatial scales. We tested three different techniques for the first interpolation of GCM output: bilinear, bicubic, and kriging. The results were evaluated by comparing downscaled temperature and precipitation at local sites with paleoclimate reconstructions based on paleoclimate archives (archaeozoological and palynological data). Our results show that the simulated, downscaled temperature and precipitation values are in good agreement with paleoclimate reconstructions at local sites confirming that our method for producing fine-grained paleoclimate simulations suitable for conducting paleo-anthropological research is sound. In addition, the bilinear and bicubic interpolation techniques were shown to distort either the temporal variability or the values of the response variables, while the kriging method offers the best compromise. Since climate variability is an aspect of their environment to which human populations may have responded in the past this is an important distinction.