Present and future Köppen-Geiger climate classification maps at 1-km resolution

Abstract

We present new global maps of the Köppen-Geiger climate classification at an unprecedented 1-km resolution for the present-day (1980–2016) and for projected future conditions (2071–2100) under climate change. The present-day map is derived from an ensemble of four high-resolution, topographically-corrected climatic maps. The future map is derived from an ensemble of 32 climate model projections (scenario RCP8.5), by superimposing the projected climate change anomaly on the baseline high-resolution climatic maps. For both time periods we calculate confidence levels from the ensemble spread, providing valuable indications of the reliability of the classifications. The new maps exhibit a higher classification accuracy and substantially more detail than previous maps, particularly in regions with sharp spatial or elevation gradients. We anticipate the new maps will be useful for numerous applications, including species and vegetation distribution modeling. The new maps including the associated confidence maps are freely available via www.gloh2o.org/koppen.

Design Type(s)	modeling and simulation objective
Measurement Type(s)	climate change
Technology Type(s)	computational modeling technique
Factor Type(s)
Sample Characteristic(s)	Earth (Planet) • climate system

Machine-accessible metadata file describing the reported data (ISA-Tab format)

Background & Summary

The Köppen-Geiger system classifies climate into five main classes and 30 sub-types. The classification is based on threshold values and seasonality of monthly air temperature and precipitation. Considering vegetation as “crystallized, visible climate”¹, this classification aims to empirically map biome distributions around the world: different regions in a similar class share common vegetation characteristics. The first version of this classification was developed in the late 19th century²; it is still widely used today, for many applications and studies conditioned on differences in climatic regimes, such as ecological modeling or climate change impact assessments^3–8. This wide use reflects the fact that climate has since long been recognized as the major driver of global vegetation distribution^9–11. In species distribution models¹², climate variables are considered the primary driver to explain species ranges at larger spatial extents, while habitat and topography are considered to only be modifiers of plant species distributions at smaller extents^13–15. The Köppen-Geiger climate classification is a highly suitable means to aggregate complex climate gradients into a simple but ecologically meaningful classification scheme. It is therefore often used as input when analyzing the distribution^4,16,17 or growth behavior¹⁸ of species, or to set-up dynamic global vegetation models¹⁹.

Three recent versions of the world maps of the Köppen-Geiger climate classification exist^20–22. Kottek et al.²⁰ produced a map (0.5° resolution) based on CRU TS 2.1²³ for temperature and VASClimO V1.1²⁴ for precipitation. CRU was based on approximately 7000–17,000 stations (depending on the year) and VASClimO on 9343 stations. The Peel et al. map²¹ (0.1° resolution) was derived from 4844 air temperature stations and 12,396 precipitation stations. Kriticos et al.²² produced a map (0.083° resolution) based on WorldClim V1 temperature and precipitation datasets²⁵, which are based on 24,542 and 47,554 stations, respectively.

All maps have a relatively low resolution (≥0.1°) and the map of Peel et al.²¹ has not been explicitly corrected for topographic effects, which influences air temperature²⁶ and precipitation²⁷ in mountainous regions. In addition, the maps of Kottek et al.²⁰ and Peel et al.²¹ are based on a relatively small number of stations. This can lead to widespread misclassifications, particularly in regions with a low station density and/or strong climatic gradients such as mountain ranges²⁸. Moreover, since these maps do not include corresponding uncertainty estimates, they may provide users a false sense of confidence.

Here, we present a new and improved Köppen-Geiger climate classification map for the present (1980–2016) with an unprecedented 0.0083° resolution (approximately 1 km at the equator), providing more accurate representation of highly heterogeneous regions (Fig. 1a). To maximize the accuracy and assess uncertainties in map classifications, we combine climatic air temperature and precipitation data from multiple independent sources, including WorldClim V1 and V2, CHELSA V1.2, and CHPclim V1 (Table 1). These datasets have all been explicitly corrected for topographic effects and, with the exception of the CHELSA V1.2 temperature dataset, been based on a large number of stations (≥34,542 for precipitation and ≥20,268 for temperature). The use of multiple data sources allows us to provide an estimate of uncertainty in the derived classes. Further, we combine climate change projections from 32 Coupled Model Intercomparison Project phase 5 (CMIP5²⁹) models to map future (2071–2100) climate classes at the same spatial resolution (Fig. 1b).

Figure 1: New and improved Köppen-Geiger classifications.

Part (a) shows the present-day map (1980–2016) and panel (b) the future map (2071–2100). The color scheme was adopted from Peel et al.²¹.

Table 1 Global monthly datasets used for deriving the Köppen-Geiger maps.

Methods

Köppen-Geiger climate classification

We follow the Köppen-Geiger climate classification as described in Peel et al.²¹, which was also used by Kriticos et al.²² (Table 2). This classification is identical to that presented by Köppen in 1936¹ with three differences. First, temperate (C) and cold (D) climates are distinguished using a 0 °C threshold instead of a 3 °C threshold, following the suggestion of Russell³⁰. Second, the arid (B) sub-climates W (desert) and S (steppe) were identified depending on whether 70% of precipitation occurred in summer or winter. Third, the sub-climates s (dry summer) and w (dry winter) within the C and D climates were made mutually exclusive by assigning s when more precipitation falls in winter than in summer and assigning w otherwise. Note that the tropical (A), temperate (C), cold (D), and polar (E) climates are mutually exclusive but may intersect with the arid (B) class. To account for this, climate type B was given precedence over the other classes.

Table 2 Overview of the Köppen-Geiger climate classes including the defining criteria.

Climate data

The present Köppen-Geiger classification map was derived from three climatic datasets for air temperature (WorldClim V1 and V2, and CHELSA V1.2) and four climatic datasets for precipitation (WorldClim V1 and V2, CHELSA V1.2, and CHPclim V1; Table 1). All datasets have a 0.0083° resolution with the exception of CHPclim V1.2, which has a 0.05° resolution. For consistency CHPclim V1.2 was downscaled to 0.0083° using bilinear interpolation.

The future Köppen-Geiger classification was produced using monthly historical and future air temperature and precipitation data from the CMIP5 archive²⁹. For the future scenario, we used Representative Concentration Pathway 8.5 (RCP8.5³¹). All climate models with data during the 1980–2016 and 2071–2100 periods were used. Data for 1980–2016 was derived by concatenating historical runs (which end in 2005) and future runs (which begin in 2006). For each model, we only considered a single initialization ensemble. In total 32 models had sufficient data and hence were used for deriving the future map (Table 1).

Present-day Köppen-Geiger map

The present-day Köppen-Geiger map (Fig. 1a) was derived from an ensemble of high-resolution climatic datasets (Table 1) using the criteria listed in Table 2. Since the climatic datasets have inconsistent temporal coverages, we first adjusted them to reflect the period 1980–2016. To this end, we calculated, for each climatic dataset, monthly 0.5° climatologies for temperature using CRU TS V4.01 and for precipitation using GPCC FDR V7, for both the 1980–2016 period and the temporal span of the climatic dataset. Next, for each month we calculated climate change offsets (for temperature) or factors (for precipitation) between the two periods, and resampled these offsets or factors from 0.5° to 0.0083° resolution using bilinear interpolation, and adjusted the climatic maps by addition (for temperature) or multiplication (for precipitation).

For each adjusted temperature and precipitation climatic dataset combination, we derived a Köppen-Geiger map at 0.0083° resolution. From this ensemble of 4×3=12 maps we derived a final map by selecting, for each grid-cell, the most common class (Fig. 1a). A corresponding confidence map was derived by dividing the frequency of occurrence of the most common class by the ensemble size and converting these fractions to percentages (Fig. 2a). For example, if Csa is the most common class for a particular grid-cell, and it has been assigned eight times out of 12, the resulting confidence level is $100\times \frac{8}{12}=66.6$%. This confidence level should be interpreted as the degree of trust we place in our final present-day classification. Confidence levels are generally lower in the vicinity of borders between climate zones, in particular at high latitudes where the climatic data show more uncertainty.

Figure 2: The confidence levels (%) associated with the new Köppen-Geiger classifications.

Part (a) shows the present-day confidence map (1980–2016) and panel (b) the future confidence map (2071–2100). These maps provide an indication of classification accuracy.

Future Köppen-Geiger map

The future Köppen-Geiger map (Fig. 1b) was derived by the so-called “anomaly method”³² based on an ensemble of climate projections from the 32 CMIP5 models (Table 1). First, observed monthly present-day reference temperature and precipitation climatologies (0.0083° resolution) were derived, by simple averaging of the ensemble of temporally-homogenized, high-resolution climatic maps. Then, for each climate model and each month, we subsequently calculated climate change offsets (for temperature) or factors (for precipitation) between 1980–2016 and 2071–2100 and resampled these offsets or factors from the native model resolution to 0.0083° using bilinear interpolation (Fig. 3). Finally, future high-resolution climatic temperature and precipitation maps were derived from the present-day, observed reference climatologies by addition of the offsets (for temperature) or multiplication by the factors (for precipitation). We want to emphasize that the change factors are never excessively high (i.e., >5; Fig. 3), because (i) model simulations tend to overestimate the precipitation frequency³³ (resulting in the near-absence of areas with close to zero precipitation), and (ii) over the majority of arid regions the future projections tend toward drying rather than wetting³⁴ (resulting in factors <1).

Figure 3: Projected changes in mean air temperature (°C) and precipitation (unitless) between 1980–2016 and 2071–2100 derived from climate model outputs.

Part (a) presents air temperature change offsets and part (b) precipitation change factors. The values represent the mean over all models and months.

For each climate model, we derived a future Köppen-Geiger map at 0.0083° resolution from the downscaled future temperature and precipitation data. From this ensemble of 32 maps we derived a final map by selecting, for each grid-cell, the most common class (Fig. 1b). A corresponding confidence map was derived by dividing the frequency of occurrence of the most common class by the ensemble size and converting these fractions to percentages (Fig. 2b). For example, if Cfa is the most common class for particular grid-cell, and it has been assigned 24 times out of 32, the corresponding confidence level is $100\times \frac{24}{32}=75.0$%. This confidence level should be interpreted as the degree of trust we have in our final future classification based on the uncertainties in climate change projections. Thus, uncertainties are larger than for the present-day map. In particular, they are larger at high latitudes because of the greater model spread in projected warming in those regions.

Code availability

The new Köppen-Geiger classifications have been produced using MathWorks MATLAB version R2017a. The function used to classify the temperature and precipitation data according to the criteria listed in Table 2 (KoppenGeiger.m) is freely available via (Data Citation 1) and www.gloh2o.org/koppen. The other codes are available upon request from the first author.

Data Records

The present and future Köppen-Geiger classification maps and the corresponding confidence maps are freely available for download at (Data Citation 1) and www.gloh2o.org/koppen. The maps are stored in GeoTIFF format as unsigned 8-bit integers. We also provide a legend file (legend.txt) linking the numeric values in the maps to the Köppen-Geiger climate symbols and providing the color scheme used for displaying the maps in the current study (adapted from Peel et al.²¹). The maps are referenced to the World Geodetic Reference System 1984 (WGS 84) ellipsoid and made available at three resolutions (0.0083°, 0.083°, and 0.5°; approximately 1 km, 10 km, and 50 km at the equator, respectively). The classifications are upscaled from 0.0083° to 0.083° and 0.5° using majority resampling and the confidence levels using bilinear averaging. Table 3 presents the file naming convention. The maps can be visualized and analyzed using most Geographic Information Systems (GIS) software (e.g., QGIS, ArcGIS, and GRASS).

Table 3 File naming convention.

Technical Validation

We validated the new high-resolution present-day Köppen-Geiger classification (Fig. 1a), and previous maps from Kottek et al.²⁰, Peel et al.²¹, and Kriticos et al.²², by calculating the classification accuracy (defined as the percentage of correct classes) using station observations as reference. An initial database was compiled from the Global Historical Climatology Network-Daily (GHCN-D) database³⁵ and the Global Summary Of the Day (GSOD) database (https://data.noaa.gov). For each station, we calculated monthly mean temperature and precipitation time series (discarding months with <25 daily values), and subsequently monthly climatologies by averaging the monthly means (if ≥10 values were present). Stations with gaps in the climatologies or missing data for one of the four maps were discarded, resulting in a final dataset comprising 22,078 stations which we used to calculate the classification accuracy of each map.

The newly derived high-resolution present-day Köppen-Geiger classification (Fig. 1a) exhibited a classification accuracy of 80.0%, while the maps of Kottek et al.²⁰, Peel et al.²¹, and Kriticos et al.²² exhibited classification accuracies of 66.1, 70.9, and 73.4%, respectively. These results confirm that the new map is more accurate, which is primarily due to its high (1 km) resolution and use of an ensemble of topographically-corrected climatic datasets. The map of Kottek et al.²⁰ showed the lowest classification accuracy, due to its low (0.5°) resolution. The map of Peel et al.²¹ also performed less well, due to a lack of topographic corrections and the use of a relatively small number of stations.

We also tested the usefulness of the confidence map associated with the new present-day classification (Fig. 2a) using station observations. We obtained a mean confidence level of 92.6% for the correctly classified stations (n=17,667) and 77.4% for the misclassified stations (n=4411). The mean confidence level was thus substantially lower for the misclassified stations, confirming that the confidence map provides a useful indication of the classification accuracy.

Figures 4 and 5 show historic Köppen-Geiger classification maps from all three previous studies and our present-day map for the Alps (Europe) and the central Rocky Mountains (North America), respectively, illustrating the enhanced detail in our map. The other maps sometimes fail to depict important topographic features; the map of Peel et al.²¹, for example, does not represent the Apennine mountains (Italy), due to a lack of topographic corrections (Fig. 4e). The new map (Figs. 4a and 5a) also exhibits better agreement with a Landsat-based forest cover map³⁶ (30-m resolution; Figs. 4c and 5c). The spatial extent of the polar (E) climate, for example, corresponds closely with treelines in the forest cover maps. Additionally, the new present-day and future Köppen-Geiger maps (Fig. 4a and g, respectively) agree well with equivalent high-resolution maps derived for the Alps⁸ (their Figs. 1 and 2, respectively).

Figure 4: Köppen-Geiger classifications, and associated maps, for the European Alps.

Part (a) present-day results from our study (1980–2016); (b) confidence levels associated with our present-day map; and (c) forest cover map³⁶ (2000). Historic Köppen-Geiger classification maps for the three previous studies are provided as: (d) Kottek et al.²⁰ (1951–2000); (e) Peel et al.²¹ (1916–1992); and (f) Kriticos et al.²² (1960–1990). Our future Köppen-Geiger map (2071–2100) is presented in (g) and the corresponding confidence map in (h). The representative period of each map is listed in parentheses. Thin black lines are country borders and the unmapped white area is part of the Mediterranean Sea.

Figure 5: Köppen-Geiger classifications, and associated maps, for the central Rocky Mountains (North America).

Part (a) present-day results from our study (1980–2016); (b) confidence levels associated with our present-day map; and (c) forest cover map³⁶ (2000). Historic Köppen-Geiger classification maps for the three previous studies are provided as: (d) Kottek et al.²⁰ (1951–2000); (e) Peel et al.²¹ (1916–1992); and (f) Kriticos et al.²² (1960–1990). Our future Köppen-Geiger map (2071–2100) is presented in (g) and the corresponding confidence map in (h). The representative period of each map is listed in parentheses. The thick black line at 49° latitude represents the Canada-US border.

Usage Notes

The future Köppen-Geiger classification (Fig. 1b) should be viewed as providing insights into potential spatial changes in regional climatic zones under climate change. However, caution should be exerted not to equate those changes directly with changes in actual biomes. First, vegetation changes by 2100 may lag the change in climate zones. Secondly, factors not accounted for in the Köppen-Geiger classification, such as higher atmospheric CO₂ levels, may alter the relationship between climate classes and vegetation. It is thus advised to interpret the future Köppen-Geiger classification first and foremost from a ‘climatic conditions’ perspective.

The rationale for using the anomaly method to build future maps using climate model projections, instead of directly computing present and future maps from model outputs, is that superimposing future modeled anomalies onto the observed climate removes mean biases from climate model outputs. This is a widely used method in climate change impact assessments³². However, an unavoidable limitation of this approach is that because of model spatial biases, modeled climate change anomalies may not be fully geographically consistent with the baseline observed climatology to which they are added³⁷ (e.g., if the climate of one region in a given model is spatially shifted relative to reality).

Another irredeemable limitation is that because of their coarser resolution (typically 1–2°), climate model outputs do not resolve future climate change at the same scale as our baseline climatology. Thus, in cases where there could be significant heterogeneities in precipitation change and/or warming below the model resolution (e.g., along coastlines and/or in regions with strong land-cover differences and/or elevation gradients), future changes at the 1-km scale might be under- or over-estimated, because only the model-scale mean anomalies are used to compute future changes. High-elevation mountainous regions are a prime example of this because they are expected to experience considerably more warming than adjacent valleys³⁸.