Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Easy Heading Macro
navigationExpandOptioncollapse-all-but-headings-1

Show If
spacePermission@self:edit

Major modification in progress for this article - please DO NOT UPDATE THIS PAGE or your edits will get overwritten when the new modified article will be published. Please edit the working copy of this page in the C3S Modified CKB articles or CAMS Modified CKB articles section.

Introduction

Here, we document the ERA5 dataset, which covers the period from January 1940 to the present and continues to be extended forward in near real time. For up to date information on ERA5, please consult the C3S Announcements on the Copernicus user forum.

...

The original ERA5 release contained data from 1979 onwards. The final ERA5 back extension for 1940-1978 has been produced and is available alongside the original/main release. The final ERA5 back extension differs from the preliminary version in several respects. 

An ERA5 back extension 1950-1978 (Preliminary version) was produced. Although in many other respects the quality is relatively good, this preliminary data does di suffer from excessively intense tropical cyclones. This dataset is available as a separate entry in the CDS catalogue (and in MARS), though soon it will be now deprecated.

Data update frequency

...

Expand
titleData organisation on the CDS disks

ERA5 data on the CDS disks can be downloaded either from the relevant CDS download page or, for larger data volumes in particular, using the CDS API. Subdivisions of the data are labelled using dataset and product_type.

Datasets reanalysis-era5-single-levels and reanalysis-era5-pressure-levels contain the following (sub-daily) product types:

  • reanalysis
  • ensemble_mean
  • ensemble_spread
  • ensemble_members

Datasets reanalysis-era5-single-levels-monthly-means and reanalysis-era5-pressure-levels-monthly-means contain the following (monthly) product types:

  • monthly_averaged_reanalysis
  • monthly_averaged_reanalysis_by_hour_of_day
  • monthly_averaged_ensemble_members
  • monthly_averaged_ensemble_members_by_hour_of_day

    particular, using the CDS API. Subdivisions of the data are labelled using dataset and product_type.

    Datasets reanalysis-era5-single-levels-preliminary-back-extension and reanalysis-era5-pressure-levels-preliminary-back-extension contain the following (sub-daily) product types:

    • reanalysis
    • ensemble_mean
    • ensemble_spread
    • ensemble_members

    Datasets reanalysis-era5-single-levels-monthly-means-preliminary-back-extension and reanalysis-era5-pressure-levels-monthly-means-preliminary-back-extension contain the following (monthly) product types:

    • monthly_averaged_reanalysis
    • monthly_averaged_reanalysis_by_hour_of_day
    • monthly_averaged_ensemble_members
    • monthly_averaged_ensemble_members_by_hour_of_day

    ...

    Expand
    titleData organisation in MARS

    ERA5 data in MARS can be accessed with the CDS API by specifying dataset whereas member state users can access data in MARS by specifying class and expver, according to the following table:

    ERA5 back extension 1950-1978

    CDS API access to MARS

    (specify the dataset)

    Member state access to MARS

    (specify class and expver)

    ERA5reanalysis-era5-completeclass=ea, expver=0001
    ERA5.1

    reanalysis-era5.1-complete

    class=ea, expver=0051
    ERA5Treanalysis-era5-complete1class=ea, expver=0005

    (Preliminary version)

    reanalysis-era5-complete-preliminary-back-extensionclass=ea, expver=0098

    1ERA5T data for a month is overwritten with the final ERA5 data about two months after the month in question.

    Subdivisions of the data are labelled using the keywords stream, type and levtype:

    Stream:

    • oper (HRES sub-daily)
    • wave (HRES sub-daily, for waves)
    • mnth (HRES synoptic monthly means, ie by hour of day)
    • moda (HRES monthly means of daily means)
    • wamo (HRES synoptic monthly means, ie by hour of day, for waves)
    • wamd (HRES monthly means of daily means, for waves)
    • enda (EDA sub-daily)
    • ewda (EDA sub-daily, for waves)
    • edmm (EDA synoptic monthly means, ie by hour of day)
    • edmo (EDA monthly means of daily means)
    • ewmm (EDA synoptic monthly means, ie by hour of day, for waves)
    • ewmo (EDA monthly means of daily means, for waves)

    Type:

    • an: analyses
    • fc: forecasts
    • em: ensemble mean
    • es: ensemble standard deviation

    Levtype:

    • sfc: surface or single level
    • pl: pressure levels
    • pt: potential temperature levels
    • pv: potential vorticity level
    • ml: model levels

    ...

    In order to define the surface geopotential in ERA5, the IFS uses surface elevation data interpolated from a combination of SRTM30 and other surface elevation datasets. For more details please see the IFS documentation, Cycle 41r2, Part IV. Physical processes, section 11.2.2 Surface elevation data at 30 arc seconds.

    Spatial reference systems and Earth model

    The IFS assumes that the underlying shape of the Earth is a perfect sphere, with of radius 6371.229 km, with the surface elevation specified relative to that sphere. The geodetic latitude/longitude of the surface elevation datasets are used as if they were the spherical latitude/longitude of the IFS.

    ERA5 data is referenced in the horizontal with respect to the WGS84 ellipse (which defines the major/minor axes) and in the vertical it is referenced to the EGM96 geoid over land but over ocean it is referenced to mean sea level, with the approximation that this is assumed to be coincident with the geoid. For more information on the relationship between mean sea level and the geoid, see for example Gregory et al. (2019).

    ...

    For data in GRIB1 format the earth model is a sphere with radius = 6367.47 km (note, this is inconsistent with what is actually used in the IFS),, as defined in the the WMO GRIB Edition 1 specifications, Table 7, GDS Octet 17.

    For data in GRIB2 format the earth model is a sphere with radius = 6371.2290 km2229 km (note, this is consistent with what is actually used in the IFS), as defined in the the WMO GRIB2 specifications, section 2.2.1, Code Table 3.2, Code figure 6.

    ...

    count

    name

    units

    Variable name in CDS

    shortName

    paramId

    an

    fc

    1

    Convective inhibition

    J kg**-1

    convective_inhibition

    cin

    228001


    x

    2

    Friction velocity

    m s**-1

    friction_velocity

    zust

    228003


    x

    3

    Lake mix-layer temperature

    K

    lake_mix_layer_temperature

    lmlt

    228008

    x

    x

    4

    Lake mix-layer depth

    m

    lake_mix_layer_depth

    lmld

    228009

    x

    x

    5

    Lake bottom temperature

    K

    lake_bottom_temperature

    lblt

    228010

    x

    x

    6

    Lake total layer temperature

    K

    lake_total_layer_temperature

    ltlt

    228011

    x

    x

    7

    Lake shape factor

    dimensionless

    lake_shape_factor

    lshf

    228012

    x

    x

    8

    Lake ice temperature

    K

    lake_ice_temperature

    lict

    228013

    x

    x

    9

    Lake ice depth

    m

    lake_ice_depth

    licd

    228014

    x

    x

    10

    UV visible albedo for direct radiation

    (0 - 1)

    uv_visible_albedo_for_direct_radiation

    aluvp

    15

    x

    x

    11

    Minimum vertical gradient of refractivity inside trapping layer

    m**-1

    minimum_vertical_gradient_of_refractivity_inside_trapping_layer

    dndzn

    228015


    x

    12

    UV visible albedo for diffuse radiation

    (0 - 1)

    uv_visible_albedo_for_diffuse_radiation

    aluvd

    16

    x

    x

    13

    Mean vertical gradient of refractivity inside trapping layer

    m**-1

    mean_vertical_gradient_of_refractivity_inside_trapping_layer

    dndza

    228016


    x

    14

    Near IR albedo for direct radiation

    (0 - 1)

    near_ir_albedo_for_direct_radiation

    alnip

    17

    x

    x

    15

    Duct base height

    m

    duct_base_height

    dctb

    228017


    x

    16

    Near IR albedo for diffuse radiation

    (0 - 1)

    near_ir_albedo_for_diffuse_radiation

    alnid

    18

    x

    x

    17

    Trapping layer base height

    m

    trapping_layer_base_height

    tplb

    228018


    x

    18

    Trapping layer top height

    m

    trapping_layer_top_height

    tplt

    228019


    x

    19

    Cloud base height

    m

    cloud_base_height

    cbh

    228023


    x

    20

    Zero degree level

    m

    zero_degree_level

    deg0l

    228024


    x

    21

    Instantaneous 10 metre wind gust

    m s**-1

    instantaneous_10m_wind_gust

    i10fg

    228029


    x

    22

    Sea ice area fraction

    (0 - 1)

    sea-ice_cover

    ci

    31

    x

    x

    23

    Snow albedo

    (0 - 1)

    snow_albedo

    asn

    32

    x

    x

    24

    Snow density

    kg m**-3

    snow_density

    rsn

    33

    x

    x

    25

    Sea surface temperature

    K

    sea_surface_temperature

    sst

    34

    x

    x

    26

    Ice temperature layer 1

    K

    ice_temperature_layer_1

    istl1

    35

    x

    x

    27

    Ice temperature layer 2

    K

    ice_temperature_layer_2

    istl2

    36

    x

    x

    28

    Ice temperature layer 3

    K

    ice_temperature_layer_3

    istl3

    37

    x

    x

    29

    Ice temperature layer 4

    K

    ice_temperature_layer_4

    istl4

    38

    x

    x

    30

    Volumetric soil water layer 11

    m**3 m**-3

    volumetric_soil_water_layer_1

    swvl1

    39

    x

    x

    31

    Volumetric soil water layer 21

    m**3 m**-3

    volumetric_soil_water_layer_2

    swvl2

    40

    x

    x

    32

    Volumetric soil water layer 31

    m**3 m**-3

    volumetric_soil_water_layer_3

    swvl3

    41

    x

    x

    33

    Volumetric soil water layer 41

    m**3 m**-3

    volumetric_soil_water_layer_4

    swvl4

    42

    x

    x

    34

    Convective available potential energy

    J kg**-1

    convective_available_potential_energy

    cape

    59

    x

    x

    35

    Leaf area index, low vegetation3

    m**2 m**-2

    leaf_area_index_low_vegetation

    lai_lv

    66

    x

    x

    36

    Leaf area index, high vegetation3

    m**2 m**-2

    leaf_area_index_high_vegetation

    lai_hv

    67

    x

    x

    37

    Neutral wind at 10 m u-component

    m s**-1

    10m_u-component_of_neutral_wind

    u10n

    228131

    x

    x

    38

    Neutral wind at 10 m v-component

    m s**-1

    10m_v-component_of_neutral_wind

    v10n

    228132

    x

    x

    39

    Surface pressure

    Pa

    surface_pressure

    sp

    134

    x

    x

    40

    Soil temperature level 11

    K

    soil_temperature_level_1

    stl1

    139

    x

    x

    41

    Snow depth

    m of water equivalent

    snow_depth

    sd

    141

    x

    x

    42

    Charnock

    ~

    charnock

    chnk

    148

    x

    x

    43

    Mean sea level pressure

    Pa

    mean_sea_level_pressure

    msl

    151

    x

    x

    44

    Boundary layer height

    m

    boundary_layer_height

    blh

    159

    x

    x

    45

    Total cloud cover

    (0 - 1)

    total_cloud_cover

    tcc

    164

    x

    x

    46

    10 metre U wind component

    m s**-1

    10m_u-_component_of_wind

    10u

    165

    x

    x

    47

    10 metre V wind component

    m s**-1

    10m_v-_component_of_wind

    10v

    166

    x

    x

    48

    2 metre temperature

    K

    2m_temperature

    2t

    167

    x

    x

    49

    2 metre dewpoint temperature

    K

    2m_dewpoint_temperature

    2d

    168

    x

    x

    50

    Soil temperature level 21

    K

    soil_temperature_level_2

    stl2

    170

    x

    x

    51

    Soil temperature level 31

    K

    soil_temperature_level_3

    stl3

    183

    x

    x

    52

    Low cloud cover

    (0 - 1)

    low_cloud_cover

    lcc

    186

    x

    x

    53

    Medium cloud cover

    (0 - 1)

    medium_cloud_cover

    mcc

    187

    x

    x

    54

    High cloud cover

    (0 - 1)

    high_cloud_cover

    hcc

    188

    x

    x

    55

    Skin reservoir content

    m of water equivalent

    skin_reservoir_content

    src

    198

    x

    x

    56

    Instantaneous large-scale surface precipitation fraction

    (0 - 1)

    instantaneous_large_scale_surface_precipitation_fraction

    ilspf

    228217


    x

    57

    Convective rain rate

    kg m**-2 s**-1

    convective_rain_rate

    crr

    228218


    x

    58

    Large scale rain rate

    kg m**-2 s**-1

    large_scale_rain_rate

    lsrr

    228219


    x

    59

    Convective snowfall rate water equivalent

    kg m**-2 s**-1

    convective_snowfall_rate_water_equivalent

    csfr

    228220


    x

    60

    Large scale snowfall rate water equivalent

    kg m**-2 s**-1

    large_scale_snowfall_rate_water_equivalent

    lssfr

    228221


    x

    61

    Instantaneous eastward turbulent surface stress

    N m**-2

    instantaneous_eastward_turbulent_surface_stress

    iews

    229

    x

    x

    62

    Instantaneous northward turbulent surface stress

    N m**-2

    instantaneous_northward_turbulent_surface_stress

    inss

    230

    x

    x

    63

    Instantaneous surface sensible heat flux

    W m**-2

    instantaneous_surface_sensible_heat_flux

    ishf

    231

    x

    x

    64

    Instantaneous moisture flux

    kg m**-2 s**-1

    instantaneous_moisture_flux

    ie

    232

    x

    x

    65

    Skin temperature

    K

    skin_temperature

    skt

    235

    x

    x

    66

    Soil temperature level 41

    K

    soil_temperature_level_4

    stl4

    236

    x

    x

    67

    Temperature of snow layer

    K

    temperature_of_snow_layer

    tsn

    238

    x

    x

    68

    Forecast albedo

    (0 - 1)

    forecast_albedo

    fal

    243

    x

    x

    69

    Forecast surface roughness

    m

    forecast_surface_roughness

    fsr

    244

    x

    x

    70

    Forecast logarithm of surface roughness for heat

    ~

    forecast_logarithm_of_surface_roughness_for_heat

    flsr

    245

    x

    x

    71

    100 metre U wind component

    m s**-1

    100m_u-component_of_wind

    100u

    228246

    x

    x

    72

    100 metre V wind component

    m s**-1

    100m_v-component_of_wind

    100v

    228247

    x

    x

    73

    Precipitation type2

    code table (4.201)

    precipitation_type

    ptype

    260015


    x

    74

    K index2

    K

    k_index

    kx

    260121


    x

    75

    Total totals index2

    K

    total_totals_index

    totalx

    260123


    x

    ...

    count

    name

    units

    Variable name in CDS

    shortName

    paramId

    an

    fc

    1

    Significant wave height of first swell partition

    m

    significant_wave_height_of_first_swell_partition

    swh1

    140121

    x

    x

    2

    Mean wave direction of first swell partition

    degrees

    mean_wave_direction_of_first_swell_partition

    mwd1

    140122

    x

    x

    3

    Mean wave period of first swell partition

    s

    mean_wave_period_of_first_swell_partition

    mwp1

    140123

    x

    x

    4

    Significant wave height of second swell partition

    m

    significant_wave_height_of_second_swell_partition

    swh2

    140124

    x

    x

    5

    Mean wave direction of second swell partition

    degrees

    mean_wave_period_of_second_swell_partition

    mwd2

    140125

    x

    x

    6

    Mean wave period of second swell partition

    s

    mean_wave_period_of_second_swell_partition

    mwp2

    140126

    x

    x

    7

    Significant wave height of third swell partition

    m

    significant_wave_height_of_third_swell_partition

    swh3

    140127

    x

    x

    8

    Mean wave direction of third swell partition

    degrees

    mean_wave_direction_of_third_swell_partition

    mwd3

    140128

    x

    x

    9

    Mean wave period of third swell partition

    s

    mean_wave_period_of_third_swell_partition

    mwp3

    140129

    x

    x

    10

    Wave Spectral Skewness

    dimensionless

    wave_spectral_skewness

    wss

    140207

    x

    x

    11

    Free convective velocity over the oceans

    m s**-1

    free_convective_velocity_over_the_oceans

    wstar

    140208

    x

    x

    12

    Air density over the oceans

    kg m**-3

    air_density_over_the_oceans

    rhoao

    140209

    x

    x

    13

    Normalized energy flux into waves

    dimensionless

    normalized_energy_flux_into_waves

    phiaw

    140211

    x

    x

    14

    Normalized energy flux into ocean

    dimensionless

    normalized_energy_flux_into_ocean

    phioc

    140212

    x

    x

    15

    Normalized stress into ocean

    dimensionless

    normalized_stress_into_ocean

    tauoc

    140214

    x

    x

    16

    U-component stokes drift

    m s**-1

    u_component_stokes_drift

    ust

    140215

    x

    x

    17

    V-component stokes drift

    m s**-1

    v_component_stokes_drift

    vst

    140216

    x

    x

    18

    Period corresponding to maximum individual wave height

    s

    period_corresponding_to_maximum_individual_wave_height

    tmax

    140217

    x

    x

    19

    Maximum individual wave height

    m

    maximum_individual_wave_height

    hmax

    140218

    x

    x

    20

    Model bathymetry

    m

    model_bathymetry

    wmb

    140219

    x

    x

    21

    Mean wave period based on first moment

    s

    mean_wave_period_based_on_first_moment

    mp1

    140220

    x

    x

    22

    Mean zero-crossing wave period

    s

    mean_zero_crossing_wave_period

    mp2

    140221

    x

    x

    23

    Wave spectral directional width

    dimensionlessRadians

    wave_spectral_directional_width

    wdw

    140222

    x

    x

    24

    Mean wave period based on first moment for wind waves

    s

    mean_wave_period_based_on_first_moment_for_wind_waves

    p1ww

    140223

    x

    x

    25

    Mean wave period based on second moment for wind waves

    s

    mean_wave_period_based_on_second_moment_for_wind_waves

    p2ww

    140224

    x

    x

    26

    Wave spectral directional width for wind waves

    dimensionlessRadians

    wave_spectral_directional_width_for_wind_waves

    dwww

    140225

    x

    x

    27

    Mean wave period based on first moment for swell

    s

    mean_wave_period_based_on_first_moment_for_swell

    p1ps

    140226

    x

    x

    28

    Mean wave period based on second moment for swell

    s

    mean_wave_period_based_on_second_moment_for_wind_waves

    p2ps

    140227

    x

    x

    29

    Wave spectral directional width for swell

    dimensionlessRadians

    wave_spectral_directional_width_for_swell

    dwps

    140228

    x

    x

    30

    Significant height of combined wind waves and swell

    m

    significant_height_of_combined_wind_waves_and_swell

    swh

    140229

    x

    x

    31

    Mean wave direction

    degrees

    mean_wave_direction

    mwd

    140230

    x

    x

    32

    Peak wave period

    s

    peak_wave_period

    pp1d

    140231

    x

    x

    33

    Mean wave period

    s

    mean_wave_period

    mwp

    140232

    x

    x

    34

    Coefficient of drag with waves

    dimensionless

    coefficient_of_drag_with_waves

    cdww

    140233

    x

    x

    35

    Significant height of wind waves

    m

    significant_height_of_wind_waves

    shww

    140234

    x

    x

    36

    Mean direction of wind waves

    degrees

    mean_direction_of_wind_waves

    mdww

    140235

    x

    x

    37

    Mean period of wind waves

    s

    mean_period_of_wind_waves

    mpww

    140236

    x

    x

    38

    Significant height of total swell

    m

    significant_height_of_total_swell

    shts

    140237

    x

    x

    39

    Mean direction of total swell

    degrees

    mean_direction_of_total_swell

    mdts

    140238

    x

    x

    40

    Mean period of total swell

    s

    mean_period_of_total_swell

    mpts

    140239

    x

    x

    41

    Mean square slope of waves

    dimensionless

    mean_square_slope_of_waves

    msqs

    140244

    x

    x

    42


    Expand
    title10 metre wind speed

    This 10m wind parameter is the wind speed that has been used by the wave model, which is coupled to the atmospheric model.

    For this reason:

    •  it is archived on the wave model's native grid, with the same land-sea mask as that model. Therefore, this parameter is not defined over land and wherever else the wave model is not defined, where it is encoded as missing data. Improper decoding of the missing value usually results in very large values being given for these land points.
    • the wave model resets all values below 2 m/s to 2m/s. The reason for this is that as the winds become weak, the long waves (swell) try to drive the wind from below but this is not modelled in the IFS, as it assumes that the wind profile should be logarithmic (+- stability correction). To account for this effect, the whole of the boundary layer scheme would need to be revised. A simple trick to avoid the problem is to boost the weak winds to 2m/s, which is outside the range where the waves can potentially drive the wind.
    • this parameter is actually the 10m neutral wind speed as determined from the atmospheric surface stress (see documentation on Ocean Wave model output parameters). 
    • If wave altimeter data were assimilated, the analysis of this parameter also contains wind speed updates that come directly out of the wave height updates.

    This parameter should not be used for looking at the quality of reanalysis surface wind - the u and v components of the 10m wind (atmospheric parameters 165 and 166) should be used instead.


    m s**-1

    ocean_surface_stress_equivalent_10m_neutral_wind_speed

    wind

    140245

    x

    x

    43

    10 metre wind direction

    degrees

    ocean_surface_stress_equivalent_10m_neutral_wind_direction

    dwi

    140249

    x

    x

    44

    Wave spectral kurtosis

    dimensionless

    wave_spectral_kurtosis

    wsk

    140252

    x

    x

    45

    Benjamin-Feir index

    dimensionless

    benjamin_feir_index

    bfi

    140253

    x

    x

    46

    Wave spectral peakedness

    dimensionless

    wave_spectral_peakedness

    wsp

    140254

    x

    x

    47

    Altimeter wave height

    m

    Not available from the CDS disks

    awh

    140246

    x


    48

    Altimeter corrected wave height

    m

    Not available from the CDS disks

    acwh

    140247

    x


    49

    Altimeter range relative correction

    ~

    Not available from the CDS disks

    arrc

    140248

    x


    50

    2D wave spectra (single)1

    m**2 s radian**-1

    Not available from the CDS disks

    2dfd

    140251

    x


    ...

    1. In general, we recommend that the hourly (analysed) "2 metre temperature" be used to construct the minimum and maximum over longer periods, such as a day, rather than using the forecast parameters "Maximum temperature at 2 metres since previous post-processing" and "Minimum temperature at 2 metres since previous post-processing".
    2. ERA5: compute pressure and geopotential on model levels, geopotential height and geometric height
    3. ERA5: How to calculate wind speed and wind direction from u and v components of the wind?
    4. Sea surface temperature and sea-ice cover (sea ice area fraction), see Table 2 above, are available at the usual times, eg hourly for the HRES, but their content is only updated once daily. However, for inland water bodies (lakes, reservoirs, rivers and coastal waters) the FLake model calculates the surface temperature (ie the lake mixed-layer temperature or lake ice temperature) and does include diurnal variations.
    5. Mean rates/fluxes and accumulations at step=0 have values of zero because the length of the processing period is zero.
    6. Convective Inhibition (CIN). A missing value is assigned to CIN for values of CIN > 1000 or where there is no cloud base. This can occur where convective available potential energy (CAPE) is low.

    7. Expand
      titleERA5: mixing CDS and MARS data

      In the ECMWF data archive (MARS), ERA5 data is archived on various native grids. For the CDS disks, ERA5 data have been interpolated and are stored on regular latitude/longitude grids. For more information, see ERA5: data documentation#Spatialgrid Spatialgrid.

      Storing the data on these different grids can cause incompatibilities, particularly when comparing native spherical harmonic, pressure level, MARS data with CDS disk data on a third, coarse grid.

      Native spherical harmonic, pressure level parameters are comprised of: Geopotential, Temperature, U component of wind, V component of wind, Vertical velocity, Vorticity, Divergence and Relative humidity. When these parameters are retrieved from MARS and a coarse output grid is specified, the default behaviour is that the spherical harmonics are truncated to prevent aliasing on the output grid. The coarser the output grid, the more severe the truncation. This truncation removes the higher wavenumbers, making the data smoother. However, the CDS disk data has been simply interpolated to the third grid, without smoothing.

      This incompatibility is particularly relevant when comparing ERA5.1 data (which are only available from MARS - see DataorganisationandhowtodownloadERA5 - and only for 2000-2006) with ERA5 data on the CDS disks.

      The simplest means of minimising such incompatibilities is to retrieve the MARS data on the same grid as that used to store the ERA5 CDS disk data.



    8. Expand
      titleERA5: Land-sea mask for wave variables

      The land-sea mask in ERA5 is an invariant field.

      This parameter is the proportion of land, as opposed to ocean or inland waters (lakes, reservoirs, rivers and coastal waters), in a grid box.

      This parameter has values ranging between zero and one and is dimensionless.

      In cycles of the ECMWF Integrated Forecasting System (IFS) from CY41R1 (introduced in May 2015) onwards, grid boxes where this parameter has a value above 0.5 can be comprised of a mixture of land and inland water but not ocean. Grid boxes with a value of 0.5 and below can only be comprised of a water surface. In the latter case, the lake cover is used to determine how much of the water surface is ocean or inland water. 

      The ERA5 land-sea mask provided is not suitable for direct use with wave parameters, as the time variability of the sea-ice cover needs to be taken into account and wave parameters are undefined for non-sea points.

      In order to produce a land-sea mask for use with wave parameters, users need to download the following ERA5 data (for the required period):

      1. the model bathymetry (Model bathymetry. Fig 1)
      2. the sea-ice cover (Sea ice area fraction, Fig 2)

      and combine these data to produce the land-sea mask (Fig 3). See attached pictures:

      Model bathymetry fieldSea ice cover fieldCombined mask

      Fig 1: Model bathymetry                                                 Fig 2: Sea-ice cover                                                          Fig 3: Combined mask


      Note

      Please note that sea-ice cover is only updated once daily.

      Please see the Toolbox workflow below to see a possible way to proceed. The results is a carousel of land-sea mask for each time step requested:

      Code Block
      titleToolbox workflow
      collapsetrue
      import cdstoolbox as ct
      
      @ct.application(title='Download data')
      @ct.output.download()
      @ct.output.carousel()
      
      def download_application():
          count = 0
          years=['1980']
          months = [
                  '01', #'02', '03',
              #    '04', '05', '06',
              #    '07', '08', '09',
              #    '10', '11', '12'
          ]
      # For hourly data hourly=True
      # For monthly data monthly=True
          hourly = True
          monthly = False
          for yr in years:
              for mn in months:
                  if hourly == True:
                      mb,si = get_hourly_data(yr, mn)
                  elif monthly == True:
                      mb,si = get_monthly_data(yr, mn)                
                  print(mb)
      # Check values are >= 0.0 in the model bathymetry mask
                  compare_ge_mb = ct.operator.ge(mb, 0.0)
                  print(si)
      # Check values are > 0.5 in the sea ice mask
                  compare_ge_si = ct.operator.gt(si, 0.500)
      
      # Invert model bathymetry mask
                  new =  ct.operator.add(compare_ge_mb, -1.0)
                  new1 =  ct.operator.mul(new, -1.0)
      # Add the Bathymetry Mask to the Sea Ice Mask
                  new_all = ct.operator.add(compare_ge_si,new1)
      # Reset scale to land=1, ocean=0
                  new_all_final = ct.operator.ge(new_all, 1.0)
                  print(new_all_final)
      
                  if count == 0:
                     combined_mask = new_all_final
                  else:
                     combined_mask = ct.cube.concat([combined_mask, new_all_final], dim = 'time')
                  count =  count + 1
      
          renamed_data = ct.cdm.rename(combined_mask, "wavemask")  
          new_data = ct.cdm.update_attributes(renamed_data, attrs={'long_name': 'Wave Land Sea Mask'})
          combined_mask = new_data
          print("combined_mask")  
          print(combined_mask)    
      
      # Plot mask for first timestep
      
          fig_list = ct.cdsplot.geoseries(combined_mask)
          return combined_mask, fig_list
      
      def get_monthly_data(y,m):
          m,s = ct.catalogue.retrieve(
              'reanalysis-era5-single-levels-monthly-means',
              {
                  'product_type': 'monthly_averaged_reanalysis',
                  'variable': [
                      'model_bathymetry', 'sea_ice_cover',
                  ],
                  'year': y,
                  'month': m,
                  'time': '00:00',
              }
          )
          return m, s
          
      def get_hourly_data(y,m):
          m,s = ct.catalogue.retrieve(
              'reanalysis-era5-single-levels',
              {
                  'product_type': 'reanalysis',
                  'variable': [
                      'model_bathymetry', 'sea_ice_cover',
                  ],
                  'year': y,
                  'month': m,
                  'day': [
                  '01', '02', '03',
                  '04', '05', '06',
                  '07', '08', '09',
                  '10', '11', '12',
                  '13', '14', '15',
                  '16', '17', '18',
                  '19', '20', '21',
                  '22', '23', '24',
                  '25', '26', '27',
                  '28', '29', '30',
                  '31',
                  ],
                  'time': [
                  '00:00', '01:00', '02:00',
                  '03:00', '04:00', '05:00',
                  '06:00', '07:00', '08:00',
                  '09:00', '10:00', '11:00',
                  '12:00', '13:00', '14:00',
                  '15:00', '16:00', '17:00',
                  '18:00', '19:00', '20:00',
                  '21:00', '22:00', '23:00',
                  ],
      
                  }
                  )
          return m, s
      
      




    9. Expand
      titleAltimeter wave parameters

      The following wave parameters are sparse observations, or quantities derived from the observations, that have been interpolated to the wave model grid and contain many missing values:

      • altimeter_wave_height (140246)
      • altimeter_corrected_wave_height (140247)
      • altimeter_range_relative_correction (140248)

      These parameters are not available from the CDS disks but can be retrieved from MARS using the CDS API. For further guidelines, please see: Altimeter wave height in the Climate Data Store (CDS)



    10. Expand
      titleComputation of near-surface humidity

      Near-surface humidity is not archived directly in ERA datasets, but the archive contains near-surface (2m from the surface) temperature (T) and dew point temperature (Td), and also surface pressure (sp), from which you can calculate specific and relative humidity at 2m.

      • Specific humidity can be calculated using equations 7.4 and 7.5 from Part IV, Physical processes section (Chapter 7, section 7.2.1b) in the documentation of the IFS for CY41R2. Use the 2m dew point temperature and surface pressure (which is approximately equal to the pressure at 2m) in these equations. The constants in 7.4 are to be found in Chapter 12 (of Part IV: Physical processes) and the parameters in 7.5 should be set for saturation over water because the dew point temperature is being used.
      • Relative humidity should be calculated from: RH = 100 * es(Td)/es(T)

       Relative humidity can be calculated with respect to saturation over water, ice or mixed phase by defining es(T) with respect to saturation over water, ice or mixed phase (water and ice). The usual practice is to define near-surface relative humidity with respect to saturation over water. Note that in ERA5, the relative humidity on pressure levels has been calculated with respect to saturation over mixed phase.



    11. Expand
      titleComputation of snow cover

      In the ECMWF model (IFS), snow is represented by an additional layer on top of the uppermost soil level. The whole grid box may not be covered in snow. The snow cover gives the fraction of the grid box that is covered in snow.

      For ERA5, the snow cover (SC) is computed using snow water equivalent (ie parameter SD (141.128)) as follows:

      Panel
      titleERA5 Snow cover formula

      snow_cover (SC) = min(1, (RW*SD/RSN) / 0.1 )

      where RW is density of water equal to 1000 and RSN is density of snow (parameter 33.128).


      ERA5 physical depth of snow where there is snow cover is equal to RW*SD/(RSN*SC).



    12. Expand
      title"Forecast albedo" is only for diffuse radiation

      The parameter "Forecast albedo" is only for diffuse radiation and assuming a fixed spectrum of downward short-wave radiation at the surface. The true broadband, all-sky, surface albedo can be calculated from accumulated parameters:

      (SSRD-SSR)/SSRD

      where SSRD is parameter 169.128 and SSR is 176.128. This true surface albedo cannot be calculated at night when SSRD is zero. For more information, see Radiation quantities in the ECMWF model and MARS.



    13. Expand
      titleActual and potential evapotranspiration

      Actual evapotranspiration in the ERA5 single levels datasets is called "Evaporation" (param ID 182) and is the sum of the following four evaporation components (which are not available separately in ERA5 but only for ERA5-Land):

      1. Evaporation from bare soil
      2. Evaporation from open water surfaces excluding oceans
      3. Evaporation from the top of canopy
      4. Evaporation from vegetation transpiration

      For the ERA5 single levels datasets, actual evapotranspiration can be downloaded from the C3S Climate Data Store (CDS) under the category heading "Evaporation and Runoff", in the "Download data" tab.

      For details about the computation of actual evapotranspiration, please see Chapter 8 of Part IV : Physical processes, of the IFS documentation:

      ERA5 IFS cycle 41r2

      The potential evapotranspiration in the ERA5 single levels CDS dataset is given by the parameter potential evaporation (pev)

      Pev data can be downloaded from the CDS under the category heading "Evaporation and Runoff", in the "Download data" tab for the ERA5 single levels datasets.

      Note

      The definitions of potential and reference evapotranspiration may vary according to the scientific application and can have the same definition in some cases. Users should therefore ensure that the definition of this parameter is suitable for their application.


      Note

      Please note that based on ERA5 atmospheric forcing, other independent (offline) methods such us "Priesley-Taylor1 (1972) , Schmidt2 (1915) or de Bruin3 (2000)" can also be used to estimate Potential evapotranspiration.

      1PRIESTLEY, C. H. B., & TAYLOR, R. J. (1972). On the Assessment of Surface Heat Flux and Evaporation Using Large-Scale Parameters, Monthly Weather Review, 100(2), 81-92. Retrieved Aug 27, 2021, from https://journals.ametsoc.org/view/journals/mwre/100/2/1520-0493_1972_100_0081_otaosh_2_3_co_2.xml 

      2Schmidt, W., 1915: Strahlung und Verdunstung an freien Wasserflächen; ein Beitrag zum Wärmehaushalt des Weltmeers und zum Wasserhaushalt der Erde (Radiation and evaporation over open water surfaces; a contribution to the heat budget of the world ocean and to the water budget of the earth). Ann. Hydro. Maritimen Meteor., 43, 111–124, 169–178.

      3de Bruin, H. A. R., , and Stricker J. N. M. , 2000: Evaporation of grass under non-restricted soil moisture conditions. Hydrol. Sci. J., 45, 391406, doi:10.1080/02626660009492337.




    14. Expand
      title"Evaporation" and "Instantaneous moisture flux"

      The "Instantaneous moisture flux" (units: kg m-2 s-1; paramId=232) incorporates the same processes as "Evaporation" (units: m of water equivalent; paramId=182), but the latter is accumulated over a particular time period (during the hour preceeding the validity date/time, in the ERA5 HRES), whereas the former is an instantaneous parameter. Note, the different units of these two parameters.

      For the atmosphere, these two parameters only involve water vapour. Cloud liquid does not sediment and the cloud ice sedimentation flux is included in the snowfall flux.

      Here are some further details about the processes in the "Instantaneous moisture flux" and "Evaporation":

      Surface characteristics

      Process from surface to atmosphere

      (defined to be negative)

      Process from atmosphere to surface

      (defined to be positive)

      Warm surfaceEvaporation from liquid water to water vapourDew deposition from water vapour
      Cold vegetation surfaceEvaporation from liquid water to water vapourDew deposition from water vapour
      Ice surfaceSublimation from ice to water vapourIce deposition from water vapour
      Snow surfaceSublimation from snow to water vapourSnow deposition from water vapour



    ...

    1. ERA5T: from 1 September to 13 December 2021, the final ERA5 product is different to ERA5T due to the correction of the assimilation of incorrect snow observations in central Asia. Although the differences are mostly limited to that region and mainly to surface parameters, in particular snow depth and soil moisture and to a lesser extent 2m temperature and 2m dewpoint temperature, all the resulting reanalysis fields can differ over the whole globe but should be within their range of uncertainty (which is estimated by the ensemble spread and which can be large for some parameters). On the CDS disks, the initial, ERA5T, fields have been overwritten (with the usual 2-3 month delay), i.e., for these months, access to the original CDS disk, ERA5T product is not possible after it has been overwritten. Potentially incorrect snow observations have been assimilated in ERA5 up to this time, when the effects became noticeable. The quality control of snow observations has been improved in ERA5 from September 2021 and from 15 November 2021 in ERA5T.
    2. ERA5 uncertainty: although small values of ensemble spread correctly mark more confident estimates than large values, numerical values are over confident. The spread does give an indication of the relative, random uncertainty in space and time.
    3. ERA5 suffers from an overly strong equatorial mesospheric jet, particularly in the transition seasons.
    4. From 2000 to 2006, ERA5 has a poor fit to radiosonde temperatures in the stratosphere, with a cold bias in the lower stratosphere. In addition, a warm bias higher up persists for much of the ERA5 period. The lower stratospheric cold bias was rectified in a re-run for the years 2000 to 2006, called ERA5.1, see "Resolved issues" below.
    5. Discontinuities in ERA5: The historic ERA5 data was produced by running several parallel experiments, each for a different period, which were then spliced together to create the final product. This can create discontinuities at the transition points.
    6. The analysed "2 metre temperature" can be larger than the forecast "Maximum temperature at 2 metres since previous post-processing".
    7. The analysed 10 metre wind speed (derived from the 10 metre wind components) can be larger than the forecast "10 metre wind gust since previous post-processing".
    8. ERA5 diurnal cycle for near surface winds: the hourly data reveals a mismatch in the analysed near surface wind speed between the end of one assimilation cycle and the beginning of the next (which occurs at 9:00 - 10:00 and 21:00 - 22:00 UTC). This problem mostly occurs in low latitude oceanic regions, though it can also be seen over Europe and the USA. We cannot rectify this problem in the analyses. The forecast near surface winds show much better agreement between the assimilation cycles, at least on average, so if this mismatch is problematic for a particular application, our advice would be to use the forecast winds. The forecast near surface winds are available from MARS, see the section, Data organisation and how to download ERA5.
    9. ERA5 diurnal cycle for near surface temperature and humidity: some locations do suffer from a mismatch in the analysed values between the end of one assimilation cycle and the beginning of the next, in a similar fashion to that for the near surface winds (see above), but this problem is thought not to be so widespread as that for the near surface winds. The forecast values for near surface temperature and humidity are usually smoother than the analyses, but the forecast low level temperatures suffer from a cold bias over most parts of the globe. The forecast near surface temperature and humidity are available from MARS, see the section Data organisation and how to download ERA5.
    10. ERA5: large 10m winds: up to a few times per year, the analysed low level winds, eg 10m winds, become very large in a particular location, which varies amongst a few apparently preferred locations. The largest values seen so far are about 300 ms-1.
    11. ERA5 rain bombs: up to a few times per year, the rainfall (precipitation) can become extremely large in small areas. This problem occurs mostly over Africa, in regions of high orography.
    12. Large values of CAPE: occasionally, the Convective available potential energy in ERA5 is unrealistically large.
    13. Ship tracks in the SST: prior to September 2007, in the period when HadISST2 was used, ship tracks can be visible in the SST.
    14. Prior to 2014, the SST was not used over the Great Lakes to nudge the lake model. Consequently, the 2 metre temperature has an annual cycle that is too strong, with temperatures being too cold in winter and too warm in summer.
    15. The Potential Evaporation field (pev, parameter Id 228251) is largely underestimated over deserts and high-forested areas. This is due to a bug in the code that does not allow transpiration to occur in the situation where there is no low vegetation.
    16. Wave parameters (Table 7 above) for the three swell partitions: these parameters have been calculated incorrectly. The problem is most evident in the swell partition parameters involving the mean wave period: Mean wave period of first swell partition, Mean wave period of second swell partition and Mean wave period of third swell partition, where the periods are far too long.
    17. Surface photosynthetically available radiation (PAR) is too low in the version (CY41R2) of the ECMWF Integrated Forecasting System (IFS) used to produce ERA5, so PAR and clear sky PAR have not been published in ERA5. There is a bug in the calculation of PAR, with it being taken from the wrong parts of the spectrum. The shortwave bands include 0.442-0.625 micron, 0.625-0.778 micron and 0.778-1.24 micron. PAR should be coded to be the sum of the radiation in the first of these bands and 0.42 of the second (to account for the fact that PAR is normally defined to stop at 0.7 microns). However, in CY41R2, PAR is in fact calculated from the sum of the second band plus 0.42 of the third. We will try to fix this in a future cycle.

    18. Expand
      titleThe instantaneous turbulent surface stress components (eastward and northward) and friction velocity tend to be too small

      The ERA5 analysed and forecast step=0, instantaneous surface stress components and surface roughness and the forecast step=0, friction velocity (friction velocity is not available from the analyses in ERA5) tend to suffer from values that are too low over the oceans.

      The analysis for such parameters is obtained by running the surface module to connect the surface with the model level analysed variables.

      However, at that stage, the surface aero-dynamical roughness length scale (z0) over the oceans is not initialised from its actual value but a constant value of 0.0001 is used instead.

      This initial value of z0 is needed to determine the initial value of u* and the surface stress based on solving for a simple logarithmic wind profile between the surface and the lowest model level. This initial u* is in turn used to determine an updated value of z0 based on the input Charnock parameter and then the value of the exchange coefficients needed to determine the output 10m winds (normal and neutral) and u* (see (3.91) to (3.94) with (3.26) in the IFS documentation). The surface stress is output as initialised.

      This initial value for z0 is generally too low ( by one order of magnitude or more):

      Over the oceans, for winds above few m/s, z0 is modelled using the Charnock relation:

      z0 ~ (alpha/g) u*2

      where alpha is the Charnock parameter, g is gravity, and u* is the friction velocity

      with typical values of

      alpha ~ 0.018

      g=9.81

      u*2 = Cd U102

      where Cd is the drag coefficient

      Cd ~ 0.008 + 0.0008 U10

      for U10=10m/s =>  z0 ~ 0.003


      As a consequence, the analysed instantaneous surface stress components will tend to be too low and even the updated value of z0 (surface roughness) will also tend to be too low.

      For forecast, instantaneous surface stress components, surface roughness and friction velocity, the same problem affects step 0. However, this problem will not affect the accumulated surface stress parameters (recall the accumulated parameters are produced by running short range forecasts), because the accumulation starts from the first time step (i.e. at time step 0 all accumulated variables are initialised to 0).

      This problem can easily be fixed, by using the initial value of Charnock that is available at the initial time.

      Note, in ERA5 the parameter for surface roughness is called "forecast surface roughness", even when it's analysed.


    19. ERA5 forecast parameters are missing for the validity times of 1st January 1940 from 00 UTC to 06 UTC (except for forecast step=0). This problem occurs because the first forecast in ERA5 was initiated from 1st January 1940 at 06 UTC.


    20. Maximum temperature at 2 metres since previous post-processing: in a small region over Peru, at 19 UTC, 2 August 2013, this forecast parameter exhibited erroneous values, which were greater than 50C. This occurrence is under investigation. Note, in general, we recommend that the hourly (analysed) "2 metre temperature" be used to construct the minimum and maximum over longer periods, such as a day.


    21. Expand
      titleFour reasons why hourly data might not be consistent with their monthly mean

      The ERA5 monthly means are calculated from the hourly (3 hourly for the EDA) data, on the native grid (including spherical harmonics) from the GRIB data, in each production "stream" or experiment. This can give rise to inconsistencies between the sub-daily data and their monthly mean, particularly in the CDS. In general, the inconsistencies will be small.

      • In the CDS, the ERA5 data (sub-daily and monthly mean) has been interpolated to a regular latitude/longitude grid. This interpolated sub-daily data will be slightly different to the native sub-daily data used in the production of the ERA5 monthly means.
      • The netCDF data available in the CDS has been packed, see What are NetCDF files and how can I read them, which states "unpacked_data_value = (packed_data_value * scale_factor) + add_offset" and "packed_data_value = nint((unpacked_data_value - add_offset) / scale_factor)". This netCDF packing will change the sub-daily values slightly, compared with the native sub-daily data used in the production of the ERA5 monthly means.
      • The GRIB data in the ERA5 monthly means (and sub-daily data) has been packed using a binning algorithm (which is different to the netCDF packing algorithm). Monthly means produced in other formats, such as netCDF, will differ from the ERA5 monthly means because of this packing.
      • Finally, there is a further reason why monthly mean values might be different to the mean of the sub-daily values, which even occurs in MARS. This cause only affects forecast parameters (the CDS provides analysed parameters unless the parameter is only available from the forecasts), such as the Total precipitation, and only occurs sporadically. In order to speed up production, ERA5 is produced in several parallel "streams" or experiments, which are then spliced together to produce the final product. Consider, the "stream" change at the beginning of 2015. The ERA5 forecast monthly means for January 2015 have been produced from the sub-daily data from that "stream", the first few hours of which (up until 06 UTC on 1st January 2015) come from the 18 UTC forecast on 31 December 2014. However, the sub-daily forecast data published in ERA5, is based on the date of the start of the forecast, so these first few hours of 2015 originate from the "stream" that produced December 2014. These two "streams" are different experiments, with different data values. The resulting inconsistencies might be larger than for the other three causes, above, depending on how consistent the two streams are.



    22. ERA5 sea-ice cover and 2 metre temperature: in the period 1979-1989, in a region just to the north of Greenland, the sea-ice cover outside of the melt season is too low and hence the 2 metre temperature is too high. For more information, see Section 3.5.4 of Low frequency variability and trends in surface air temperature and humidity from ERA5 and other datasets
    23. ERA5 sea-ice cover is missing in the Caspian Sea from late 2007 to 2013, inclusive.
    24. ERA5 sea-ice surface temperature (skin temperature) in the Arctic, during winter, can have a warm bias of 5K or more. This issue is most pronounced over thick snow-covered sea ice under cold clear-sky conditions, when the modelled conductive heat flux from the warm ocean underneath the ice and snow layer is too high. More information can be found in Batrak and Müller (2019) and Zampieri et al., (2023), the latter of which, also describes a method to improve on this bias.
    25. Altimeter wave height observations have not been available for ERA5 in the following periods (since coverage began in mid-1991): early February 2021 to mid-January 2022; mid-October 2023 onwards.
    26. ERA5 CDS: wind values are far too low on pressure levels at the poles in the CDS
    27. ERA5 back extension 1950-1978 (Preliminary version): tropical cyclones are too intense (Dataset deprecated in August 2023)
    28. ERA5 back extension 1950-1978 (Preliminary version): large bias in surface analysis over Australia prior to 1970 (Dataset deprecated in August 2023)
    29. ERA5 back extension 1950-1978 (Preliminary version): the deep soil moisture tends to be too dry (Dataset deprecated in August 2023)

    Resolved issues

    1. ERA5.1 is a re-run of ERA5, for the years 2000 to 2006 only, and was produced to improve upon the cold bias in the lower stratosphere seen in ERA5

      Expand
      titleMore information and details for downloading ERA5.1

      ERA5.1 is a re-run of ERA5 for the years 2000 to 2006 only. ERA5.1 was produced to improve upon the cold bias in the lower stratosphere exhibited by ERA5 during this period. Moreover, ERA5.1 analyses have a better representation of the following features:

      • upper stratospheric temperature
      • stratospheric humidity

      The lower and middle troposphere in ERA5.1 are similar to those in ERA5, as is the synoptic evolution in the extratropical stratosphere.

      For access to ERA5.1 data read Data organisation and how to download ERA5. The dataset is 'reanalysis-era5.1-complete' in the CDS API.


    2. ERA5.1 CDS: If you retrieved ERA5.1 using the CDS API anytime before 20/05/2020 08:00 UTC, for any stream other than oper (i.e. streams: wave, enda, edmo, ewmo, edmm, ewmm, ewda, moda, wamd, mnth, wamo), you will need to request the data again. Prior to this date, stream oper would be delivered regardless of which stream was requested.
    3. ERA5 CDS: incorrect values of U/V on pressure levels in the CDS
    4. ERA5 CDS: Data corruption

    ...