Global and regional reanalyses


Improved reanalyses of past weather can be obtained by retrospectively assimilating reprocessed observational datasets ranging from surface stations and satellites with a up-to-date Numerical Weather Prediction (NWP) model. The resulting time series of the atmospheric state is both dynamically consistent and close to observations. The interest in extracting climate information from reanalysis is rising and creating a request for reanalysis uncertainty estimation at various temporal-spatial scales.
These research questions have been addressed in EU-funded research projects (e.g.ERA-CLIM, EURO4M and UERRA). Regional reanalyses are now available for Europe and specific sub-domains, e.g. produced by national meteorological services. Global and regional reanalyses are also an important element of the Copernicus Climate Change Services.

This session invites papers that:
• Explore and demonstrate the capability of global and regional reanalysis data for climate applications
• Compare different reanalysis (global, regional) with each other and/or observations
• Improve recovery, quality control and uncertainty estimation of related observations
• Analyse the uncertainty budget of the reanalyses and relate to user applications

Convener: Frank Kaspar | Co-conveners: Eric Bazile, Jan Keller
Lightning talks
| Mon, 06 Sep, 09:00–12:30 (CEST)

Lightning talks: Mon, 6 Sep

Chairpersons: Frank Kaspar, Jan Keller
Overview of European activities in developing regional reanalyses
Andreas Hense and the The Bonn/Cologne group of Hans Ertel Center for Weather Research

In recent years, the reanalysis efforts have been focused on increasing the spatio-temporal resolution to better capture the variability of the atmospheric state between the mesoβ up to synoptic scales. An approach to gain a substantial increase in horizontal resolution are regional reanalyses which utilize a limited area model instead of a global one and thus focus on a specific region of the world. For the European continent several regional reanalyses are now available, e.g., HARMONIE (11km) or COSMO-REA6 (6km).

This work presents ongoing work of a variety of applications and evaluations made possible through the hourly 6km data of the COSMO-REA6 reanalysis product. 
Snow and precipitation analysis and events are going to represent better the smaller scales. Here we will show the representation of the scaling properties of the precipitation extremes and an evaluation of the 2005 extreme early winter snow event which happened in Western Germany.  With this reanalysis we also investigated the potential of season identification, finding that gridpoint-wise vertical temperature and wind information provides efficient ways to classify probabilistically a single day into a specific season.  Moreover, the calculation of the bio-climatic variables from COSMO-REA6 allowed the comparison with the Worldclim 2.0 dataset. The main advantage of the newly derived bio-climatic variables is the internal consistency between the temperature and precipitation-related variables. Furthermore, from the increment analysis, we investigate the soil moisture biases and how they relate to the missing representation of irrigation processes in the COSMO model. Results show that there is a similarity in the specific humidity increments and both the irrigation annual cycle and spatial distribution of irrigated areas. With physically consistent spatio-temporal fields it is possible to investigate human discomfort measurements and their relation to temperature indexes. Lastly, extended use of the COSMO-REA6 data was made for renewable energy applications both for wind turbines and photovoltaic cells. The examples given in this presentation will highlight the general potential of high-resolution reanalyses to be part of transdisciplinary problems related to climate and weather variability.  

How to cite: Hense, A. and the The Bonn/Cologne group of Hans Ertel Center for Weather Research: The  regional reanalysis COSMO-REA6 and its transdisciplinary applications, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-464, https://doi.org/10.5194/ems2021-464, 2021.

Semjon Schimanke, Ludvig Isaksson, Lisette Edvinsson, Martin Ridal, Lars Berggren, Susanna Hopsch, Adam El-Said, Michael Glinton, Eric Bazile, Patrick Le Moigne, Antoine Verrelle, Per Dahlgren, and Roger Randriamampianina

The Copernicus European regional reanalysis (https://climate.copernicus.eu/regional-reanalysis-europe) is produced as part of the Copernicus Climate Change Service (C3S). The presentation will introduce the service and its main objectives as well as it will give and overview of available data. Data quality will be demonstrated by comparison with ERA5 and other gridded datasets.

In the first phase of the service, systems inherited from the FP7 project UERRA (Uncertainties in Ensembles of Regional ReAnalyses, http://www.uerra.eu) were applied extending the UERRA-HARMONIE as well as the MESCAN-SURFEX datasets. These datasets contain analyses of the atmosphere, the surface and the soil. UERRA-HARMONIE is a full model system including a 3D-Var data assimilation scheme for upper air observations and an OI-scheme for surface observations. MESCAN-SURFEX is a complementary 2D surface analysis system interfaced to a land surface model. Data is available for entire Europe at a horizontal resolution of 11 km for UERRA-HARMONIE and at 5.5 km for MESCAN-SURFEX. The systems provide four analyses per day – at 0 UTC, 6 UTC, 12 UTC, and 18 UTC. Between the analyses ranges, forecasts of the systems are available with hourly resolution. More than fifty parameters are available on various level types. Data are available for the period 1961 – July 2019 through Copernicus Climate Data Store (CDS).

In spring 2020, the service started the production of the next generation regional reanalysis. The successor comprises three components:
- CERRA (5.5 km horizontal resolution)
- CERRA-EDA (10-member ensemble at 11 km resolution)
- CERRA-Land (5.5 km horizontal resolution)

In addition to the higher resolution, CERRA is more sophisticated than UERRA. For instance, more observations are assimilated into CERRA, in particular remote sensing data. CERRA is produced with 3-hourly cycling and a flow depending part of the B-matrix is derived from CERRA-EDA. The production of CERRA, CERRA-EDA and CERRA-Land will complete in September/October 2021 and data will become available in the CDS shortly thereafter.

The quality of the regional reanalysis in comparison to ERA5 will be shown with results of the standard HARMONIE-verification package as well as based on certain case studies. For instance, the winter storm Gudrun (January 2005, southern Sweden) will be investigated.

How to cite: Schimanke, S., Isaksson, L., Edvinsson, L., Ridal, M., Berggren, L., Hopsch, S., El-Said, A., Glinton, M., Bazile, E., Le Moigne, P., Verrelle, A., Dahlgren, P., and Randriamampianina, R.: Copernicus European regional reanalysis, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-389, https://doi.org/10.5194/ems2021-389, 2021.

Antoine Verrelle, Michael Glinton, Eric Bazile, and Patrick Le Moigne

The Copernicus Climate Change Service of the European Commission aims to produce and deliver a land regional reanalysis for Europe covering the period from 1984 to the present at a resolution of 5.5 km.

Accurate spatial and temporal datasets of precipitation and surface variables are essential for water resource management and climate change studies. The need for higher spatial and temporal resolution of precipitation and surface variables is a demand that leads regional reanalysis.

CERRA-Land (Copernicus European Regional Re-Analysis) is a regional land surface reanalysis dataset which describes the evolution of soil moisture, soil temperature and snowpack. CERRA-Land is the result of a single stand alone integration of the SURFEX V8.1 land surface model driven by meteorological forcing from the CERRA atmospheric reanalysis and a daily accumulated surface precipitation analysis produced with the MESCAN system (Soci et al., 2016), which performs an optimal interpolation between CERRA forecast precipitation and in situ rain gauges.

No downscaling method was used to build up the forcing data because the system has the same grid as the CERRA atmospheric reanalysis. To better represent the permafrost active layer, the effect of soil organic carbon on hydraulic and thermal soil properties was taken into account in the northern part of the domain. To solve both heat and water transfer equations in the soil, a discretization of the soil into 14 layers was used. A 15 year spin-up of the model was required to reach an equilibrium.

The quality of CERRA-Land is assessed by comparisons to ground-based observations, such as of snow depth and turbulent, latent and sensible fluxes. A comparison with ERA5-Land will be done and the added value of this regional dataset will be discussed. The entire CERRA-Land dataset from 1984 to present is expected to be available for public release later in 2021.

How to cite: Verrelle, A., Glinton, M., Bazile, E., and Le Moigne, P.: CERRA-Land : A new land surface reanalysis at 5.5 km resolution over Europe, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-492, https://doi.org/10.5194/ems2021-492, 2021.

Sabrina Wahl, Clarissa Figura, and Jan D. Keller

Reanalysis is a procedure to merge numerical model integrations and observations to obtain a synergetic representation of the past climatological state of a system, e.g., of the atmosphere. An alternative to running a full reanalysis scheme is a so-called surface reanalysis. Here, an existing reanalysis is used as prior information (for the near-surface state). This first guess is then corrected in a data assimilation step preferrably by applying observations not used in the original assimilation. In such a scheme, an additional downscaling is often performed to enhance the spatial representation of the surface reanalysis.

We present here the development of a new approach aiming to establish such a data set based on the COSMO-REA6 regional reanalysis of the Hans-Ertel-Centre and Deutscher Wetterdienst (DWD). The data assimilation step is based on the operational Local Ensemble Transform Kalman Filter (LETKF) of DWD. While the data assimilation is often performed univariately in such surface reanalysis schemes, here we apply it to various parameters at once thus conserving the covariances among the parameters and allowing for a consistent multivariate utilization of the data. Further, this reanalysis will not be restricted to the ground level and near-surface parameters. Instead, it will be extended to the lower part of the boundary layer aiming at an improved representation of wind speeds in wind turbine hub heights especially relevant for renewable energy applications. The envisaged resolution is approximately 1km and therefore enables an enhanced representation of spatial variability and heterogeneity on small scales. In addition, the LETKF is an ensemble-based data assimilation scheme which also provides uncertainty estimates through an ensemble of the re-analyzed parameters which can also be used as input for downstream applications.

How to cite: Wahl, S., Figura, C., and Keller, J. D.: A novel approach to surface reanalysis, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-135, https://doi.org/10.5194/ems2021-135, 2021.

Adam El-Said, Pierre Brousseau, Roger Randriamampianina, and Martin Ridal

A new augmented Ensemble of Data Assimilations (EDA) technique, which estimates background error covariances (B-matrix), has been developed for the new Copernicus European Regional Re-Analysis (CERRA-EDA). CERRA-EDA has 10 members with two main pools of forecast differences: seasonal and daily. The seasonal component is pre-prepared (`offline') at reanalysis-resolution (5.5km). The new augmentation governs the time-dependent mixture of winter and summer differences of this seasonal component with respect to the time of year. The daily component is (`online') and averaged in moving succession over 2.5 days with subsequent B-matrix computation every 2 days. This daily component runs at 11km and the forecasts are interpolated to 5.5km prior to use. The seasonal-daily split is set to a fixed value of 80-20\% for CERRA production. The EDA is cycled 6-hourly while CERRA has a 3-hour analysis cycle. The B-matrix is modelled on a bi-Fourier limited area weather model, where dependence of vertical correlations on horizontal scale (non-separability), horizontal homogeneity and isotropy are assumed. The mass-wind and specific humidity fields are related via vorticity and geopotential and the relationships are estimated via multiple linear regressions enforcing simplified analogues of flow-dependence. 

We demonstrate the potential of CERRA-EDA to estimate rapid changes in weather regime change over Europe by assessing B-matrix statistics and forecast skill scores in a case study. The case study assesses two like-periods bearing different weather regimes, Mar-03 (blocking regime) and Mar-18 (NAO- regime). The aptitude of the B-matrix to reflect weather regime change is shown to be mostly dependent on the observation network in a given year. We also illustrate the impact of: change in observation networks over time, and varying the seasonal-daily split. This is shown through analysing the spatio-temporal evolution of background standard deviations. Finally, analysis and forecast skill scores up to 24-hours are also shown to offer improvements worth considering.

How to cite: El-Said, A., Brousseau, P., Randriamampianina, R., and Ridal, M.: B-matrix estimation in the Copernicus European Regional Re-Analysis (CERRA), EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-161, https://doi.org/10.5194/ems2021-161, 2021.

Stéphane Van Hyfte, Patrick Le Moigne, Eric Bazile, and Antoine Verrelle

Within the UERRA project, a daily precipitation reanalysis at a 5,5km resolution has been realized from 1961 to 2015. The reanalysis was obtained by the MESCAN analysis system which combines an a priori estimate of the atmosphere – called background – and observations using an optimum interpolation (OI) scheme. Such method requires the specification of observations and background errors. In general, constant standard deviation errors are used but more errors are made when high precipitation are observed. Then, to take this effect into account and to avoid a model over-estimation in case of light precipitation, a variable formula of the observation standard deviation error was purposed with a small value for null precipitation and greater values when precipitation are higher, following a linear equation.

Desroziers et al proposed a method to determine observations and background errors called a posteriori diagnosis. To use this iterative method, the analysis has to be ran several times until it converged. In this study, the a posteriori diagnosis is used per precipitation class to determine the observation standard deviation error formula. MESCAN was tested using the French operational model AROME at 1,3km resolution and the atmopsheric UERRA analysis downscaled to 5,5km background and combined to the French observational network over the 2016-2018 period. The observation standard deviation error formula obtained by the a posteriori diagnosis is then used in the MESCAN analysis system to produce precipitation analysis over the 2016-2018 period. Results are compared to UERRA precipitation reanalysis over independant observations by comparing bias, RMSE and scores per precipitation class.

How to cite: Van Hyfte, S., Le Moigne, P., Bazile, E., and Verrelle, A.: Estimation of the observation standard deviation error formula thanks to the a posteriori diagnosis, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-491, https://doi.org/10.5194/ems2021-491, 2021.

Marie Doutriaux-Boucher, Roger Huckle, Alessio Lattanzio, Olivier Sus, Jaap Onderwaater, Mike Grant, and Joerg Schulz

This presentation provides an overview of the different upper-air wind data records available at EUMETSAT for usage in global and regional reanalysis. The assimilation of Atmospheric Motion Vectors (AMV) is recognised to be important to reduce the forecast errors in NWP model runs. In support of the Copernicus Climate Change Service (C3S), EUMETSAT produced several AMV Climate Data Records (CDR) from geostationary and low-earth orbit satellites for assimilation into ECMWF’s next global reanalysis ERA6.

Since the launch of its first generation of geostationary satellites, EUMETSAT has developed its own unique algorithms to derive atmospheric motion vectors (AMVs). These algorithms are used to provide real time AMVs using images acquired from instruments on-board both polar and geostationary satellites. These AMVs are routinely assimilated into weather forecast models. EUMETSAT archived all image data from its instruments (MVIRI and SEVIRI) in geostationary orbit and the global record of Advanced Very High Resolution Radiometer (AVHRR) data back to the late 1970s providing a suitable data source for climate research allowing the production of consistent AMV CDRs over the entire period.

Two long AMV data records are available now from the geostationary sensors on Meteosat-2 to Meteosat-10 covering 1981-2017 over Africa and Europe and from AVHRR Global Area Coverage (GAC) data from 16 AVHRR instruments starting with the TIROS-N satellite and covering polar AMVs over the Northern and Southern hemisphere from 1978-2019. In addition, full resolution AVHRR images (Local Area Coverage (LAC)) from the AVHRR aboard the polar orbiting Metop-A and -B satellites were used to generate a CDR containing polar AMVs from single satellite retrievals and global AMVs from the combined Metop-A/B dual satellite retrieval starting in 2007 and 2013, respectively.

For all data records, the EUMETSAT AMV algorithm adapted for climate purposes was used and extensive validation of the data records were performed. It shows that the CDR are homogeneous and very stable over the period. They are suitable for usage in model reanalysis and climate analysis. The CDR are in agreement with ground based radiosonde and model data. For the polar AMVs, a remarkable agreement with MODIS AMVs has been found.

To better serve closer to real time needs for reanalysis, EUMETSAT is experimenting with the continuous production of an Interim Climate Data Record (ICDR) with a timeliness close to real-time. With a still not completely operational low-cost approach, a timeliness of 83% within 18 hours at similar quality was achieved.

In addition to the existing data records the presentation provides the plan for future improvements and new CDR releases for AMV data records in the coming years. In particular, the use of better information on multi-layer cloud objects in AMV retrievals is a central part for the improvements of the AMVs from geostationary orbit.  

How to cite: Doutriaux-Boucher, M., Huckle, R., Lattanzio, A., Sus, O., Onderwaater, J., Grant, M., and Schulz, J.: Climate data record of atmospheric motion vectors at EUMETSAT for usage in reanalysis, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-445, https://doi.org/10.5194/ems2021-445, 2021.

Chairpersons: Eric Bazile, Frank Kaspar
Evaluation and applications of reanalysis datasets.
María Ortega, Enrique Sánchez, Claudia Gutiérrez, and María Ofelia Molina

Regional winds are caused by small-scale pressure differences in a way that important air flows can arise in a very small and specific region. Sometimes an orographic feature, such as a channel like the Ebro Valley or the Strait of Gibraltar, lead the wind, due to mass conservation, to acquire a certain specific range of directions and considerable speed. In the regions where they are observed, the wind is of great importance not only for the climatology and meteorology of these areas but also for their culture and identity. However, it is difficult to analyze them using the most common reanalysis products, since their spatial resolutions are not high enough to properly describe the orographic characteristics that lead to the regional winds in specific locations. Here, we will explore the application of the COSMO-REA6 high resolution reanalysis system for the assessment of the main regional winds in the Iberian Peninsula: the cierzo wind in the Ebro Valley and the levante and poniente winds in the Strait of Gibraltar, for the 2000-2018 period. COSMO-REA6 uses a spatial resolution of 6 km (0.055º), which is much larger than previous reanalysis and regional modelling databases, so it can better capture the orography of the areas and therefore the regional winds we intend to study. The cierzo, levante and poniente winds are very relevant in the Iberian Peninsula due to their intensity and their frequency. Defined with a 5 m/s threshold for each hour, and their specific direction range, around 95, 150 and 110 wind days per year are obtained, respectively. Their study may also be important for other reasons, such as the production of renewable energy in these areas. First, we conduct a preliminary assessment of wind speed and direction with hourly data from weather stations, which have been obtained from the HadISD global sub-daily dataset. Then, we compare data from stations with COSMO-REA6 reanalysis in each location and produce a spatial description of the reanalysis in the Peninsula. We also study the atmospheric patterns associated with the regional winds characterized above. Due to the few studies that have been carried out on regional winds in the Iberian Peninsula, these results can be of great interest for various fields, such as meteorology, climatology and the generation of renewable energy.

How to cite: Ortega, M., Sánchez, E., Gutiérrez, C., and Molina, M. O.: Regional winds over the Iberian Peninsula with COSMO-REA6 high resolution regional reanalysis: cierzo, levante and poniente, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-140, https://doi.org/10.5194/ems2021-140, 2021.

Thomas Spangehl, Michael Borsche, Deborah Niermann, Frank Kaspar, and Birger Tinz

The exploitation of offshore wind energy is an essential part of the German energy transition (Energiewende). The planning of new offshore wind farms demands detailed information on wind conditions at turbine hub heights in the North Sea and Baltic Sea. High-resolution reanalyses which are based on state-of-the-art numerical weather prediction (NWP) models combined with data assimilation systems offer the required meteorological data which are suitable for climatological assessment.

The regional reanalysis COSMO-REA6 operated by Germany’s national meteorological service (Deutscher Wetterdienst, DWD) provides hourly data of 6 km horizontal resolution for 1995-2019/08 (Kaspar et al., 2020). Moreover, hourly data of 31 km horizontal resolution for 1950 to present are available from the global reanalysis ERA5 produced by the European Centre for Medium-Range Weather Forecasts (ECMWF). DWD delivers reanalysis data and statistical evaluation results to Bundesamt für Seeschifffahrt und Hydrographie (BSH) in order to facilitate offshore site tenders. Data and a report were recently published as part of the tenders for 2021 (https://pinta.bsh.de/).

Here we present an evaluation of the 100 m wind speed and direction from COSMO-REA6 and ERA5 based on a comprehensive statistical analysis. On the reference side the FINO measurements (Research platforms in the North Sea and Baltic Sea, https://www.fino-offshore.de/en/index.html) from FINO1 and FINO2 are used. The FINO measurements are not used by the data assimilation schemes of the two reanalyses and therefore constitute independent reference data. The focus is on episodes prior to the installation of wind farms in the direct vicinity of the FINO platforms to avoid wake effects. The quality of the two reanalyses is compared to other state-of-the-art reanalyses and wind atlas data.


Kaspar et al. (2020): Regional atmospheric reanalysis activities at Deutscher Wetterdienst: review of evaluation results and application examples with a focus on renewable energy, Adv. Sci. Res., 17, 115–128, https://doi.org/10.5194/asr-17-115-2020.

How to cite: Spangehl, T., Borsche, M., Niermann, D., Kaspar, F., and Tinz, B.: Evaluation of high-resolution reanalyses for wind energy application in the North Sea and Baltic Sea, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-241, https://doi.org/10.5194/ems2021-241, 2021.

Sebastian Brune, Jan D. Keller, and Sabrina Wahl

The correct spatio-temporal representation of wind speed is of large interest for the wind energy sector. Therefore, this study compares wind measurements in different heights from several locations in Central Europe with two global (ERA5, MERRA-2) and one regional reanalysis (COSMO-REA6). Employing a two-parameter Weibull distribution, the shape and scale parameters as well as mean, standard deviation and RMSE are investigated at and around common wind turbine hub height. We find that COSMO-REA6 best describes wind fields closer to the surface possibly due to its high horizontal resolution. Here, it also exhibits a good alignment with the diurnal cycle. However, for common wind turbine hub heights and above, ERA5 outperforms the other two reanalyses possibly due to its higher vertical resolution. MERRA-2 overestimates wind speed in the lower boundary layer at nearly all sites.

In the next step, a diagnostic and mass-consistent wind model is applied to the COSMO-REA6 wind field. The resolution of the wind field will be increased by a factor of 8 from originally 6 km to approximately 800 m. In addition to the vertical stability of the lower atmosphere, the orography on the finer grid and the corresponding effects are taken into account. We expect that especially in complex terrain the wind field will be corrected and thus should fit better to the observations. Channeling effects, shadowing and increased wind speed in exposed locations can be better represented. The new high-resolution wind field forms the basis for a statistical wind model to obtain post-processed wind estimates in the lower boundary layer. This approach will utilize generalized linear model and/or an artificial neural network techniques.

How to cite: Brune, S., Keller, J. D., and Wahl, S.: Evaluation and Post-Processing of Reanalysis Wind Speeds for Renewable Energy Applications, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-449, https://doi.org/10.5194/ems2021-449, 2021.

Vladimir Platonov and Anna Shestakova

The number of severe weather events at the Arctic region increased significantly. Its formation related generally to the mesoscale processes including downslope windstorms over Novaya Zemlya, Svalbard, Tiksi bay accompanied by strong winds. Therefore, its investigation required detailed hydrometeorological and climatic information with a horizontal resolution of at least several kilometers. This work aims to investigate extreme wind speeds statistics associated with downslope windstorms and evaluate it according to the COSMO-CLM Russian Arctic hindcast, ASR reanalysis, stations and satellite data.

COSMO-CLM Russian Arctic hindcast created in 2020 covers the 1980–2016 period with grid size ~12 km and 1-hour output step, containing approximately a hundred hydrometeorological characteristics, as well at surface, as on the 50 model levels. The primary assessments of the surface wind speed and temperature fields showed good agreement with ERA-Interim reanalysis in large-scale patterns and many added values in the regional mesoscale features reproduction according to the coastlines, mountains, large lakes, and other surface properties.

Mean values, absolute and daily maxima of wind speed, high wind speed frequencies were estimated for the COSMO-CLM Russian Arctic hindcast and the well-known Arctic System Reanalysis (ASRv2) for a 2000-2016 period. COSMO-CLM showed higher mean and daily maximal wind speed areas concerned to coastal regions of Svalbard and Scandinavia, over the northern areas of Taymyr peninsula. At the same time, the absolute wind speed maxima are significantly higher according to ASRv2, specially over the Barents Sea, near the Novaya Zemlya coast (differences are up to 15-20 m/s). The same pattern observed by a number of days with wind speed above the 30 m/s threshold. Compared with station data, the ASRv2 reproduced mean wind speeds better at most coastal and inland station, MAE are within 3 m/s. For absolute wind speed maxima differences between two datasets get lower, the COSMO-CLM hindcast is quite better for inland stations.

Model capability to reproduce strong downslope windstorms evaluated according to the observations timeseries over Novaya Zemlya, Svalbard and Tiksi stations during bora conditions. Generally, the ASRv2 reproduced the wind direction closer to observations and the wind speed worser than COSMO-CLM. The extreme wind speed frequencies during bora cases have less errors according to COSMO-CLM hindcast (up to ~5%) compared to the ASRv2 data (up to 10%). At the same time, moderate wind speed frequencies are reproduced by ASRv2 better.

Five specific Novaya Zemlya bora cases were evaluated according to SAR satellite wind speed data. Both ASRv2 and COSMO-CLM overestimated mean wind speed (MAE 0.5-6 m/s), maximal wind speed bias has different signs, however, the COSMO-CLM is better in most cases. Extreme percentiles biases (99 and 99.9%), correlation, structure and amplitude (according to the SAL method) are closer to observations by the COSMO-CLM hindcast.

How to cite: Platonov, V. and Shestakova, A.: Downslope windstorms and associated extreme wind speed statistics evaluation according to the COSMO-CLM Russian Arctic hindcast, ASR reanalysis and observations, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-300, https://doi.org/10.5194/ems2021-300, 2021.

Monika Goeldi, Stefanie Gubler, Christian Steger, Simon C. Scherrer, and Sven Kotlarski

Snow cover is a key component of alpine environments and knowledge of its spatiotemporal variability, including long-term trends, is vital for a range of dependent systems like winter tourism, hydropower production, etc. Snow cover retreat during the past decades is considered as an important and illustrative indicator of ongoing climate change. As such, the monitoring of surface snow cover and the projection of its future changes play a key role for climate services in alpine regions.

In Switzerland, a spatially and temporally consistent snow cover climatology that can serve as a reference for both climate monitoring and for future snow cover projections is currently missing. To assess the value and the potential of currently available long term spatial snow data we compare a range of different gridded snow water equivalent (SWE) datasets for the area of Switzerland, including three reanalysis-based products (COSMO-REA6, ERA5, ERA5-Land). The gridded data sets have a horizontal resolution between 1 and 30 km. The performance of the data sets is assessed by comparing them against three reference data sets with different characteristics (station data, a high-resolution 1km snow model that assimilates snow observations, and an optical remote sensing data set). Four different snow indicators are considered (mean SWE, number of snow days, date of maximum SWE, and snow cover extent) in nine different regions of Switzerland and six elevation classes.

The results reveal high temporal correlations between the individual datasets and, in general, a good performance regarding both countrywide and regional estimates of mean SWE. In individual regions, however, larger biases appear. All data sets qualitatively agree on a decreasing trend of mean SWE during the previous decades particularly at low elevations, but substantial differences can exist. Furthermore, all data sets overestimate the snow cover fraction as provided by the remote sensing reference. In general, reanalysis products capture the general characteristics of the Swiss snow climatology but indicate some distinctive deviations – e.g. like a systematic under- respectively overestimation of the mean snow water equivalent.

How to cite: Goeldi, M., Gubler, S., Steger, C., Scherrer, S. C., and Kotlarski, S.: The Swiss snow climatology as seen by reanalyses and further gridded datasets, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-282, https://doi.org/10.5194/ems2021-282, 2021.

Vojtěch Bližňák, Lucie Pokorná, and Zuzana Rulfová

Regional reanalysis constitute of potentially attractive new data source for many applications. They can offer added value benefiting from their higher spatial and temporal resolution. On the other hand, similarly as other data sources in regular network, reanalysis comes with uncertainties especially in the case of extremes.

The monitoring capability of reanalysis is essential for their usage as the reference datasets for climate model validation as well as for hydrological models. In this contribution, we evaluate the agreement of precipitation between modern reanalysis products (Era5, Era5 Land, Harmonie and Mescan-Surfex; with resolution between 5.5 and 32 km) and observed data at different time scales, from annual to subdaily. Studied characteristics of precipitation are for instance annual cycle of precipitation amount and the number of wet days, diurnal cycle of precipitation, and extremes. The common period for all datasets is 2002 – 2018.

Observed data used in this study are represented by adjusted radar-derived precipitation totals in 1km raster over the Czech Republic. The adjusted radar-derived precipitation totals are gained as follows. First, the radar-derived rain field is spatially adjusted to measurements from the rain gauges as a whole. For each day, the ratio between the mean 1-day precipitation total calculated from all rain gauges and the mean 1-day precipitation total estimated from the corresponding radar pixels is determined and used for a multiplication of radar-derived precipitation in every pixel of the radar domain. Second, the spatially adjusted radar-derived rain rates are locally adjusted in individual pixels based on the distance from the closest rain gauge, whereas the weight of the observed precipitation is decreasing with increasing distance to the given pixel. Adjusted daily precipitation total is then divided according to 10 min radar-derived estimates, from which the precipitation accumulations of longer duration are calculated.

How to cite: Bližňák, V., Pokorná, L., and Rulfová, Z.: Uncertainties in ERA-5 and UERRA reanalyses detected by adjusted radar precipitation - a regional study for the Czech Republic, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-60, https://doi.org/10.5194/ems2021-60, 2021.

Platon Patlakas, Christos Stathopoulos, Nikolaos S. Bartsotas, Helena Flocas, and George Kallos

The Arabian Peninsula is an area of diverse climatic conditions due to its location and geomorphological characteristics. It is affected by the Indian monsoon in the South and the Mediterranean synoptic scale systems in the North. It has mostly a desert-type climate with extreme heat during the daytime and very low annual rainfall. Extreme precipitation events related to convective activities are not that rare though. Such events often lead to flooding and may pose threat to human life and activities. At the same time, their local nature makes the recording and forecasting very difficult but essential. Towards a better understanding of the spatiotemporal features, the links and the feedback mechanisms associated with precipitation, a thirty year regional climatic analysis has been prepared with the aid of the state-of-the-art modeling system RAMS/ICLAMS. Its two-way interactive nesting capabilities, explicit cloud microphysical schemes with seven categories of hydrometeors and the ability to handle dust aerosols as predictive quantities make it suitable in an area where dust is a dominant factor. Special focus is given to densely populated regions and in several cities in order to cover the climatic features of the Peninsula to a satisfactory extent. An extended evaluation based on in-situ measurements and satellite records is performed. A monthly climatic analysis is performed alongside with a trend analysis to better assess the rainfall patterns throughout the thirty year period. The extremes are studied under the principles of the extreme value theory focusing not only on the duration but also on the intensity of the events. To enrich the analysis of the precipitation climate an examination of droughts is also performed. Apart from the strict scientific interest, the outcome has an added value in civil protection and several industries such as constructions and reinsurance.

How to cite: Patlakas, P., Stathopoulos, C., Bartsotas, N. S., Flocas, H., and Kallos, G.: A precipitation climatology for the Arabian Peninsula: from droughts to extremes, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-23, https://doi.org/10.5194/ems2021-23, 2021.

Iason Markantonis, Diamando Vlachogiannis, Thanasis Sfetsos, Ioannis Kioutsioukis, and Nadia Politi

Climate change is set to affect extreme climate and meteorological events. The combination of interacting physical processes (climate drivers) across various spatial and temporal scales resulting to an extreme event is referred to as compound event. So far, climate change impacts on compound events in Greece such as daily cold-wet events have not been explored. The complex geography and topography of Greece forms a variety of regions with different local climate and a great range in daily minimum temperature and precipitation distributions. This leads to the assumption that there we will also observe a variety in the distribution of cold-wet events depending on the region. Aim of our study in this work is first to identify the cold-wet events based on observational data and then to examine the predictive capability of regional different climate models and ERA-Interim against observations from the Hellenic National Meteorological Service (HNMS) stations for the occurrence of cold-wet compound events in the present climate. The study will focus on the colder and wetter period of the year (November-April) to determine the extremes for this period. Specifically, the datasets employed are from two EURO-CORDEX Regional Climate Models (RCMs) with 0.11° horizontal resolution and validated ERA-Interim Reanalysis downscaled with the Weather Research and Forecasting (WRF) model at 5km horizontal resolution, for the historical period 1980-2004. In particular, the RCM datasets analyses have been produced from SMHI-RCA4 driven by MPI-M-MPI-ESM-LR Global Climate Model (GCM) and CLMcom-CLM-CCLM4-8-17 driven by MOHC-HadGEM2-ES GCM. After the comparison with the observations, the gridded data from the models will give us the ability to observe the spatial distribution of the compound events.

How to cite: Markantonis, I., Vlachogiannis, D., Sfetsos, T., Kioutsioukis, I., and Politi, N.: An Investigation of cold-wet Compound Events in Greece, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-188, https://doi.org/10.5194/ems2021-188, 2021.

Annika Vogel, Ghazi Alessa, Robert Scheele, Lisa Weber, and Stephanie Fiedler

Aerosols are known to affect atmospheric processes on a wide range of spatio-temporal scales, from dust storms reducing incoming solar radiation to aerosol-climate feedbacks. Although plenty of studies address aerosol radiative forcing, there are persistent differences in current aerosol estimates from both, observations and models. Global reanalyses are able to provide consistent estimates of aerosol distributions by combining these two data sources. However, continuous assimilation of single satellite products forces the analyses towards the satellites climatology including possible inaccuracies. This study investigates differences between current estimates of aerosol optical depth (AOD) by addressing two questions: (1.) How well do we know the large-scale spatio-temporal pattern of present-day AOD across state-of-the-art data? (2.) How does current global aerosol reanalyses perform in comparison to other model- and observation-based data sets? To answer these questions, AOD from the global CAMS and MERRA-2 reanalyses is compared to 8 satellite products, 1 established climatology and 4 multi-model ensembles. The comprehensive data set used in this study allows to evaluate the performance of individual products concerning different spatial and temporal aspects. The evaluation covers results from 1998 to 2019, including most recently available products like the climate model inter-comparison project CMIP6.

Spatially and temporally averaged AOD from MERRA-2 agrees well with the mean satellite climatology, while the CAMS climatology is higher than most other products. With relative standard deviations of about 11%, temporal variations of CAMS and MERRA-2 agree well with the mean satellite variation. However, averaged AOD from the individual satellites show large differences, ranging from 0.124 for MISR to 0.164 for MODIS. In addition to average differences, spatial patterns vary significantly between the individual data sets. Because the CAMS reanalysis only assimilates AOD from MODIS, it remains close to the MODIS climatology which overestimates AOD in most regions in comparison to other products. This overestimation is considerably increased over eastern China were CAMS simulates regional values of more than 1.2 during summer. By assimilating both, MODIS and MISR data, the MERRA-2 reanalysis is closer to the satellite mean under most conditions. Although annual deviations remain small compared to other models, MERRA-2 tends to underestimate AOD at the equator and overestimates AOD at higher latitudes especially during the winter-season. The spatio-temporal differences between individual aerosol data sets underline the need for further research on both satellite retrievals and model simulations for aerosols. For example, integrating multiple observations in a reanalysis system would allow to compensate for inaccuracies of the individual products. Further developing the multi-scale coupled ICON-ART system at the German Weather Service provides a promising environment to achieve accurate aerosol climatologies on high spatial resolution.

How to cite: Vogel, A., Alessa, G., Scheele, R., Weber, L., and Fiedler, S.: Assessment of present-day estimates of AOD from global reanalyses against different satellite products and multi-model ensembles, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-228, https://doi.org/10.5194/ems2021-228, 2021.

Pedro Salvador, Jorge Pey, Noemí Pérez, Xavier Querol, and Begoña Artíñano

All African dust outbreaks (ADO) that increased the regional-background levels of PM at the SE, SW, E, Central, NW, N and NE regions of the Iberian Peninsula and the Balearic Islands have been identified from 2001 to 2020 using a well-known procedure (https://www.miteco.gob.es/images/es/metodologiaparaepisodiosnaturales-revabril2013_tcm30-186522.pdf). However, the meteorological tools (air mass back-trajectories, satellite imagery and numerical models for prediction of dust levels) and data bases (time series of PM levels registered at air quality monitoring stations) needed to perform the identification procedure are scarce before the year 2000. For this reason the occurrence of ADO in the western Mediterranean basin during the former decades has not been well addressed so far.

In this study we have used the NCEP/NCAR global reanalysis dataset fields of meteorological parameters to characterize the main distinctive synoptic dynamic and thermodynamic features that were associated to the development of ADO and analyze their time evolution in the 1948-2020 period.

First, the main synoptic circulation types that favored the occurrence of the ADO identified in 2001-2020 were obtained. With this aim, a circulation classification methodology was applied for classifying the daily fields of geopotential height at the 850 hPa level at 12 UTC into prevalent atmospheric circulation types for SW Europe and NW Africa.

Next, the daily mean values of the 1000-500 hPa layer geopotential thickness (GT), the mean 925-700 layer potential temperature (TPOT) and the anomalies of temperature at 850 hPa (TANOM) were computed over all the regions of study. High values of these thermodynamic variables are associated with the presence of warm, stable and dry air masses. In fact, significantly higher values of GT, TPOT and TANOM were obtained in all regions during days under ADO circulation types than during other days in 2001-2020.

Finally, we analyzed the time evolution of all the days under ADO circulation types and their associated daily mean values of GT, TPOT and TANOM over the regions of study using the Theil-Sen method from 1948 to 2020.

Our results show that the monthly number of days under ADO circulation types display an upward trend of 0.06 monthly days per year at the 99.9% confidence level. Statistically significant upward trends for the monthly mean values of GT, TPOT and TANOM were also obtained over all the regions during days under prevailing ADO circulation types in summer, spring and winter.

In summary, the frequency of the dynamic and thermodynamic synoptic conditions favouring the development of ADO over the western Mediterranean basin has increased over the last 70 years. These results are in line with the exacerbation of warm conditions registered in southern Europe during the last decades.


This study was funded by research project POSAHPI (Agencia Estatal de Investigación, PID2019-108101RB-I00).

How to cite: Salvador, P., Pey, J., Pérez, N., Querol, X., and Artíñano, B.: Characterization of dynamic and thermodynamic features of African dust outbreaks over the western Mediterranean basin: trend analysis for the 1948-2020 period, EMS Annual Meeting 2021, online, 6–10 Sep 2021, EMS2021-169, https://doi.org/10.5194/ems2021-169, 2021.


Comments on the session UP3.6

to access the discussion

Supporters & sponsors