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.

Due to the major role of the sun in heating the earth's surface, the atmospheric planetary boundary layer over land is inherently marked by a diurnal cycle. The afternoon transition, the period of the day that connects the daytime dry convective boundary layer to the night-time stable boundary layer, still has a number of unanswered scientific questions. This phase of the diurnal cycle is challenging from both modelling and observational perspectives: it is transitory, most of the forcings are small or null and the turbulence regime changes from fully convective, close to homogeneous and isotropic, toward a more heterogeneous and intermittent state. These issues motivated the BLLAST (Boundary-Layer Late Afternoon and Sunset Turbulence) field campaign that was conducted from 14 June to 8 July 2011 in southern France, in an area of complex and heterogeneous terrain. A wide range of instrumented platforms including full-size aircraft, remotely piloted aircraft systems, remote-sensing instruments, radiosoundings, tethered balloons, surface flux stations and various meteorological towers were deployed over different surface types. The boundary layer, from the earth's surface to the free troposphere, was probed during the entire day, with a focus and intense observation periods that were conducted from midday until sunset. The BLLAST field campaign also provided an opportunity to test innovative measurement systems, such as new miniaturized sensors, and a new technique for frequent radiosoundings of the low troposphere. Twelve fair weather days displaying various meteorological conditions were extensively documented during the field experiment. The boundary-layer growth varied from one day to another depending on many contributions including stability, advection, subsidence, the state of the previous day's residual layer, as well as local, meso- or synoptic scale conditions. Ground-based measurements combined with tethered-balloon and airborne observations captured the turbulence decay from the surface throughout the whole boundary layer and documented the evolution of the turbulence characteristic length scales during the transition period. Closely integrated with the field experiment, numerical studies are now underway with a complete hierarchy of models to support the data interpretation and improve the model representations.

A six‐year climatology of radiation fog has been compiled at two sites: the Research Centre for the Lower Atmosphere (CIBA, Spain) and the Cabauw Experimental Site for Atmospheric Research (CESAR, The Netherlands). These sites are contrasted in terms of geographical situation, climate zone, altitude, humidity and soil water availability. Therefore, several climatological differences in fog abundance, onset, dissipation and duration have been quantified between the two sites. The more humid site (CESAR) is characterised by relatively short radiation fog events distributed throughout the year. However, radiation fog at the drier site (CIBA) is more persistent and appears during late autumn/winter months. In general, its formation requires more time after sunset (∼2 h more), since further cooling is required to reach saturation. The forecast of these fog events has been evaluated through two different approaches. First, we extend the statistical method presented by Menut et al. (2014) (M14). This method uses statistics to define threshold values on key variables for fog formation (pre‐fog) and verifies its predictability using observations and numerical model output. We present some of the most appropriate threshold values for the forecasting of pre‐fog periods at both sites, which differ from those presented in M14 and depend on the optimisation of the hit rate or the false‐alarm rate. Additionally, we also extend M14 by suggesting other variables as potential predictors for fog formation (friction velocity and visibility tendency). Finally, we focus on fog simulation by the Weather Research and Forecasting (WRF) model in terms of liquid water content. The WRF model was able to simulate radiation fog when configured with sophisticated physical options and high resolution. However it failed in simulating the onset, dissipation and the vertical extent of fog (which was overestimated). The model results were extremely sensitive to the spin‐up time.

The prediction of fogs is one of the processes not well reproduced by the Numerical Weather Prediction (NWP) models. In particular, the role of turbulence in the formation or dissipation of fogs is one of the physical processed not well understood, and therefore, not well parameterized by the NWP models. Observational analysis of three different periods with fogs at the Spanish Northern Plateau has been carried out. These periods have also been simulated with the Weather Research and Forecasting (WRF) numerical model and their results have been compared to observations. The study includes a comparison of the skill of different planetary boundary layer (PBL) parameterizations, surface layer schemes and a test of the gravitational settling of clouds/fogs droplets option. A statistical analysis of this comparison has been evaluated in order to study differences between the periods and between the various parameterizations used. The model results for each PBL parameterization were different, depending on the studied period, due to differences in the features of each fog. This fact made it difficult to obtain generalized conclusions, but allowed us to determine which parameterization performed better for each case. In general, judging from the models results of liquid water content (LWC), none of the PBL schemes were able to correctly simulate the fogs, being Mellor-Yamada Nakanishi and Niino (MYNN) 2.5 level PBL scheme the best one in most of the cases. This conclusion is also supported by the root mean square error (RMSE) calculated for different meteorological variables.