...
Info | ||||
---|---|---|---|---|
| ||||
|
Easy Heading Macro |
---|
History of modifications to this User Guide
...
Expand | ||||||||||||||||||||||||||||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||
|
...
In this section we will briefly mention the main features of the system. In general, CERRA consists of three components: CERRA (i) CERRA and CERRA-EDA, three-dimensional reanalysis systems, and (ii) CERRA-Land, a two-dimensional surface reanalysis system.
...
The CERRA system is based on the HARMONIE-ALADIN data assimilation system which is has been developed and used within the HIRLAM and ALADIN consortia. The CERRA system ACCORD consortium. It is implemented and optimized for the entire European area with surrounding sea areas (see Fig. 2) with a horizontal resolution of 5.5 km and 106 vertical levels. The system uses lateral boundaries conditions obtained from the ERA5 global reanalysis (Fig. 5). Also, the large scales in the regional system are constrained by data from the global reanalysis. The increase of resolution from the global reanalysis (RA) as well as the precursor regional reanalysis (RRA) is depicted in Figure 6.
...
Figure 5: Three different stages of RA: the global reanalyses from ERA5 are used as lateral boundary conditions for CERRA and CERRA-EDA reanalyses. Subsequently, short-forecasts forecast data from CERRA are used as background fields for the CERRA-Land surface reanalysis. As indicated by the vertical arrows, the amount of assimilated observations per area unit increases, in principle, from the global to the regional reanalysis as indicated by the arrows.
The CERRA system employs the 3D variational analysis (3D-VAR) method depicted schematically in Figure 4. At fixed points in time the model state is adjusted based on the observed state, taking into account the error statistics of both model and observations. The CERRA high-resolution system is has been running with eight assimilation cycles per day performing analyses at 00 UTC, 03 UTC, 06 UTC, 09 UTC, 12 UTC, 15 UTC, 18 UTC and 21 UTC. The forecasts lengths vary between 6 and 30 hours (see section 24.1.3 for more information) depending on the starting hour.
A flow-dependent background error covariance matrix (B) matrix is estimated using a 10-member ensemble of data assimilation (EDA) system, which is briefly described in the next section.
Anchor | ||||
---|---|---|---|---|
|
...
The main purpose of the CERRA-EDA system is to create an 'online', continuously serially updated flow-dependent B background error covariance matrix for use in the CERRA high-resolution CERRA system. CERRA-EDA is set up for the same geographical region as CERRA but at high-resolution, has the same number of vertical levels, but a horizontal resolution of 11 km. The EDA system utilises lateral boundary conditions from ERA5 and the same type of observations, including satellite observations , as CERRA high-resolution. CERRA-EDA is a ten member EDA system, where nine members differ through the perturbation of observations. But, one of the members is running with unperturbed observations forming the so called control member.
...
More information on the CERRA-EDA system and especially regarding the construction of the B background error covariance matrix can be found in El-Said et al. (2021).
...
The CERRA-Land analysis system uses the 2D-analysis system MESCAN and the land surface platform SURFEX to generate a coherent surface and soil analysis. The system combines CERRA forecast fields and additional surface observation (e.g. precipitation), to generate high-resolution (5.5 km) 2-dimensional analyses over Europe. MESCAN is a surface analysis system which uses an optimal interpolation algorithm for the analysis of 2-m 2m temperature and relative humidity and 24-h total accumulated precipitation (Soci et al., 2016). SURFEX is a land surface platform, which is driven by temperature, humidity, precipitation, wind and radiative fluxes (Bazile et al., 2017).
...
As depicted in Figure 3, all parameters are computed in grid points. Having a horizontal grid spacing of 5.5 km, as for the CERRA system, implies that each grid point characterizes values for an area of roughly 30 km2 (5.5km*5.5km). This needs to be considered when for instance CERRA data are compared with observations. The CERRA-EDA system has a horizontal resolution of 11km, meaning that each grid point characterizes values for an area of 121 km2 (11km*11km).
Vertical resolution
As already mentioned in Section 13.3.1, the CERRA system has 106 levels in the vertical levels (also called model levels) from the surface up to 1 hPa. However, only a very restricted number of parameters is stored on model levels. The main reason for that is the amount of needed storage space, when all parameters would be stored for all levels. Moreover, the vertical model grid is on hybrid-sigma coordinates, which do not match with any standard pressure level. Hence, this makes it quiet complex to use. Information about the model levels is here.
Therefore, the major part of the data is post-processed and stored on 29 selected pressure levels between 1000 - 1 hPa with a higher number of levels at lower altitudes. In addition, some parameters are also stored on 11 height levels between 15-500m. One reason to provide atmospheric variables on height levels is for applications in the wind energy sector. The exact levels both for pressure and height levels are given in Section 5 in the corresponding tables.
The CERRA system also contains a soil model which has 3 layers in the vertical levels. The three levels layers represent approximately the surface, the soil at root depth and the deep soil. Due to the used force-restore scheme in the soil model it is not possible to relate the levels layers with a certain depth in metre. Users interested in soil parameters should consider to use the data from the CERRA-Land system.
The CERRA-Land soil model has 14 layers in the vertical levels, which range from the surface to a depth of 12m. The edges between different levels are at 0.01m, 0.04m, 0.1m, 0.2m, 0.4m, 0.6m, 0.8m, 1.0m, 1.5m, 2m, 3m, 5m, 8m , and 12m. Values for a certain level reflect the mean value over the level layer thickness.
Time resolution
In general, data are stored with hourly resolution for the CERRA dataset. However, for all time steps, users have different options to select from and this is not always an easy choice. The preferred selection might vary for different parameters and the application of the user, respectively.
Figure 7 gives an overview on available analysis and forecast times. First, there are the eight analyses at 00 UTC, 03 UTC, 06 UTC, 09 UTC, 12 UTC, 15 UTC, 18 UTC and 21 UTC highlighted in red. Analysis data at these hours are assumed to be of higher quality than the forecasts valid at the same hours as they are in general closer to observations. They are available only every third hour and not all parameters are available for the analyses (the available analysis parameters are listed in Section 35). The forecast model is then started from the analysis and the output is saved hourly for the first six hours as indicated in dark blue in Fig. 7. Whereas the forecasts initiated at 00 UTC and 12 UTC continue until the forecast hour 30, all forecasts initiated at other hours (e.g. 15 UTC) stop after 6 hours. However, note that the output frequency of the forecasts initiated at 00 and 12 UTC is three hourly after the first six forecast hours (see blue boxes in Fig. 7) are completed. Due to the forecast lengths, the forecasts are overlapping and for every hour of the day data might be chosen from forecasts initiated at different forecasts hours and eventually from the analysis, respectively. For instance, at 12 UTC, the users can choose between the analysis and four different forecasts.
Anchor | ||||
---|---|---|---|---|
|
Figure 7: This table illustrates how different forecasts overlap and which options users have at a certain hour of the day. The availability of data is illustrated for the example date 2009/12/10. The color colour coding reflects analysis (red) and forecasts (blue). Moreover, different shades of blue correspond to frequency of the saved forecasts – hourly in dark blue and 3 hourly forecasts in blue.
...
Similar to above, results in complex terrain, such as mountainous regions or coastal areas, are generally less reliable than results over a more homogeneous terrain. The models cannot can hardly represent the strong local gradients that sometimes are caused by the variable (e.g. temperature gradients) that are usually driven by the complex terrain.
Figure 8 illustrates this behaviorbehaviour. Here, we show locations in Sweden having the best (blue) and the worst (red) correlations between the UERRA 2m - temperature and observational sites. A total of 853 measurement sites have been investigated. The Figure shows 50 locations each of highest and lowest correlation. Clearly, correlations are the lowest in the Swedish mountains and along the (east) coast.
Users need to be aware that the reanalysis provides gridded data where each grid point value describes an entire grid box area. That is in contrast to observations, which are usually point measurements. In case users need information with a higher horizontal resolution than provided by the CERRA systems, further downscaling (statistically or dynamically) needs to be considered. For instance, the correlations indicated in Fig. 8 increased when a linear interpolation from the four closest grid points to the observational site was applied than purely taken the values from the closest grid point.
...
Figure 8: A validation of the UERRA 2m - temperature data with Swedish observations. 50 places each with highest (blue) and lowest (red) correlation are shown out of 853 measurement sites included in the investigation.
...
Parameters on single level
Metadata for CERRA surface parameters | |
Horizontal coverage | The model domain spans from northern Africa beyond the northern tip of Scandinavia. In the west it ranges far into the Atlantic Ocean and in the east it reaches to the Ural Mountains. Herewith, it covers entire Europe. |
Horizontal resolution | 5.5 km x 5.5 km for CERRA high-resolution reanalysis 11 km x 11 km for CERRA ensemble members |
Vertical coverage | Each surface parameter is valid for one level in the vertical. There are four different (near) surface levels:
The cloud cover is provided for 3 atmospheric layers. |
Vertical resolution | single level |
Temporal coverage | 1984-09-01 00 UTC – 2021-06-30 21 UTC |
Temporal resolution | CERRA high-resolution reanalysis:
CERRA ensemble members:
|
Data type and format | Gridded data in GRIB2 |
Grid | Lambert conformal conic grid; 1069x1069 grid points for CERRA high-resolution reanalysis; 565x565 grid points for CERRA-EDA |
Table 1: Overview of the surface parameters Anchor table1 table1
Name | Unit | GRIB code | Analysis 3 hourly | Forecast | Height | |
1. | 10m wind speed | m/s | 207 | yes | yes | 10m |
2. | 10m wind gust since previous post-processing | m/s | 49 | - | yes | 10m |
3. | 10m wind direction | degree of true North | 260260 | yes | yes | 10m |
4. | 2m relative humidity | % | 260242 | yes | yes | 2m |
5. | 2m temperature | K | 167 | yes | yes | 2m |
6. | Albedo | % | 260509 | yes | yes | surface |
7. | Evaporation | kg/m2 | 260259 | - | yes | surface |
8. | Total column integrated water vapour | kg/m2 | 260057 | yes | yes | vertically integrated above the surface |
9. | Total precipitation | kg/m2 | 228228 | - | yes | surface |
10. | Maximum 2m temperature since previous post-processing | K | 201 | - | yes | 2m |
11. | Minimum 2m temperature since previous post-processing | K | 202 | - | yes | 2m |
12. | Skin temperature | K | 235 | yes | yes | surface |
13. | Surface latent heat flux | J/m2 | 147 | - | yes | surface |
14. | Surface sensible heat flux | J/m2 | 146 | - | yes | surface |
15. | Time-integrated surface direct short-wave radiation | J/m2 | 260264 | - | yes | surface |
16. | Surface net solar radiation | J/m2 | 176 | - | yes | surface |
17. | Surface solar radiation downwards | J/m2 | 169 | - | yes | surface |
18. | Surface net thermal radiation | J/m2 | 177 | - | yes | surface |
19. | Surface thermal radiation downwards | J/m2 | 175 | - | yes | surface |
20. | Surface net solar radiation, clear sky | J/m2 | 210 | - | yes | surface |
21. | Surface net thermal radiation, clear sky | J/m2 | 211 | - | yes | surface |
22. | Momentum flux at the surface u-component | N/m2 | 235017 | - | yes | surface |
23. | Momentum flux at the surface v-component | N/m2 | 235018 | - | yes | surface |
24. | Mean sea level pressure | Pa | 151 | yes | yes | surface |
25. | Surface pressure | Pa | 134 | yes | yes | surface |
26. | High cloud cover | % | 3075 | yes | yes | above 5000m |
27. | Low cloud cover | % | 3073 | yes | yes | surface-2500m |
28. | Medium cloud cover | % | 3074 | yes | yes | 2500m-5000m |
29. | Total cloud cover | % | 228164 | yes | yes | above ground |
30. | Snow density | kg/m3 | 33 | yes | yes | surface |
31. | Snow depth | m | 3066 | yes | yes | surface |
32. | Snow depth water equivalent | kg/m2 | 228141 | yes | yes | surface |
33. | Snowfall water equivalent | kg/m2 | 228144 | - | yes | surface |
34. | Land-sea mask | dimensionless | 172 | yes | - | surface |
35. | Orography | m2/s2 | 228002 | yes | - | surface |
36. | Surface roughness | m | 173 | yes | yes | surface |
37. | Soil temperature | K | 260360 | yes | yes | top layer of soil |
38. | Liquid Volumetric soil moisture (non-frozen) | m3/m3 | 260199 | yes | yes | top layer of soil |
39. | Volumetric soil moisture | m3/m3 | 260210 | yes | yes | top layer of soil |
10m wind speed
The 10-metre (10m) wind speed is the wind speed valid for the grid area determined for a height of 10m above the surface. The parameter is given in m/s. It is computed from both the zonal (u) and the meridional (v) wind components by wind speed=u2+v2. The 10m wind speed is available for the analysis and the forecast time steps. For the forecast, the value is instantaneous meaning that it is valid for the last time step of the integration at the issued time step.
...
Parameters on height levels
Metadata | |||||||
Horizontal coverage | The model domain spans from northern Africa beyond the northern tip of Scandinavia. In the west it ranges far into the Atlantic Ocean and in the east it reaches to the Ural Mountains. Herewith, it covers entire Europe. | ||||||
Horizontal resolution | 5.5 km x 5.5 km for CERRA high-resolution reanalysis 11 km x 11 km for CERRA ensemble members | ||||||
Vertical coverage | 11 height levels (from 15m up to 500m) | ||||||
Vertical levels | 15, 30, 50, 75, 100, 150, 200, 250, 300, 400 and 500m | ||||||
Temporal coverage | 1984-09-01 00 UTC – 2021-06-30 21 UTC | ||||||
Temporal resolution |
CERRA high-resolution reanalysis:
CERRA ensemble members:
| ||||||
Data type and format | Gridded data in GRIB2 | ||||||
Grid | Lambert conformal conic grid; 1069x1069 grid points for CERRA high-resolution reanalysis; 565x565 grid points for CERRA-EDA |
Table 2: Overview of the parameters on height levels Anchor table2 table2
Parameter | Unit | GRIB code | Analysis 3 hourly | forecast | |
1. | Wind speed | m/s | 10 | yes | yes |
2. | Wind direction | degree of true North | 3031 | yes | yes |
3. | Pressure | Pa | 54 | yes | yes |
4. | Relative humidity | % | 157 | yes | yes |
5. | Temperature | K | 130 | yes | yes |
6. | Specific cloud liquid water content | kg/kg | 246 | - | yes |
7. | Specific cloud ice water content | kg/kg | 247 | - | yes |
8. | Specific rain water content | kg/kg | 75 | - | yes |
9. | Specific snow water content | kg/kg | 76 | - | yes |
10. | Turbulent kinetic energy | J/kg | 260155 | - | yes |
Wind speed
Wind speed is the wind speed valid for the grid area determined for a certain height (15m-500m) above the surface. The parameter is given in m/s. It is computed from both the zonal (u) and the meridional (v) wind components by wind speed=u2+v2. The wind speed is available for the analysis and the forecast time steps. The value is instantaneous meaning that it is valid for the last time step of the integration at the issued time step.
...
Parameters on pressure levels
Metadata | |
Horizontal coverage | The model domain spans from northern Africa beyond the northern tip of Scandinavia. In the west it ranges far into the Atlantic Ocean and in the east it reaches to the Ural Mountains. Herewith, it covers entire Europe. |
Horizontal resolution | 5.5 km x 5.5 km for CERRA high-resolution reanalysis 11 km x 11 km for CERRA ensemble members |
Vertical coverage | From 1000 hPa to 1 hPa |
Vertical levels | 29 pressure levels (1000, 975, 950, 925, 900, 875, 850, 825, 800, 750, 700, 600, 500, 400, 300, 250, 200, 150, 100, 70, 50, 30, 20, 10, 7, 5, 3, 2, 1) |
Temporal coverage | 1984-09-01 00 UTC – 2021-06-30 21 UTC |
Temporal resolution | CERRA high-resolution reanalysis:
CERRA ensemble members:
|
Data type and format | Gridded data in GRIB2 |
Grid | Lambert conformal conic grid; 1069x1069 grid points for CERRA high-resolution reanalysis; 565x565 grid points for CERRA-EDA |
Table 3: Overview of the parameters on pressure levels Anchor table3 table3
Parameter | Unit | GRIB code | Analysis 3 hourly | forecast | |
1. | Cloud cover | % | 260257 | - | yes |
2. | Specific cloud liquid water content | kg/kg | 246 | - | yes |
3. | Specific cloud ice water content | kg/kg | 247 | - | yes |
4. | Specific rain water content | kg/kg | 75 | - | yes |
5. | Specific snow water content | kg/kg | 76 | - | yes |
6. | Turbulent kinetic energy | J/kg | 260155 | - | yes |
7. | Relative humidity | % | 157 | yes | yes |
8. | Temperature | K | 130 | yes | yes |
9. | U-component of wind | m/s | 131 | yes | yes |
10. | V-component of wind | m/s | 132 | yes | yes |
11. | Geopotential | m2/s2 | 129 | yes | yes |
Cloud cover
Cloud cover is the percentage of sky covert with clouds. It is valid for the grid at the corresponding height. The parameter is given in %. Total cloud cover is only available for the forecast time steps. The value is instantaneous meaning that it is valid for the last time step of the integration at the issued time step.
...
Table 4: Overview of parameters on model levels
Metadata | |
Horizontal coverage | The model domain spans from northern Africa beyond the northern tip of Scandinavia. In the west it ranges far into the Atlantic Ocean and in the east it reaches to the Ural Mountains. Herewith, it covers entire Europe. |
Horizontal resolution | 5.5 km x 5.5 km for CERRA high-resolution reanalysis 11 km x 11 km for CERRA ensemble members |
Vertical coverage | From approximately 10m (model level 106) above the surface to a height of 1 hPa (model level 1) |
Vertical levels | 106 hybrid atmospheric model levels (106, 105, 104 ... 3, 2, 1) |
Temporal coverage | 1984-09-01 00 UTC – 2021-06-30 21 UTC |
Temporal resolution | CERRA high-resolution reanalysis: 3-hourly analyses at 00, 03, 06, 09, 12, 15, 18 and 21 UTC CERRA ensemble members: 6-hourly analyses at 00, 06, 12 and 18 UTC Note: forecast data are not saved for the parameters on model levels |
Data type and format | Gridded data in GRIB2 |
Grid | Lambert conformal conic grid; 1069x1069 grid points for CERRA high-resolution reanalysis; 565x565 grid points for CERRA-EDA |
Table 4: Overview of parameters on model levels Anchor table4 table4
Parameter | Unit | GRIB code | Analysis 3 hourly | forecast | |
1. | Specific humidity | kg/kg | 133 | yes | - |
2. | Temperature | K | 130 | yes | - |
3. | U-velocity | m/s | 131 | yes | - |
4. | V-velocity | m/s | 132 | yes | - |
Specific humidity
The specific humidity is the mass of water vapour per unit mass of air valid for the grid area at the corresponding model level. Only analyses are stored for parameters on model levels. The parameter is given in kg/kg.
...
Below are the essential parameters describing the grid and the used Lambert Conformal Conic projection. More information about the grid and coordinates can be found in the FAQ.
Number of points along x-axis: 1069
Number of points along y-axis: 1069
X-direction grid length: 5500 m
Y-direction grid length: 5500 m
Projection: Lambert Conformal Conic
Central meridian: 8
Standard parallel 1: 50
Standard parallel 2: 50
Latitude of origin: 50
Earth assumed spherical with radius: 6371229 m
Latitude and longitude of the corner grid points in decimal degrees | ||
Grid point | Latitude | Longitude |
Upper-left | 63.7695 | -58.1051 |
Upper-right | 63.7695 | 74.1051 |
Lower-right | 20.2923 | 33.4859 |
Lower-left | 20.2923 | -17.4859 |
CERRA-EDA
CERRA-EDA comprises the same domain as CERRA and has exactly the same set of parameters as the high-resolution CERRA dataset. It differs only in the horizontal resolution, which is 11km as well as in the number of available time steps. CERRA-EDA has four analyses per day, at 00, 06, 12 and 18 UTC. Starting from the analyses, forecasts are run for six hours. Forecast fields are saved with hourly resolution.
Metadata for CERRA-EDA parameters | |
Horizontal coverage | Same as CERRA. |
Horizontal resolution | 11 km x 11 km |
Vertical coverage | Same as for CERRA parameters. |
Vertical resolution | Same as for CERRA-parameters. |
Temporal coverage | Same as for CERRA. |
Temporal resolution | Analyses are available at 00, 06, 12, and 18 UTC. |
Data type and format | Gridded data in GRIB2 |
Grid | Lambert conformal conic grid with 565x565 grid points |
The CERRA-EDA grid description
Below are the essential parameters describing the CERRA-EDA grid and the used Lambert Conformal Conic projection. More information about the grid and coordinates can be found in the FAQ in section 2.3.
Number of points along x-axis: 565
Number of points along y-axis: 565
X-direction grid length: 11000 m
Y-direction grid length: 11000 m
Projection: Lambert Conformal Conic
Central meridian: 8
Standard parallel 1: 48
Standard parallel 2: 48
Latitude of origin: 48
Earth assumed spherical with radius: 6371229 m
Latitude and longitude of the corner grid points in decimal degrees | ||
Grid point | Latitude | Longitude |
Upper-left | 63.4028 | -60.4047 |
Upper-right | 63.4028 | 76.4047 |
Lower-right | 17.6121 | 34.3203 |
Lower-left | 17.6121 | -18.3203 |
Known issues
Wrong metadata for the 2m maximum and minimum temperature
The metadata for the forecast step range is incorrect in the GRIB2 files for the maximum 2m temperature since previous post-processing and the minimum 2m temperature since previous post-processing. The correct step range is given in the table below.
Step range in the metadata | Correct step range |
---|---|
0-1 | 0-1 |
0-2 | 0-2 |
0-3 | 0-3 |
0-4 | 3-4 |
0-5 | 3-5 |
0-6 | 3-6 |
0-9 | 6-9 |
0-12 | 9-12 |
0-15 | 12-15 |
0-18 | 15-18 |
0-21 | 18-21 |
0-24 | 21-24 |
0-27 | 24-27 |
0-30 | 27-30 |
Minor data assimilation issues (issues when some observations were not available)
...
For completeness, a list of the occurrences of missing observational data and affected periods are listed below.
CERRA | |
---|---|
Affected period | Description |
2020-10-01 to 2021-06-30 | Fewer SYNOP data than usual. About 1850 instead of 2100 stations. |
2021-02-01 to 2021-03-31 | Only few AMV data assimilated in this period. |
2020-04-01 to 2021-03-31 2019-01-01 to 2019-04-09 | Only very few ocean buoy observations. |
2019-01-01 to 2019-03-09 | Data assimilation used a slightly degraded B-matrix. The climatological part of the B-matrix was shifted by two months, i.e. the November climatology was used instead of the January climatology. |
CERRA-EDA | |
---|---|
Affected period | Description |
2021-04-01 to 2021-04-30 | IASI missing for both METOP-A and METOP-B |
2019-01-01 to 2019-05-31 | No AMV included. Only very few ocean buoy observations. |
2019-01-01 to 2019-01-02 | No additional local observations included for Greenland, Iceland, Norway, Sweden, Finland, and France. |
2016-10-03 | The 18UTC cycle was ran without TEMP and PILOT data. |
1984-09-01 to present | No MSU data. |
References
- Bazile E, R. Abida, A. Verelle, P. Le Moigne and C. Szczypta (2017): MESCAN-SURFEX surface analysis, deliverable D2.8 of the UERRA project, http://www.uerra.eu/publications/deliverable-reports.html
- El-Said A., P. Brousseau, M. Ridal and R. Randriamampianina (2021): A new temporally flow-dependent EDA estimating background errors in the new Copernicus European Regional Re-Analysis (CERRA), Earth and Space Science Open Archive, pp. 28, doi 10.1002/essoar.10507207.1, https://doi.org/10.1002/essoar.10507207.1
- Niermann D. et al. (2017): Scientific report on assessment of regional analysis against independent data sets, deliverable D3.6 of the UERRA project, http://www.uerra.eu/publications/deliverable-reports.html
- Ridal M., S. Schimanke and S. Hopsch (2018): Documentation of the RRA system: UERRA (C3S deliverable D322_Lot1.1.1.2, Documenting the UERRA system)
- Soci C., E. Bazile, F. Besson and T. Landelius (2016). High-resolution precipitation re-analysis system for climatological purposes. Tellus A, Dynamic Meteorology and Oceanography, 68:1, DOI: 10.3402/tellusa.v68.29879
- Verver Gé (2017): User Guidance, deliverable D8.4 of the UERRA project, http://www.uerra.eu/publications/deliverable-reports.html
...