Motivation: In order to understand how species evolutionarily responded to Plio-Pleistocene climate oscillations (e.g. in terms of speciation, extinction, migration and adaptation), it is first important to have a good understanding of those past climate changes per se. This, however, is currently limited due to the lack of global-scale climatic datasets with high temporal resolution spanning the Plio-Pleistocene. To fill this gap, I here present Oscillayers, a global-scale and region-specific bioclim dataset, facilitating the study of climatic oscillations during the last 5.4 million years at high spatial (2.5 arc-minutes) and temporal (10 kyr time periods) resolution. This data set builds upon interpolated anomalies (Δ layers) between bioclim layers of the present and the Last Glacial Maximum (LGM) that are scaled relative to the Plio-Pleistocene global mean temperature curve, derived from benthic stable oxygen isotope ratios, to generate bioclim variables for 539 time periods. Evaluation of the scaled, interpolated estimates of palaeo-climates generated for the Holocene, Last Interglacial and Pliocene showed good agreement with independent General Circulation Models (GCMs) for respective time periods in terms of pattern correlation and absolute differences. Oscillayers thus provides a new tool for studying spatial-temporal patterns of evolutionary and ecological processes at high temporal and spatial resolution. Main types of variable contained: 19 bioclim variables for time periods throughout the Plio-Pleistocene. Input data and R script to recreate all 19 bioclim variables. Spatial location and grain: Global at 2.5 arc-minutes (4.65 x 4.65 = 21.62 km2 at the equator). Time period and grain: The last 5.4 million years. The grain is 10 kyr (= 539 time periods). Level of measurement: Data are for terrestrial climates (excluding Antarctica) taking sea level changes into account. Software format: All data are available as ASCII (ESRI) grid files.