...

The Copernicus Climate Change Service (C3S) of the European Commission aims to provide authoritative information about the past, present and future climate in Europe and the rest of the World. It is in this context that the pan-European (Fig. 1) Copernicus European Regional ReAnalysis (CERRA) has been developed, under the contract C3S_322_Lot1. This The land surface reanalysis (CERRA-Land) covers the period 1984-2021 and has a horizontal resolution of 5.5km. The datasets dataset can be used in support of adaptation action and policy development as well as contribute to climate monitoring and research, but also provide valuable information to climate services.

...

Figure1: CERRA-Land domain. Orography [m] is presented in color. The system used the Lambert Conformal Conic projection.

 The The need for precipitation and surface variables at an ever-increasing spatial and temporal resolution is a recurrent demand. These variables allow, among other things, to address water resource management issues and to carry out climate change impact studies. Regional surface reanalyses are a way to reconstruct these variables for past periods covering several decades using state-of-the-art models.

...

 CERRA-Land is the result of a unique standalone integration of the SURFEX V8.1 land surface model (Le Moigne et al. 2020) driven by meteorological forcing data from the CERRA atmospheric reanalysis and an offline analysis of daily accumulated surface precipitation (Soci et al., 2016). The accumulated precipitation is created by an optimal interpolation between an initial estimate (first guess) based on CERRA predicted precipitation and daily rain gauge data. Temperature and relative humidity forcing data at 2 metres are from the CERRA surface analysis and the short- and long-wave downwelling radiation, wind speed at 10 metres, and surface pressure are from CERRA forecast data. A brief description of the CERRA-Land system is available (C3S_D322_Lot1.1.1.12_202110_documentation_CERRA-Land, add link)in Documentation tab.

Details about the data fields

The outputs were archived at hourly time step and the precipitation analysis was archived at daily time step. The ouputs data are available in the same way as the CERRA dataset as if analysis were done every 3 hours ( at 00, 03, 06, 09, 12, 15, 18 and 21 UTC). The date of the analysis is related to the two metere metre temperature and humidity forcing fields. The forecast fields are available at three time step ( +1, 2  and 3 hours).

...

The surface pressure, 10-metre wind speed and direction, two metre temperature and relative humidity analysis were not archived into the CERRA-Land dataset, but are available through the CDS from the CERRA on single levels dataset (add link).

The available CERRA-Land variables in the CDS are eitherinstantaneousinstantaneous, accumulated or static. The natural, urban and inland water fraction will be available in netCDF format only. This is specified for each of the variables listed in the tables below.

Table1Table 1: Overview of Fields variables calculated as a mean value of a grib box.

Table2: Overview of Fields available for the natural land fraction.

...

Name

...

ShortName

...

Unit

...

GRIB2 CODE

...

Analysis

0, 3 …, 21

(or daily)

...

Forecast range

1, 2, 3

...

Height

Snow depth

(Instantaneous)

...

sde

...

m

...

3066

...

no

...

yes

...

surface

Snow depth water equivalent

(Instantaneous)

...

sd

...

kg.m^-2

228141

...

no

...

yes

...

surface

Snow density

(Instantaneous)

...

rsn

...

kg.m^-3

...

33

...

no

...

yes

...

surface

Fraction of snow cover

(Instantaneous)

...

fscov

...

[0-1]

...

260289

...

no

...

yes

...

surface

...

Snow albedo

(Instantaneous)

...

asn

...

[0-1]

...

228032

...

no

...

yes

...

surface

Temperature of snow layer

(Instantaneous)

...

tsn

...

K

...

238

...

no

...

yes

...

surface

...

Snow melt

(accumulated)

...

snom

...

kg m^-2

...

3099

...

no

...

yes

...

surface

...

Percolation

(accumulated)

...

perc

...

kg m^-2

...

260430

...

no

...

yes

...

surface

...

Surface runoff

(accumulated)

...

sro

...

kg m^-2

...

174008

...

no

...

yes

...

surface

Albedo

The albedo [0-100%] is the total reflectance of downward solar radiation at the surface over the grid box. The albedo is the ratio of one hour time-integrated surface solar radiation upward by the one hour time-integrated surface solar radiation downward. Multiplying the albedo with the one hour accumulated downward solar radiation gives the one hour accumulated upward solar radiation.

Evaporation

Evaporation is the amount of water that has evaporated from the earth’s surface from the initial time of the forecast to the forecast time step. It is given as a mean for the grid area between the three surface types in the grid - inland water, natural land and urban. Hence, evaporation is represented by negative values and positive values represent condensation. By model convention downward fluxes are positive. It is an accumulated variable.

Total precipitation

The total precipitation is the amount of precipitation falling at the surface during the last 24-hours. It includes all kind of precipitation forms as convective precipitation, large scale precipitation, liquid and solid precipitation. The total precipitation is available only for the analyses at 06 UTC. It is an accumulated field from the previous day at 06 UTC to the present day at 06 UTC. The date in the metadata GRIB2 file represents the end of the accumulation period.

Skin temperature 

It represents the average air temperature at the surface of each grid box. The skin temperature is an average of temperatures given by the three surface types in the grid - inland water, natural land and urban. It is an instantaneous variable.

Surface latent heat flux

The surface latent heat flux is the accumulated exchange of latent heat (due to phase transitions - evaporation, condensation) with the surface (ground and water) through turbulent diffusion from the initial time of the forecast to the forecast time step. It is given as a mean for the grid area between the three surface types in the grid - inland water, natural land and urban. By model convention downward fluxes are positive. It is an accumulated variable.

Surface sensible heat flux

The surface sensible heat flux is the accumulated exchange of heat (no phase transition) with the surface (ground and water) through turbulent diffusion from the initial time of the forecast to the forecast time step. It is given as a mean for the grid area between the three surface types in the grid - inland water, natural land and urban. By model convention downward fluxes are positive. It is an accumulated variable.

Surface net solar radiation

The surface net solar radiation is the accumulated solar short-wave radiation that is absorbed at the surface from the initial time of the forecast to the forecast time step. It is calculated as the difference between the downward solar energy and the upward solar energy at the surface. By model convention downward fluxes are positive. It is an accumulated variable.

Surface solar radiation downwards

The surface solar radiation downward is the accumulated total (direct and diffuse) solar short-wave radiation reaching the surface from the initial time of the forecast to the forecast time step. By model convention downward fluxes are positive. It is an accumulated variable.

Surface net thermal radiation

The net thermal radiation at the surface is accumulated from the initial time of the forecast to the forecast time step. It is calculated as the difference between the thermal radiation downwards and the upward thermal radiation at the surface. By model convention downward fluxes are positive. It is an accumulated variable.

Surface thermal radiation downwards

The surface thermal radiation downward is the amount of thermal (long-wave) radiation reaching the surface accumulated from the initial time of the forecast to the forecast time step. By model convention downward fluxes are positive. It is an accumulated variable.

Soil heat flux 

The soil heat flux is the energy receive by the soil to heat it per unit of surface and time. The Soil heat flux is positive when the soil receives energy (warms) and negative when the soil loses energy (cools). It is an instantaneous variable.

Surface roughness 

The surface roughness describes the aerodynamic roughness length. It is given as a mean for the grid area between the three surface types in the grid (inland water, natural land and urban) and has missing values over the ocean. The roughness length of the surface is the height above the surface at which the wind profile is assumed to become zero. Each grid point has one value representing the mean over the grid point. The effective surface roughness is depending on the orographic component (constant part), the snow depth, the evolution of the Leaf Area Index and the fraction of vegetation, which is different for each month. It is an instantaneous variable.

Snow variables

The snow variables are related to natural land only. The fraction of snow cover represents the fraction of natural land which has snow on the ground. Snow on urban or lake fraction is not available.

Table 2: Overview of variables available for the natural land fraction.

 Fraction of snow cover

It represents the fraction (0-1) occupied by snow for nature fraction only. It is an instantaneous variable.

 Snow depth

Snow thickness on the ground. It is an instantaneous variable.

 Snow depth water equivalent

The mass of liquid water obtained from melting the snow per unit area. This is equivalent to the depth of the liquid water in units of mm. It is an instantaneous variable.

 Snow density

The mean snow density is calculated as the ratio of snow depth water equivalent by the snow depth. It is an instantaneous variable.

Snow albedo

It is defined as the fraction of solar (short-wave) radiation reflected by the snow, across the solar spectrum, for both direct and diffuse radiation. Values vary between 0 and 1. It is an instantaneous variable.

 Temperature of snow layer

It is the mean temperature of the 12 snow layers. It is an instantaneous variable.

 Snow melt

Melting of snow. This variable is accumulated from the beginning of the forecast time to the end of the forecast step. It is an accumulated variable.

Percolation

The mass per unit area of water that drains below the deepest soil level in the model. The percolation (or drainage) is accumulated from the initial time of the forecast to the forecast time step. This variable is calculated for the natural land, including soil, vegetation and snow. (not for urban and water bodies fraction). It is an accumulated variable calculated for the natural land, including soil, vegetation and snow whereas.

Surface runoff

The mass per unit area of water at the surface when saturation occurs.  It is an accumulated variable calculated for the natural land, including soil, vegetation and snow.

Soil temperature

The soil temperature is provided for each soil layer. The SURFEXmodel has 14 soil layers. The vertical discretization (bottom depth of each layer in metres) is as follows: 0.01, 0.04, 0.1, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 3, 5, 8, and 12 m. Soil temperature is an instantaneous variable.

Volumetric soil moisture 

The volume concentration of liquid and ice water. The vertical discretization (bottom depth of each layer in metres) is as follows: 0.01, 0.04, 0.1, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 3, 5, 8, and 12 m. The volumetric soil moisture is available for each soil layer. It is an instantaneous variable.

Liquid volumetric soil moisture (non-frozen)

It is the volume concentration of liquid water only. The vertical discretization (bottom depth of each layer in metres) is as follows: 0.01, 0.04, 0.1, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 3, 5, 8, and 12 m. A liquid volumetric soil moisture is available for each soil layer. It is an instantaneous variable.

Table 3: Overview of variables available for the inland water fraction.

Lake bottom temperature

Temperature of water at the bottom of inland water bodies (lakes). The model keeps lake depth and surface area (or fractional cover) constant in time. It is an instantaneous variable.

Lake ice depth

The thickness of ice on inland water bodies (lakes). A single ice layer is represented. This parameter is the thickness of that ice layer. It is an instantaneous variable.

Lake ice temperature

The temperature of the uppermost surface of ice on inland water bodies (lakes). A single ice layer is represented. It is an instantaneous variable.

Lake mix-layer depth

The thickness of the upper most layer of an inland water body (lake) that is well mixed and has a near constant temperature with depth (uniform distribution of temperature). The Flake model  represents inland water bodies with two layers in the vertical, the mixed layer above and the thermocline below. Thermoclines upper boundary is located at the mixed layer bottom, and the lower boundary at the lake bottom. Mixing within the mixed layer can occur when the density of the surface (and near-surface) water is greater than that of the water below. Mixing can also occur through the action of wind on the surface of the lake. It is an instantaneous variable.

Lake mix-layer temperature

The temperature of the upper most layer of inland water bodies (lakes) that is well mixed. The Flake model represents inland water bodies with two layers in the vertical, the mixed layer above and the thermocline below. Thermoclines upper boundary is located at the mixed layer bottom, and the lower boundary at the lake bottom. Mixing within the mixed layer can occur when the density of the surface (and near-surface) water is greater than that of the water below. Mixing can also occur through the action of wind on the surface of the lake. It is an instantaneous variable.

Lake shape factor

This parameter describes the way that temperature changes with depth in the thermocline layer of inland water bodies (lakes). It is used to calculate the lake bottom temperature and other lake-related parameters. The Flake model represents inland water bodies with two layers in the vertical, the mixed layer above and the thermocline below where temperature changes with depth. It is an instantaneous variable.

Lake total layer temperature

The mean temperature of total water column in inland water bodies (lakes). The Flake model represents inland water bodies with two layers in the vertical, the mixed layer above and the thermocline below where temperature changes with depth. This parameter is the mean over the two layers. It is an instantaneous variable.

Static variables

Static variable do not change depending on the model initial time or the forecast length (in other words they are time-independent). These include the land-sea mask, that is the fraction of land in a given model grid box of 5.5 x 5.5 km² in units of %, and the orography in units of m. For each model grid box in CERRA-Land, 3 tile fractions are defined in units of %: (1) the fraction of inland water (lakes and rivers), (2) the fraction of urban areas, and (3) the fraction of nature, i.e. land areas that are not inland water or urban. There is no data available for the fraction of sea. The fraction data will be available directly from the climate data store website in netCDF format.

Table4: Static variable. The variables marked with * and labelled with TBD ("To Be Determined") do not have yet short

names in the WMO GRIB2 code definitions and have not been uploaded to the CDS in the first batch of released data.

Lake depth

Depth of inland water. It is defined for positive fractions only. It is a static field.

Volumetric wilting point - if in the model is a function of soil type and volumetric transpiration stress-onset shouldn't this be mentioned?

When the soil moisture reaches the wilting point, the vegetation is not able to extract the soil water. The soil moisture content is too low to be absorbed by the vegetation.

Volumetric transpiration stress-onset (soil moisture)

The soil moisture is the water content of a soil after gravitational drainage. When the water content of the soil reaches this value, the water cannot drain any more by gravity.

Land-sea mask

The land-sea mask is a field that contains, for every grid point, the proportion of land (including inland water) in the grid box. It is the sum of the three fractions - natural land, urban and inland water. The variable is dimensionless and takes values between 0 (sea) and 1 (land). It is a static variable.

Orography

The height above sea level of the land surface. This variable does not change with snow cover. It is a static variable.

Inland water tile fraction

This variable represents the fraction of water (e.g lakes, rivers) in the grid-box. It takes values between 0 and 1.

Urban tile fraction

This variable represents the urban (e.g town) in the grid-box. It takes values between 0 and 1.

Nature tile fraction

This variable represents the fraction of natural land areas that are neither inland water nor urban in the grid-box. It takes values between 0 and 1.

Soil variables

The prognostic variables of soil temperature and soil moisture are represented in the model by a diffusive approach. Such a method proposes a discretization of the soil into 14 layers, resulting in a total depth of 12 m, with a fine description of the subsurface layers to capture the diurnal cycle. The vertical discretization (bottom depth of each layer in metres) is as follows: 0.01, 0.04, 0.1, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 3, 5, 8, and 12 m. Heat transfer is resolved over the total depth, while moisture transfer is resolved only over the depth of the roots, which depends on the type of vegetation and its geographical location.

Accumulated surface fluxes

All energy fluxes at the surface are accumulated variables from the initial time of the forecast to the forecast hour in question with the unit J/m². They are considered positive downward to the surface. Average hourly energy fluxes in W/m² can be computed by subtracting two successive hourly accumulated variables and dividing by 3600 s. The albedo in units of % is defined by the ratio of average hourly upward radiation upward by the average hourly solar radiation downward. The net solar radiation is the difference between the solar radiation downward and the solar radiation upward. The solar radiation upward can be calculated by subtracting the net solar radiation from the downward solar radiation.

The net thermal radiation is the difference between the downward thermal radiation and the upward thermal radiation. The upward thermal radiation can be calculated by subtracting the net thermal radiation from the downward thermal radiation. The accumulated surface sensible heat flux is the conductive energy from the atmosphere to the surface. If this is going from the surface to the atmosphere it has negative values. The accumulated latent heat flux is the sum of all latent energy fluxes that are due to the phase transitions of water. Here condensation causes a positive latent heat flux to the surface and evaporation causes a negative heat flux from the surface.

References

Albergel C., Dorigo W., Reichle R. H., Balsamo G., de Rosnay P., Muñoz-Sabater J., Isaksen L., de Jeu R., and Wagner

Soil temperature

(Instantaneous)

...

sot

...

K

...

260360

...

no

...

yes

...

Soil

(14 layers)

Volumetric soil moisture

(Instantaneous)

...

wsw

...

m³m^-3

...

260199

...

no

...

yes

...

Soil

(14 layers)

Liquid volumetric soil   moisture

(Instantaneous)

...

liqvsm

...

 m³m^-3

260210

...

no

...

yes

...

Soil

(14 layers)

Table3: Overview of Fields available for the inland water fraction.

...

Name

...

ShortName

...

Unit

...

GRIB2 CODE

...

Analysis
0, 3 …, 21
(or daily)

...

Forecast range
1, 2, 3

...

Height

...

Lake bottom temperature
(Instantaneous)

...

lblt

...

K

...

228010

...

yes

...

no

...

surface

...

Lake ice depth
(Instantaneous)

...

licd

...

m

...

228014

...

no

...

yes

...

surface

...

Lake ice surface temperature
(Instantaneous)

...

lict

...

K

...

228013

...

no

...

yes

...

surface

...

Lake mix-layer depth
(Instantaneous)

...

lmld

...

m

...

228009

...

no

...

yes

...

surface

...

Lake mix-layer temperature
(Instantaneous)

...

lmlt

...

K

...

228008

...

no

...

yes

...

surface

...

Lake shape factor
(Instantaneous)

...

lshf

...

dimensionless

...

228012

...

no

...

yes

...

surface

...

Lake total layer temperature
(Instantaneous)

...

ltlt

...

K

...

228011

...

no

...

yes

...

surface

Table4: Static fields. Parameters labelled with TBD ("To Be Determined") do not yet have short name and GRIB2 code definitions

...

Name

...

ShortName

...

Unit

...

GRIB2 CODE

...

Analysis

0, 3 …, 21

(or daily)

...

Forecast range

1, 2, 3

...

Height

...

Lake total depth

...

dl

...

m

...

228007

...

no

...

no

...

surface

Volumetric wilting point

...

vwiltm

...

m³m^-3

...

260200

...

no

...

no

...

Soil

(14 layers)

Volumetric field capacity

...

voltso

...

m³m^-3

...

260211

...

no

...

no

...

Soil

(14 layers)

Land-sea mask

...

lsm

...

[0-1]

...

172

...

no

...

no

...

surface

Orography

...

orog

...

m

...

228002

...

no

...

no

...

surface

...

Inland water tile fraction

...

TBD

...

[0-1]

...

netCDF

...

no

...

no

...

surface

...

Urban tile fraction

...

TBD

...

[0-1]

...

netCDF

...

no

...

no

...

surface

...

Nature tile fraction

...

TBD

...

[0-1]

...

netCDF

...

no

...

no

...

surface

Soil parameters

The prognostic variables of soil temperature and soil moisture are represented in the model by a diffusive approach. Such a method proposes a discretization of the soil into 14 layers, resulting in a total depth of 12 m, with a fine description of the subsurface layers to capture the diurnal cycle. The vertical discretization (bottom depth of each layer in metres) is as follows: 0.01, 0.04, 0.1, 0.2, 0.4, 0.6, 0.8, 1, 1.5, 2, 3, 5, 8, and 12 m. Heat transfer is resolved over the total depth, while moisture transfer is resolved only over the depth of the roots, which depends on the type of vegetation and its geographical location.

Static fields

Static fields are variables that do not change depending on the model initial time or the forecast length (in other words they are time-independent). These include the land-sea mask, that is the fraction of land in a given model grid box of 5.5 x 5.5 km² in units of %, and the orography in units of m. For each model grid box in CERRA-Land 3 tile fractions are defined in units of %:
the fraction of inland water (lakes and rivers), the fraction of urban areas, and the fraction of nature, i.e. land areas that are not inland water or urban. There is no data available for the fraction of sea. The fraction data will be available directly from the climate data store website in netCDF format.

Accumulated surface fluxes

All energy fluxes at the surface are accumulated variables from the initial time of the forecast to the forecast hour in question with the unit J/m². They are considered positive downward to the surface. Average hourly energy fluxes in W/m² can be computed by subtracting two successive hourly accumulated variables and dividing by 3600 s. The albedo in units of % is defined by the ratio of average hourly upward radiation upward by the average hourly solar radiation downward. The net solar radiation is the difference between the solar radiation downward and the solar radiation upward. The solar radiation upward can be calculated by subtracting the net solar radiation from the downward solar radiation.

The net thermal radiation is the difference between the downward thermal radiation and the upward thermal radiation. The upward thermal radiation can be calculated by subtracting the net thermal radiation from the downward thermal radiation. The accumulated surface sensible heat flux is the conductive energy from the atmosphere to the surface. If this is going from the surface to the atmosphere it has negative values. The accumulated latent heat flux is the sum of all latent energy fluxes that are due to the phase transitions of water. Here condensation causes a positive latent heat flux to the surface, and evaporation causes a negative heat flux from the surface.

Precipitation and water fluxes

The daily surface precipitation analysis is the accumulation of total water (liquid and solid) that fall at ground level during the last 24 hours in kg/m² . The date represents the end of the accumulation period. It is available only at 0600 UTC.

The percolation (or drainage) is the mass per unit area of water that drains below the deepest soil level in the model whereas the surface runoff is the mass per unit area of water at the surface when saturation occurs. Percolation and surface runoff are calculated for the natural land, including soil, vegetation and snow whereas evaporation is available only as a mean value of a grid box. Those two variables can be used as input for hydrological model. Evaporation over the natural land fraction is not available.

Snow variables

The snow variables are related to natural land only. The fraction of snow cover represents the fraction of natural land which has snow on the ground. Snow on urban or lake fraction is not available.

References

Albergel, C., Dorigo, W., Reichle, R. H., Balsamo, G., de Rosnay, P., Muñoz-Sabater, J., Isaksen, L., de Jeu, R., and Wagner, W.: Skill and Global Trend Analysis of Soil Moisture from Reanalyses and Microwave Remote Sensing, Journal of Hydrometeorology, 14, 1259–1277, 825
, https://doi.org/10.1175/JHM-D-12-0161.1, 2013

Bazile E., E; Abida R. Abida , Verelle A. Verelle, Le Moigne P. Le Moigne, Szczypta C. Szczypta , (2017): MESCAN-SURFEX surface analysis, deliverable D2.8 of the UERRA project, http://www.uerra.eu/publications/deliverable-reports.html

Le Moigne P., Bazile E., Glinton M. and Verrelle A. (2022) : Documentation of the CERRA-Land system ( C3S delivrable C3S_D322_Lot1.1.1.12_202110_documentation_CERRA-Land)

Le Moigne , P., Besson , F., Martin , E., Boé , J., Boone , A., Decharme , B., Etchevers , P., Faroux , S.,
Habets , F., Lafaysse , M., Leroux , D., & and Rousset-Regimbeau , F. (2020). : The latest improvements with SURFEX v8.0 of the Safran–Isba–Modcou hydrometeorological model for France. Geoscientific Model Development, 13(9), 3925–3946France, Geosci. Model Dev., 13, 3925–3946, https://doi.org/10.5194/gmd-13-3925-2020, 2020.

Muñoz-Sabater J., Dutra , E., Agustí-Panareda , A., Albergel , C., Arduini , G., Balsamo , G., Boussetta , S., Choulga , M., Harrigan , S., Hersbach , H., Martens , B., Miralles , D. G., Piles , M., Rodríguez-Fernández , N. J., Zsoter , E., Buontempo , C., and Thépaut , J.-N.: ERA5-Land: A state-of-the-art global reanalysis dataset for land applications, Earth Syst.  SciSci. Data,13, 4349–4383, 2021.   https://doi.org/10.5194/essd-13-4349-2021 

Schimanke S., Isaksson L., and  Edvinsson L.: CERRA data user guide (C3S deliverable C3S_322_Lot1.4.1.3_CERRA_ data_user_guide)

Soci , C., Bazile , E., Besson , F., & and Landelius , T. (2016). High-resolution precipitation re-analysis system for climatological purposes. Tellus A, Dynamic Meteorology and Oceanography, 68:1, DOI: 10.3402/tellusa.v68.29879

...