Contributors: Marco Cucchi (B-Open), Vivien Mavel Alessandro Amici (B-Open), Chiara Cagnazzo ), Graham Weedon (Met Office), Nicolas Bellouin (University of Reading), Stefan Lange (PIK), Hannes Müller Schmied (SBIK-F), Hans Hersbach (ECMWF), Marie-Louise John Carlo Buontempo (ECMWF), Niclas Rieger Chiara Cagnazzo (ECMWF)

Issued by: B-Open / Marco Cucchi

Issued Date: 12 15/1204/20252021

Ref:   C3S2_420a_Provision_of_WFDE5_extension    C3S_322_Lot1.4.1.3_CERRA_data_user_guide – version 1HPC ECMWF

Official reference number service contract:2024 2017/C3S2C3S_322_420a_B-Open Lot1_SMHI/SC2

Executive summary

The W5E5 dataset [1] was compiled to support the bias adjustment of climate input data for the impact assessments carried out in phase 3b of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3b). It consists of 14 present dataset, also known as WATCH Forcing Data methodology applied to ERA5 (WFDE5), is a meteorological forcing dataset for land surface and hydrological models. It consists of eleven variables (see Table 1) with a daily an hourly temporal resolution on a regular longitude-latitude half-degree grid, with global spatial coverage . W5E5 is a merged dataset, which combines WFDE5 [2] data over land with ERA5 reanalysis [3] data over the ocean. In addition, it makes use of the Global Precipitation Climatology Project (GPCP) data [3] for precipitation data rescaling.

W5E5 v3.0 is distributed through the C3S Climate Data Store as monthly files in netCDF format. It uses a full half-degree grid (720 × 360 grid boxes) covering the whole globe. See Table 2 for a summary of all the input datasets used for W5E5 v3.0, and Table 3 for a description of of the dataset's general attributes.

...

Variable name

...

Description

...

Units

...

Time coverage

...

mask

...

WFDE5-ERA5 mask

...

1

...

-

...

orog

...

Surface Altitude

...

m

...

-

...

sfcWind

...

Near-Surface Wind Speed

...

m s-1

...

1979-01-01 to 2024-12-31

...

tas

...

Near-Surface Air Temperature

...

K

...

1979-01-01 to 2024-12-31

...

tasmax

...

Daily Maximum Near-Surface Air Temperature

...

K

...

1979-01-01 to 2024-12-31

...

tasmin

...

Daily Minimum Near-Surface Air Temperature

...

K

...

1979-01-01 to 2024-12-31

...

ps

...

Surface Air Pressure

...

Pa

...

1979-01-01 to 2024-12-31

...

psl

...

Sea Level Pressure

...

Pa

...

1979-01-01 to 2024-12-31

...

huss

...

Near-Surface Specific Humidity

...

kg kg-1

...

1979-01-01 to 2024-12-31

...

hurs

...

Near-Surface Relative Humidity

...

kg kg-1

...

1979-01-01 to 2024-12-31

...

rlds

...

Surface Downwelling Longwave Radiation

...

W m-2

...

1979-01-01 to 2024-12-31

...

rsds

...

Surface Downwelling Shortwave Radiation

...

W m-2

...

1979-01-01 to 2024-12-31

...

pr

...

Precipitation

...

kg m-2 s-1

...

1979-01-01 to 2020-12-31

...

prsn

...

Snowfall Flux

...

kg m-2 s-1

...

1979-01-01 to 2020-12-31

and values defined only for land and lake points. The dataset was derived applying sequential elevation and monthly bias correction methods described in [1] and [2], and briefly outlined in Section 2, to half-degree aggregated ERA5 reanalysis products [3].

The monthly observational datasets used for bias correction are CRU TS [4] for all variables and the GPCC full data monthly product [5] [6] [7] for precipitation variables only, with specific versions detailed in Table 2. As a result, two different datasets for each precipitation variable (Rainf and Snowf) are obtained: one corrected using only the CRU TS dataset and one corrected using both the CRU TS and GPCC datasets. Other input data used for the generation of WFDE5 dataset are: a) catch corrections for rain- and snow-precipitation gauges and b) shortwave downward radiation corrections due to atmospheric aerosol loading effects. See Table 3 for a summary of all the input datasets used.

WFDE5 dataset is distributed through the C3S Climate Data Store as monthly files in netCDF format. It uses a full half-degree grid (720 × 360 grid boxes) with sea and large lakes grid-points flagged as missing data, comprising a total of 92889 land points (Antarctica included). General dataset attributes are described in Table 4.

A detailed description of the dataset can be found in [8].

Algorithms Description

 The algorithms described here correspond to those of [1] and [2] developed for the WATCH forcing data. The WFDE5 dataset is computed using a Python software which take the input variables shown in Table 3 and automatically apply the following key processing steps:

1. ERA5 reanalysis aggregation from the quarter- to the CRU half-degree grid;
2. Sequential “elevation correction” of Tair, PSurf, Qair, and LWdown to account for differences in surface heights between quarter- and half-degree grids and to ensure consistency between the different corrected variables. Tair is bias-adjusted to observed monthly averaged near-surface air temperature and observed monthly averaged diurnal temperature range before the processing of PSurf;
3. Adjustment of SWdown, Rainf and Snowf at the monthly scale via CRU TS and, for precipitation variables only, GPCC observational datasets.

Aggregation to the half-degree grid is applied on sea-level adjusted values of the different ERA5 variables, which are obtained as part of the elevation-correction procedure.

Elevation-correction procedures are physically-based sequential adjustments applied to account for differences in surface heights between input and output grid and ensure consistency among the different variables. The variables to which elevation-correction is applied are i) Tair, ii) PSurf, iii) Qair and iv) SWdown, which have to be processed in the cited order: each corrected variable is indeed needed to process the following one.

Further adjustments based on the monthly observational datasets listed in Table 2 and Table 3 are applied to variables Tair, SWdown, Rainf and Snowf.

Two corrections are applied to Tair values, both driven by differences between pre-processed ERA5 values and CRU monthly observations: first, a monthly scale correction based on the differences between temperature values is applied; then, an hourly scale correction taking into account differences in diurnal temperature ranges (DTR) leads to the final values. It is worth noticing that, before applying the aforementioned processing steps, CRU DTR values are corrected for known anomalies over a few specific grid boxes

. Also, to avoid the occurrence of anomalously high Tair values noticed in WFDE5v1.0 and v1.1, the hourly scale correction algorithm based on DTR has been update in WFDE5v2.0 with respect to the one detailed in [1] and [2], and it is briefly described in Appendix A.

A further modification to the original algorithm detailed in [1] and [2] has been applied to generate WFDE5v2.1, in order to fix the issue affecting Qair (and, with a minor impact, LWdown) values described in the "Known issues" section. The implemented change is described in Appendix B.

SWdown values are first adjusted to be consistent with CRU observed cloud-cover fractions: this is done using local linear relationship between anomalies in monthly short-wave radiation and cloud cover in aggregated ERA5 data along with CRU cloud cover anomalies to reconstruct the associated short-wave radiation anomalies. Then, corrections due to changes in direct and indirect effects of atmospheric aerosol loading on surface short-wave radiation fluxes are applied.

Finally, the following corrections are applied to Rainf and Snowf values: a) adjustment of precipitation fluxes values to ensure matching of the monthly number of wet days between interpolated ERA5 and CRU data; b) correction to monthly precipitation totals to match CRU/GPCC values; c) correction to rainfall and snowfall fluxes to account for catch corrections of precipitation gauges [9]; d) precipitation phase (snowfall/rainfall) correction when adjusted Tair lays beyond thresholds derived from ERA5 temperature extremes at fixed precipitation phase within each grid box and calendar month.

For all meteorological variables, values over Antarctica are obtained simply via aggregation and elevation-correction, given the absence of observations over that area.

ASurf values, corresponding to WFDE5 dataset's grid-points altitude, are obtained applying ERA5 sea- land and lake cover mask to CRU grid-points heights: in this way, only those which are identified as land or lake points by both ERA5 and CRU datasets are retained, resulting in a total of 92889 grid- points.

Data Description

File naming convention

File names adhere to the following convention:

<var>_WFDE5_<reference_dataset>_<YYYYMM>_v<version>.nc, 

where

• <var>: variable name, as in Table 1

• <reference_dataset>: one between CRU (all variables) and CRU+GPCC (Rainf and Snowf only)
• <YYYYMM>: year and month
• <version>: WFDE5 dataset version

As an example, possible file names are:

• Tair_WFDE5_CRU_201801_v1.0.nc
• Rainf_WFDE5_CRU_201801_v1.0.nc
• Rain_WFDE5_CRU+GPCC_201801_v1.0.nc

File content

Actual file content and dimensions size varies from month to month and depending on the variable, but the general file structure is constant and analogous to the following example:

netcdf Tair_WFDE5_CRU_201801_v1.0 { 
	dimensions:
		lon = 720 ;
		lat = 360 ;
		time = 744 ;
	variables:
		double lon(lon) ;
			lon:_FillValue = NaN ;
			lon:standard_name = "longitude" ;
			lon:units = "degrees_east" ; 
			lon:axis = "X" ;
			lon:long_name = "Longitude" ;
		double lat(lat) ;
			lat:_FillValue = NaN ;
			lat:standard_name = "latitude" ;
			lat:units = "degrees_north" ;
			lat:axis = "Y" ;
			lat:long_name = "Latitude" ;
		double time(time) ;
			time:_FillValue = NaN ;
			time:standard_name = "time" ;
			time:long_name = "Time" ; time:axis = "T" ;
			time:units = "hours since 1900-01-01" ;
			time:calendar = "gregorian" ;
		float Tair(time, lat, lon) ;
			Tair:_FillValue = 1.e+20f ;
			Tair:long_name = "Near-Surface Air Temperature" ;
			Tair:standard_name = "air_temperature" ;
			Tair:units = "K" ;
// global attributes:
:title = "WATCH Forcing Data methodology applied to ERA5 data" ;
:institution = "Copernicus Climate Change Service" ;
:contact = "http://copernicus-support.ecmwf.int" ;
:summary = "ERA5 data regridded to half degree regular lat-lon; Genuine land points from CRU grid and ERA5 land-sea mask only; Tair elevation & bias-corrected using CRU TS4.03 mean monthly temperature and mean diurnal temperature range" ;
:Conventions = "CF-1.7" ;
:comment = "Methodology implementation for ERA5 and dataset production by B-Open Solutions for the Copernicus Climate Change Service in the context of contract C3S_25c" ;
:reference = "Weedon et al. 2014 Water Resources Res. 50, 7505-7514, doi:10.1002/2014WR015638; Harris et al. 2014 Int. J. Climatol. 34, 623-642, doi:10.1002/joc.3711; Cucchi et al. 2020 Earth Syst. Sci.
Data Discuss., doi:10.5194/essd-2020-28, 2020" ;
:licence = "The dataset is distributed under the Licence to Use Copernicus Products. The corrections applied are based upon CRU TS4.03, distributed under the Open Database License (ODbL)" ;
}

Known issues

1. 
   

   Expand
   title Error for Near Surface specific humidity values
   Error: An error has been reported for the variable near-surface specific humidity (kg kg-1), in regions/times where/when T < 0. For this variable, given ERA5 near-surface specific and relative humidity are not available in the CDS, the following procedure is used: ERA5 vapour pressure and saturation vapour pressure at the surface, VP and SVP, respectively, are computed following Buck (1981) [10] VP and SVP are used to compute ERA5 relative humidity (RHE5) at surface as Mathinline RH_{E5} = \frac{VP}{SVP}*100 final bias corrected specific humidity is calculated by following the methodology described in [1] See Cucchi et al (2020) [8] for further details. The problems are due to the calculation of RHE5, because in WFDE5 this is done with different coefficients for ice and water phases, whereas the WMO standard is to compute RH with respect to water at all temperatures. Details: The WFDE5 near-surface specific humidity (QairW5) is calculated in four steps. Step 1 Computation of ERA5 near-surface relative humidity (RHE5), by using: ERA5 near-surface air temperature (TairE5), near-surface dewpoint temperature (TdewE5) and surface pressure (PSurfE5) using the formula: Mathdisplay RH = 100 \ast VP_{E5}(Tdew{E5}, PSurf{E5}) / SVP_{E5}(Tair{E5}, PSurf{E5}) where both VPE5 (ERA5 near-surface vapour pressure) and SVPE5 (ERA5 near-surface saturated vapour pressure) have been computed using eq 4a (with sets of coefficients ew4 and ei3 respectively for T(C)>=0 and T(C)<0) and corrections factor fw4 and fi4 all defined in Buck (1981) [10] for both TairE5 and TdewE5, i.e. with different coefficients for both TairE5 >= 0 / TairE5 < 0 and TdewE5 >= 0 / TdewE5 < 0. (See eq 12 and 13 from Weedon et al (2010) [1]. Step 2 Computation of WFDE5 near-surface saturated water content (SWCW5) from WFDE5 near-surface air temperature (TairW5) and surface pressure (PSurfW5), using the formula (as equation 13 of Weedon et al (2010) [1]): Mathdisplay SWC_{W5} = \frac{SVP_{W5}}{PSurf_{W5} - SVP_{W5} * 0.37802} * 0.62198 where SVPW5 (WFDE5 near-surface saturated vapour pressure) has been computed using eq 4a (with sets of coefficients ew4 and ei3 respectively for T(C)>=0 and T(C)<0) and corrections factor fw4 and fi4 all defined in Buck (1981) [10]. Step 3 Computation of WFDE5 near-surface specific humidity (QairW5) using WFDE5 near-surface saturated water content (SWCW5) and implied near-surface relative humidity (i.e. ERA5 near-surface relative humidity RHE5), using the formula (as Eq 16 of Weedon et al (2010) [1]): Mathdisplay Qair{W5} = \frac{SWC_{W5} \ast RH_{E5}}{100} Why is it calculated in this way and which are the differences w.r.t. WFDEI? The bias correction procedure for WFDE5 near-surface specific humidity (QairW5) and WFDEI near-surface specific humidity (QairWI) is the same and follows Weedon et al (2010) [1]. However, since ERA5 near-surface specific humidity and relative humidity are not available in the Copernicus Climate Change Service (C3S) Climate Data Store (CDS) (from which, for data provenance and algorithmic reasons all non-observational data for the generation of the WFDE5 dataset are taken), they had to be derived from other available variables (see above). In the case of WFDEI (based on ERA-Interim), instead, QairWI was estimated by taking directly as input data the ERA-Interim lowest model level specific humidity (at model level=137, close to 10m). Implications: Users have reported an unexpected behaviour when estimating the relative humidity against temperature. As an example, Figure 1 below shows, starting from the top-left corner and proceeding clockwise, observed and ERA5-, WFDE5- and WFDEI-derived values of near-surface relative humidity against surface air temperature on a grid-point (68.75N, 133.75W) close to a station located in Canada, for 2013: in the figure, WFDE5-derived values show an anomalous behaviour by never reaching ice saturation. Image Added  Figure 1: Unexpected behaviour when estimating relative humidity against temperature close to a station located in Canada at 68.75N, 133.75W.  Solution: The issue here described has been solved in WFDE5v2.1 (see Appendix B for a description of the required update to the dataset generation algorithm). In Figure 2 a comparison of near-surface relative humidity (RH) against near-surface air temperature (Tair) for WFDEI, WFDE5v2.0 and the new WFDE5v2.1 at 68.75N, 133.75W, including data for the whole 2013 year, is shown. Image Added  Figure 2: Comparison of RH vs Tair for WFDEI (blue), WFDE5v2.0 (red) and WFDE5v2.1 (yellow) at 68.75N, 133.75W.  Figure 2 shows how, in the updated WFDE5v2.1, relative humidity values at low air temperatures realign with the expected behaviour, i.e. reach the ice saturation curve as it happens for WFDEI. A more general analysis of the impact of the new algorithm has been performed. Here the key findings: As expected, the impact on Qair values (changes with respect to the previous version) is higher on cold regions/seasons (see Figure 3 below); The impact on LWdown values computation, which uses bias-corrected Qair values as input data, is negligible. Image Added Image Added Figure 3: relative difference of Qair in WFDE5v2.1 and WFDE5v2.0 for two timesteps, in January and June 2013.

   
   

Licence and Citation

The described dataset is distributed under the Licence to Use Copernicus Products. The corrections applied are based upon CRU TS 4.03, distributed under the Open Database License (ODbL), CRU TS 4.04 and 4.09, distributed under the Open Government Licence (OGL), and GPCCv2018/v2020/v2025, distributed under the Creative Commons Attribution 4.0 International Licence (CC BY 4.0).

If publishing using this dataset please cite [8], where a more detailed description of the dataset can be found.

Appendix A – WFDE5v2.0: DTR-based correction algorithm update

As mentioned in section 2, the Tair correction algorithm based on CRU DTR has been updated in WFDE5v2.0 with respect to the one used for WFDE5v1.0 and WFDE5v1.1. While details on the latter can be found in [1] and [2], a brief description of the new algorithm is reported below:

1. For each month, interpolated ERA5 hourly near-surface air temperatures are used to compute the ERA5 interpolated monthly average diurnal temperature range (DTR_E5,mm);

2. DTR correction factors A and R are computed in the following way: 

   Mathdisplay
   A = DTR_{CRU} - DTR_{E5,mm}

   
   

   Mathdisplay
   R = \frac{DTR_{CRU}}{DTR_{E5,mm}}

   
   

3. ERA5 interpolated diurnal temperature range values are corrected (DTR_E5,corr), depending on the value of R at each grid-point, in the following way: 

   
   

   Mathdisplay
   DTR_{E5,corr} = \begin{cases} DTR_{E5} * R & A \ge 0 \\ DTR_{E5} + A & A < 0 \\ \end{cases}

   
   

4. A multiplicative correction factor F is computed as: 

   Mathdisplay
   F = \frac{DTR_{E5,corr}}{DTR_{E5}}

   
   

5. After applying the correction for monthly mean biases, deviations of corrected hourly ERA5 near-surface air temperature values (Tair_E5,corr) from their respective daily means (Tair_E5,corr,dm) are multiplied by F , and results are added to ERA5 near-surface air temperature daily means (Tair_E5,dm) to obtain final corrected hourly WFDE5 near-surface air temperature values (Tair_W5): 

   Mathdisplay
   Tair_{W5} = Tair_{E5,corr,dm} + (Tair_{E5,corr} - Tair_{E5,corr,dm}) \ast F

   
   

Appendix B – WFDE5v2.1: Qair derivation algorithm update

As mentioned in section 2, the Qair derivation algorithm has been updated in WFDE5v2.1 with respect to the one used for previous versions of the dataset. While details on the latter can be found in [1] and [2], a brief description of the new algorithm is reported below.

1) Computing ERA5 implied relative humidity

The first step needed for the derivation of bias corrected specific humidity at surface (Qair) is the derivation of ERA5 relative humidity at surface (RH_E5), which is not available from the Climate Data Store (CDS) catalogue. To do so, the following steps are carried out.

Step 1

Compute ERA5 saturated vapour pressure (SVP_E5) and vapour pressure (VP_E5) using ERA5 surface pressure (PSurf_E5) and, respectively, ERA5 surface air temperature (Tair_E5) and dewpoint temperature (Tdew_E5), with the following formulas [10]:

\begin{array}{ll}
SVP_{E5}(Tair_{E5}, PSurf_{E5}) =
EA_{w/i}*\exp((EB_{w/i}-\frac{Tair_{E5}}{ED_{w/i}})*\frac{Tair_{E5}}{Tair_{E5}+EC_{w/i}})*F_{w/i}(Tair_{E5}, PSurf_{E5}) \\
VP_{E5}(Tdew_{E5}, PSurf_{E5}) =
EA_{w}*\exp((EB_{w}-\frac{Tdew_{E5}}{ED_{w}})*\frac{Tdew_{E5}}{Tdew_{E5}+EC_{w}})*F_{w}(Tdew_{E5}, PSurf_{E5})
\end{array}

with enhancement factors F_w and F_i computed as:

\begin{array}{ll}
F_{w}(T=[Tair_{E5}, Tdew_{E5}], PSurf_{E5}) = PSurf_{E5} * (FC_{w} * T^{2} + FB_w) + FA_{w} \\ 
F_{i}(T=[Tair_{E5}, Tdew_{E5}], PSurf_{E5}) = PSurf_{E5} * (FC_{i} * T^{2} + FB_i) + FA_{i} 
\end{array}

where:

• Tair_E5: ERA5 surface air temperature (C)
• Tdew_E5: ERA5 surface dewpoint temperature (C)
• PSurf_E5: ERA5 surface pressure (hPa)
• EA_w = 6.1121
• EB_w = 18.729
• EC_w = 257.87
• ED_w = 227.3
• EA_i = 6.1115
• EB_i = 23.036
• EC_i = 279.82
• ED_i = 333.7
• FA_w = 1.00072
• FB_w = 3.2 * 10^-6
• FC_w = 5.9 * 10^-10
• FA_i = 1.00022
• FB_i = 3.83 * 10^-6
• FC_i = 6.4 * 10^-10
• For formulas involving Tair_E5 only, constants with subscript “w” or “i” are used when respectively Tair_E5 >= 0 and Tair_E5 < 0

NOTE: These corresponds to formulas 4a (with sets of coefficients ew4 and ei3 respectively for Tair_E5 >= 0 and Tair_E5 < 0) with enhancement factors fw4 and fi4 from Buck (1981) [10]

Step 2

Compute ERA5 saturated water content (SWC_E5) and water content (WC_E5) using ERA5 surface pressure (PSurf_E5) and, respectively, SVP_E5 and VP_E5 with the following formulas:

\begin{array}{ll}
SWC_{E5} = \frac{SVP_{E5}}{PSurf_{E5} - SVP_{E5} * 0.37802} * 0.62198 \\ 
WC_{E5} = \frac{VP_{E5}}{PSurf_{E5} - VP_{E5} * 0.37802} * 0.62198
\end{array}

Step 3

Compute ERA5 implied relative humidity (RH_E5) as:

RH_{E5} = \frac{WC_{E5}}{SWC_{E5}} * 100

2) Computing WFDE5 saturated vapour pressure

Here, WFDE5 saturated vapour pressure (SVP_W5) is derived from WFDE5 surface pressure (PSurf_W5) and WFDE5 surface air temperature (Tair_W5) using the same formulas applied in Step 1 above, i.e.:

SVP_{W5}(Tair_{W5}, PSurf_{W5}) =
EA_{w/i}*\exp((EB_{w/i}-\frac{Tair_{W5}}{ED_{w/i}})*\frac{Tair_{W5}}{Tair_{W5}+EC_{w/i}})*F_{w/i}(Tair_{W5}, PSurf_{W5})

with:

\begin{array}{ll}
F_{w}(Tair_{W5}, PSurf_{W5}) = PSurf_{W5} * (FC_{w} * Tair_{W5}^{2} + FB_w) + FA_{w} \\ 
F_{i}(Tair_{W5}, PSurf_{W5}) = PSurf_{W5} * (FC_{i} * Tair_{W5}^{2} + FB_i) + FA_{i} 
\end{array}

where:

• Tair_W5: WFDE5 surface air temperature (C)
• PSurf_W5: WFDE5 surface pressure (hPa)
• Coefficient values are the same reported in Step 1 above

3) Computing WFDE5 saturated water content

Here, WFDE5 saturated water content (SWC_W5) is derived from WFDE5 surface pressure (PSurf_W5) and SVP_W5,Tair using the following formula:

SWC_{W5} = \frac{SVP_{W5}}{PSurf_{W5} - SVP_{W5} * 0.37802} * 0.62198

4) Computing WFDE5 specific humidity at surface

Finally, WFDE5 specific humidity at surface (Qair) is derived from SWC_W5 and RH_E5 using the formula:

Qair = \frac{SWC_{W5} * RH_{E5}}{100}

Anchor
references references

References

...

Dataset

...

Summary

...

Variables used

...

ERA5

...

ECMWF reanalysis product

...

• Geopotential
• 10 m u-component of wind
• 10 m v-component of wind
• 2 m temperature
• Surface pressure
• 2 m dewpoint temperature
• Mean surface downward longwave radiation flux
• Mean surface downward shortwave radiation flux
• Total precipitation
• Snowfall

...

WFDE5v3.0

...

Near surface meteorological variables from 1979 to 2024 derived from bias-corrected reanalysis

...

• Grid-points altitude
• Near-surface wind speed
• Near-surface air temperature
• Surface air pressure
• Near-surface specific humidity
• Surface downwelling longwave radiation
• Surface downwelling shortwave radiation
• Rainfall flux (CRU+GPCC corrected)
• Snowfall flux (CRU+GPCC corrected)

...

GPCPv2.3

...

Global Precipitation Climatology Project data

...

• Precipitation

...

Dataset attribute

...

Details

...

Horizontal coverage

...

Global

...

Horizontal resolution

...

0.5° x 0.5°

...

Vertical coverage

...

Surface

...

Temporal coverage

...

Depends on variable (see Table 2)

...

Temporal resolution

...

Daily

...

File format

...

netCDF

...

Data type

...

Grid

...

Versions

...

v3.0

Algorithms Description

The algorithms used for the generation of W5E5 v3.0 are identical to the ones used for the generation of previous versions v1.0 and v2.0, and are described at [1]. For the generation of W5E5 v3.0, a Python implementation of these algorithm has been developed and used, and is available for download in the "Documentation" tab of the dataset's CDS catalogue entry.

Data Description

File naming convention

File names adhere to the following convention:

<var>_W5E5_<YYYYMM>_v<version>.nc, 

where

• <var>: variable name, as in Table 1

• <YYYYMM>: year and month (where applicable)
• <version>: W5E5 dataset version

As an example, possible file names are:

• tas_W5E5_201801_v3.0.nc
• pr_W5E5_201801_v3.0.nc
• mask_W5E5_v3.0.nc

File content

Actual file content and dimensions size varies from month to month and depending on the variable, but the general file structure is constant and analogous to the following example:

netcdf tas_W5E5_201801_v3.0 {
dimensions:
	lat = 360 ;
	lon = 720 ;
	time = 31 ;
variables:
	double lat(lat) ;
		lat:_FillValue = NaN ;
		lat:long_name = "Latitude" ;
		lat:standard_name = "latitude" ;
		lat:units = "degrees_north" ;
		lat:axis = "Y" ;
	double lon(lon) ;
		lon:_FillValue = NaN ;
		lon:long_name = "Longitude" ;
		lon:standard_name = "longitude" ;
		lon:units = "degrees_east" ;
		lon:axis = "X" ;
	int64 time(time) ;
		time:long_name = "Time" ;
		time:standard_name = "time" ;
		time:axis = "T" ;
		time:units = "hours since 1900-01-01" ;
		time:calendar = "proleptic_gregorian" ;
	float tas(time, lat, lon) ;
		tas:_FillValue = 1.e+20f ;
		tas:units = "K" ;
		tas:long_name = "Near-Surface Air Temperature" ;
		tas:standard_name = "air_temperature" ;

// global attributes:
		:title = "WFDE5 over land merged with ERA5 over the ocean (W5E5) version 3.0" ;
		:institution = "Potsdam Institute for Climate Impact Research (PIK) & European Center for Medium Range Weather Forecast (ECMWF) under Copernicus Climate Change Service" ;
		:project = "Inter-Sectoral Impact Model Intercomparison Project phase 3 (ISIMIP3)" ;
		:contact = "ISIMIP cross-sectoral science team <info@isimip.org> <https://www.isimip.org>; ECMWF Support Portal <http://copernicus-support.ecmwf.int>" ;
		:summary = "WFDE5 (with GPCC precipitation correction) over land merged with ERA5 (with GPCP precipitation correction) over the ocean for ISIMIP3" ;
		:references = "Cucchi et al. (2020) <https://doi.org/10.5194/essd-12-2097-2020> and Lange et al. (2021) <https://doi.org/10.48364/ISIMIP.342217> for W5E5; Weedon et al. (2014) <https://dx.doi.org/10.1002/2014WR015638> and Cucchi et al. (2020) <https://doi.org/10.5194/essd-12-2097-2020> for WFDE5; Hersbach et al. (2020) <https://dx.doi.org/10.1002/qj.3803> and Copernicus Climate Change Service Climate Data Store (2017) <https://cds.climate.copernicus.eu/#!/home> for ERA5" ;
		:version = "3.0" ;
		:Conventions = "CF-1.7" ;
}

Licence and Citation (???)

The described dataset is distributed under the Licence to Use Copernicus Products. The corrections applied are based upon CRU TS 4.03, distributed under the Open Database License (ODbL), CRU TS 4.04, distributed under the Open Government Licence (OGL), and GPCCv2018/v2020, distributed under the Creative Commons Attribution 4.0 International Licence (CC BY 4.0).

If publishing using this dataset please cite [2], where a more detailed description of the dataset can be found.

...

Anchorreferences1references1[1]    Lange, Stefan (2019): WFDE5 over land merged with ERA5 over the ocean (W5E5). V. 1.0. GFZ Data Services. https://doi.org/10.5880/pik.2019.023 Anchorreferences2references2[2]    M. Cucchi, G. P. Weedon, A. Amici, N. Bellouin, S. Lange, H. M. Schmied, H. Hersbach, C. Buontempo, "WFDE5: bias-adjusted ERA5 reanalysis data for impact studies", Earth System Science Data, https://essd.copernicus.org/articles/12/2097/2020, 2020. Anchor. Lange, H. Müller Schmied, H. Hersbach and C. Buontempo, "WFDE5: bias-adjusted ERA5 reanalysis data for impact studies", Earth Syst. Sci. Data, vol. 12, no. 3, pp. 2097–2120, doi: 10.5194/essd-12-2097-2020, 2020.

...