the Creative Commons Attribution 4.0 License.

the Creative Commons Attribution 4.0 License.

# Improving mesoscale wind speed forecasts using lidar-based observation nudging for airborne wind energy systems

### Martin Dörenkämper

### Gerald Steinfeld

### Curran Crawford

Airborne wind energy systems (AWESs) aim to operate at altitudes above conventional wind turbines where reliable high-resolution wind data are scarce. Wind light detection and ranging (lidar) measurements and mesoscale models both have their advantages and disadvantages when assessing the wind resource at such heights. This study investigates whether assimilating measurements into the mesoscale Weather Research and Forecasting (WRF) model using observation nudging generates a more accurate, complete data set. The impact of continuous observation nudging at multiple altitudes on simulated wind conditions is compared to an unnudged reference run and to the lidar measurements themselves. We compare the impact on wind speed and direction for individual days, average diurnal variability and long-term statistics. Finally, wind speed data are used to estimate the optimal traction power and operating altitudes of AWES. Observation nudging improves the WRF accuracy at the measurement location. Close to the surface the impact of nudging is limited as effects of the air–surface interaction dominate but becomes more prominent at mid-altitudes and decreases towards high altitudes. The wind speed frequency distribution shows a multi-modality caused by changing atmospheric stability conditions. Therefore, wind speed profiles are categorized into various stability conditions. Based on a simplified AWES model, the most probable optimal altitude is between 200 and 600 m. This wide range of heights emphasizes the benefit of such systems to dynamically adjust their operating altitude.

- Article
(7986 KB) - Full-text XML
- BibTeX
- EndNote

The prospects of higher energy potential, more consistent strong winds and less turbulence in comparison to near-surface winds has sparked interest in mid-altitude wind energy systems, with “mid-altitude” defined here as heights above 100 m and below 1500 m,. Airborne wind energy systems are a novel class of renewable energy technology that harvest stronger winds at altitudes which are unreachable by current wind turbines, at potentially significantly reduced capital cost (Lunney et al., 2017; Fagiano and Milanese, 2012). For practical and economical reasons we focus on resource assessment within the lower part of the atmosphere, an altitude range spanned by the highly variable boundary layer. Unlike conventional wind energy, which has converged to a single concept with three blades and a conical tower, several different AWES designs are under investigation by numerous companies and research institutes worldwide (Cherubini et al., 2015). Various concepts are competing for entry into the market, ranging from ring-shaped aerostats to rigid wings and soft kites with different sizes, rated power and altitude ranges. Since this technology is still in an early stage, none are currently commercially available.

Developers and operators of large conventional wind turbines, AWES and drones require accurate wind data to estimate power output and mechanical loads. They currently rely on oversimplified approximations such as the logarithmic wind profile (Optis et al., 2016) or coarsely resolved reanalysis data sets (Archer and Caldeira, 2009; Bechtle et al., 2019) as the applicability of conventional spectral wind models (Burton, 2011) has not been verified for these altitudes. The first investigations (Fechner and Schmehl, 2018) resorting to the Mann model (Mann, 1994; IEC, 2005) have been conducted.

Recent advancements in wind lidar technology enable measurements at higher altitudes. This measurement technique, however, suffers from reduced data availability with increasing altitude caused by a decrease in aerosol density, which is needed for the backscattering of the lidar signal (Peña et al., 2015). No mid-altitude measurement device can reliably gather long-term, high-frequency data. The temporal and spatial resolution of lidar devices is insufficient to precisely measure high-frequency fluctuations, but estimated turbulence intensity correlates with sonic turbulence measurements for lower altitudes (Sathe et al., 2011). Balloon-mounted sonic anemometers are in early development (Canut et al., 2016). The expensive and time-consuming nature of measurements motivates the usage of numerical weather prediction models, such as the mesoscale Weather Research and Forecasting (WRF) model, as adequate tools to assess the synoptic characteristics of the atmospheric boundary layer (ABL) (Al-Yahyai et al., 2010). These models typically have a spatial resolution that ranges from 1 km to tens of kilometers and a temporal resolution of the order of minutes. Subgrid-scale high-frequency variations in resolved quantities are parameterized. Mesoscale models can be used to produced long-term reference data sets up to higher altitudes such as the New European Wind Atlas (Witha et al., 2019).

This work is a continuation of a previous investigation of mid-altitude wind lidar measurements (Sommerfeld et al., 2019). The measurements used in these studies were gathered as part of the OnKites II project (Gambier et al., 2017) at the Fraunhofer Institute for Wind Energy Systems (IWES) with the goal of evaluating the potential of AWES. This paper makes use of various statistical tools to describe the relationship between the mesoscale WRF model and lidar measurements to determine the impact of wind speed observation nudging (Mylonas-Dirdiris et al., 2016).

Section 2 describes the measurement campaign. Section 3 introduces the mesoscale model and observation nudging methodology used in this article. Section 4 quantifies the impact of observation nudging and summarizes the statistical differences between WRF and lidar. Results are applied to estimate optimal operating altitude and power output based on a simplified AWES model in Sect. 4.7. Section 5 concludes the article with an outlook and motivation for future work.

The lidar data used in this study (Bastigkeit et al., 2017) were collected between 1 September 2015 and 29 February 2016 at the Pritzwalk Sommersberg airport (lat: 53^{∘}10^{′}47.00^{′′} N, lon: 12^{∘}11^{′}20.98^{′′} E) in northern Germany (see white X in Fig. 1). The area surrounding the airport mostly consists of flat agricultural land with the town of Pritzwalk to the south.
A Galion4000 single-beam pulsed wind lidar from SgurrEnergy was used (Gottschall et al., 2009). Wind speed data were collected using the Doppler beam swinging (DBS) method (opening angle of 62^{∘}), which averaged multiple line-of-sight measurements at a constant elevation angle and four azimuth angles to calculate the 10 min mean wind speed at 40 range gates up to an altitude of about 1100 m. Reference measurement found the mean lidar error to be around 1 % with a standard deviation of 5 % (Gottschall, 2013).
The resulting wind speed is inherently spatially and temporally averaged.
At an altitude of 1100 m the radius of the averaging disk defined by the four azimuth positions with 90^{∘} increments is about 585 m. For the reconstruction of 10 min mean wind speed it is thus assumed that the wind vector does not change over this area, a valid assumption for these heights over flat terrain.

Lidar data availability highly depends on the applied carrier-to-noise ratio (CNR) filter and the aerosol content of the air as the wind speed is calculated based on the backscatter of the emitted laser beam. Most aerosols originate from the surface and are transported aloft. Particle density decreases with height and drops to almost zero within the free atmosphere above the ABL (Matthias and Bösenberg, 2002). Data quality quantified by the CNR dropped on average by approximately 5 dB over the course of 1000 m. A fixed CNR threshold of CNR${}_{\mathrm{dB}}\phantom{\rule{0.25em}{0ex}}>-\mathrm{25}$ dB, combined with additional self-defined filters (Sommerfeld et al., 2019), was applied and insufficient data were discarded. As a result, data availability dropped from about 81 % at 100 m to about 24 % at 1000 m. Low data availability caused by weather effects (e.g., strong precipitation) further emphasizes the importance of simulations for mid-altitude wind resource assessment as no measurement technique with sufficient spatial and temporal resolution is available at this point.

To complement the 6-month lidar data set, two WRF 3.6.1 simulations using the Advanced Research Weather Research and Forecasting (ARW) model (Skamarock and Klemp, 2008) were carried out. The “baseline run”, which is hereinafter referred to as NoOBS, is a 12-month study of the area around the measurement location (see Fig. 1) from 1 September 2015 used to derive annual statistics. Lidar measurements (Sommerfeld et al., 2019) were incorporated into the 6-month test model between September 2015 and February 2016 using OBSGRID (Wang et al., 2015), which is hereinafter referred to as OBS.

This methodology uses the difference between the model and measurements to calculate a nonphysical forcing term which is added to the governing conservation equations of the simulation to gradually nudge the model towards the observation (see Eq. 1) (Stauffer et al., 1991; Deng et al., 2007). Each simulation is composed of three nested domains with 27, 9 and 3 km grid spacing and horizontal grid dimensions of about 120×120 elements at 60 heights along the terrain-following vertical hybrid pressure coordinate *η*. Differences between the simulation runs (see Sect. 3.1) are compared within the innermost domain of the simulation. Output data were stored in 10 min intervals. Figure 1 shows the topography map of the simulation. Initial and boundary conditions of both simulations are based on the ERA-Interim (Dee et al., 2011) reanalysis data set by the European Centre for Medium-Range Weather Forecasts, which consists of 6-hourly atmospheric fields with a spatial resolution of roughly 80 km horizontally and 60*η* levels. Turbulent kinetic energy (TKE) closure within the ABL was achieved by using the Mellor–Yamada–Nakanishi–Niino (MYNN) 2.5 scheme, which predicts subgrid TKE as a prognostic variable (Nakanishi and Niino, 2004; Lee and Lundquist, 2017). The Noah-MP land-surface model and MYNN surface layer scheme were used. The Rapid Radiative Transfer Model (RRTM) longwave radiation and Dudhia shortwave radiation scheme were used (see Appendix A). In addition to observation nudging (see Sect. 3.1) analysis nudging was performed on every domain of each simulation. Analysis nudging nudges each grid point towards a time-interpolated value from gridded analyses of synoptic observations (Stauffer et al., 1991), whereas observation nudging directly drives the simulation towards the additional observations. Within the planetary boundary layer (PBL) of the inner domain analysis nudging was switched off (see nudging settings in Appendix A). All simulations were run on the EDDY^{1} high-performance computing clusters at the University of Oldenburg.

## 3.1 Observation nudging

Observation nudging, also referred to as “dynamic analysis”, is a form of four-dimensional data assimilation (FDDA) whereby each grid point within the radius of influence and time window is nudged towards observations using a weighted average of differences between the model (*q*_{m} interpolated at the observation location) and observations (*q*_{o}) (Dudhia, 2012; Reen, 2016). In this study horizontal wind speed *U* and direction Φ were nudged towards measurements with a time interval of 6 h between an altitude of 66 and 1100 m in order to not overly constrain the simulation. Nudging could not be performed at times and altitudes at which lidar data were not available. The nonphysical forcing term is implemented in the form of prognostic equations (Deng et al., 2007):

Here *q* refers to the quantity that is nudged, *μ* is the dry hydrostatic pressure, *F*_{q}(*x*, *y*, *z*, *t*) is the physical tendency term of *q*, *G*_{q} is the nudging strength of *q*, *N* is the total number of assimilated observations, *i* is the index of the current observation, and *W*_{q} is the weighting-function-based temporal and spatial separation between grid cells and observations (Dudhia, 2012). The four weighting functions *G*_{q}, *W*_{t}(*x*, *y*, *z*, *t*), *W*_{z}(*x*, *y*, *z*, *t*) and *W*_{xy}(*x*, *y*, *z*, *t*) describe the temporal and spatial nudging strength. Values used in this study can be found in Appendix A. The inverse of *G*_{q} can be interpreted as a nudging timescale as it dictates how quickly the model approaches the observation.

*W*_{xy} and *W*_{z} define the spatial nudging weight, while the temporal weighting function *W*_{t} defines the duration and weighting strength in time. *W*_{t} ramps from 0 to 1 and back to 0 (Reen, 2016). The nudging time window and the time between implemented observations was chosen to be 6 h so that the implemented observations do not overlap each other. This ensures all time steps are nudged while not excessively limiting the model.

Vertical influence was set very small so that observations only affect their own *η* level (Dudhia, 2012). The horizontal weighting factor *W*_{xy} (see Eq. 2) is calculated based on the radius of influence *R* and the distance between the observation and the grid location *D*. We used the Cressman scheme as the horizontal nudging weighting function with a radius of influence of *R*=180 km, thereby affecting the whole inner domain.

It is important to keep the differences in temporal and spatial resolution between lidar measurements and WRF simulations in mind. Furthermore, data availability highly influences the ability to nudge the simulation (see Sect. 2) and compare wind speed statistics.

To quantify the local effect of observation nudging, we investigate the cell closest to the lidar measurement location and compare measured and modeled horizontal wind speeds *U* and direction. Additionally, we investigate several sections at different locations and altitudes within the inner domain to quantify the spatial and temporal impact of single-location observation nudging on the entire domain. Vector values of each WRF cell are calculated on the faces of each cell, linearly interpolated to the cell center and rotated from the grid projection to the earth coordinate system.

## 4.1 Impact of nudging on wind statistics

Figure 2 shows the scatter plots of measured and simulated horizontal wind speed at various altitudes for times at which lidar data are available. The continuous line represents the linear regression of the data (the regression coefficient is displayed in the legend), while the dotted line shows an ideal correlation. The color of the scatter points corresponds to the frequency of occurrence. Multiple wind speed clusters caused by stratification can be identified. While there is a trend towards higher wind speeds with increasing altitude, low wind speeds (*U*<6 m s^{−1}) still occur at high altitudes. Both simulations overpredict horizontal wind speeds at low altitudes, which is a known problem of WRF and could be attributed to the model not resolving subgrid-scale roughness elements properly (e.g., modeling a strongly simplified parameterization of forests and/or cities) or flaws in the planetary boundary layer model; this could lead to overly geostrophic winds over land (Mass and Ovens, 2011). Observation nudging improves the overall correlation with measurements at the measurement location as the surface influence decays. Both models approach similar values at higher altitudes, which could be caused by the lack of observations, and therefore observation nudging due to reduced data availability, or is indicative of WRF generally being better at modeling more geostrophic winds.

The statistical analysis of the absolute difference between the WRF-simulated quantities at the measurement location and the lidar observations ($\mathrm{\Delta}U={U}_{\mathrm{WRF}}-{U}_{\mathrm{lidar}}$; $\mathrm{\Delta}\mathrm{\Phi}={\mathrm{\Phi}}_{\mathrm{WRF}}-{\mathrm{\Phi}}_{\mathrm{lidar}}$ wrapped on an interval [−*π*, *π*]) is shown in Fig. 3 in the form of a box plot. The circle corresponds to the median, the colored box indicates the 25th and 75th percentile, and the whiskers to both sides mark ±2.7 times the standard deviation (*σ*). Outliers beyond ±2.7*σ* are hidden to maintain clarity and readability. The continuous line in Fig. 3a represents the root mean square error between the measured *U*_{lidar} and simulated wind speed *U*_{WRF}.

The simulation with observation nudging generally outperforms the unnudged simulation and is in better agreement with the measurements, particularly at the altitudes of interest to high-altitude wind energy systems. It furthermore reduces the spread of the bias, illustrated by the smaller whiskers and boxes. The root mean square error (RMSE) Δ*U* shows similar results for both simulations below 100 m and above 700 m. The largest improvement or smallest error can be found between 300 and 600 m. This could be explained by a better performance of the mesoscale model at these altitudes due to a reduced impact of the air–surface interaction, which is strongly parameterized.

The NoOBS shows an almost constant wind direction bias at all altitudes. Observation nudging substantially reduces the directional bias ΔΦ up to high altitudes as can be seen in Fig. 3b. Similar to the wind speed bias, wind direction bias at 1100 m is almost the same for both simulations. The negative wind direction bias represents an anticlockwise deviation. Other studies (Carvalho et al., 2014; Giannakopoulou and Nhili, 2014) have found similar wind direction biases. A possible reason for this systematic error is that WRF does not adequately resolve surface roughness, resulting in lower surface friction and leading to overly geostrophic winds (Mass and Ovens, 2011). The almost constant median wind direction bias indicates that WRF is able to capture the clockwise rotation of the “Ekman spiral” in the Northern Hemisphere.

## 4.2 Representative nudging results

We compare 10 min mean horizontal wind speed for 24 h on 21 September 2015 in Fig. 4 to visualize the impact of observation nudging on the mesoscale model output. The white spaces in the lidar measurements (see Fig. 4a) are data points that have been filtered out due to insufficient data quality. The dashed line is the WRF-modeled surface heat flux (SHF) used to estimate atmospheric stability (see Sect. 4.5). The color of the profiles indicates the wind direction, and lidar-measured profiles are shown in grey for comparison. The black dot in each profile marks the altitude of highest wind speed, while the black circle indicates the optimal altitude for the operation of an airborne wind energy system based on a simplified power approximation (see Sect. 4.7). However, the single-point representation is only a rough measure of operational altitude since an AWES generally sweeps a range of altitudes.

Even though observation nudging leads to statistical improvements in wind speed and wind direction prediction over the entire period (compare Sect. 4.1 and 4.4), individual days can still show a decline in model accuracy. The low-level jet (LLJ) and high wind speeds at higher altitudes, which the NoOBS model captures fairly well, are significantly weaker in the OBS model. Implementing additional measurements at a higher frequency might yield results closer to measurements, but adding too many unphysical forcing terms might overly restrict the simulation.

The planetary boundary layer height (PBLH) (black line), which in the MYNN scheme is calculated from the profile of virtual potential temperature and from the profile of the TKE (Brunner et al., 2015; Nakanishi and Niino, 2004), is directly affected by wind speed observation nudging. During the investigated day, observation nudging leads to a lower daytime PBLH.

## 4.3 Spatial influence

Single-location observation nudging influences the area within the radius of influence (*R*_{xy}=180 km; see Appendix A), which here includes the entire inner domain (150 km × 150 km).
Figure 5 shows the mean absolute difference of horizontal wind speed ($\mathrm{\Delta}U=\left|{U}_{\mathrm{OBS}}\right|-\left|{U}_{\mathrm{NoOBS}}\right|$) between the OBS and NoOBS model along lines of constant longitude and latitude for the entire simulation period. The grid cell in which observations were assimilated is indicated by the vertical line and highlighted by the square marker. The four colors indicate different altitudes. As the outer domains remain unnudged, the boundary conditions of the inner domain remain the same, which leads to a rapid decline in absolute difference towards the outside of the domain. The difference in wind speed does not go to exact zero because the results are interpolated to the center of each grid cell. Near-surface results close to the measurement location, which is highlighted by the black vertical line, experience the largest change in wind speed (red line, *z*=12 m). The asymmetry could be caused by the downstream transportation of nudging effects (dominant wind direction: west).

## 4.4 Diurnal variability

Average diurnal variation indicates typical wind speed variations for a given location and period. It further reinforces the benefit of dynamically adapting operating altitudes of AWES. The hourly average lidar wind speed depends on data availability, as described in Sect. 2. Lidar availability below 100 m on average decreases by about 10 percentage points during the noon hours, while it remains fairly constant at altitudes between 100 and 300 m. Above this altitude, data availability increases in the afternoon by up to about 15 percentage points (Sommerfeld et al., 2019).

Figure 6a shows the lidar-measured and mesoscale modeled diurnal wind speed variation at the measurement location filtered by lidar availability; i.e., times when no lidar data were available were disregarded. A clear diurnal wind speed variation resulting from the cycle of stable and unstable stratification can be identified. On average OBS shows lower hourly wind speeds than NoOBS and is closer to measurements. The diurnal variation of the unfiltered 12-month NoOBS, 6-month OBS and 6-month NoOBS data sets (Fig. 6, right) deviates significantly from the measurements. Observation nudging leads to overall lower wind speeds and wind shear throughout the day in the unfiltered data set. Due to the large difference in average measured and unfiltered modeled diurnal wind speeds, it seems that lidar measurements alone cannot appropriately represent average wind conditions aloft due to availability bias, which has also been observed at other locations (Gryning and Floors, 2019). Therefore, we believe that the nudged data set yields more representative results than the unnudged model or the measurements alone.

## 4.5 Wind speed probability distribution

The common way to approximate the probability distribution of the horizontal wind speed *f*(*U*) is the Weibull distribution fit (Eq. 3), which describes the statistical distribution as a function of the scale parameter *A* and the shape parameter *k* (Troen and Lundtang Petersen, 1989).

Previous investigation of the lidar measurements showed a multi-modality in the wind speed frequency of occurrence caused by different atmospheric stability (Sommerfeld et al., 2019). Figure 7a visualizes the entire measured and simulated wind speed frequency distribution. Its corresponding Weibull fit is shown in Fig. 7b, and the difference between the two can be found in Fig. 7c. Each row summarizes the various data sets: first 6-month lidar, then 6-month OBS and 6-month NoOBS, followed by 12-month NoOBS.

All 6-month data sets show a high occurrence of low and high wind speeds, which indicates a multi-modal frequency distribution. This effect is most pronounced in the lidar data set. The comparison of wind speed frequency with the Weibull fit further emphasizes the multi-modality as a simple Weibull fit is not able to capture the higher probability at low and high wind speeds. These distinct flow situations further drift apart with increasing surface distance. As a result the Weibull distribution overestimates the occurrence of wind speeds between the two peaks. Both OBS and NoOBS slightly overestimate low-altitude wind speed (see Fig. 3) compared to lidar measurements. Both models and the lidar measurements show a broadening of the frequency distribution towards higher altitudes. High wind speeds become more likely, while low wind speeds still occur. Therefore, an AWES needs to be able to operate in a wide range of wind speeds or be controlled in a way that avoids extreme conditions. The 12-month NoOBS simulation shows lower wind speeds than the 6-month simulations as the included summer months generally have lower wind speeds due to the lower synoptic pressure gradients. The Weibull fit of this simulation tends to overestimate higher wind speeds and underestimate low wind speeds at all altitudes.

Using the sign of the WRF-calculated SHF as a simple proxy to differentiate stable and unstable wind conditions similar to Sommerfeld et al. (2019), the wind speed distribution follows the expected trends of low wind shear during unstable stratification and higher wind shear and wind speeds during stable stratification (Arya and Holton, 2001). Observation nudging reduces the occurrence of high wind speeds at high altitudes in comparison to NoOBS and leads to an increase in the probability of wind speeds around 5 m s^{−1} during times of positive SHF. The Weibull distribution fit of these sub-states is generally better at representing the modeled wind conditions.

Figure 8 shows the scale parameter *A*, shape parameter *k* and Hellinger distance *H* (Upton and Cook, 2008) between the wind speed frequency distribution and the corresponding Weibull distribution fit for lidar (first row), 6-month OBS (second row), 6-month NoOBS (third row) and 12-month NoOBS (fourth row).

The different trends under positive and negative SHF of both Weibull parameters visualize the existence of entirely different flow regimes.
The Hellinger distance between the Weibull fit and frequency distribution (negative SHF: blue, positive SHF: red), the total data and a simple fit (black), and between the total data and the weighted sum of both Weibull fits (green) is shown. All WRF models show an overall smaller *H* than a similar analysis of the lidar data set (Sommerfeld et al., 2019). The sharp bend in both *A* and *k* of the lidar data above 750 m is likely caused by insufficient data availability.
NoOBS results show a sharp increase in *A* up to 250 m and a slight reduction above, while OBS shows a trend close to the surface, and *A* values remain almost constant above 500 m. No data set shows a convergence of *A* at higher altitudes, indicating that these wind conditions are driven by different conditions in the free atmosphere. The 12-month NoOBS simulations show lower-scale parameter values as they include generally slower winds during summer. While *A* trends are quite different for lidar and WRF, *k* trends are more similar. They peak between 150 and 250 m and are especially high during stable stratification (Monahan et al., 2011). OBS trends of *k* are generally closer to measurement results than NoOBS.

Even though the Hellinger distance of individual Weibull fits for times of positive or negative SHF is generally higher than the Weibull fit of the entire data set, the weighted sum of both individual fits yields the best result at all altitudes. The 12-month Weibull fit using the entire data set performs comparably to a weighted sum up to an altitude of about 250 m.

## 4.6 Effect of stability on average wind shear

Atmospheric stability highly influences the shape of wind speed profiles, which is important for determining optimal operating conditions for an AWES (see Sect. 4.7). Obukhov length *L* (Obukhov, 1971; Sempreviva and Gryning, 1996) is commonly used to categorize the stability of the boundary layer. Here the application is extended to mid-altitudes. *L* is defined by the simulated friction velocity *u*_{*}, virtual potential temperature *θ*_{v}, potential temperature *θ*, kinematic virtual sensible surface heat flux *Q*_{S}, kinematic virtual latent heat flux *Q*_{L}, the von Kármán constant *k* and gravitational acceleration *g*. Table 1 summarizes the frequency of occurrence of each stability class.

In comparison with the unnudged simulation, OBS shows an increase in unstable and nearly unstable situations. Stable and nearly stable stratification seem almost unaffected by OBS nudging, while neutral and very stable stratification occur slightly less often. This might improve the overall predicting capabilities of WRF as the MYNN 2.5 boundary layer scheme overestimates the frequency of very stable conditions with an error of up to 9 % (Krogsæter and Reuder, 2015). Neutral conditions, still commonly used in many wind energy siting applications, only occur about 30 % of the time during the measurement period and only about 20 % of the time during the 1-year reference NoOBS simulation.

Figure 9 shows the frequency distribution of the different stability categories with the mean highlighted by white squares. All categories show distinct trends and distributions that are consistent between data sets, which contribute to the multi-modality of the overall wind speed frequency distribution. The difference in high-altitude wind speeds between stratifications indicates the influence of different geostrophic wind conditions. The categorization by *L* is based on surface data and seems to be valid within the lower part of the atmosphere where the spread of the corresponding frequency distribution is relatively small in comparison to high altitudes. This is particularly true for stable and neutral stratification whereby wind speeds above approximately 200 m spread widely. Unstable conditions are probably more consistent because of increased mixing from the surface up to high altitudes. The divergence of wind speeds towards higher altitudes indicate inhomogeneous atmospheric stability and suggests that surface-based stability categorization is insufficient for higher altitudes. Wind speed extrapolation based on low-altitude measurements can lead to a misestimation of mid-altitude wind conditions, especially during neutral and stable conditions close to the surface (Konow, 2015).

Altitudes below 200 m are least affected by observation nudging as OBS remains almost unchanged from NoOBS (see Sect. 4.1). Stable profiles show a peak at around 300 m, which is indicative of a characteristic low-level jet. Comparing OBS and NoOBS for 6 months, observation nudging seems to reduce the spread at higher altitudes within each category except very stable. The impact of observation nudging on wind profiles during unstable stratification is relatively low, while wind speed profiles under neutral and stable stratification are more affected.

## 4.7 Optimal operating altitude and power production

We estimate optimal operating altitude and traction power of a ground-generator AWES using a simple ground-generator (pumping-mode) AWES point-mass model adapted from Schmehl et al. (2013). We focus on 6-month OBS, as we previously proved increased accuracy, and use 12-month NoOBS to estimate annual values. The estimated optimal power per unit of lifting area of the wing *p*_{opt} is described by

Air density *ρ*_{air} is calculated by a linear approximation of the standard atmosphere (, ) (*ρ*_{air}(*z*)=1.225–0.00011*z*; kg m^{−3}). Losses associated with the mispositioning of the aircraft relative to the wind direction, expressed by azimuth angle *ϕ* and elevation angle *ε* relative to the ground station, are included in the model. Additional losses caused by gravity, tether sagging and tether drag are neglected. As a result, lift *F*_{L} and drag *F*_{D} force and therefore the lift (*c*_{L}=1.7) and drag coefficients (*c*_{D}=0.06), which are assumed to be constant, are geometrically related to the apparent wind velocity. Assuming an optimal tether speed and a quasi-steady state with the wing moving directly cross-wind with a zero azimuth angle (*ϕ*=0) relative to the wind direction, we can estimate the optimal traction power. The optimal elevation angle (*ε*_{opt}) and operating altitude (*z*_{opt}) are geometrically related to the assumed-constant tether length (*l*_{tether}) ($\mathrm{sin}{\mathit{\epsilon}}_{\mathrm{opt}}=\frac{{z}_{\mathrm{opt}}}{{l}_{\mathrm{tether}}}$).

Figure 10 summarizes the frequency of optimal operating altitude and optimal power assuming a constant tether length of 1500 m. The white solid line shows the cumulative frequency of optimal operating altitude. Both simulations for this particular location and time period show similar trends, with the most probable optimal altitude between approximately 200 and 400 m. Times of very high traction power are fairly rare and likely associated with low-level jets. Lower power at higher altitudes is caused by misalignment losses.

Figure 11 estimates the optimal traction power and operating altitude as a function of tether length based on the mean wind speed profile of atmospheric stability conditions (Fig. 9). The tether length of each estimation is assumed to be constant and used to calculate the optimal elevation angle. The axis limits of different atmospheric conditions had to be adjusted as the calculated power varied by orders of magnitude. All estimates show diminishing benefits of a longer tether. These incremental gains would probably be negated by additional drag- and weight-associated losses. Winds during times of very stable and unstable stratification lead to a clear optimal altitude independent of tether length between 200 and 400 m, while weakly stable and shear-driven wind speed profiles lead to higher optimal operating altitudes and a broader range of optimal altitudes as a function of tether length.

A full 6 months of lidar measurements up to 1100 m were assimilated into a mesoscale model using observation nudging. An unnudged reference model (NoOBS), the nudged model (OBS) outputs and lidar measurements were compared in terms of wind speed and direction statistics, as well as wind profile shape at the measurement site, and spatial differences were quantified. Observation nudging only has a marginal impact on simulated surface layer wind speeds as ground effects dominate the WRF model. Wind speeds between 300 and 500 m were most affected by observation nudging. Modeled wind speeds at these altitudes are statistically closest to the measurements, making this an adequate approach for resource assessment at mid-altitudes as measurement availability decreases. The impact of nudging weakens above these altitudes. Whether this is caused by lower measurement data availability or a generally better performance of the mesoscale model above the surface layer could not be determined. Observation nudging reduced the seemingly systematic wind direction bias between simulations and measurements at all altitudes. Due to the lack of high-resolution measurements at high altitudes, unnudged mesoscale model data are the best we have in terms of preliminary resource assessment.

Filtering the mesoscale model data according to lidar data availability yields similar diurnal variation, with OBS being closer to measurements. Comparing the diurnal variation of the unfiltered model wind speeds to measurements shows a significant deviation, which is likely caused by insufficient lidar data availability at higher altitudes. The bias between real and lidar-measured wind speed, which depends on the applied CNR threshold and data availability, can result in a misrepresentation of the actual wind conditions, especially at higher altitudes. Mesoscale models, particularly with observation nudging, can be used to account for this error. Lidar measurements seem to be biased towards high wind speeds as measured winds are generally higher than the unfiltered mesoscale model data. The impact of observation nudging on the wind profiles in the case of an unstably stratified boundary layer is relatively low, while wind speed profiles under stable stratification are significantly affected. At the measurement location OBS is overall closer to measurements, especially between 200 and 600 m. Variations of stratification, primarily those associated with the diurnal cycle, lead to a multi-modal wind speed frequency distribution which is better represented by the weighted sum of two Weibull fits than by a single Weibull fit. Obukhov-length-categorized wind speed profiles, especially during neutral and stable conditions close to the surface, show a divergence with height. This indicates inhomogeneous atmospheric stability and suggests that surface-based stability categorization is insufficient for higher altitudes.

Optimal AWES operating altitudes and power output per wing area were estimated based on a simplified model for 6 months of OBS and 12 months of NoOBS. The model neglects kite and tether weight as well as tether drag. Accounting for these losses, which are proportional to tether length, will reduce the performance of the AWES. Results for both wind speed data sets show the highest potential at an altitude between 200 and 600 m above which the losses associated with the elevation angle are too high. A comparison of different tether lengths under average wind speeds associated with different atmospheric stability conditions shows diminishing returns in terms of power output for tether lengths longer than 1500 m. While higher altitudes can potentially be reached, the optimal operating altitude remains almost unchanged. The highest energy potential and operating altitude is associated with neutral and stable stratification. Unstable conditions result in significantly lower energy potential due to lower, almost altitude-independent average wind speeds.

Future studies will include using the enhanced mesoscale model output to drive large-eddy simulations to provide better insight into mid-altitude turbulence. The resulting data set will lead to the development of a mid-altitude engineering wind model which can be used for the design, load estimation, control and optimization of airborne wind energy systems. Mesoscale model data will be implemented into an AWES optimization framework to quantify the impact of various wind speed profiles on power production, optimal trajectory and system size. Furthermore, the possibility of merging the mesoscale output with lidar measurements to fill gaps in the measurement data set to reduce the wind speed bias introduced by lidar availability is being investigated.

The data are not publicly available because they are subject to an NDA with the Fraunhofer IWES. Furthermore, the overall WRF data size makes up multiple terabytes.

WRF input parameter |
Value |

grid_fdda | 1, 1, 1, |

gfdda_inname | “wrffdda_d<domain>”, |

gfdda_end_h | 99 999, 99 999, 99 999, |

gfdda_interval_m | 360, 360, 360, |

fgdt | 0, 0, 0, |

if_no_pbl_nudging_uv | 0, 0, 1, |

if_no_pbl_nudging_t | 0, 0,1, |

if_no_pbl_nudging_q | 0, 0, 1, |

if_zfac_uv | 0, 0, 0, |

k_zfac_uv | 0, 0, 30, |

if_zfac_t | 0, 0, 0, |

k_zfac_t | 0, 0, 30, |

if_zfac_q | 0, 0, 0, |

k_zfac_q | 0, 0, 30, |

guv | 0.0003, 0.0003, 0.0003, |

gt | 0.0003, 0.0003, 0.0003, |

gq | 0.0003, 0.0003, 0.0003, |

if_ramping | 1, |

dtramp_min | 60.0, |

io_form_gfdda | 2, |

obs_nudge_opt | 0,0,1 |

Cressman scheme | 1 |

time_step | 60 |

obs_rinxy | 240, 240, 180 |

obs_rinsig | 0.1 |

obs_twindo | 3, 3, 3 |

auxinput11_interval_s | 360, 360, 360 |

obs_dtramp | 40 |

obs_nudge_wind | 1, 1, 1 |

obs_coef_wind | $\mathrm{6}\times {\mathrm{10}}^{-\mathrm{4}}$, $\mathrm{6}\times {\mathrm{10}}^{-\mathrm{4}}$, $\mathrm{6}\times {\mathrm{10}}^{-\mathrm{4}}$ |

iobs_onf | 2, 2, 2 |

auxinput11_interval_s | 360, 360, 360 |

auxinput11_end_h | 6, 6, 6 |

if_no_pbl_nudging_uv | 0, 0, 1 |

if_zfac_uv (max_dom) | 0, 0, 30 |

sf_sfclay_physics | 5, 5, 5 |

sf_surface_physics | 4, 4, 4 |

bl_pbl_physics (max_dom) | 5, 5, 5 |

bl_mynn_tkeadvect | .true., .true., .true. |

ra_lw_physics | 1, 1, 1 |

ra_sw_physics | 1, 1, 1 |

mp_physics | 5, 5, 5 |

MD and GS helped set up the numerical simulation and contributed to the meteorological evaluation of the data. MS evaluated the data and wrote the paper in consultation with and under the supervision of CC.

The corresponding author (Markus Sommerfeld) confirms on behalf of all authors that there have been no involvements that might raise the question of bias in the work reported or in the conclusions, implications or opinions stated.

The authors thank the Federal Ministry for Economic Affairs and Energy (BMWi) for funding of the OnKites I and OnKites II project (grant number 0325394A) on the basis of a decision by the German Bundestag and project management by Projektträger Jülich. The simulations were performed at the HPC Cluster EDDY, located at the University of Oldenburg (Germany), and funded by the Federal Ministry for Economic Affairs and Energy (BMWi) under grant number 0324005. We thank the PICS and the DAAD for their funding. We further thank all the technicians and staff at the Fraunhofer Institute for Wind Energy Systems (IWES) for carrying out the measurement campaign at Pritzwalk and their support in evaluating the data.

This research has been supported by the Federal Ministry of Economics and Technology (Germany) (grant no. 0325394A), the Federal Ministry for Economic Affairs and Energy (BMWi) (grant no. 0324005), and the Pacific Institute for Climate Solutions (PICS) (student scholarship project).

This paper was edited by Jakob Mann and reviewed by Roland Schmehl and Rogier Floors.

Al-Yahyai, S., Charabi, Y., and Gastli, A.: Review of the use of numerical weather prediction (NWP) models for wind energy assessment, Renew. Sustain. Energy Rev., 14, 3192–3198, https://doi.org/10.1016/j.rser.2010.07.001, 2010. a

Archer, C. L. and Caldeira, K.: Global Assessment of High-Altitude Wind Power, Energies, 2, 307–319, https://doi.org/10.3390/en20200307, 2009. a

Arya, P. and Holton, J.: Introduction to Micrometeorology, in: International Geophysics, Elsevier Science, available at: https://www.elsevier.com/books/introduction-to-micrometeorology/arya/978-0-12-059354-5 (last access: 3 October 2019), 2001. a

Bastigkeit, I., Gottschall, J., Gambier, A., Sommerfeld, M., Wolken-Möhlmann, G., and Rudolph, C.: Abschlussbericht-OnKites-Juni 2017_Final-5, detailled report AP1-AP2-AP5, Fraunhofer-Institut für Windenergie und Energiesystemtechnik IWES Nordwest, Bremerhaven, 2017. a

Bechtle, P., Schelbergen, M., Schmehl, R., Zillmann, U., and Watson, S.: Airborne wind energy resource analysis, Renew. Energy, 141, 1103–1116, https://doi.org/10.1016/j.renene.2019.03.118, 2019. a

Brunner, D., Savage, N., Jorba, O., Eder, B., Giordano, L., Badia, A., Balzarini, A., Baró, R., Bianconi, R., Chemel, C., Curci, G., Forkel, R., Jiménez-Guerrero, P., Hirtl, M., Hodzic, A., Honzak, L., Im, U., Knote, C., Makar, P., Manders-Groot, A., van Meijgaard, E., Neal, L., Pérez, J. L., Pirovano, G., San Jose, R., Schröder, W., Sokhi, R. S., Syrakov, D., Torian, A., Tuccella, P., Werhahn, J., Wolke, R., Yahya, K., Zabkar, R., Zhang, Y., Hogrefe, C., and Galmarini, S.: Comparative analysis of meteorological performance of coupled chemistry-meteorology models in the context of AQMEII phase 2, Atmos. Environ., 115, 470–498, https://doi.org/10.1016/j.atmosenv.2014.12.032, 2015. a

Burton, T. (Ed.): Wind energy handbook, 2nd Edn., Wiley, Chichester, West Sussex,, https://doi.org/10.1002/9781119992714, 2011. a

Canut, G., Couvreux, F., Lothon, M., Legain, D., Piguet, B., Lampert, A., Maurel, W., and Moulin, E.: Turbulence fluxes and variances measured with a sonic anemometer mounted on a tethered balloon, Atmos. Meas. Tech., 9, 4375–4386, https://doi.org/10.5194/amt-9-4375-2016, 2016. a

Carvalho, D., Rocha, A., Gómez-Gesteira, M., and Silva Santos, C.: WRF wind simulation and wind energy production estimates forced by different reanalyses: Comparison with observed data for Portugal, Appl. Energy, 117, 116–126, https://doi.org/10.1016/j.apenergy.2013.12.001, 2014. a

Cherubini, A., Papini, A., Vertechy, R., and Fontana, M.: Airborne Wind Energy Systems: A review of the technologies, Renew. Sustain. Energy Rev., 51, 1461–1476, https://doi.org/10.1016/j.rser.2015.07.053, 2015. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597, https://doi.org/10.1002/qj.828, 2011. a

Deng, A., Stauffer, D. R., Dudhia, J., Otte, T., and Hunter, G. K.: Update on analysis nudging FDDA in WRF-ARW, in: Proceedings of the 8th WRF Users' Workshop, 11–15 June 2007, National Center for Atmospheric Research, Boulder, Colorado, p. 35, 2007. a, b

Dudhia, J.: WRF Four-Dimensional Data Assimilation (FDDA), available at: http://www2.mmm.ucar.edu/wrf/users/tutorial/200801/WRF_FDDA_Dudhia.pdf (last access: 3 October 2019), 2012. a, b, c

Fagiano, L. and Milanese, M.: Airborne Wind Energy: An overview, in: IEEE 2012 American Control Conference (ACC), 27–29 June 2012, Montreal, QC, Canada, 3132–3143, https://doi.org/10.1109/ACC.2012.6314801, 2012. a

Fechner, U. and Schmehl, R.: Flight Path Planning in a Turbulent Wind Environment, in: Airborne Wind Energy: Advances in Technology Development and Research, edited by: Schmehl, R., Springer, Singapore, 361–390, https://doi.org/10.1007/978-981-10-1947-0_15, 2018. a

Floors, R., Batchvarova, E., Gryning, S.-E., Hahmann, A. N., Peña, A., and Mikkelsen, T.: Atmospheric boundary layer wind profile at a flat coastal site – wind speed lidar measurements and mesoscale modeling results, Adv. Sci. Res., 6, 155–159, https://doi.org/10.5194/asr-6-155-2011, 2011. a

Gambier, A., Bastigkeit, I., and Nippold, E.: Projekt OnKites II: Untersuchung zu den Potentialen von Flugwindenergieanlagen (FWEA) Phase II: Abschlussbericht (ausführliche Darstellung), Fraunhofer Institut für Windenergie und Energiesystemtechnik, Bremerhaven, https://doi.org/10.2314/GBV:1009915452, 2017. a

Giannakopoulou, E.-M. and Nhili, R.: WRF Model Methodology for Offshore Wind Energy Applications, Adv. Meteorol., 2014, 319819, https://doi.org/10.1155/2014/319819, 2014. a

Gottschall, J.: Galion Lidar Performance Verification, technical report, Fraunhofer-Institut für Windenergie und Energiesystemtechnik IWES Nordwest, Bremerhaven, available at: http://halo-photonics.com/Galion_LiDAR_system.htm (last access: 18 November 2016), 2013. a

Gottschall, J., Lindelöw-Marsden, P., and Courtney, M.: Executive summary of key test results for SgurrEnergy Galion, Executive summary, Technical University of Denmark DTU, Roskilde, available at: https://www.woodgroup.com/__data/assets/pdf_file/0024/15693/Riso_ExecutiveSummary_-Galion_lidar.pdf (last access: 18 November 2016), 2009. a

Gryning, S.-E. and Floors, R.: Carrier-to-Noise-Threshold Filtering on Off-Shore Wind Lidar Measurements, Sensors, 19, 592, 2019. a

IEC: 61400-1: Wind turbines part 1: Design requirements, International Electrotechnical Commission, International Electrotechnical Commission (IEC), Geneva, Switzerland, p. 177, 2005. a

ISO 2533:1975: Standard Atmosphere, Standard, International Organization for Standardization, Geneva, Switzerland, 1975. a

Konow, H.: Tall wind profiles in heteorogeneous terrain, PhD thesis, Universität Hamburg Hamburg, 2015. a

Krogsæter, O. and Reuder, J.: Validation of boundary layer parameterization schemes in the Weather Research and Forecasting model (WRF) under the aspect of offshore wind energy applications – Part II: Boundary layer height and atmospheric stability, Wind Energy, 18, 1291–1302, https://doi.org/10.1002/we.1765, 2015. a

Lee, J. C. Y. and Lundquist, J. K.: Observing and Simulating Wind-Turbine Wakes During the Evening Transition, Bound.-Lay. Meteorol., 164, 449–474, https://doi.org/10.1007/s10546-017-0257-y, 2017. a

Lunney, E., Ban, M., Duic, N., and Foley, A.: A state-of-the-art review and feasibility analysis of high altitude wind power in Northern Ireland, Renew. Sustain. Energy Rev., 68, 899–911, https://doi.org/10.1016/j.rser.2016.08.014, 2017. a

Mann, J.: The spatial structure of neutral atmospheric surface-layer turbulence, J. Fluid Mech., 273, 141–168, https://doi.org/10.1017/S0022112094001886, 1994. a

Mass, C. and Ovens, D.: Fixing WRF's High Speed Wind Bias: A New Subgrid Scale Drag Parameterization and the Role of Detailed Verification, American Meteorological Society, 24th Conference on Weather and Forecasting/20th Conference on Numerical Weather Prediction, 22–27 January 2011, Seattle, WA, available at: https://ams.confex.com/ams/91Annual/webprogram/Paper180011.html (last access: 3 October 2019), 2011. a, b

Matthias, V. and Bösenberg, J.: Aerosol climatology for the planetary boundary layer derived from regular lidar measurements, Atmos. Res., 63, 221–245, https://doi.org/10.1016/S0169-8095(02)00043-1, 2002. a

Monahan, A. H., He, Y., McFarlane, N., and Dai, A.: The Probability Distribution of Land Surface Wind Speeds, J. Climate, 24, 3892–3909, https://doi.org/10.1175/2011JCLI4106.1, 2011. a

Mylonas-Dirdiris, M., Barbouchi, S., and Herrmann, H.: Mesoscale modelling methodology based on nudging to reduce the error of wind resource assessment, in: Conference: European Geosciences Union General Assembly at, Vienna, Austria, 2016. a

Nakanishi, M. and Niino, H.: An Improved Mellor–Yamada Level-3 Model with Condensation Physics: Its Design and Verification, Bound.-Lay. Meteorol., 112, 1–31, https://doi.org/10.1023/B:BOUN.0000020164.04146.98, 2004. a, b

Obukhov, A. M.: Turbulence in an atmosphere with a non-uniform temperature, Bound.-Layer Meteorol., 2, 7–29, https://doi.org/10.1007/BF00718085, 1971. a

Optis, M., Monahan, A., and Bosveld, F. C.: Limitations and breakdown of Monin–Obukhov similarity theory for wind profile extrapolation under stable stratification, Wind Energy, 19, 1053–1072, https://doi.org/10.1002/we.1883, 2016. a

Peña, A., Gryning, S.-E., and Floors, R.: Lidar observations of marine boundary-layer winds and heights: a preliminary study, Meteorol. Z., 24, 581–589, https://doi.org/10.1127/metz/2015/0636, 2015. a

Reen, B.: A Brief Guide to Observation Nudging in WRF, Tech. rep., Army Research Laboratory, available at: http://www2.mmm.ucar.edu/wrf/users/docs/ObsNudgingGuide.pdf (last access: 3 October 2019), 2016. a, b

Sathe, A., Mann, J., Gottschall, J., and Courtney, M. S.: Can Wind Lidars Measure Turbulence?, J. Atmos. Ocean. Tech., 28, 853–868, https://doi.org/10.1175/JTECH-D-10-05004.1, 2011. a

Schmehl, R., Noom, M., and van der Vlugt, R.: Traction Power Generation with Tethered Wings, in: Airborne Wind Energy, chap. 2, Springer, Berlin, Heidelberg, 23–45, https://doi.org/10.1007/978-3-642-39965-7_2, 2013. a

Sempreviva, A. M. and Gryning, S.-E.: Humidity fluctuations in the marine boundary layer measured at a coastal site with an infrared humidity sensor, Bound.-Lay. Meteorol., 77, 331–352, https://doi.org/10.1007/BF00123531, 1996. a

Skamarock, W. C. and Klemp, J. B.: A time-split nonhydrostatic atmospheric model for weather research and forecasting applications, J. Comput. Phys., 227, 3465–3485, https://doi.org/10.1016/j.jcp.2007.01.037, 2008. a

Sommerfeld, M., Crawford, C., Monahan, A., and Bastigkeit, I.: LiDAR-based characterization of mid-altitude wind conditions for airborne wind energy systems, Wind Energy, https://doi.org/10.1002/we.2343, in press, 2019. a, b, c, d, e, f, g

Stauffer, D. R., Seaman, N. L., and Binkowski, F. S.: Use of Four-Dimensional Data Assimilation in a Limited-Area Mesoscale Model Part II: Effects of Data Assimilation within the Planetary Boundary Layer, Mon. Weather Rev., 119, 734–754, https://doi.org/10.1175/1520-0493(1991)119<0734:UOFDDA>2.0.CO;2, 1991. a, b

Troen, I. and Lundtang Petersen, E.: European Wind Atlas, Risø National Laboratory, Risø, 1989. a

Upton, G. and Cook, I.: A Dictionary of Statistics, Oxford University Press, Oxford, 2008. a

Wang, W., Bruyère, C., Duda, M., Dudhia, J., Gill, D., Kavulich, M., Keene, K., Lin, H.-C., Michalakes, J., Rizvi, S., Zhang, X., Berner, J., and Smith, K.: ARW Version 3.6 User's Guide, in: Chapter 7: Objective Analysis (OBSGRID), available at: http://www2.mmm.ucar.edu/wrf/users/docs/user_guide_V3.6/ARWUsersGuideV3.6.1.pdf (last access: 3 October 2019), 2015. a

Witha, B., Hahmann, A., Sīle, T., Dörenkämper, M., Ezber, Y., García-Bustamante, E., González-Rouco, J. F., Leroy, G., and Navarro, J.: WRF model sensitivity studies and specifications for the NEWA mesoscale wind atlas production runs, Technical report, The NEWA consortium, 73 pp., https://doi.org/10.5281/zenodo.2682604, 2019. a

^{1}

EDDY: HPC cluster at the Carl von Ossietzky Universität Oldenburg; see https://www.uni-oldenburg.de/fk5/wr/hochleistungsrechnen/hpc-facilities/eddy/ (last access: 3 October 2019).