This is the second of two papers that document the creation of the New European Wind Atlas (NEWA). In Part 1, we described the sensitivity experiments and accompanying evaluation done to arrive at the final mesoscale model setup used to produce the mesoscale wind atlas. In this paper, Part 2, we document how we made the final wind atlas product, covering both the production of the mesoscale climatology generated with the Weather Research and Forecasting (WRF) model and the microscale climatology generated with the Wind Atlas Analysis and Applications Program (WAsP). The paper includes a detailed description of the technical and practical aspects that went into running the mesoscale simulations and the downscaling using WAsP. We show the main results from the final wind atlas and present a comprehensive evaluation of each component of the NEWA model chain using observations from a large set of tall masts located all over Europe. The added value of the WRF and WAsP downscaling of wind climatologies is evaluated relative to the performance of the driving ERA5 reanalysis and shows that the WRF downscaling reduces the mean wind speed bias and spread relative to that of ERA5 from $-\mathrm{1.50}\pm \mathrm{1.30}$ to 0.02±0.78 m s⁻¹. The WAsP downscaling has an added positive impact relative to that of the WRF model in simple terrain. In complex terrain, where the assumptions of the linearized flow model break down, both the mean bias and spread in wind speed are worse than those from the raw mesoscale results.

This is the first of two papers that document the creation of the New European Wind Atlas (NEWA). It describes the sensitivity analysis and evaluation procedures that formed the basis for choosing the final setup of the mesoscale model simulations of the wind atlas. The suitable combination of model setup and parameterizations, bound by practical constraints, was found for simulating the climatology of the wind field at turbine-relevant heights with the Weather Research and Forecasting (WRF) model. Initial WRF model sensitivity experiments compared the wind climate generated by using two commonly used planetary boundary layer schemes and were carried out over several regions in Europe. They confirmed that the most significant differences in annual mean wind speed at 100 m a.g.l. (above ground level) mostly coincide with areas of high surface roughness length and not with the location of the domains or maximum wind speed. Then an ensemble of more than 50 simulations with different setups for a single year was carried out for one domain covering northern Europe for which tall mast observations were available. We varied many different parameters across the simulations, e.g. model version, forcing data, various physical parameterizations, and the size of the model domain. These simulations showed that although virtually every parameter change affects the results in some way, significant changes in the wind climate in the boundary layer are mostly due to using different physical parameterizations, especially the planetary boundary layer scheme, the representation of the land surface, and the prescribed surface roughness length. Also, the setup of the simulations, such as the integration length and the domain size, can considerably influence the results. We assessed the degree of similarity between winds simulated by the WRF ensemble members and the observations using a suite of metrics, including the Earth Mover's Distance (EMD), a statistic that measures the distance between two probability distributions. The EMD was used to diagnose the performance of each ensemble member using the full wind speed and direction distribution, which is essential for wind resource assessment. We identified the most realistic ensemble members to determine the most suitable configuration to be used in the final production run, which is fully described and evaluated in the second part of this study (Dörenkämper et al., 2020).

Europe's offshore wind resource mapping is part of the New European Wind Atlas (NEWA) international consortium effort. This study presents the results of analysis of synthetic aperture radar (SAR) ocean wind maps based on Envisat and Sentinel-1 with a brief description of the wind retrieval process and Advanced Scatterometer (ASCAT) ocean wind maps. The wind statistics at 10 and 100 m above mean sea level (a.m.s.l.) height using an extrapolation procedure involving simulated long-term stability over oceans are presented for both SAR and ASCAT. Furthermore, the Weather Research and Forecasting (WRF) offshore wind atlas of NEWA is presented. This has 3 km grid spacing with data every 30 min for 30 years from 1989 to 2018, while ASCAT has 12.5 km and SAR has 2 km grid spacing. Offshore mean wind speed maps at 100 m a.m.s.l. height from ASCAT, SAR, WRF and ERA5 at a European scale are compared. A case study on offshore winds near Crete compares SAR and WRF for flow from the north, west and all directions. The paper highlights the ability of the WRF model to simulate the overall European wind climatology and the near-coastal winds constrained by the resolution of the coastal topography in the WRF model simulations.

Climate change has been shown to impact aspects of agriculture and phenology. This study aims to quantify changes in the timing of garden strawberry blooms and harvests in the Baltic States using Regional Climate Models (RCMs). First, parameters for a strawberry phenology model based on the growing degree day (GDD) methodology were determined. Growing degree days were calculated using a modified sine wave method that estimates the diurnal temperature cycle from the daily maximum and minimum temperature. Model parameters include the base temperature and the required cumulative GDD sum, estimated from phenological and meteorological observations in Latvia for the years 2010–2013 via iterative calibration. Then an ensemble of bias-corrected RCM results (ENSEMBLES project) was used as input to the phenological model to estimate the timing of strawberry phenological processes for the years 1951–2099. The results clearly show that strawberry phenological processes can be expected to occur earlier in the future, with a significant change in regional patterns. Differences between coastal and inland regions are expected to decrease over time. The uncertainty of the results was estimated using the RCM ensemble spread, with northern coastal locations showing the largest spread.

Wind power forecasting greatly relies on wind speed forecasts. Numerical Weather Prediction (NWP) models are a reliable source of meteorological forecasts and they can also be used in wind resource assessment. In this work we carry out the verification of wind speed results from the NWP model Weather Research and Forecast (WRF), grid resolution - 3km. Results from 172 model runs in May and November 2013 are compared with meteorological observations in 24 stations in Latvia.The model usually predicts wind speed values that are larger than the observed and the diurnal cycle has a large impact on verification results. Verification results obtained by interpolating model results between gridpoints to the exact station location are compared with results acquired using values from the nearest gridpoint. This distinction is significant in coastal stations. The influence of wind speed and wind direction on model error is analyzed

This is the first of two papers that documents the creation of the New European Wind Atlas (NEWA). It describes the sensitivity analysis and evaluation procedures that formed the basis for choosing the final setup of the mesoscale model simulations of the wind atlas. An optimal combination of model setup and parameterisations was found for simulating the climatology of the wind field at turbine-relevant heights with the Weather Research and Forecasting (WRF) model. Initial WRF model sensitivity experiments compared the wind climate generated by using two commonly used planetary boundary layer schemes and were carried out over several regions in Europe. They confirmed that the largest differences in annual mean wind speed at 100 m above ground level mostly coincide with areas of high surface roughness length and not with the location of the domains or maximum wind speed. Then an ensemble of more than 50 simulations with different setups for a single year was carried out for one domain covering Northern Europe, for which tall mast observations were available. Many different parameters were varied across the simulations, for example, model version, forcing data, various physical parameterisations and the size of the model domain. These simulations showed that although virtually every parameter change affects the results in some way, significant changes on the wind climate in the boundary layer are mostly due to using different physical parameterisations, especially the planetary boundary layer scheme, the representation of the land surface, and the prescribed surface roughness length. Also, the setup of the simulations, such as the integration length and the domain size can considerably influence the results. The degree of similarity between winds simulated by the WRF ensemble members and the observations was assessed using a suite of metrics, including the Earth Mover's Distance (EMD), a statistic that measures the distance between two probability distributions. The EMD was used to diagnose the performance of each ensemble member using the full wind speed distribution, which is important for wind resource assessment. The most realistic ensemble members were identified to determine the most suitable configuration to be used in the final production run, which is fully described and evaluated in the second part of this study.