...
Contributors: K. P. Nielsen (DMI), X. Yang (DMI), S. Agersten (MET Norway), P. Dahlgren (MET Norway), M. A. Ødegaard Køltzow (MET Norway), H. Schyberg (MET Norway), E. Støylen (MET Norway), J. Bojarova (SMHI)
...
Info | ||||
---|---|---|---|---|
| ||||
|
Introduction
This user guide describes the datasets released from the Arctic Regional Reanalysis service, which is part of the Copernicus Climate Change Service (C3S). The datasets will include the actual grid point reanalysis information on different levels (atmospheric vertical levels, surface including soil). This version also provides details on the uncertainty information provided. We will refer to the dataset as the CARRA (Copernicus Arctic Regional ReAnalysis) dataset.
First in this user guide you will find sections on the data availability, formats and data content followed by the information on the uncertainty of the data. For a short description of the principles and methods behind this reanalysis, numerical weather prediction and data assimilation, see the Annex at the end of this document.
Data availability
The Arctic regional reanalysis data can be downloaded from the Copernicus Climate Data Store (CDS) at the following locations:
...
The scripts in section 2.1 are also available for download via the open GitHub athttps://github.com/metno/carra_cds
CDS API example, C3S Arctic reanalysis
The example below shows a Python script for retrieving a subset of the C3S Arctic Reanalysis data from CDS for the month of July 2012 for the West domain into an output file named: "CARRA-West_2012-07.grib". Such script can be derived from the download form (see "Show API request" at the bottom of the form).
...
In the example above single level data are retrieved from the western domain. If instead pressure level data from the eastern domain is required the following command should be used in the script: c.retrieve('reanalysis-carra-pressure-levels', 'domain': 'east_domain'...
)
Data formats
The reanalysis data is available in two formats: GRIB and NetCDF.
GRIB
The GRIB data format is a standard binary format for modelled meteorological data governed by the World Meteorological Organisation (WMO). This data format stored the data in a very compact way and the metadata is contained within the files. It can be read with the tool eccodes that is available from the European Centre for Medium-range Weather Forecasts (ECMWF): https://confluence.ecmwf.int/display/ECC/GRIB+tools. National weather services and weather companies use a wide array of graphical tools to display GRIB data. The GRIB data format is recommended for users comfortable and familiar with this format.
NetCDF
NetCDF (Network Common Data Form, https://www.unidata.ucar.edu/software/netcdf/docs/index.html ) is an open standard data format developed and supported by the American University Corporation for Atmospheric Research (UCAR). It is a widely used format for geoscientific data. A large array of access libraries and applications for reading and plotting NetCDF data exist: https://www.unidata.ucar.edu/software/netcdf/software.html.
Data content
Single level variables
Diagnostic single level output is available in three hourly intervals at 00, 03, 06, 09, 12, 15, 18 and 21 UTC.
Long forecasts are available from the forecasts initiated at 00 and 12 UTC. Long forecasts include forecast lengths of 1, 2, 3, 4, 5, 6, 9, 12, 15, 18, 21, 24, 27 and 30 hours.
Short forecasts of 1, 2 and 3 hours are made for the forecasts initiated at 03, 06, 09, 15, 18 and 21 UTC.
For most of the variables the shortest forecast data are recommended to use. In general the data quality decreases with forecast length. On the other hand, for variables that are affected by spin-up effects - that is the model needs to run for a certain number of hours before these variables have an optimal quality - the longer forecasts can be better to use. The cloud and precipitation variables are directly affected by spin-up. For time integrated quantities such as theseprecipitation, accumulation over 12 hours between +6 and +18 h forecasts are recommended. Likewise 24 hour accumulation can be obtained as the difference between +30 hour and +6 hour forecasts. However, if accuracy in timing of precipitation events is of very high importance in the application, an option could be to combine hourly forecasts from each of the 8 analysis times (00, 03, 06, 09, 12, 15, 18, 21 UTC) for lead times 1, 2 and 3, which then will be slightly affected by spin-up. It is not possible to make general recommendations on this issue. Users , therefore users are advised to check this for themselves.make their own choices based on the general guidelines described here. (On spin-up of precipitation, see also here.)
Table 1: Overview of single level variables. Some variables have not yet been uploaded to CDS, these are marked with * and unfortunately not be included in the first batch of released reanalysis data. Parameters labelled with TBD ("To Be Determined") do not yet have short name and GRIB code definitions.Most static fields (except land-sea mask and orography) marked with * are available only as NetCDF-files below. Anchor table1 table1
Precipitation, cloud water and humidity | Name | |||||
Precipitation, cloud water and humidity | ||||||
Name | ||||||
Short Name | Unit | |||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
2r | % | 260242 | yes | yes | 2m | |
2sh | kg/kg | 174096 | yes | yes | 2m | |
tciwv | kg/m2 | 260057 | yes | yes | vertically integrated above the surface | |
tclw | kg/m2 | 78 | ||||
no | yes | vertically integrated above the surface | ||||
tciw | kg/m2 | 79 | ||||
no | yes | vertically integrated above the surface | ||||
tcolg | kg/m2 | 260001 | yes | yes | vertically integrated above the surface | |
tp | kg/m2 | 228228 | no | yes | surface | |
tirf | kg/m2 | 235015 | no | yes | surface | |
titspf | kg/m2 | 260645 | no | yes | surface | |
ptype | integer code | 260015 | no | yes | surface | |
sro | kg/m2 | 174008 | no | yes | surface | |
Percolation (drainage) | perc | kg/m2 | 260430 | no | yes | sub-surface |
Temperature and wind speed | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
10si | m/s | 207 | yes | yes | 10m | |
10wdir | degrees | 260260 | yes | yes | 10m | |
10m u-component of wind (defined relative to the rotated model grid) | 10u | m/s | 165 | yes | yes | 10m |
10m v component of wind | 10v | m/s | 166 | yes | yes | 10m |
| 10efg | m/s | 260646 | no | yes | 10m |
| 10nfg | m/s | 260647 | no | yes | 10m |
10fg | m/s | 49 | no | yes | 10m | |
mx2t | K | 201 | no | yes | 2m | |
mn2t | K | 202 | no | yes | 2m | |
2t | K | 167 | yes | yes | 2m | 2m temperature for the sea tile* | TBD | K | TBD | yes
2m temperature for the inland water tile* | TBD | K | TBD | yes | yes | 2m |
2m temperature for the nature tile* | TBD | K | TBD | yes | yes | 2m |
2m temperature for the urban tile* | TBD | K | TBD | yes | yes | 2m |
yes | 2m | |||||
skt | K | 235 | yes | yes | Surface |
Accumulated fluxes | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
al | % | 260509 | ||||
yes** | yes | surface | ||||
eva | kg/m2 | 260259 | no | yes | surface | |
tisef | kg/m2 | 235072 | no | yes | surface | |
sshf | J/m2 | 146 | no | yes | surface | |
slhf | J/m2 | 147 | no | yes | surface | |
tislhef | J/m2 | 235019 | no | yes | surface | |
tislhsf | J/m2 | 235071 | no | yes | surface | |
radiation | ||||||
dsrp | J/m2 | 47 | no | yes | surface | |
tidirswrf | J/m2 | 260264 | no | yes | surface | |
ssr | J/m2 | 176 | no | yes | surface | |
ssrd | J/m2 | 169 | no | yes | surface | |
ssrc | J/m2 | 210 | no | yes | surface | |
str | J/m2 | 177 | no | yes | surface | |
strd | J/m2 | 175 | no | yes | surface | |
strc | J/m2 | 211 | no | yes | surface | |
Top net solar radiation | ||||||
tsr | J/m2 | 178 | no | yes | surface | |
ttr | J/m2 | 179 | no | yes | surface | |
| tisemf | kg⋅m/s | 235017 | no | yes | surface |
tisnmf | kg⋅m/s | 235018 | no | yes | surface | |
Pressure | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
msl | Pa | 151 | yes | yes | surface (scaled to sea level) | |
sp | Pa | 134 | yes | yes | surface | |
Geometric cloud properties | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
hcc | % | 3075 | yes | yes | above 5000m | |
mcc | % | 3074 | yes | yes | 2500m - 5000m | |
lcc | % | 3073 | yes | yes | surface - 2500m | |
tcc | % | 228164 | yes | yes | above ground | |
fog | % | 260648 | ||||
no | yes | lowest model level | ||||
vis | m | 3020 | yes | yes | lowest model level | |
cdcb | m | 260107 | yes | yes | - | |
cdct | m | 260108 | yes | yes | - | |
Snow | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
rsn | kg/m3 | 33 | yes | yes | surface | |
sd | kg/m2 | 228141 | yes | yes | surface | |
fscov | 0-1 | 260289 | yes | |||
no | surface | |||||
asn | % | 228032 | yes | no | ||
surface | ||||||
Surface roughness lengths | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
sr | m | 173 | yes | no | surface | |
srlh | m | 260651 | yes | no | surface | |
Sea states | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
sst | K | 34 | yes | no | surface | |
ci | 0-1 | 31 | yes | no | surface | |
sist | K | 260649 | yes | |||
yes | surface | |||||
sithick | m | 174098 | yes | no | surface | |
sitd | m | 260650 | yes | yes | surface | |
Static fields | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | Height | |||
lsm | ||||||
0-1 | 172 | no | no | surface | ||
Sea tile fraction* | ||||||
NA | ||||||
0-1 | ||||||
NA | no | no | surface | |||
Inland water tile fraction* | ||||||
NA | ||||||
0-1 | ||||||
NA | no | no | surface | |||
Urban tile fraction* | ||||||
NA | ||||||
0-1 | ||||||
NA | no | no | surface | |||
Nature tile fraction* | ||||||
NA | ||||||
0-1 | ||||||
NA | no | no | surface | |||
Glacier fraction* | NA | 0-1 | NA | no | no | surface |
Subgrid orography average slope* | TBD | % | ||||
NA | 0-1 | NA | no | no | surface | |
Subgrid orography standard deviation* | NA | m | NA | |||
no | no | surface | ||||
orog | m | 228002 | no | no | surface |
Soil level variables
The static fields marked as * above are available as NetCDF files for the West and East domain respectively here: fractions.west.nc and fractions.east.nc
** Albedo is available at analysis time, but as a static climatological value which we discourage for use. Only the forecast time albedos should be used, see end of Section 5.5 below.
Soil level variables
Soil level variables are given for two model depths, where the first depth is the soil surface and the Soil level variables are given for two model depths, where the first depth is the soil surface and the second depth is the so-called root depth. The root depth varies with the cover type climatology. Please note that the soil level variables are accommodated in the Arctic Regional Reanalysis single level variables catalogue entry.
...
Long forecasts are available from the forecasts initiated at 00 and 12 UTC. Long forecasts include forecast lengths of 1, 2, 3, 4, 5, 6, 9, 12, 15, 18, 21, 24, 27 and 30 hours.
Short forecasts of 1, 2 and 3 hours are made for the forecasts initiated at 03, 06, 09, 15, 18 and 21 UTC.
Anchor | ||||
---|---|---|---|---|
|
Soil level variables | |||||
Name | Short Name | Unit | |||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | |||
vsi | m³/m³ | 260644 | yes | yes | |
vsw | m³/m³ | 260199 | yes | yes |
Model level variables
Model level variables are output at 65 hybrid model levels of the HARMONIE-AROME model. These follow the surface at the lowest levels and are gradually evolved into pure pressure levels at the highest levels. These are the levels at which the model computations are done. The height level and pressure level variables are interpolated from these data.
...
Table 3: Overview of model level variables Anchor table3 table3
Model level variables | |||||
Name | Short Name | Unit |
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2 | |||
q | kg/kg | 133 | yes | yes | |
t | K | 130 | yes | yes | |
u-component of wind (defined relative to the rotated model grid) | u | m/s | 131 | yes | yes |
v-component of wind (defined relative to the rotated model grid) | v | m/s | 132 | yes | yes |
ccl | % | 260257 | yes | yes | |
clwc | kg/kg | 246 | yes | yes | |
ciwc | kg/kg | 247 | yes | yes | |
crwc | kg/kg | 75 | yes | yes | |
cswc | kg/kg | 76 | yes | yes | |
grle | kg/kg | 260028 | yes | yes | |
tke | J/kg | 260155 | yes | yes |
Pressure level variables
Pressure level variables are interpolated to 23 specific pressure levels: 1000, 950, 925, 875, 850, 800, 750, 700, 600, 500, 400, 300, 200, 100, 70, 50, 30, 20 and 10 hPa. Thus, they are on isobaric surfaces.
...
Long forecasts are available from the forecasts initiated at 00 and 12 UTC. Long forecasts include forecast lengths of 1, 2, 3, 4, 5, 6, 9, 12, 15, 18, 21, 24 and 30 hours.
Short forecasts of 1, 2 and 3 hours are made for the forecasts initiated at 03, 06, 09, 15, 18 and 21 UTC.
Anchor | ||||
---|---|---|---|---|
|
Pressure level variables | ||||||
Name | Short Name | Unit | ||||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | ||||
r | % | 157 | yes | yes | ||
t | K | 130 | yes | yes | ||
u-component of wind (defined relative to the rotated model grid) | u | m/s | 131 | yes | yes | |
v-component of wind (Component defined relative to the rotated model grid) | v | m/s | 132 | yes | yes | |
wz | m/s | 260238 | yes | yes | ||
ccl | % | 260257 | yes | yes | ||
clwc | kg/kg | 246 | yes | yes | ||
ciwc | kg/kg | 247 | yes | yes | ||
crwc | kg/kg | 75 | yes | yes | ||
cswc | kg/kg | 76 | yes | yes | ||
grle | kg/kg | 260028 | yes | yes | ||
papt | K | 3014 | yes | yes | ||
z | m²/s² | 129 | yes | yes | ||
pv | K·m²/ (kg·s) | 60 | yes | yes |
Height level variables
Height level variables are interpolated to 11 specific height levels: 15, 30, 50, 75, 100, 150, 200, 250, 300, 400 and 500 metres above the surface.
...
Anchor | ||||
---|---|---|---|---|
|
Height level variables | |||||
Name | Short Name | Unit | |||
Param ID | Analysis: 0,3,...,21 | Forecast: 1,2,3,… | |||
r | % | 157 | yes | yes | |
t | K | 130 | yes | yes | |
ws | m/s | 10 | yes | yes | |
wdir | deg | 3031 | yes | yes | |
clwc | kg/kg | 246 | yes | yes | |
ciwc | kg/kg | 247 | yes | yes | |
Pressure | pres | Pa | 54 | yes | yes |
Details about the data fields
All data fields are model grid box (2.5 km X 2.5 km = 6.25 km2) averages. When output data are compared to local data near the coast or close to a glacier boundary, it should be remembered that the model variables represent such averages. For instance, to compute the surface forcing at specific glacial sites near a glacier boundary, it will be most representative to choose output from a model grid box that is fully on the glacier, rather than the closest model grid box if this only partially covers the glacier. Also, very local rain shower intensity can be higher than the modelled average grid box intensity.
Output variables can be either instantaneous, accumulated, or maxima/minima from a given period. This is specified for each of the variables listed below.
Precipitation and water fluxes
Precipitation at the surface is output as two separate types: Rain and total solid precipitation. These have the unit kg/m2, which for rain with a density of 1000 kg/m3 is the same unit as mm. For the solid precipitation, using this unit removes confusion between mm water equivalent and mm thickness. Total solid precipitation is the sum of the model variables snow and graupel. Total precipitation is the sum of all three precipitation types. The precipitation species include both convective and stratiform precipitation and are available only for the forecast time steps. They are accumulated variables meaning that they are accumulated from the beginning of the forecast. For instance, the 24h-forecast includes accumulated precipitation over 24 hours. Hourly precipitation can be retrieved by subtracting two accumulated precipitation fields with one-hour separation.
...
Water fluxes at the surface are output as surface runoff, percolation (drainage), water evaporation and snow sublimation. All these variables have units of kg/m2. Surface runoff occurs, when the model soil is saturated with water and more precipitation comes. Percolation is drainage of water below the deep soil level in the model. The water evaporation can be both positive and negative, where negative values signify condensation. This is also true for the snow sublimation.
Cloud and water vapour integrated variables
All cloud and water vapour variables are instantaneous, i.e. they are given for the time step at which they are output. Vertically integrated water vapour is given in units of kg/m2. It is vertically integrated from the surface to the top of the atmosphere. In practice it is computed from the specific water vapour on the 65 model levels (see section 4.4). Likewise, integrated cloud liquid water, integrated cloud ice water, and integrated graupel are computed from the specific cloud liquid water, cloud ice and graupel on the 65 model levels.
...
Visibility is given with the unit m and is calculated from the cloud water, cloud ice, rain, snow and graupel present at the lowest model level. For cloud and precipitation free conditions mist is calculated from the lowest model level relative humidity and the concentration of cloud condensation nuclei. Empirical relations are used (pers. comm. Esbjörn Olsson, SMHI, 2013). Cloud base height and cloud top height are output in units of m above the surface. They are defined for the highest and lowest model level with more than 4/8 cloud cover.
Variables at 10-metre height
All 10-metre height wind variables are instantaneous, i.e. they are given for the time step at which they are output. Diagnosed u and v wind components (u and v) are output at the WMO standard height of 10 m above the surface with the unit m/s. The u- and v-components follow the direction of the Lambertian model grid with the u-component being directed 90 degrees clockwise relative to the v-component. The diagnosed winds are computed from the winds at the lowest full model level (see section 4.3). The 10-metre u wind component (u10m) in the case of a stable or neutral boundary layer is calculated as
...
The wind direction D clockwise from North can be calculated as
Anchor | ||||
---|---|---|---|---|
|
LaTeX Formatting |
---|
$$ D = atan tan_{2}^{-1} \frac({u/U},{v/V}) + 180^{\circ} + \alpha, (4) $$ |
where 𝛼 where 𝛼 is the local rotation of the model grid relative to North, and tan2-1 atan2 is the very specific 2-argument arcus tangens function atan2, which is included in most programming languages. With atan2 both the sign in the numerator and the denominator are independently important for the resulting angle. Take care For definition of the atan2 function, see for instance at https://en.wikipedia.org/wiki/Atan2. Take care to check if the atan2 result is in radians, in which case it should be converted to degrees with the factor 180°/π. Also take care to check if the resulting direction is between 0° and 360°. Note that the wind direction is the direction from which the wind comes! The grid rotation angle 𝛼 can be computed with this script: https://github.com/metno/NWPdocs/wiki/From/Examples/#wind-direction-obtained-from-x-y-wind-to-wind-direction.
10-metre u and v wind gust components are also output. These are computed from the diagnosed 10-metre winds and the turbulent kinetic energy (pers. comm. Gwenaëlle Hello, Meteo France, 2007).
Variables at 2-metre height
2-metre temperature (T2m) is diagnosed from the so-called skin temperature at 0 metre height (T0m) Since the 2-metre temperature is often referred to as the surface temperature we here, in order to avoid misunderstandings, explicitly use 2-metre and 0-metre subscripts. and the temperature at the lowest model level (Tlml). For stable and neutral boundary layers it is calculated as
Anchor | ||||
---|---|---|---|---|
|
...
2-metre specific humidity in units of kg/kg is diagnosed in the same way as the 2-metre temperature, only from the specific humidity at 0 metre height and at the lowest model level rather than the temperatures at these levels. From this and the saturation specific humidity the 2-metre relative humidity in % is computed.
Variables at 0-metre height
The surface pressure is given in units of Pascal (Pa). From this the mean sea level pressure [Pa] is computed by reducing the surface pressure to the mean sea level. To avoid confusion with the 2-metre “surface” temperature, the 0-metre temperature is referred to as the skin temperature. It is given in units of K. Both the pressure and temperature variables are instantaneous, i.e. they are given for the time step at which they are output. The surface is the lowest model half-level. Thus, the variables at this level are explicitly calculated on each model iteration.
All energy fluxes at the surface are output as accumulated variables from the initial time of the forecast to the forecast hour in question with the unit J/m2. They . Note that a different accumulation applies to the albedo, see below. The other accumulated variables are considered positive downward to the surface. Energy fluxes are not output variables at the analysis times. Average hourly energy fluxes in W/m2 can be computed by subtracting two successive hourly accumulated variables and dividing by 3600 s. The solar radiation variables at the surface are accumulated downward, direct, direct normal and net solar radiation. Direct normal solar radiation is considered on a plane perpendicular to the direction to the sun, while the other solar variables are considered on a horizontal surface. The albedo in units of % is also given. Multiplying this with the accumulated downward solar radiation gives the accumulated upward solar radiation. The net solar radiation is the difference between the downward solar radiation and the upward solar radiation. Note that the albedo should only be used from forecast data, as the analysis time albedo are is incorrect – and not used in the model. The variable accumulated net clear sky solar radiation is the net solar radiation of a cloud free atmosphere. Dividing this with one minus the albedo gives the accumulated downward clear sky solar radiation. The thermal radiance variables at the surface are accumulated downward, net, and net clear sky thermal radiation. The thermal radiation variables are all considered on a horizontal surface. 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. The latent heat due to evaporation and sublimation are given as individual output variables.
The u and v components of the accumulated surface momentum flux are given as output variables in units of kg m/s. The momentum roughness length and the heat roughness length as used in the model are given as output in units of m. Note that the roughness lengths should only be used from forecast data, as the analysis time values are incorrect – and not used in the model.
Variables at the top of the atmosphere
The albedo is accumulated with a special procedure, noting that it changes over time depending on e.g. snowfall. It is derived from the formula: albedo = 1 - SWnet/SW↓, where SW↓ is the downward solar flux at the surface, and SWnet is the net (downward minus upward) solar flux at the surface. The solar fluxes are accumulated over a given time interval. This time interval is one hour until the +6h forecast range and three hours afterwards. So these albedos represent averages over hourly or 3-hourly periods prior to the forecast time. Instantaneous albedos at the output times would be less precise than using these accumulated variables, see alsosection 4 in the note by Hogan, 2015. The problem with this procedure is that in some cases unphysical albedo values (>100%) may occur due to rounding errors. There is also an analysis time albedo provided in the data set, which can be quite different from the actual albedo at forecast times ad therefore it is not recommended to be used.
Variables at the top of the atmosphere
At the top of the atmosphere (TOA) the accumulated solar net TOA radiation and the accumulated thermal net TOA radiation are output variables in units of J/m2. These are both considered on a horizontal surface and are both positive in the downward direction. Since the downward solar TOA radiation is always larger than the upward solar TOA radiation, the solar net TOA radiation is always positive. Since there is virtually no downward thermal TOA radiation, the thermal net TOA radiation is always negative. The TOA is the highest model half-level. Thus, the variables at this level are At the top of the atmosphere (TOA) the accumulated solar net TOA radiation and the accumulated thermal net TOA radiation are output variables in units of J/m2. These are both considered on a horizontal surface and are both positive in the downward direction. Since the downward solar TOA radiation is always larger than the upward solar TOA radiation, the solar net TOA radiation is always positive. Since there is virtually no downward thermal TOA radiation, the thermal net TOA radiation is always negative. The TOA is the highest model half-level. Thus, the variables at this level are explicitly calculated.
Snow variables
The snow variables are given as instantaneous values from the most recent model time step relative to the output time. The snow density output unit is kg/m3, the snow water equivalent (SWE) output is in units of kg/m2, and the snow fraction output has fractional units in the range 0-1. The snow depth can be derived from these variables.
Sea and sea
...
ice variables
For the sea the sea surface temperature is is output in units of K. For areas partially or completely covered with sea - ice, the following variables are output: Sea - ice area fraction [-], upper layer sea ice temperature [K], sea - ice thickness [m], and sea - ice snow thickness [m]. Here the For sea -ice fraction and sea surface temperatures are in fact interpolated input data and are only updated once every day. During the ice thickness please note that the routine that computes this variable does not reproduce the evolution of ice thickness with all its complexity. Rather this variable should be treated as a rough estimate in order to get reasonable estimations for the energy fluxes. The sea ice fraction and sea surface temperatures are in fact interpolated input data and are only updated once every day. During the course of a forecast they are kept constant. All other sea and sea - ice variables are given as instantaneous values from the most recent model time step relative to the output time.
Soil level variables
The soil level variables are computed at the skin depth immediately at the surface and at the "root depth", that is the average depth from which the vegetation retrieves its water. This level varies depending on the surface type and is also defined for surface types without vegetation - for instance rock. The volumetric soil moisture content is the volume concentration of liquid water at root depth. The volumetric soil frozen water content is the volume concentration of ice at root depth. Please note that the soil level variables are listed in the "single level" catalogue entry. The root depth variables are used to account for effect of the evapotranspiration on the temperature and humidity at and immediately above the surface. They are a means to an end rather than a detailed model of what the humidity is at specific depths in the soil. For anyone who wishes to perform detailed modelling of processes below the surface we recommend to use surface fluxes from our reanalysis dataset as forcing for more advanced sub-surface models.
Model level variables
The model level variables are computed at the full model levels, and are given as instantaneous values from the most recent model time step relative to the output time. There are 65 vertical model levels in HARMONIE-AROME. These full model levels are hybrid-sigma coordinates that are counted from the model top toward the surface. They go from being pure pressure levels, i.e. levels with constant pressure starting at 10 hPa, 30 hPa, etc. to being relative to the surface topography in height. Level 64 is at approximately 30 m height and level 65 is at approximately 12 m height above the surface. For a more detailed description of the vertical model layers, see Annex 8.3. The following thermodynamic variables are the output variables at model levels: Temperature [K], u-component of wind [m/s], v-component of wind [m/s], turbulent kinetic energy [J/kg]. Here turbulent kinetic energy is the mean kinetic energy per unit mass from eddies in turbulent flow. Note that the HARMONIE-AROME weather forecasting model with 2.5 x 2.5 km2 resolution does not explicitly resolve this turbulent energy. The u- and v- wind components follow the direction of the Lambertian model grid with the u-component being directed 90 degrees clockwise relative to the v-component. From these model level wind components, the model level wind speed and wind direction relative north can be calculated with Equations 3 and 4. The grid rotation angle 𝛼 can be computed with this script: https://github.com/metno/NWPdocs/wiki/From-x-y-wind-to-wind-direction.
The following moisture, cloud and precipitation variables are output at model levels: Specific humidity, specific cloud liquid water, specific cloud ice water content, specific cloud rain water content, specific cloud snow water content, and graupel. All of these variables are given in units of kg/kg. Additionally, the cloud cover in % is output on model levels. The model level cloud covers can be used to compute the total cloud cover based on other cloud cover overlap assumptions than the one described in section 5.2.
Pressure level variables
Pressure level variables are output at the pressure levels: 1000, 950, 925, 900, 875, 850, 825, 800, 750, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, and 10 hPa, and are given as instantaneous values from the most recent model time step relative to the output time. The model level variables described in section 4.3 are vertically interpolated to fixed pressure levels. The following thermodynamic variables are output at pressure levels: Temperature [K], u-component of wind [m/s], v-component of wind [m/s], geometric vertical velocity [m/s], geopotential [m2/s2], pseudo-adiabatic potential temperature [K], and potential vorticity [(K m2)/(kg s)]. Here the geopotential is the work required to lift an air parcel of unit mass from mean sea level to the given pressure level. The pseudo-adiabatic potential temperature is the temperature that an air parcel would have if it were first expanded through a pseudo-adiabatic process to 0 hPa pressure and thereafter compressed to a pressure of 1000 hPa through a dry-adiabatic process. The potential vorticity is a measure of the capacity for air to rotate in the atmosphere. The u- and v- wind components follow the direction of the Lambertian model grid with the u-component being directed 90 degrees clockwise relative to the v-component. From these pressure level wind components, the pressure level wind speed and wind direction relative north can be calculated from Equations 3 and 4. The grid rotation angle 𝛼 can be computed with this script: https://github.com/metno/NWPdocs/wiki/From-x-y-wind-to-wind-direction.
The following moisture, cloud and precipitation variables are output at pressure levels: Relative humidity [%], cloud cover [%], specific cloud liquid water content [kg/kg], specific cloud ice water content [kg/kg], specific cloud rain water content [kg/kg], specific cloud snow water content [kg/kg], and graupel [kg/kg]. These pressure level variables are interpolated from the full model level variables described in section 4.3.
Height level variables
Height level variables are output at the following fixed heights above the surface: 15 m, 30 m, 50 m, 75 m, 100 m, 150 m, 200 m, 250 m, 300 m, 400 m and 500 m. They are interpolated to these heights from the full model level variables described in section 4.3. The following thermodynamic variables are output at height levels: Temperature [K], pressure [Pa], wind speed [m/s], and wind direction [degrees from North]. Note here that the wind direction is defined as the direction from which the wind comes!
The following moisture, cloud and precipitation variables are output at height levels: Relative humidity [%], specific cloud liquid water [kg/kg], and specific cloud ice water content [kg/kg]. These height level variables are interpolated from the full model level variables described in section 4.3.
Static fields
Static fields are output 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 2.5 x 2.5 km2 in units of %, and the orography in units of m. There are two more orography-related static parameters: subgrid orography average slope and subgrid orography standard deviation. For each model grid box in HARMONIE-AROME 4 tile fractions are defined in units of %fraction. These are: The fraction of sea, 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. The fraction of glaciers is also output. This is assumed to be a constant field with glacier extents representative of the middle of the full reanalysis period (19971991-2021present). Glacier extent in remote Arctic locations is not available as accurately mapped yearly datasets. The official maps are outdated due to major calving events in the recent decades. HARMONIE-AROME has not yet been designed to deal with changing land-sea masks or other surface classifications. Thus, these are static fields.
What are the uncertainties of the data fields?
A user of the reanalysis fields may be interested in knowing the accuracy or uncertainty of the provided data. Uncertainty is an expression of the expected or average deviations from what is considered the "true" value, and several approaches to estimate uncertainties are possible.
...
Anchor | ||||
---|---|---|---|---|
|
Name | Levels |
Pressure | Surface as well as mean sea level |
U-component of wind | Pressure levels (50 - 1000 hPa) |
V-component of windy | Pressure levels (50 - 1000 hPa) |
Temperature | Pressure levels (50 - 1000 hPa) |
Geopotential (East domain only) | Pressure levels (50 - 1000 hPa) |
Relative humidity (East domain)/Specific humidity (West domain) | Pressure levels (50 - 1000 hPa) |
In section 6.5 we present uncertainties measured as statistics of actual deviations from observations (also known as verification statistics). Note that observations are a reference not identical to the actual truth, as they also will have uncertainties as well as representativeness issues. This verification statistics is provided for a set of near-surface quantities which are covered by the meteorological observation network, including 2m temperature, 10m winds and precipitation.
Field based atmospheric uncertainty estimation
The Earth system observations are unevenly distributed in space and time and the largest part of the Arctic Regional Reanalysis (CARRA) domains has sparse data coverage. Additional information in the form of numerical simulations is needed in order to construct the gridded dataset. Reanalysis employs a data assimilation system for this purpose, the same device that is used to initialize a numerical weather prediction forecasting system with the latest observations. The analysis gridded dataset is obtained as an optimally weighted average of observations and short-range weather forecast (numerical simulations) where weights are dependent on the distance to observations, the forecast error and observation error statistics and the statistical relationships valid between meteorological quantities (model state variables), such as horizontal wind, temperature, specific humidity, surface pressure, geopotential.
...
The methodology does not account for systematic errors or for uncertainties associated with the forecasting model used for numerical simulations, such as for example parameterization of unresolved processes. In addition to the climatological upper bound analysis error standard deviation that we have discussed above, we provide a refined climatological analysis error standard deviation. The refined climatological analysis error standard deviation is the climatological upper bound analysis error standard deviation multiplied by a rescaling factor. The rescaling factor is obtained by confronting the climatological upper bound forecast error variance against observation error variance in observation space. The multiplication by rescaling factor allows compensating for the limitations of the methodology mentioned above to some extent. The refined climatological analysis error standard deviation is a more reliable estimate of the uncertainty of the CARRA datasets than the climatological upper bound analysis error standard deviation in the areas with dense observation coverage. At the same time, the short-range forecast error variance is always lower in the areas covered by observations than in the areas of sparse observation coverage. Thus, the refined climatological analysis error standard deviation can underestimate the actual uncertainty of the Arctic reanalysis data in the areas with poor observation coverage. We propose to use the refined analysis error standard deviation in data rich areas and to use the upper bound standard deviation in data sparse areas as a measure of uncertainty.
Uncertainties for the CARRA-East model domain
In this section we provide the actual estimates of the uncertainty in the form of the "Upper Bound" standard deviation and the "Refined" standard deviation for analysis errors over the CARRA-East domain. We provide both summer and winter statistics. The uncertainty estimation is given both in form of tables for surface pressure and mean sea level pressure (Table 7), u-component of wind (Table 8), v-component of wind (Table 9), temperature (Table 10), geopotential (Table 11) and relative humidity (Table 12) and in form of figures for u-wind component (Figure 2), v-wind component (Figure 3), temperature (Figure 4) geopotential (Figure 5) and relative humidity (Figure 6). All tables contain the name of the meteorological variable, short abbreviation, units and vertical levels. Uncertainty of the surface pressure field is provided at 0 m above ground and at 0 m above mean sea level. The uncertainty of the u- and v- wind components, temperature, geopotential and relative humidity is provided at standard pressure levels: 50hPa, 100hPa, 150 hPa, 200hPa, 250 hPa, 300hPa, 400hPa, 500hPa, 600hPa, 700hPa, 800hPa, 850hPa, 900hPa, 925hPa, 950hPa and 1000hPa. The figures show uncertainty estimates as a function of standard pressure levels. The Upper Bound standard deviation (STDV) of the analysis errors is shown as a solid curve and the Refined Standard Deviation of the analysis error is shown as a dashed curve. Summer statistics are presented in the plots to the left and winter statistics are presented in the plots to the right.
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Summer statistics | Winter statistics | ||
Upper Bound STDV | Refined STDV | Upper Bound STDV | Refined STDV | ||||
Surface pressure | sp | Pa | 0m above ground | 37.20 | 26.78 | 41.14 | 38.67 |
Mean sea level pressure | msl | Pa | 0m above sea level | 38.18 | 27.49 | 42.15 | 39.62 |
Anchor | ||||
---|---|---|---|---|
|
Figure 2: Climatological analysis error standard deviation for u-component of wind: summer statistics (left), winter statistics (right) as function of standard vertical pressure levels for the CARRA-East domain.
...
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Summer statistics | Winter statistics | ||
Upper Bound STDV | Refined STDV | Upper Bound STDV | Refined STDV | ||||
u-component of wind | u | m/s | 50hPa | 0.545 | 0.223 | 0.556 | 0.228 |
u-component of wind | u | m/s | 100hPa | 0.372 | 0.153 | 0.397 | 0.169 |
u-component of wind | u | m/s | 150hPa | 0.584 | 0.239 | 0.619 | 0.254 |
u-component of wind | u | m/s | 200hPa | 0.865 | 0.355 | 0.923 | 0.378 |
u-component of wind | u | m/s | 250hPa | 1.077 | 0.442 | 1.089 | 0.447 |
u-component of wind | u | m/s | 300hPa | 1.290 | 0.530 | 1.240 | 0.508 |
u-component of wind | u | m/s | 400hPa | 1.282 | 0.526 | 1.292 | 0.530 |
u-component of wind | u | m/s | 500hPa | 1.284 | 0.526 | 1.364 | 0.559 |
u-component of wind | u | m/s | 600hPa | 1.322 | 0.648 | 1.392 | 0.738 |
u-component of wind | u | m/s | 700hPa | 1.346 | 0.740 | 1.404 | 0.927 |
u-component of wind | u | m/s | 800hPa | 1.344 | 0.833 | 1.324 | 1.006 |
u-component of wind | u | m/s | 850hPa | 1.331 | 0.865 | 1.274 | 1.045 |
u-component of wind | u | m/s | 900hPa | 1.311 | 0.905 | 1.275 | 1.122 |
u-component of wind | u | m/s | 925hPa | 1.315 | 0.921 | 1.295 | 1.178 |
u-component of wind | u | m/s | 950hPa | 1.330 | 0.958 | 1.305 | 1.227 |
u-component of wind | u | m/s | 1000hPa | 1.302 | 0.937 | 1.172 | 1.102 |
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Summer statistics | Winter statistics | ||
Upper Bound STDV | Refined STDV | Upper Bound STDV | Refined STDV | ||||
v-component of wind | v | m/s | 50hPa | 0.541 | 0.222 | 0.541 | 0.222 |
v-component of wind | v | m/s | 100hPa | 0.376 | 0.154 | 0.384 | 0.157 |
v-component of wind | v | m/s | 150hPa | 0.597 | 0.245 | 0.599 | 0.246 |
v-component of wind | v | m/s | 200hPa | 0.884 | 0.362 | 0.898 | 0.368 |
v-component of wind | v | m/s | 250hPa | 1.105 | 0.453 | 1.061 | 0.435 |
v-component of wind | v | m/s | 300hPa | 1.337 | 0.548 | 1.209 | 0.496 |
v-component of wind | v | m/s | 400hPa | 1.313 | 0.538 | 1.270 | 0.521 |
v-component of wind | v | m/s | 500hPa | 1.328 | 0.544 | 1.351 | 0.554 |
v-component of wind | v | m/s | 600hPa | 1.363 | 0.668 | 1.391 | 0.737 |
v-component of wind | v | m/s | 700hPa | 1.383 | 0.761 | 1.402 | 0.925 |
v-component of wind | v | m/s | 800hPa | 1.382 | 0.857 | 1.332 | 1.012 |
v-component of wind | v | m/s | 850hPa | 1.363 | 0.886 | 1.282 | 1.051 |
v-component of wind | v | m/s | 900hPa | 1.342 | 0.926 | 1.277 | 1.124 |
v-component of wind | v | m/s | 925hPa | 1.344 | 0.941 | 1.293 | 1.177 |
v-component of wind | v | m/s | 950hPa | 1.361 | 0.980 | 1.299 | 1.221 |
v-component of wind | v | m/s | 1000hPa | 1.318 | 0.949 | 1.167 | 1.097 |
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Summer statistics | Winter statistics | ||
Upper Bound STDV | Refined STDV | Upper Bound STDV | Refined STDV | ||||
temperature | t | K | 50hPa | 0.141 | 0.058 | 0.142 | 0.058 |
temperature | t | K | 100hPa | 0.122 | 0.050 | 0.122 | 0.050 |
temperature | t | K | 150hPa | 0.230 | 0.094 | 0.224 | 0.092 |
temperature | t | K | 200hPa | 0.466 | 0.191 | 0.441 | 0.181 |
temperature | t | K | 250hPa | 0.490 | 0.201 | 0.438 | 0.180 |
temperature | t | K | 300hPa | 0.392 | 0.161 | 0.406 | 0.167 |
temperature | t | K | 400hPa | 0.314 | 0.129 | 0.341 | 0.140 |
temperature | t | K | 500hPa | 0.346 | 0.142 | 0.359 | 0.147 |
temperature | t | K | 600hPa | 0.399 | 0.196 | 0.399 | 0.211 |
temperature | t | K | 700hPa | 0.465 | 0.256 | 0.465 | 0.307 |
temperature | t | K | 800hPa | 0.537 | 0.333 | 0.485 | 0.369 |
temperature | t | K | 850hPa | 0.586 | 0.381 | 0.508 | 0.417 |
temperature | t | K | 900hPa | 0.614 | 0.424 | 0.520 | 0.458 |
temperature | t | K | 925hPa | 0.604 | 0.423 | 0.531 | 0.483 |
temperature | t | K | 950hPa | 0.612 | 0.441 | 0.546 | 0.513 |
temperature | t | K | 1000hPa | 0.582 | 0.419 | 0.576 | 0.504 |
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Summer statistics | Winter statistics | ||
Upper Bound STDV | Refined STDV | Upper Bound STDV | Refined STDV | ||||
geopotential | z | m²/s² | 50hPa | 42.17 | 17.29 | 41.54 | 17.03 |
geopotential | z | m²/s² | 100hPa | 41.29 | 16.93 | 40.39 | 16.56 |
geopotential | z | m²/s² | 150hPa | 39.06 | 16.01 | 38.22 | 15.67 |
geopotential | z | m²/s² | 200hPa | 30.34 | 12.44 | 30.41 | 12.47 |
geopotential | z | m²/s² | 250hPa | 29.65 | 12.16 | 29.87 | 12.25 |
geopotential | z | m²/s² | 300hPa | 31.28 | 12.82 | 30.92 | 12.68 |
geopotential | z | m²/s² | 400hPa | 29.68 | 12.17 | 29.56 | 12.12 |
geopotential | z | m²/s² | 500hPa | 27.60 | 11.32 | 27.56 | 11.30 |
geopotential | z | m²/s² | 600hPa | 26.57 | 13.02 | 26.47 | 14.03 |
geopotential | z | m²/s² | 700hPa | 25.99 | 14.29 | 25.93 | 17.11 |
geopotential | z | m²/s² | 800hPa | 25.40 | 15.75 | 25.48 | 19.36 |
geopotential | z | m²/s² | 850hPa | 25.13 | 16.33 | 25.68 | 21.06 |
geopotential | z | m²/s² | 900hPa | 25.40 | 17.53 | 26.80 | 23.58 |
geopotential | z | m²/s² | 925hPa | 25.89 | 18.12 | 27.77 | 25.27 |
geopotential | z | m²/s² | 950hPa | 26.74 | 19.25 | 29.02 | 27.28 |
geopotential | z | m²/s² | 1000hPa | 29.47 | 21.22 | 32.43 | 30.48 |
Table 12: Climatological analysis error standard deviation relative humidity for the CARRA-East domain. Anchor table12 table12
Name | Short name | Unit | Level | Summer statistics | Winter statistics | ||
Upper Bound STDV | Refined STDV | Upper Bound STDV | Refined STDV | ||||
relative humidity | r | % | 50hPa | 0.024 | 0.010 | 1.256 | 0.515 |
relative humidity | r | % | 100hPa | 0.010 | 0.004 | 0.194 | 0.080 |
relative humidity | r | % | 150hPa | 0.040 | 0.016 | 0.425 | 0.174 |
relative humidity | r | % | 200hPa | 0.534 | 0.219 | 1.913 | 0.784 |
relative humidity | r | % | 250hPa | 3.223 | 1.321 | 4.900 | 2.009 |
relative humidity | r | % | 300hPa | 6.153 | 2.523 | 8.861 | 3.633 |
relative humidity | r | % | 400hPa | 8.236 | 3.377 | 12.634 | 5.180 |
relative humidity | r | % | 500hPa | 7.181 | 2.944 | 12.562 | 5.150 |
relative humidity | r | % | 600hPa | 6.819 | 3.341 | 11.824 | 6.267 |
relative humidity | r | % | 700hPa | 6.740 | 3.707 | 11.522 | 7.605 |
relative humidity | r | % | 800hPa | 6.544 | 4.057 | 9.638 | 7.325 |
relative humidity | r | % | 850hPa | 6.118 | 3.977 | 8.194 | 6.719 |
relative humidity | r | % | 900hPa | 5.215 | 3.598 | 6.820 | 6.002 |
relative humidity | r | % | 925hPa | 4.674 | 3.272 | 6.303 | 5.736 |
relative humidity | r | % | 950hPa | 4.285 | 3.085 | 5.767 | 5.421 |
relative humidity | r | % | 1000hPa | 3.477 | 2.503 | 4.277 | 4.020 |
Uncertainties for the CARRA-West model domain
In this section we provide the actual estimates of the uncertainty in the form of the Upper Bound Standard Deviation and the Refined Standard Deviation for analysis errors over CARRA-West domain. For this domain, we provide all-year-around averaged statistics. The uncertainty estimation is given both in form of tables for surface pressure and mean sea level pressure (Table 13), u-wind component (Table 14), v-wind component (Table 15), temperature (Table 16) and specific humidity (Table 17), and in form of figures for the u- and v- wind components (Figure 7, left and right panel, respectively), temperature (Figure 8 , left panel) and specific humidity (Figure 8, right panel). All tables contain the name of the meteorological variable, short abbreviation, units and vertical levels. Uncertainty of the surface pressure field is provided at 0 m above ground and at 0 m above mean sea level. The uncertainty of the u- and v- components of wind, temperature and specific humidity is provided at standard pressure levels : 50 hPa, 100 hPa, 150 hPa, 200 hPa, 250 hPa, 300 hPa, 400 hPa, 500 hPa, 600hPa, 700 hPa, 800 hPa, 850 hPa, 900 hPa, 925 hPa, 950 hPa and 1000 hPa. The figures show uncertainty estimates as a function of standard pressure levels. Note that the uncertainty estimates for the CARRA-West domain are not given for exactly the same variables as for the CARRA-East domain due to technical processing challenges. The Upper Bound Standard Deviation of the analysis errors is shown as a solid curve and the Refined standard deviation (STDV) of the analysis error is shown as a dashed curve.
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Upper Bound STDV | Refined STDV |
Surface pressure | sp | Pa | 0m above ground | 26.11 | 7.04 |
Mean sea level pressure | msl | Pa | 0m above sea level | 39.44 | 10.64 |
Anchor | ||||
---|---|---|---|---|
|
Figure 7: Climatological analysis error standard deviation for u-component of wind (left plot) and for v-component of wind (right plot) as function of standard vertical pressure levels for the CARRA-West domain.
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Upper Bound STDV | Refined STDV |
u-component of wind | u | m/s | 50hPa | 0.540 | 0.146 |
u-component of wind | u | m/s | 100hPa | 0.370 | 0.100 |
u-component of wind | u | m/s | 150hPa | 0.447 | 0.121 |
u-component of wind | u | m/s | 200hPa | 0.611 | 0.165 |
u-component of wind | u | m/s | 250hPa | 0.732 | 0.198 |
u-component of wind | u | m/s | 300hPa | 0.871 | 0.235 |
u-component of wind | u | m/s | 400hPa | 0.869 | 0.235 |
u-component of wind | u | m/s | 500hPa | 0.862 | 0.233 |
u-component of wind | u | m/s | 600hPa | 0.887 | 0.240 |
u-component of wind | u | m/s | 700hPa | 0.893 | 0.241 |
u-component of wind | u | m/s | 800hPa | 0.856 | 0.231 |
u-component of wind | u | m/s | 850hPa | 0.831 | 0.224 |
u-component of wind | u | m/s | 900hPa | 0.810 | 0.219 |
u-component of wind | u | m/s | 925hPa | 0.803 | 0.217 |
u-component of wind | u | m/s | 950hPa | 0.799 | 0.216 |
u-component of wind | u | m/s | 1000hPa | 0.773 | 0.209 |
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Upper Bound STDV | Refined STDV |
v-component of wind | v | m/s | 50hPa | 0.535 | 0.144 |
v-component of wind | v | m/s | 100hPa | 0.360 | 0.097 |
v-component of wind | v | m/s | 150hPa | 0.449 | 0.135 |
v-component of wind | v | m/s | 200hPa | 0.603 | 0.162 |
v-component of wind | v | m/s | 250hPa | 0.714 | 0.193 |
v-component of wind | v | m/s | 300hPa | 0.851 | 0.230 |
v-component of wind | v | m/s | 400hPa | 0.853 | 0.230 |
v-component of wind | v | m/s | 500hPa | 0.852 | 0.230 |
v-component of wind | v | m/s | 600hPa | 0.877 | 0.237 |
v-component of wind | v | m/s | 700hPa | 0.884 | 0.239 |
v-component of wind | v | m/s | 800hPa | 0.852 | 0.230 |
v-component of wind | v | m/s | 850hPa | 0.826 | 0.223 |
v-component of wind | v | m/s | 900hPa | 0.801 | 0.216 |
v-component of wind | v | m/s | 925hPa | 0.794 | 0.214 |
v-component of wind | v | m/s | 950hPa | 0.792 | 0.214 |
v-component of wind | v | m/s | 1000hPa | 0.763 | 0.206 |
Anchor | ||||
---|---|---|---|---|
|
Figure 8: Climatological analysis error standard deviation for temperature (left plot) and specific humidity (right plot) as a function of standard pressure levels for the CARRA-West domain.
Anchor | ||||
---|---|---|---|---|
|
Name | Short name | Unit | Level | Upper Bound STDV | Refined STDV |
temperature | t | K | 50hPa | 0.130 | 0.035 |
temperature | t | K | 100hPa | 0.112 | 0.030 |
temperature | t | K | 150hPa | 0.169 | 0.045 |
temperature | t | K | 200hPa | 0.276 | 0.075 |
temperature | t | K | 250hPa | 0.285 | 0.077 |
temperature | t | K | 300hPa | 0.263 | 0.071 |
temperature | t | K | 400hPa | 0.227 | 0.061 |
temperature | t | K | 500hPa | 0.235 | 0.063 |
temperature | t | K | 600hPa | 0.259 | 0.070 |
temperature | t | K | 700hPa | 0.350 | 0.095 |
temperature | t | K | 800hPa | 0.473 | 0.128 |
temperature | t | K | 850hPa | 0.500 | 0.135 |
temperature | t | K | 900hPa | 0.508 | 0.137 |
temperature | t | K | 925hPa | 0.509 | 0.137 |
temperature | t | K | 950hPa | 0.520 | 0.140 |
temperature | t | K | 1000hPa | 0.542 | 0.146 |
Table 17: Climatological analysis error standard deviation for specific humidity for the CARRA-West domain. Anchor table17 table17
Name | Short name | Unit | Level | Upper Bound STDV | Refined STDV |
specific humidity | r | g/kg | 50hPa | 0. | 0. |
specific humidity | r | g/kg | 100hPa | 0. | 0. |
specific humidity | r | g/kg | 150hPa | 0. | 0. |
specific humidity | r | g/kg | 200hPa | 0.001 | 0.0002 |
specific humidity | r | g/kg | 250hPa | 0.002 | 0.0004 |
specific humidity | r | g/kg | 300hPa | 0.006 | 0.001 |
specific humidity | r | g/kg | 400hPa | 0.020 | 0.005 |
specific humidity | r | g/kg | 500hPa | 0.048 | 0.013 |
specific humidity | r | g/kg | 600hPa | 0.082 | 0.022 |
specific humidity | r | g/kg | 700hPa | 0.113 | 0.031 |
specific humidity | r | g/kg | 800hPa | 0.127 | 0.034 |
specific humidity | r | g/kg | 850hPa | 0.126 | 0.034 |
specific humidity | r | g/kg | 900hPa | 0.119 | 0.032 |
specific humidity | r | g/kg | 925hPa | 0.115 | 0.031 |
specific humidity | r | g/kg | 950hPa | 0.112 | 0.030 |
specific humidity | r | g/kg | 1000hPa | 0.099 | 0.027 |
Questions & Answers on field uncertainty estimates
For information and "questions and answers" (Q&A) on the uncertainty estimates for the ERA5 host reanalysis used on the lateral boundaries of this Arctic Regional Reanalysis, see this link:
ERA5: uncertainty estimation
...
For advanced users needing more detailed uncertainty information, uncertainties from ERA5 for the actual variable, time and position, might provide an additional indication (upper bound as concerns the near-surface variables) of the actual uncertainty of CARRA.
Uncertainties relative to observations: Verification statistics
Any observation or analysis of a meteorological variable is uncertain. For an observation, the uncertainty may be assessed by comparing measurements performed with similar instruments. In average, these instruments should observe the same values, but they will not get exactly the same results. Often the distribution of such random and independent deviations relative to the average measurement is a Gaussian distribution. This is defined by having a mean value 𝛍 and a standard deviation 𝛔 (STDV) as shown in Figure 9. The numbers on the distribution shown here are arbitrary and are just meant as an illustration.
...
For further details see the test and verification report for CARRA: https://datastore.copernicus-climate.eu/documents/reanalysis-carra/CARRATestVerificationFinal.pdf and the full system documentation: https://datastore.copernicus-climate.eu/documents/reanalysis-carra/CARRAFullSystemDocumentationFinal.pdf
Known Issues
Since the regional reanalysis is run nested into the ERA5 global reanalysis, it is affected by the known issues of ERA5. In addition to those issues, we have found that ERA5 uses incorrect glacier masks for most of the glaciers in the regional Arctic reanalysis domain, and the glaciers in ERA5 always have an analysis albedo of 0.85. This is wrong, since for instance exposed glacier ice albedos during summer are unaccounted for. These areas affect the general circulation and thermodynamic state in ERA5 and can affect the quality of the Arctic regional reanalysis.
Additionally, the Arctic reanalysis has the following known issues:
Annex: Brief outline of the Arctic reanalysis system
The details about the Numerical Weather Prediction model system used can be found in the full system documentation: https://datastore.copernicus-climate.eu/documents/reanalysis-carra/CARRAFullSystemDocumentationFinal.pdf
To supplement this we provide in the following a short description of the service as well as a description of basic principles of data assimilation in numerical weather prediction aimed at a non-specialist audience.
The service
The C3S_322_Lot2 contract of the Copernicus Climate Change Service produces and delivers a regional reanalysis (RRA) for the Arctic including long-term datasets of Essential Climate Variables (ECVs) for the 24 year period from 1997 to 2021. The model domains are shown in Figure 1. The modelling system has a resolution of 2.5x2.5 km2 and 65 vertical levels in the atmosphere. The produced datasets are freely available and can be used by anyone who wishes detailed long-term atmospheric data, for instance to study typical ranges of meteorological variables, as a reference for climate model runs, or to investigate highly resolved changes in the ECVs during the reanalysis period. An extended model version for the entire Arctic (pan-Arctic) area is run for a period of 1 year as a proof-of-concept.
The reanalysis model is the weather forecasting model HARMONIE-AROME cy40h1.1.1, which has been enhanced with more Arctic input data, and more extensive surface and atmospheric data assimilation. Data assimilation is explained further in section 7.2. Additionally, the model formulation has been improved with a specific focus on processes essential in the Arctic.
Principles of reanalysis systems
Atmospheric reanalysis is a method to reconstruct the atmospheric states by using historical observations (in situ, surface and satellite remote sensing) together with a weather forecasting model. It mainly provides a physically and dynamically coherent description of the state of the atmosphere. Surface variables are included mainly to the extent that they affect the atmosphere. This synthesis is accomplished by assimilating the observational data into a meteorological model and thereby forcing the model to reproduce the observations as closely as possible and make a gridded dataset that covers also locations from which no observations are available. The advantage of a reanalysis is that it provides a spatially and temporally complete, and consistent record of the atmospheric state. This cannot be achieved with an observational dataset alone, since interpolation methods cannot adequately represent the non-linear variability of the atmosphere.
The main advantages of reanalysis dataset are:
...
The Arctic reanalysis system applies the so-called 3D variational data assimilation (3D-VAR) reanalysis method. The 3D-VAR method is depicted schematically in Figure 19. At fixed points in time the model state is adjusted based on the observed state, taking into account the statistics of model and observation errors. The Arctic reanalysis system is run with eight cycles per day performing analyses at 00 UTC, 03 UTC, 06 UTC, 09 UTC, 12 UTC, 15 UTC, 18 UTC and 21 UTC. The forecast lengths vary between 3 and 30 hours.
Info | ||
---|---|---|
| ||
This document has been produced in the context of the Copernicus Climate Change Service (C3S). The activities leading to these results have been contracted by the European Centre for Medium-Range Weather Forecasts, operator of C3S on behalf of the European Union (Delegation agreementAgreement signed on 11/11/2014 and Contribution Agreement signed on 22/07/2021). All information in this document is provided "as is" and no guarantee or warranty is given that the information is fit for any particular purpose. The users thereof use the information at their sole risk and liability. For the avoidance of all doubt , the European Commission and the European Centre for Medium - Range Weather Forecasts have no liability in respect of this document, which is merely representing the author's view. |
Definition of the 65 vertical layer structure in HARMONIE
The CARRA vertical coordinate system is a terrain-following hybrid vertical coordinate, which means that it is terrain following at the bottom and pressure based on the top of the atmosphere. It has the advantage to describe the surface terrain properly, but also benefitting the advantage of having the pressure coordinate at the top of the atmosphere.
CARRA uses 65 model levels (level 65 is the surface and level 1 is the top of the atmosphere), which is further splitted into the so called half levels. CARRA has 66 half levels and the pressure of each half level can be obtained by the following formula:
P (k+1/2) = A (k+1/2) + B (k+1/2) * Ps
where k=0.... 65, Ps is surface pressure and the A and B coefficients (see below) valid at each half level.
The full model level pressure [1, 2, ...65] is defined as the mean of the pressure of each pair of neighbouring half levels [0,5, 1.5, ...65.5]. The model variables are defined in the full model levels.
The A and B coefficients are listed below (from the top to the bottom).
'AHALF'=>'0.00000000, 2000.00000000, 4000.21287319, 6002.09662113, 7911.25838577, 9633.01049417,11169.37146237, 12522.57753978, 13695.00149653, 14689.11546998, 15507.49052823, 16154.69697732, 16632.12471208, 16940.14949960, 17082.34869816, 17065.28164099, 16898.18367797, 16592.58939571, 16161.90395878, 15620.94340550, 14985.46502362, 14271.70773051, 13495.95994372, 12674.16909910, 11821.60314859, 10952.57042620, 10080.20053763, 9216.28565403, 8371.17893039, 7553.74479607, 6771.35457397, 6029.92021691, 5333.95880836, 4686.68074804, 4090.09511346, 3545.12645110, 3051.73811264, 2609.05813936, 2215.50455766, 1868.90774223, 1566.62821060, 1305.66882073, 1081.85503306, 890.47596795, 727.74548529, 590.17748096, 474.58767980, 378.08857614, 298.07947335, 232.23312781, 178.48015386, 134.99207440, 100.16369201, 72.59529482, 51.07508967, 34.56216490, 22.17022046, 13.15225964, 6.88641310, 2.86306141, 0.67344356, 0.00000000, 0.00000000, 0.00000000, 0.00000000, 0.00000000
'BHALF'=>'0.00000000, 0.00000000, 0.00000000, 0.00000000, 0.00095468, 0.00382570, 0.00862327, 0.01535782, 0.02404046, 0.03468314, 0.04729839, 0.06195102, 0.07868187, 0.09744325, 0.11815586, 0.14071098, 0.16497348, 0.19078554, 0.21797086, 0.24633925, 0.27569119, 0.30582244, 0.33652825, 0.36760726, 0.39886479, 0.43011564, 0.46118624, 0.49191624, 0.52215946, 0.55178443, 0.58067442, 0.60872709, 0.63585388, 0.66197911, 0.68703898, 0.71098036, 0.73375964, 0.5534143, 0.77569737, 0.79480486, 0.81264598, 0.82920633, 0.84454000, 0.85875505, 0.87191802, 0.88409276, 0.89534045, 0.90571965, 0.91528643, 0.92409452, 0.93219549, 0.93963895, 0.94647277, 0.95274328, 0.95849551, 0.96377340, 0.96862008, 0.97307803, 0.97718944, 0.98099640, 0.98454132, 0.98786727, 0.99102462, 0.99406510, 0.99703923, 1.00000000
Related articles
Content by Label | ||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
|