...
This document describes the "Copernicus Interactive Climate Atlas gridded monthly dataset" which constitutes an evolution of the IPCC Atlas gridded monthly dataset recently published in the CDS1. This dataset has been entirely produced using a workflow ensuring total reproducibility within the CDS. The climate information used to construct the dataset is fully available in C3S CDS. The dataset has a corresponding visualisation tool, which is the Copernicus Interactive Climate Atlas, available at https://atlas.climate.copernicus.eu/atlas.
In addition to the description of the dataset, we are summarising the fit-for-purpose quality control applied to the datasets, documenting and curating the dataset based on the results from basic products (maps and climate stripes), and highlighting and fixing any issue encountered.
Info | ||
---|---|---|
| ||
1 Copernicus Climate Change Service (C3S) (2023): Gridded monthly climate projection dataset underpinning the IPCC AR6 Interactive Atlas. C3S Climate Data Store (CDS). https://doi.org/10.24381/cds.5292a2b0 |
...
Table 1: Descriptions of the datasets and CDS catalogue entries used. Color codes indicate the different temporal resolutions (hourly: green, daily: blue and monthly: purple) of the datasets available in the CDS catalogues (one or several for each catalogue). The column “horizontal resolution” indicates the grids used for the different datasets (products regridded to the common regular sub-grids: 2°, 1°, 0.5°, 0.25°, 0.125°, are indicated with an asterisk).
Project | Datasets | Horizontal resolution | Time period | CDS-catalogues |
CMIP6 | CMIP6 | 1° (*) | Historical: 1850 - 2014 Scenarios: 2015 - 2100 | Projections-cmip6 [daily/monthly] |
CMIP5 | CMIP5 | 2° (*) | Historical: 1850 - 2005 Scenarios: 2006 - 2100 | Projections-cmip5-daily-single-levels [daily] Projections-cmip5-monthly-single-levels [monthly] |
CORDEX | CORDEX-CORE | 0.25° (*) | Historical: 1970 - 2005 Scenarios: 2006 - 2100 | Projections-cordex-domains-single-levels [daily/monthly] |
CORDEX-EUR-11 | 0.125° (*) | Historical: 1970 - 2005 Scenarios: 2006 - 2100 | ||
ERA5 | ERA5 | 0.25° | 1940 - 2022 | Reanalysis-era5-single-levels [hourly] Reanalysis-era5-single-levels-monthly-means [monthly] |
ERA5-Land | ERA5-Land | 0.1° | 1950 - 2022 | Reanalysis-era5-land [hourly] Reanalysis-era5-land-monthly-means [monthly] |
E-OBS | E-OBS | 0.125° (*) | 1950 - 2021 | Insitu-gridded-observations-europe (0.1.deg) [daily] |
ORAS5 | ORAS5 | 0.25° (*) | 1958 - 2014 | Reanalysis-oras5 (consolidated, single levels) [daily] |
Data harmonization across the different datasets and catalogues is crucial in this contract. For instance, ERA5 and ERA5-Land do not provide directly daily data, so this have to be computed aggregating the hourly values. However, both datasets follow different conventions for coding hourly data (ERA5 provides hourly accumulated values whereas ERA5-Land provides accumulated values since the beginning of the forecast, at 00 for a period of 24 hours) and require different post-processing to obtain daily values https://confluence.ecmwf.int/pages/viewpage.action?pageId=197702790 . This needs to be taken into account when harmonizing the different datasets in order to obtain a standard common definition and units for each of the variables (Sec. 2.3).
...
Table 2: Description of the 21 climate variables and indices included in the IPCC Interactive Atlas, as well as the additional 9 new ones for the first version of the C3S Interactive Atlas dataset (indicated with "*" after the sequence number). Color codes are included for variables/indices characterizing different Climatic Impact Driver (CID) categories: heat and cold (red), wet and dry (green), wind and radiation (black), snow and ice (violet), coastal and open ocean (blue), circulation (grey). Italics indicates the indices derived from other variables (Table 4).
CID category | # | Code | Index | Units |
heat and cold | 1 | t | Monthly mean of daily mean temperature | °C |
2 | tn | Monthly mean of daily minimum temperature | °C | |
3 | tx | Monthly mean of daily maximum temperature | °C | |
4 | tnn | Monthly minimum of daily minimum temperature | °C | |
5 | txx | Monthly maximum of daily maximum temperature | °C | |
6 | tx35 | Monthly count of days with maximum temperature above 35°C | 1 | |
7 | tx35ba | Monthly count of days with bias adjusted maximum temperature above 35°C | 1 | |
8 | tx40 | Monthly count of days with maximum temperature above 40°C | 1 | |
9 | tx40ba |
Monthly count of days with bias adjusted maximum temperature above 40°C | 1 | |||
10 | fd | Monthly count of frost days | 1 | |
11 | hd | Annual heating degree-days | °C day | |
12 | cd | Annual cooling degree-days | °C day | |
wet and dry | 13 | pr | Monthly mean of daily accumulated precipitation | mm |
14 | rx1day | Monthly maximum of 1-day accumulated precipitation | mm | |
15 | rx5day | Monthly maximum of 5-day accumulated precipitation | mm | |
16 | cdd | Annual consecutive dry days | days | |
17 | spi6 | Monthly Standardized Precipitation Index (SPI) for 6 months cumulation period | 1 | |
18* | spei6 | Monthly Standardized Precipitation-Evapotranspiration Index (SPEI) for 6 months cumulation period | 1 | |
19* | huss | Monthly near surface specific humidity | 1 | |
20* | evspsbl | Monthly evaporation including sublimation and transpiration | mm | |
21* | mrsos | Monthly soil moisture in upper soil portion | kg m-2 | |
22* | mrro | Monthly total runoff | mm | |
snow and ice | 23 | prsn | Monthly mean of daily accumulated snowfall precipitation | mm |
24 | siconc | Monthly mean of sea-ice area percentage | % | |
wind and radiation | 25 | sfcwind | Monthly mean of daily mean wind speed | m s-1 |
26* | clt | Monthly fraction of cloud cover | % | |
27* | rsds | Monthly surface solar radiation downwards | W m-2 | |
28* | rlds | Monthly surface thermal radiation downwards | W m-2 | |
ocean | 29 | sst | Monthly mean of sea surface temperature | °C |
circulation | 30* | psl | Monthly sea level pressure | Pa |
Table 3 shows the variables from Table 2 which are calculated for the different datasets (see Table 1). Note that bias adjusted indices are not computed for reanalysis/observational datasets, but only for climate projections (CMIP6 and CORDEX-EUR-11 by now). In the case of CMIP5, the only variables/indices considered are those included in the IPCC Atlas gridded monthly dataset due to the complex harmonization required for some of the new variables, which makes comparison across datasets difficult and potentially misleading.
Table 3. Description of variables/indices calculated for the different datasets. Data not available is indicated by 'N/A' and data not considered is indicated by 'N/C'. (*) Indicates variables not considered due to lack of data (only available for a reduced subset).
Code | CMIP6 | CMIP5 | CORDEX- | CORDEX- | ERA5 | ERA5- | E-OBS | ORAS5 |
t | X | X | X | X | X | X | X | N/A |
tn | X | X | X | X | X | X | X | N/A |
tx | X | X | X | X | X | X | X | N/A |
tnn | X | X | X | X | X | X | X | N/A |
txx | X | X | X | X | X | X | X | N/A |
tx35 | X | X | X | X | X | X | X | N/A |
tx35ba | X | N/C | X | N/C | N/A | N/A | N/A | N/A |
tx40 | X | X | X | X | X | X | X | N/A |
tx40ba | X | N/C | X | N/C | N/A | N/A | N/A | N/A |
fd | X | N/C | X | X | X | X | X | N/A |
hd | X | N/C | N/C | N/C | X | X | X | N/A |
cd | X | N/C | N/C | N/C | X | X | X | N/A |
pr | X | X | X | X | X | X | X | N/A |
rx1day | X | X | X | X | X | X | X | N/A |
rx5day | X | X | X | X | X | X | X | N/A |
cdd | X | N/C | X | X | X | X | X | N/A |
spi6 | X | N/C | X | X | X | X | X | N/A |
spei6 | X | N/C | N/C | N/C | X | N/C | X | N/A |
huss | X | N/C | X | X | N/A | N/A | N/A | N/A |
evspsbl | X | N/C | X | X | X | X | N/A | N/A |
mrsos | X | N/C | N/A | N/A | X | X | N/A | N/A |
mrro | X | N/A | X | N/A | X | X | N/A | N/A |
prsn | X | X | N/A | N/A | X | X | N/A | N/A |
siconc | X | N/C | N/A | N/A | X | N/A | N/A | X |
sfcwind | X | X | X | X | X | N/C | X | N/A |
clt | X | N/C | X | N/A (*) | X | N/A | N/A | N/A |
rsds | X | N/C | X | X | X | X | X | N/A |
rlds | X | N/C | X | X | X | X | N/A | N/A |
sst | X | N/C | N/A | N/A | X | N/A | N/A | X |
psl | X | N/C | X | N/A (*) | X | N/A | X | N/A |
2.3. Harmonization of variables across catalogues
...
Table 4: Harmonization work across the different catalogues carried out for the variables of the dataset. Color codes indicate the temporal resolutions used for the different datasets available in the CDS catalogues, which correspond to the colors indicated in Table 1 (green for hourly, blue for daily and purple for monthly). Unit harmonization (when needed) is indicated with "==>".
Code | Standard name and (units) | From raw input data to harmonized results (standard units) |
t | near_surface_air_temperature | CMIP6: near_surface_air_temperature (K) ==> -273.15 (°C) |
tn | daily_minimum_near_surface_air_temperature | CMIP6: daily_minimum_near_surface_air_temperature (K) ==> -273.15 (°C) |
tx | daily_maximum_near_surface_air_temperature | CMIP6: daily_maximum_near_surface_air_temperature (K) ==> -273.15 (°C) |
Pr | precipitation | CMIP6: precipitation (kg m-2 s-1) ==> *86400 (mm) |
huss | near_surface_specific_humidity | CMIP6: near_surface_specific_humidity (1) |
evspsbl | evaporation_inc_sublimation_and_transpiration | CMIP6: evaporation_including_sublimation_and_transpiration (kg m-2 s-1) ==>*86400 (mm) |
mrsos | mass_content_of_water_in_soil | CMIP6: moisture_in_upper_portion_of_soil_column (kg m-2) |
mrro | runoff_amount | CMIP6: total_runoff (kg m-2 s-1) ==> *86400 (kg m-2) |
prsn | snowfall_flux | CMIP6: snowfall_flux (kg m-2 s-1) ==> *86400 (mm) |
siconc | sea_ice_area_percentage_on_ocean_grid | CMIP6: sea_ice_area_percentage_on_ocean_grid (%) |
sfcwind | near_surface_wind_speed | CMIP6: near_surface_wind_speed (m s-1) |
clt | total_cloud_cover_percentage | CMIP6: total_cloud_cover_percentage (%) |
rsds | surface_downwelling_shortwave_radiation | CMIP6: surface_downwelling_shortwace_radiation (W m-2) |
rlds | surface_downwelling_longwave_radiation | CMIP6: surface_downwelling_longwave_radiation (W m-2) |
sst | sea_surface_temperature | CMIP6: sea_surface_temperature (K) ==> -273.15 (°C) |
psl | sea_level_pressure | CMIP6: sea_level_pressure (Pa) |
Info | ||
---|---|---|
| ||
7 Mixed Gregorian/Julian calendar as defined by UDUNITS. A deprecated alternative name for this calendar is gregorian. In this calendar, date/times after (and including) 1582-10-15 0:0:0 are in the Gregorian calendar, in which a year is a leap year if either ( i) it is divisible by 4 but not by 100 or (ii) it is divisible by 400. |
...
The climate projection datasets listed in Table 1 are based on multi-model ensembles including simulations from the historical and future scenarios (RCPs or SSPs). In order to compute climate change information, simulations from both historical and future scenarios are required. Therefore, in the Atlas dataset we consider only those models providing simulations for historical period and, at least, for one future scenario. Figures 1, 2 and 3 show the final ensembles computed for the different variables (in columns) for the CMIP6, CORDEX-EUR and CMIP5 datasets.
Figure 1: Final ensemble used in the "Copernicus Interactive Climate Atlas gridded monthly dataset" with the different models (in rows) for the different variables (in columns) for the CMIP6 dataset. Colors indicate the different scenarios available in the CDS.
Figure 2: Final ensemble used in the "Copernicus Interactive Climate Atlas gridded monthly dataset" with the different models (in rows) for the different variables (in columns) for the CORDEX-EUR-11 dataset. Colors indicate the different scenarios available in the CDS.
...
Table 5: CORDEX simulations for the different domains used to build the CORDEX-CORE dataset (note that EUR-11 simulations are upscaled to the global 0.22 grid). Empty cells indicate that the default versions for the two RCMs (RegCM4_7 and REMO2015) are used. Different versions and/or backup GCMs are indicated in the table.
AFR-22 | AUS-22 | CAM-22 | EAS-22 | EUR-11 | NAM-22 | SAM -22 | SEA-22 | WAS-22 | |
MOHC_HadGEM2-ES_REMO | |||||||||
MOHC_HadGEM2-ES_RegCM | v4-4_v0 | v4-6_v1 | v4-4-rc8 | MIROC5 | |||||
MPI-M_MPI-ESM-LR_REMO | r3i1p1 | ||||||||
MPI-M_MPI-ESM-MR_RegCM | v4-4_v0 | MPI-LR | MPI-LRv4-4-rc8 | v4-7_v0 | |||||
NCC_NorESM1-M_REMO | |||||||||
NCC_NorESM1-M_RegCM | GFDL | v4-4_v0 | v4-6_v1 | GFDLv4-4-rc8 | v4-7_v0 |
Figure 5: Final ensemble used in the "Copernicus Interactive Climate Atlas gridded monthly dataset" with the different models (in rows) for the different variables (in columns) for the CORDEX-CORE dataset. Colors indicate the different scenarios available in the CDS (see Table 5 for the specific configuration of the global mosaic).
...
Most of the variables and indices included in the datasets dataset are fully defined in Table 2 and require the application of simple temporal aggregations or threshold computations, and ; therefore, we implemented them ourselves they have been implemented in Python as part of the Atlas workflow to reduce the required computational cost. The unique indices that require more complex implementations are the Standardized exceptions are detailed below.
The Standardised Precipitation Index (SPI), the Standardized Precipitation-Evapotranspiration Index (SPEI) and the bias adjusted variables (tx35ba and tx40ba). The icclim module ({+}https://pypi.org/project/icclim/) was selected to compute the SPI and the SPEI indices and the ibicus module ({+}https://pypi.org/project/ibicus/) for the bias adjusted variables.
2.7. Bias adjustment
-6) is a monthly index that compares accumulated precipitation for 6 months with the long-term distribution (reference period: 1971-2010) for the same location and accumulation period, as the number of standard deviations from the median. This index has been calculated using using Python xclim software.
The Standardised Precipitation Evapotranspiration Index (SPEI-6) is a monthly index that compares accumulated precipitation minus potential evapotranspiration (PET, Thornthwaite definition) for 6 months with the long-term distribution (reference period: 1971-2010) for the same location and accumulation period, as the number of standard deviations from the median. This index has been calculated using using Python xclim software applying the Thornthwaite (TW) version for the calculation of PET (using tasmin and tasmax as input).
Cooling degree-days (CD) and Heating degree-days (HD) have been implemented in Python as part of the Atlas workflow, following the definition used in the IPCC AR6 Atlas (IPCC, 2021: Annex VI: Climatic Impact-driver and Extreme Indices [Gutiérrez J.M., R. Ranasinghe, A.C. Ruane, R. Vautard (eds.)]. Cambridge University Press, Cambridge, pp. 2205–2214, doi:10.1017/9781009157896.020).
2.7. Bias adjustment
In this first version of the dataset we build on recent work on an intercomparison In this first version of the dataset we build on recent work on an intercomparison of different bias adjustment methods11 which shows that a simple linear scaling method is a good choice for threshold-based temperature indices; the adjustment of the mean performed by these methods produces overall good results, as compared with more complex methods such as the ISIMIP3 method used in the IPCC Interactive Atlas. As an example, Figure 6 shows the results for the variable tasmax bias adjusted (tasmaxba) using both the ISIMIP3 and linear scaling methods for a particular model. From top to bottom, the panels display the results for linear scaling, ISIMIP3b and the difference between both methods.
All calculations are done using the ibicus package12 which implements a number of bias adjustment methods (including linear scaling and ISIMIP3 among others). In principle, climate indices including absolute thresholds are most sensitive to biases. This is more pronounced as the threshold values become more and more extreme. This is the case of the extreme maximum temperature indices as also reported by Iturbide et al. 2022 (On the need of bias adjustment for more plausible climate change projections of extreme heat, https://doi.org/10.1002/asl.1072). We use the simple linear scaling method for bias adjustment. In this version the two most extreme climate indices were bias adjusted for the CMIP6 and CORDEX-EUR-11 datasets: monthly count of days with maximum temperature above 35°Cand monthly count of days with maximum temperature above 40°C. The proposal for future versions of the Atlas dataset (if computing resources permit) is expanding the results including an additional bias adjustment method (the quantile mapping ISIMIP3 method), allowing users to explore this source of uncertainty. Note that this is aligned with the recommendation of the IPCC AR6 Report13.
Figure 6: Illustrative comparison between two commonly used bias adjustment methods of different complexity: ISIMIP3b (upper left) and linear scaling (upper right) for an illustrative model (the ACCESS-CM2 model). The panel in the bottom shows the difference between them.
...
Figure 7: Illustrations of the extent of the E-OBS (lefttop) and CORDEX-EUR-11 (rightbottom) European domains
In the longer term, we intend to introduce temporal subsetting too easing the selection of the time period and allowing to choose shorter time periods.
...
Figure 10: Example of box-plots for CMIP6 near surface air temperature (t) for different IPCC AR6 regions (in different panels) for the first version of the Copernicus Atlas (left, in blue) and the IPCC version (right, in green). The gray dashed lines join results from the same model in both datasets; small red squares indicate the ensemble means and the red lines join the medians of both ensembles(in different panels) for the first version of the Copernicus Atlas (left, in blue) and the IPCC version (right, in green). The gray dashed lines join results from the same model in both datasets; small red squares indicate the ensemble means and the red lines join the medians of both ensembles.
4. Versions of the document
The following table details the changes (and dates, in chronological order) that this document has undergone.
Date | Changes description |
---|---|
C3S Atlas V1.1 |
|
C3S Atlas V1.0 | First version of the document published together with the C3S Atlas launch (first version, v1.0), before the dataset is available in the CDS. This document serves as documentation for both the dataset and the C3S Atlas. |
Info |
---|
This document has been produced in the context of the Copernicus Climate Change Service (C3S). The activities leading to these results have been contracted by the European Centre for Medium-Range Weather Forecasts, operator of C3S on behalf of the European Union (Delegation Agreement signed on 11/11/2014 and Contribution Agreement signed on 22/07/2021). All information in this document is provided "as is" and no guarantee or warranty is given that the information is fit for any particular purpose. The users thereof use the information at their sole risk and liability. For the avoidance of all doubt , the European Commission and the European Centre for Medium - Range Weather Forecasts have no liability in respect of this document, which is merely representing the author's view |
...
Content by Label | ||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|