Skip to main content
Have a personal or library account? Click to login
On the diurnal cycle and variability of winds in the lower planetary boundary layer: evaluation of regional reanalyses and hindcasts Cover

On the diurnal cycle and variability of winds in the lower planetary boundary layer: evaluation of regional reanalyses and hindcasts

Open Access
|Jan 2021

Figures & Tables

Fig. 1.

Model domains related to the different reanalyses and hindcasts (see Table 2). CCLM-oF is integrated within the red polygon, COSMO-REA6 within the blue polygon and UE-UKMO within the green polygon. The UE-SMHI domain is denoted by the orography map [m].

Table 1.

Location, boom height above ground [m], estimation of roughness [m] and data availability for 2006–2007 [%] at the sites considered for evaluation. The roughnesses for the observational sites are not extracted in a consistent way (turbulent flux measurements at different heights vs. Richardson methods). The values of the hindcasts are extracted from the grid morphology near the site location. The values from UE-UKMO are not part of the standard output due to tile-based calculation within land-surface model JULES (Best et al., 2011; Clark et al., 2011).

StationBillwerderFalkenbergCabauwKarlsruheAbbreviationBIFACAKALocation53° 31′ 09.0″ N, 10° 06′ 10.3″ E52° 10′ 01″ N, 14° 07′ 27″ E52° 58′ 18″ N, 4°55′ 37″ E49° 5′ 33″ N, 8°25′ 33″ Ez0, Obs0.04.. 0.16 and peak: W (0.75)0.01.. 0.030.03.. 0.15 and min: SW (0.01)no dataz0, UE-SMHI0.41.. 0.440.25.. 0.270.09.. 0.100.45.. 0.47z0, COSMO0.470.140.050.43Booms50, 70, 110, 175, 250, 28010, 20, 40, 60, 80, 9810, 20, 40, 80, 140, 20020, 30, 40, 50, 60, 80, 100, 130, 160, 200Data availability94.798.7100.096.7
Table 2.

The models evaluated in this study. The PBL model levels are with respect to the lowest 1500 m. The positions of the model levels are related to flat terrain and to a surface pressure of 1013.26 hPa. For interpretation of the results, the datasets are indicated with the following abbreviations: CCLM-oF-SN, CCLM-oF, COSMO-REA6, UE-SMHI and UE-UKMO.

CCLM-oFCCLM-oF-SNCOSMO-REA6UE-SMHIUE-UKMOModelCOSMOCOSMOCOSMOHARMONIEUMDynamical coreNon-hydrostatic (Baldauf et al., 2011)Non-hydrostatic (Baldauf et al., 2011)Non-hydrostatic (Baldauf et al., 2011)Hydrostatic assumption in the core of Bubnov et al. (1995)Non-hydrostatic (Wood et al., 2014)Assimilation technique–Continuous large-scale nudging above the PBLContinuous nudging/surface analysis3D-VAR/surface analysis4D-VAR/surface analysisResolution≈7 km≈7 km≈6 km≈11 km≈12 kmDriving modelMERRA2MERRA2ERA-InterimERA-InterimERA-InterimLocation of model levels (up to 350 m)10, 35, 69, 116, 179, 259, 350As CCLM-oFAs CCLM-oF12, 38, 63, 89, 117, 146, 177, 211, 247, 287, 32910, 37, 77, 130, 197, 277PBL model levels1414142414
Table 3.

The static stability classes used for investigation. The classification is shown to be dependent on the vertical gradients of (potential) temperature. Moreover, the likelihood of occurrence [%] of static stability is listed for the flux towers.

StratificationdTdz [K m−1]dθdz [K m–1]BillwerderFalkenbergCabauwUnstable(–∞, −0.011](–∞, −0.001]12.39.918.2Near neutral[−0.011, −0.009][−0.001, 0.001]13.014.710.5Stable[−0.009, 0.001][0.001, 0.011]56.341.537.8Very stable[0.001, ∞)[0.011, ∞)18.333.933.2
Fig. 2.

Characteristic power curve as given by the ENERCON E70 wind turbine (2310 kW, with a diameter of 70 m).

Fig. 3.

Hamburg, diurnal cycle of wind speed [m/s] (hourly values) for the whole year 2006–2007 for observation (a) and five different assimilation experiments (b–f). Note, that the data for COSMO-REA6 is not available above 250 m.

Fig. 4.

Same as Figure 3 but for the tower at Falkenberg and only for the months of May to September in 2006 and 2007.

Fig. 5.

Same as Figure 3 but for the tower at Karlsruhe and only for the months of May to September in 2006 and 2007.

Table 4.

Kling-Gupta efficiency for the mean diurnal cycle in ‘MJJA’ (summer) and ‘NDJFM’ (winter). The towers are listed, and the heights are grouped for the near-surface, above 100 m and in between. The bold numbers indicate that the performance exceeds 0.4, which requires, for instance, r0.75,σr=[0.75,1.25] and μr=[0.75,1.25]. Model names are shortened by omitting the first part.

oFoF-SNREA6UKMOSMHISumWinSumWinSumWinSumWinSumWinBI 10–500.70.60.70.60.80.70.70.60.40.1BI 1100.60.30.60.40.60.10.50.20.1−2.3BI 175–2500.90.80.90.70.70.70.30.60.6−0.1KA 10–300.50.60.60.70.50.20.6−0.10.8−0.1KA 40–80−0.0−0.50.10.0−0.2−0.9−0.1−1.00.3−1.6KA 110–1800.50.40.60.60.30.00.0−0.10.4−0.1CA 10–400.80.80.80.80.80.80.80.80.50.4CA 800.60.40.60.40.50.10.5−0.10.3−1.1CA 140–2000.70.40.50.70.20.10.2−0.20.6−0.1FA 10–400.80.50.80.50.90.80.70.70.40.2FA 60-980.20.20.40.30.20.40.60.10.1−1.6
Table 5.

Estimation of the shape parameter (first number in each cell) and scale parameter (second number) of the Weibull distribution with respect to the selected measurement locations. Model names are shortened by omitting the first part. The ‘star’ indicates that the uncertainty of parameter estimation is larger than 0.1. In addition, EMD is given [m s−1].

Weibull distributionEMDLocationObsoFoF-SNREA6UKMOSMHIoFoF-SNREA6UKMOSMHIBI 50*2.2/5.72.2/6.12.2/6.02.2/5.92.2/6.72.1/5.40.40.30.20.90.1BI 1102.4/7.22.4/8.02.3/7.72.3/7.62.2/7.92.4/7.00.70.50.30.60.2BI 1752.4/8.52.5/9.52.4/9.32.3/8.92.2/8.82.4/8.00.90.60.40.40.4BI 2502.3/9.82.4/11.12.2/10.42.2/10.12.1/9.62.2/8.90.90.60.30.40.9KA 60*2.1/5.31.7/5.41.7/5.21.7/5.01.6/5.21.7/5.00.40.40.30.50.3KA 110*2.0/6.71.8/6.61.8/6.31.6/5.81.6/6.01.9/6.30.40.30.20.30.2CA 802.4/8.02.4/8.32.3/8.22.2/7.82.2/7.82.3/7.30.30.20.20.20.6CA 1402.4/9.22.4/9.62.3/9.42.2/9.12.2/8.72.3/8.40.40.20.20.40.8FA 60*2.3/5.92.3/6.52.2/6.32.2/5.92.2/6.62.4/6.00.50.40.10.70.2FA 982.4/6.72.3/7.42.3/7.22.3/6.92.2/7.42.4/7.10.60.50.20.50.3
Fig. 6.

Median wind speed bias at 00, 06, 12 and 18 UTC at the towers Cabauw (left panel) and Falkenberg (right panel). Top: The UE-SMHI (solid lines, lower x-axis) and UE-UKMO reanalyses (dashed lines, upper x-axis) are shown with respect to the years 2006 and 2007 for summer. Middle: The same is shown but for CCLM-oF-SN (solid lines, lower x-axis) and COSMO-REA6 (dashed lines, upper x-axis). Bottom: The same is shown as in the middle panel but for winter.

Fig. 7.

Diurnal cycle of occurrence of 4 different types of stratification (Table 3) for summer (‘MJJAS’, left) and winter (‘NDJFM’, right) and at towers Cabauw (a, b), Falkenberg (c, d) and Billwerder (e, f). The straight lines indicate the absolute frequency measured, and the symbols indicate the absolute frequencies simulated (‘x’: UE-UKMO, rectangle: UE-SMHI, range indicator: COSMO models). For readability, the simulated frequency of ‘st’ cases is shown only in winter. The colors belong to different stabilities (legend of the figures in the top row).

Fig. 8.

Temperature profiles depending on the stratification at Cabauw (a), Falkenberg (b) and Billwerder (c). Modified box plots are shown for the observations (gray), COSMO-REA6 (dark blue), CCLM-oF-SN (light blue), UE-SMHI (green) and UE-UKMO (red). The upper part of each figure refers to the very stable stratifications, and the lower part refers to the unstable stratifications. ‘Modified box plot’ means that only the upper and lower quartiles of the data are visualized by the box, and the median is indicated by the straight line in the box and the mean by the small black rectangle. Furthermore, transparency is used with the last gray boxes at Billwerder to distinguish between the data. Moreover, the horizontal dashed lines indicate the transition from a linear scale to a logarithmic scale (see the y-axis).

Fig. 9.

Two-year quantile-quantile plot of the vertical gradient in wind speed depending on the stratification at Cabauw (a) (80 m to 20 m), Falkenberg (b) (80 m to 20 m) and Billwerder (c) (110 m to 50 m). Each second quantile (p = 0.02, 0.04, 0.06,…, 0.98, 1.0) of the distribution of simulated gradients is plotted against the quantiles of the observed gradients. The upper part of each figure refers to the unstable stratifications (upper x-axis and left y-axis), and the lower part refers to the very stable stratifications (lower x-axis and right y-axis).

Table 6.

RMSErbc at each tower for hourly data in ‘MJJA’ (summer) and ’NDJFM’ (winter). The heights are grouped for the near-surface, above 100 m and in between.

CCLM-oFCCLM-oF-SNCOSMO-REA6UE-UKMOUE-SMHISumWinSumWinSumWinSumWinSumWinKA 10–300.510.470.440.410.380.320.420.370.420.42KA 40–800.470.420.400.360.340.290.360.310.360.33KA 110–1800.470.380.390.310.320.250.320.240.330.26BI 10–500.440.360.360.310.320.260.350.300.360.32BI 1100.400.320.320.260.280.210.270.230.280.23BI 175–2500.380.340.310.280.270.230.270.220.260.22CA 10–400.420.310.330.250.250.190.250.190.280.22CA 40–800.380.280.300.220.250.170.230.160.250.19CA 140–2000.370.280.290.220.240.170.220.140.220.16FA 10–400.440.320.350.260.290.220.300.230.330.26FA 60–980.420.310.340.240.270.200.280.210.300.23
Fig. 10.

Cumulative frequency distribution of 6-hourly ramp events depending on the relative ramp rate for Billwerder at 110 m height and Karlsruhe at 140 m height (a) as well as for Cabauw at 140 m height and Falkenberg at 98 m height (b). (c) is the same as (b) but for 2-hourly ramps.

Fig. 11.

Success ratio (SUC) and probability of detection (POD) of ramps for different measurement locations and reanalysis products. The scores are shown for small ramps defined as the percentage change between 15 and 25% and regarding the analysis times 00, 06, 12 and 18 UTC (window length equals 6 hours). The color code is the same as in Figures 9 and 10.

Language: English
Page range: 1804294 - 1804294
Published on: Jan 1, 2021
Published by: Stockholm University Press
In partnership with: Paradigm Publishing Services

© 2021 Ronny Petrik, Beate Geyer, Burkhardt Rockel, published by Stockholm University Press
This work is licensed under the Creative Commons Attribution 4.0 License.