The GABLS3 case has been selected in the NEWA project as a baseline exercise for the design of mesoscale-to-microscale methodologies for wind resource assessment. The case is suitable for the development of microscale wind farm models that incorporate realistic forcing, derived from a mesoscale model, along a typical diurnal case that leads to the development of a nocturnal low-level jet. Challenges of this case include: incorporating time- and height-dependent mesoscale forcing in microscale models, turbulence modeling at varying atmospheric stability conditions, defining suitable surface boundary conditions for momentum and heat and characterization of the wind profile in (non-logarithmic) LLJ conditions.
Data is provided open-access for registered participants. For practical reasons, a copy of the GABSL3 observational dataset, made available by KNMI, is provided as well.
Wind-energy specific objectives of the benchmark include:
- Demonstrate the capability of wind energy ABL models to incorporate realistic mesoscale forcing
- Implement surface boundary conditions suitable for wind assessment studies using mesoscale simulation data and/or observations (typical of wind energy campaigns)
- Develop suitable model calibration strategies for wind energy applications or, in other words, how to best use available measurements (typical of wind energy campaigns) to correct meso-micro predictions
- Define suitable metrics for validation of ABL models based on wind energy quantities of interest
By "typical wind energy campaigns" we would like to encourage modellers to prioritize observations that are common place in wind resource assessment campaigns (80 masts with velocity and temperature measurements, lidar profilers measuring up to 400 m).
The case set-up and input data of the original GABLS3 case can be found in the KNMI website. This is usefull if you want to compare with published results of the original SCM model intercomparison. In the original GABLS3 set-up, the simulated mesoscale tendencies are adjusted to produce a better match with the surface geostrophic wind obtained from a network of synoptic stations and the wind speed at 200-m measured at the Cabauw tower. Initial profiles are based on soundings measured near the Cabauw mast.
Alternatively, you can use inputs generated entirely from a WRF simulation, as described in . Here, instead of using observed initial projiles and adjusted mesoscale forcings you can use initial profiles and forcing produced directly from a mesoscale simulation. This is more representative of a wind energy model-chain set-up, where the inputs of a microscale model are generated by a "wind atlas" methodology that doesn't normally include corrections with local measurements. Instead, local adjustments are allowed at the microscale level by incorporating onsite measurements as if these measurements were part of a typical wind resource assessment campaign.
The WRF simulation is based on a one-way nesting configuration of three concentric square domains centered at the Cabauw site, of the same size 183x183, and at 9, 3 and 1 km horizontal resolution. The vertical grid, approximately 13 km high, is based on 46 terrain-following (eta) levels with 24 levels in the first 1000 m, the first level at approximately 13 m, a uniform spacing of 25 m over the first 300 m and then stretched to a uniform resolution of 600 m in the upper part. The U.S. Geological Survey (USGS) land-use surface data, that comes by default with the WRF model, is used together with the unified Noah land-surface model to define the boundary conditions at the surface. Other physical parameterizations used are: the rapid radiative transfer model (RRTM), the Dudhia radiation scheme and the Yonsei University (YSU) first-order PBL scheme. The simulation uses input data from ERA-Interim with a spin-up time of 24 hours. The WRF set-up follows the reference configuration of Kleczek et al , who run a sensitivity analysis of WRF showing reasonably good results at reproducing the nocturnal LLJ.
A NetCDF file is provided with the following information:
- Site coordinates and Coriolis parameter
- Time-height 2D arrays of velocity components (U,V,W) and potential temperature (Th)
- Time-height 2D arrays of mesoscale forcings (tendencies): geostrophic wind (Ug, Vg), advective wind (Uadv, Vadv) and advective potential temperature (Thadv)
- Time array of surface-layer quantities: friction velocity (ust), kinematic heat flux (wt), 2-m temperature (T2), skin temperature (TSK), surface pressure (Psfc)
Units, dimmensions and variables description are all provided in the NetCDF file. Momentum tendencies (Figure 2) are provided in [m s-1] and should be multiplied by the Coriolis parameter to obtain appropriate forces in [m s-2]. For convenience, we have ommitted information about humidity since the assumption of dry-atmosphere is typically adopted by wind energy flow models.
Figure 2: Time-height contour plots of the longitudinal wind component U and momentum budget terms:
Utend = Uadv + Ucor + Upg + Upbl. 
A python script is provided to show how to read the NetCDF input file and extract these variables.
The following quantities of interest (QoI) will be evaluated as described in , using a reference rotor size of 160 m diameter at a hub-height of 120 m (~ 7 MW turbine):
- Rotor equivalent wind speed (REWS)
- Hub-height wind direction (WDhub)
- Turbulence intensity at hub-height (TIhub)
- Wind shear (power-law exponent α) and wind veer (slope of linear fit to wind direction differences ψ) accross the rotor plane
- Surface-layer quantitites: T2, ust, wt and z/L
The evaluation consists on time-series plots of these QoIs along the diurnal cycle and mean-absolute error (MAE) integrated over the whole cycle.
The benchmark is mainly developed for microscale models that make use of the input data described above. However, mesoscale or multi-scale (online meso-micro) simulations are also welcome. The following suggestions are provided to guide the model runs:
- We shall use the 2-m temperature (T2) as our most practical reference to deduce the potential temperature surface boundary conditions using Monin Obukhov similarity theory, since this variable is routinely measured in measurement campaigns and is part of the standard output of meteorological models.
- Simulations may be based entirely on the mesoscale input data or incorporate measurements from the Cabauw mast. Priority should be given to measurements that can be found in "typical wind energy campaigns" (80 masts with velocity and temperature measurements, lidar profilers measuring up to 400 m).
- Sensitivity analysis of mesoscale models can be used to quantify the input uncertainty derived from the spread of the ensemble of simulations.
- Online multi-scale simulations models can be used as a reference for microscale models that are coupled to mesoscale asyncronously through the input data. To allow this comparison multi-scale simulations should be also run with ERA-Interim input data.
If resources allow, please use a spin-up time of 24 hours as in the input data.
Data should be provided in a single NetCDF file as described in the python template. The following output variables are requested:
- Time-height 2D arrays of: velocity components, potential temperature and turbulent kinetic energy
- Time 1D array of surface-layer quantities: friction velocity (ust, at 3 m), kinematic heat flux (wt, at 3 m) and 2-m temperature (T2)
- Time in hours since 2006-07-01 12:00 UTC
- Heights in meters (please provide model levels at least up to 4000 m)
- Sanz Rodrigo J, Churchfield M, Kosovic B (2016) A wind energy benchmark for ABL modelling of a diurnal cycle with a nocturnal low-level jet: GABLS3 revisited. J. Phys.: Conf. Ser. 753, 032024, doi: 10.1088/1742-6596/753/3/032024
- Sanz Rodrigo J., Churchfield M., Kosovic B. (2016) Atmospheric boundary layer modeling based on mesoscale tendencies and data assimilation at microscale. Wind Energy Science Discussions, doi:10.5194/wes-2016-26
- Kleczek MA, Steeveneveld GL and Holtslag AAM (2014) Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact on Boundary-Layer Schemes, Boundary Conditions and Spin-Up. Boundary-Layer Meteorol 152 213-243, doi: 10.1007/s10546-014-9925-3
- NEWA results will be first evaluated during the WP3 workshop in Barcelona, 21-23 September 2016.
- Results from external participants will be collected through February 2017 and included in a paper for the Wake Conference 2017 (abstract due on the 2nd of December 2016)
- A summary report will be included in the NEWA M24 deliverables (Feb 2017).
- The benchmark results will be presented to the Wake Conference and 2nd annual meeting of Task 31 Wakebench Phase 2, Visby, May 30 - June 2 2017.
- A follow-up analysis of turbulence quantities with LES models will be formulated as another benchmark
This benchmark has been produced within a Marie Curie International Outgoing Fellowship, MesoWake project, with funding from the European Union’s Seventh Framework Programme for research, technological development and demonstration under grant agreement no 624562.