the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Very shortterm forecast of nearcoastal flow using scanning lidars
Laura Valldecabres
Alfredo Peña
Michael Courtney
Lueder von Bremen
Martin Kühn
Wind measurements can reduce the uncertainty in the prediction of wind energy production. Today, commercially available scanning lidars can scan the atmosphere up to several kilometres. Here, we use lidar measurements to forecast nearcoastal winds with lead times of 5 min. Using Taylor's frozen turbulence hypothesis together with local topographic corrections, we demonstrate that wind speeds at a downstream position can be forecast by using measurements from a scanning lidar performed upstream in a very shortterm horizon. The study covers 10 periods characterised by neutral and stable atmospheric conditions. Our methodology shows smaller forecasting errors than those of the persistence method and the autoregressive integrated moving average (ARIMA) model. We discuss the applicability of this forecasting technique with regards to the characteristics of the lidar trajectories, the sitespecific conditions and the atmospheric stability.
 Article
(2166 KB)  Fulltext XML
 BibTeX
 EndNote
Wind energy is growing worldwide as a major source of green energy. In 2017 Denmark produced a record 43.4 % of the country's electricity with wind energy (Danish Wind Industry Association(DWIA), 2018). As the share of variable energy into the grid grows, more effort is required to increase the flexibility of power systems in such a way that they can guarantee the grid stability (Holttinen et al., 2016) and the reliability of energy supply (Ibanez and Milligan, 2012). For power systems with a high penetration of intermittent renewable energy, one of the most important sources of imbalances is wind energy forecast errors (GonzalezAparicio and Zucker, 2015). On short timescales, transmission system operators maintain the balance between electricity production and demand, activating balancing reserves. Due to their high flexibility to respond to shortterm changes in power, balancing reserves mostly come from conventional power plants, which reduces the environmental and economic benefits of wind energy. In countries like Belgium, the Netherlands and Germany, the electricity market participants can submit their intraday bids until 5 min before delivery (EPEXSPOT, 2017). With potentially shorter gate closure times wind power suppliers can better match production with demand, thus minimising the costs arising from the deviation between scheduled wind energy production and real generation (Wang et al., 2016).
In very shortterm horizons, i.e. from minutes to 1 h, wind forecasts are normally based on statistical models. They are built on relationships developed among historical measurements, assuming that these relationships are also applicable in the future. Examples of statistical methods used to predict wind speed and power can be found in Hill et al. (2012) for the autoregressive (AR) model, Torres et al. (2005) for the autoregressive moving average (ARMA) model and Kavasseri and Seetharaman (2009) for the autoregressive integrated moving average (ARIMA) model. Torres et al. (2005) applied the ARMA model to predict hourly average wind speeds at five weather stations in Navarre, Spain, during different times of the year, with a forecasting horizon from 1 to 10 h. They showed smaller errors for ARMA models compared to those of the persistence method. The classical persistence model predicts that the future value will be the same as the current value. This is the simplest version of the ARMA model and is often considered to be a benchmark for other forecasting techniques (Giebel et al., 2011).
Another forecasting technique is the spatial correlation method, which uses the wind speed at upstream neighbouring points to predict the wind speed at a downstream location. The first spatial correlation models were developed by Schlueter et al. (1986). They predicted meteorological events based on crosscorrelation curves of wind speeds at two sites using a constant delay method. Alexiadis et al. (1998) tested this method in the Greek islands of Syros and Paros in a time horizon of 10 min to some hours. Although there was a high correlation between the two sites in terms of fluctuations, the errors were higher than those of the persistence method both in magnitude and phase. They later proposed a spatial correlation predictor method, which uses linear relations to correct magnitude and phase errors.
Recently, various techniques based on artificial neural networks, which are trained with large historic data from the location, were developed (Cadenas and Rivera, 2009; Monfared et al., 2009; Li and Shi, 2010). Damousis et al. (2004) implemented a wind forecasting fuzzy model in which wind data from neighbouring meteorological stations at a radius of up to 30 km were used to predict wind speed and power in horizons of 30 min to 2 h. The model results showed significant improvement in the forecasting error of wind speed and power compared to those of the persistence model, especially when applied on flat terrain.
Wind forecasting techniques can combine physical and statistical approaches. As an example, Larson and Westrick (2006) used offsite observations in the vicinity of a wind farm in northeastern Oregon as input variables in different forecast models such as neural networks and support vector machines. They showed that the integration of realtime offsite observations significantly improves the forecasting accuracy of those algorithms.
Today, remotesensing systems like lidars are intensively being deployed for wind resource assessment (Wharton et al., 2015), turbine control (Mikkelsen et al., 2013) and turbulence characterisation (Peña et al., 2017). Lidars are proven to be relevant for very shortterm forecasting (Frehlich, 2013) as the current generation of commercially available units can scan in various atmospheric conditions up to 30 km. As an example, a 4 m s^{−1} wind speed could be observed by a lidar located in a wind farm 3.6 km upstream and 15 min ahead, thus predicting the start of power generation. Remotesensing systems could also be used to better schedule maintenance of offshore wind farms (Barthelmie et al., 2008), e.g. during periods of low wind speeds.
With our study we want to (i) experimentally investigate how lidar observations can be used to forecast wind speeds in a very shortterm horizon assuming Taylor's frozen turbulence hypothesis (Taylor, 1938) and (ii) test if with the use of lidar measurements we can predict wind speeds better than with the benchmarks ARIMA and the persistence model. For this, we use lidar observations up to 6 km in a nearcoastal area in Denmark as an input for an advectionbased wind speed forecast technique. The observations are characterised by rather high wind speeds, which limit the forecasting horizon to 5 min. The lidar measurements are described in Sect. 2.1. An insight into the wind conditions is given in Sect. 2.2. The location of the lidars in the nearcoastal area made it necessary to consider the topographic local conditions, which are modelled in Sect. 3. Section 4 gives a detailed description of the methodology used to forecast wind speeds using the lidar measurements. In Sect. 5 results are presented through comparisons between the accuracy of forecasting wind speeds based on the advection models and with persistence and ARIMA models. We discuss the suitability of using longrange lidars for very shortterm forecasting and provide main conclusions in Sect. 6.
Our study is based on measurements performed during the Reducing Uncertainty of Nearshore wind resource Estimates (RUNE) campaign (Simon and Courtney, 2016; Floors et al., 2016b). The experiment was conducted at the western coast of Denmark, north of the area of Høvsøre (see Fig. 1) and close to one of the Technical University of Denmark's wind turbine test stations. A comprehensive analysis of the wind conditions at Høvsøre during a 10year period from the test station's meteorological mast, located 1.7 km east of the North Sea (see Fig. 1, position 8), is presented in Peña et al. (2016). A pronounced cliff at the coastline (see Fig. 2) is the main feature of the terrain, which is mainly covered with grass and crops.
During the RUNE campaign, which took place during November 2015 to February 2016, profiling and scanning lidars were deployed to measure nearcoastal wind conditions (see Fig. 1). Four shortrange (positions 2, 4, 5 and 6 (later 7)) and one longrange (position 2) profiling lidar measured the wind profile. One scanning lidar (position 2) was operated in plan position indicator (PPI) mode, also known as the “sectorscan” scenario. Simultaneously, two more scanning lidars (positions 1 and 3) were configured in a dual trajectory to match at positions along three horizontal virtual lines. In what follows we will refer to them as the dual setup. In Fig. 2 the positions of the dual setup and the PPI are shown. The PPI and the dualsetup trajectories were designed so that the measurements will intersect at 5000 m offshore at 50, 100 and 150 $\mathrm{m}\phantom{\rule{0.125em}{0ex}}\mathrm{a}.\mathrm{m}.\mathrm{s}.\mathrm{l}.$ Further, a directional wave buoy (position 9) was deployed to measure waves, currents and sea surface temperature (SST) (Floors et al., 2016a; Sanchez and Rørbæk, 2016). Detailed information about the campaign can be found in Floors et al. (2016b). For this work we also use data from the sonic and cup anemometers located at the height of 100 m on the Høvsøre meteorological mast (position 8). Table 1 summarises the operational availability of all systems used in this study.
2.1 Lidar data processing and filtering
2.1.1 Dualsetup measurements
The two lidars measuring in the dualsetup trajectory acquired 45 lineofsight (LOS) wind velocities (1 s per LOS) per horizontal virtual line, separated by a distance of ≈ 200 m between points from 4 km onshore to 5 km offshore. Every trajectory, i.e. three horizontal lines at 50, 100 and 150 $\mathrm{m}\phantom{\rule{0.125em}{0ex}}\mathrm{a}.\mathrm{m}.\mathrm{s}.\mathrm{l}.$, took 145 s. In total, every position was swept four times every 10 min. Data that did not fulfil a certain distance threshold between the two lidar measurement positions were discarded. Regarding data quality, a carriertonoise ratio (CNR) threshold of − 26.50 dB was set. For every 10 min period and each point, the horizontal wind speed components were reconstructed as described in Simon and Courtney (2016). Due to the low availability in the reconstruction at positions further away from the coast (Floors et al., 2016b), we only consider data up to 2950 m. Observations close to the lidar systems (range < 500 m) were also discarded since here the angle between the beams approaches 180^{∘} and, consequently, the uncertainty of the reconstructed speed becomes very high (Stawiarski et al., 2013).
2.1.2 PPI measurements
The lidar at position 2 measured 45 different azimuthal positions over three different elevations, performing a 60^{∘} sweep every 45 s, scanning in the westerly direction (240–300^{∘}). The elevation angles were 0.27, 0.84 and 1.41^{∘}. The full trajectory lasted 145 s accounting for the 10 s that the scan needed to return to its initial position. For every azimuthal position, 156 range gates from 100 to 8150 m (separated every 50 m) were measured. The horizontal wind speed was reconstructed for every single scan and range gate, resulting in a horizontal wind speed at each range gate and elevation every 145 s; thus four measurements were performed within a 10 min period. Due to the low availability of data at long ranges when using a filtering threshold (Floors et al., 2016b), a dynamic filter is applied to “rescue” LOSs as shown in Beck and Kühn (2017). For every 10 min period, the probability density function of the data is calculated using a 2D histogram. Measurement points fulfilling a lower threshold of − 26.5 dB and an upper threshold of − 5 dB are considered. LOSs below the CNR lower threshold are still considered, if their local probability density lies within 1 standard deviation of the mean probability density. LOS measurements below − 30 dB are always discarded. A final visual checking is applied to remove outliers. In Fig. 3, a comparison of the two filtering techniques is presented. As shown in the range–CNR plot, for these data the use of a dynamic filter extends the range of measurements from 4.6 to 6 km. The availability of LOS measurements for the two filtering techniques is shown in Fig. 4. For a distance of 6000 m from the coast, the use of the dynamic filter increases the data availability from 33.65 to 73.29 %.
A comparison of the wind speeds observed by the dual setup and the PPI at their matching positions 5 km offshore can be found in Floors et al. (2016b). In general, the 10 min mean reconstructed wind speeds from the PPI show a good agreement with the dualsetup ones, especially close to the coast. At further distances from the coast, higher mean differences are found. These are related to the different size of the measuring volume of the lidars, the inherent temporal and spatial variability in the wind speed, and the distinct reconstruction methods. While the reconstruction in the PPI is performed with a sinusoidal fit of 60^{∘}, the dual setup uses two LOSs from the two lidars at a similar position in space. The uncertainties arising from the nature of the two systems are not clearly addressed. For the PPI, we need to assume horizontal flow homogeneity. At distant ranges wider areas are covered, and there is a higher uncertainty in the reconstruction. In addition, the PPI trajectories are not horizontal. For the dual setup, we do not assume horizontal flow homogeneity, but the measurement ranges are longer than in the PPI and, consequently, the uncertainty in the sensing height will be higher.
2.2 Observed wind conditions
The campaign was characterised by strong southwesterly winds. For subsequent analysis we want to estimate the atmospheric stability conditions during the campaign. Since there are no measurements of this type in the offshore area, we estimate the offshore stability based on sonic anemometer measurements from the Høvsøre met mast. We select the highest sonic anemometer, at 100 m, since this is less influenced by the land effects and by internal boundary layers growing during westerly winds. The analysis based on the derived Obukhov length L at the 100 m sonic anemometer reveals that during the winter months, there were mostly stable conditions (56.5 %) followed by neutral (27.8 %) and unstable conditions (15.7 %). Three classes are used for the stability classification, with z∕L < −0.1 for unstable, −0.1 $\le z/L\le $ 0.1 for neutral and z∕L > 0.1 for stable conditions. To test if we can estimate the offshore stability based on the onshore measurements, we conduct a comparison of the gradient of the potential temperature between the sea and the air, and the L estimated from the sonic measurements. The directional wave buoy located at position 9 measured the SST every 30 min. Due to a major failure in the buoy system, only measurements until the beginning of January 2016 were recorded (see Table 1 for more details). A comparison with the SST derived from satellite images is shown in Fig. 5. The SST was computed from nighttime observations from NOAA, AVHRR, Metop AVHRR, Terra MODIS, Aqua MODIS, Aqua AMSRE, Envisat AATSR and MSG SEVIRI satellites based on the interpolation method described in Høyer and She (2007). The spatial resolution of the satellite SST is 0.02^{∘} × 0.02^{∘} and its temporal resolution is 24 h. Figure 5 shows that the differences between both SSTs are small and both SSTs do not vary as much as the air temperature measured by the temperature sensor at 100 m on the mast.
To conduct the comparison with the stability from the met mast, we calculate the potential temperature gradient between the sea surface (buoy) and the air (met mast 100 m) for every 30 min period, first clustering the data according to the wind direction and second using the 30 min averaged L (sonic at 100 m). Only westerly and easterly sectors were analysed. The sign of the mean gradient of the potential temperature between the sea surface and the air for westerly winds (see Fig. 6) is in agreement with the stability from the sonic anemometer at 100 m. For easterly winds there is no such correspondence, as expected. For westerly winds we will assume that the stability measured by the onshore metmast at 100 m is a good indicator of the stability of the offshore area. Our further analysis refers to data from westerly winds during neutral and stable conditions.
2.2.1 Coastal gradient for westerly winds
We analyse the influence of the land on the wind speed in the nearcoastal area by using the dualsetup lidar observations at offshore positions. Due to the reduced availability of measurements at distant positions, we look at 10 min periods up to 3 km offshore. Figure 7 shows the ensemble average wind speed of all 10 min mean wind speeds reconstructed from the dualsetup observations at 50, 100 and 150 $\mathrm{m}\phantom{\rule{0.125em}{0ex}}\mathrm{a}.\mathrm{m}.\mathrm{s}.\mathrm{l}.$ for periods with neutral stratification. For all heights, the flow slows down when approaching the coast.
We will use the PPI measurements further upstream of the coast to forecast winds at positions close to the coast where we also have PPI measurements. Our forecasting technique is first based on an advection component, in which it is assumed that large turbulent structures are advected with the mean wind. Second we need to vertically extrapolate the wind because the upstream PPI observations are at different heights than those closer to the coast. Last we need corrections due to the influence of the coast; as seen in Fig. 7 the wind has been observed as decreasing as it approaches the coast. Here, we will first show the method used to account for the coastal effects. This is carried out based on the dualsetup measurements as they are independent of PPI scans and are always performed at the same heights. In this section we will only use dualsetup measurements up to 3 km during neutral conditions.
For a homogeneous and stationary flow, the mean wind speed profile is given as
where U is the mean wind speed, z the height above the ground, u_{∗} the friction velocity, κ the von Kármán constant (≈ 0.4) and z_{0} the roughness length. To account for stability effects Ψ is included, which depends on the Obukhov length L. To model the effects of the orography and roughness on the wind, which depend on the distance to the coast, we assume that the observed (obs) wind speed is
where O is an orography correction that depends on the height z, the distance to the coast x and the wind direction dd. Note that we assume that z_{0} varies with the distance to the coast.
3.1 Orography effects
The orography effects are estimated using the microscale IBZ model, which is part of the Wind Atlas Analysis and Application Program (WAsP) (Troen and Lundtang Petersen, 1989). The orography correction was determined at each position measured by the dual setup and for all wind directions using a digital terrain model as an input (Geostyrelsen, 2016). In Fig. 8 the orography corrections for the positions 500 and 2950 m from the coast and for all wind directions are shown. For westerly winds, at 500 m from the coast and 50 $\mathrm{m}\phantom{\rule{0.125em}{0ex}}\mathrm{a}.\mathrm{m}.\mathrm{s}.\mathrm{l}.$, the wind speed slows down by ≈ 2 %. Conversely, for northerly and southerly winds the flow speeds up due to the presence of the cliff. The effects are reduced further from the coast and with increasing height.
3.2 Roughness effects
We model the influence of the wind on the roughness of the water using the expression of Charnock (1955),
where α_{c} is the Charnock parameter and g the acceleration due to gravity. For open ocean α_{c} = 0.011 has been reported (Smith, 1980) while for a nearcoastal area, values between 0.008 and 0.06 can be found (Kraus, 1972). To determine the roughness length dependency with distance to shore, we apply the following strategy. Once the dualsetup observations at the different range gates are corrected by using the orography corrections, these are used together with Eqs. (1) and (3) to determine both u_{*} and z_{0}, and thus α_{c}. Figure 9b left shows the dependency of the estimated roughness length with distance to the coast after applying the orography corrections for the neutral cases. The roughness length decreases with distance from the coast. Without orography corrections, the roughness length is slightly higher than the case with corrections close to the coast, as expected. Since the roughness length varies with distance to the coast, so does the Charnock parameter (see Fig. 9b).
We test the estimated Charnock parameter dependency with distance to shore by selecting 10 min periods during neutral conditions for which both PPI and dualsetup measurements were performed simultaneously. We fit Eqs. (1) and (3) to the orographycorrected PPI measurements at those positions where we estimated the α_{c} dependency on distance to the coast with the dualsetup measurements. The comparison of the estimation of the wind using Eqs. (1) and (3) is shown in Fig. 10. As shown, with increasing distance to the coast, there is an increasing deviation of the fit from Eqs. (1) and (3) to the data, especially at the lowest height observations.
As mentioned earlier, we want to forecast wind speeds in a very shortterm horizon by assuming Taylor's frozen turbulence hypothesis. For this purpose we consider two positions: the upstream position (1) and the downstream position or forecasting position (2), with the wind blowing from (1) to (2). If at a time t a considerable change in wind speed occurs at position (1), this event will appear at position (2) after some time Δt. In other words, this event can be foreseen at position (2) with a time ahead Δt. In our analysis, the downstream position is set to 500 m from the PPI lidar (position 2) in the westerly direction at z_{2} = 33.76 m, which corresponds to the height of the intermediate PPI elevation scan. Lidar measurements are performed at multiple upstream positions (range gates) from which the forecast can be originated. This can be understood as having multiple virtual met masts over several distances west from the downstream position. To keep a fixed forecast horizon, the upstream position (1) and height z_{1}, from which the wind is advected, are determined dynamically at each time stamp using the 5 min movingaverage wind speed $\stackrel{\mathrm{\u203e}}{{v}_{\mathrm{2}}}\left(t\right)$ and direction at the downstream position. But because the vector $\stackrel{\mathrm{\u203e}}{{v}_{\mathrm{2}}}\left(t\right)$ might not be parallel to the line of virtual met masts, we use the vector projection of the advected distance on the wind direction $\left{\mathit{r}}_{\mathrm{12}}\right=\left\mathrm{\Delta}t\stackrel{\mathrm{\u203e}}{{\mathit{v}}_{\mathrm{2}}}\right(\mathit{t}\left)\right\cdot \text{cos}\left(\mathit{\theta}\right)$, with θ defined as the angle between the wind direction and 270^{∘}. Because high wind speeds were observed during the measurement campaign, and the limit for highquality PPI measurements is ≈ 6 km, we establish a forecast horizon of 5 min. We assume that a change in wind speed, observed 5 min ahead at the position (1) will propagate and travel the distance r_{12} in the time Δt = 5 min.
To incorporate the local effects and the changes in height between the upstream and downstream positions, we consider corrections in the wind speed due to height, roughness and orography. To evaluate the appropriateness of those corrections, we will compare a simple advection model (A), an advection model with height corrections (AH), an advection model with height and roughness length corrections (AHR), and, finally, an advection model with corrections due to height, roughness length change and orography changes (AHRO). We evaluate our forecasting method against the wellknown persistence method and an ARIMA model. A summary of the time periods in which the very shortterm forecasting method is applied is shown in Table 2. We select periods with mean wind speeds below 18 m s^{−1} and westerly periods, with a minimum duration of 3 h and with high availability of the data. No unstable periods fulfilled this criteria; therefore we focus here on neutral and stable conditions.
4.1 Advection model (A)
For the advection model, U_{2}(t) is estimated as follows.

The upstream position (1) at (x_{1},z_{1}) from the PPI scan is determined dynamically using the 5 min movingaverage wind speed at the downstream position and the forecast time horizon k (here 5 min). The wind direction from the previous forecasted step is used to calculate the projected distance from which the forecast is originated. For the positions in time and space in which observations at the upstream position are missing, the previous observation is used.

The observed wind speed at the upstream position is therefore advected, which means that the forecasted wind speed at the downstream position is considered to be the same as the wind speed in the upstream position, ${U}_{\mathrm{2},z\mathrm{2}}\left(t\right)={U}_{\mathrm{1},z\mathrm{1}}(tk)$.
4.2 Advection model with height correction (AH)
This is similar to the A model but the wind speed is extrapolated to match the height of the downstream observation. To do so, the following steps are carried out.

Step 1 from model A is conducted.

The logarithmic profile in Eq. (1) is fit to three consecutive PPI wind speed observations at position (1). The friction velocity ${u}_{*,\mathrm{1}}$ and a roughness length z_{0,1} are thus estimated.

The roughness length z_{0,1} is used to correct the advected wind speed to the downstream height by using
$$\begin{array}{}\text{(4)}& {\displaystyle}{\displaystyle}{U}_{\mathrm{2},z\mathrm{2}}\left(t\right)={U}_{\mathrm{1},z\mathrm{1}}(tk){\displaystyle \frac{\text{ln}({z}_{\mathrm{1}}/{z}_{\mathrm{0},\mathrm{1}}(tk\left)\right)}{\text{ln}({z}_{\mathrm{2}}/{z}_{\mathrm{0},\mathrm{1}}(tk\left)\right)}}.\end{array}$$
4.3 Advection model with height and roughness correction (AHR)

Steps 1 and 2 from the AH model are conducted.

The friction velocity ${u}_{*,\mathrm{1}}$ and the roughness length z_{0,1} are used to calculate the geostrophic wind at position (1):
$$\begin{array}{}\text{(5)}& {\displaystyle}{\displaystyle}{G}_{\mathrm{1}}\left(t\right)={\displaystyle \frac{{u}_{*,\mathrm{1}}\left(t\right)}{\mathit{\kappa}}}\sqrt{{\left(\text{ln}\left({\displaystyle \frac{{u}_{*,\mathrm{1}}\left(t\right)}{f{z}_{\mathrm{0},\mathrm{1}}}}\right)A\right)}^{\mathrm{2}}+{B}^{\mathrm{2}}},\end{array}$$where f refers to the Coriolis parameter and A=1.8 and B=4.5. We assume here that the geostrophic wind at position (2) is the same as at position (1):
$$\begin{array}{}\text{(6)}& {\displaystyle}{\displaystyle}{G}_{\mathrm{2}}\left(t\right)={G}_{\mathrm{1}}(tk).\end{array}$$ 
The geostrophic wind is used to estimate the roughness length z_{0,2} and the friction velocity ${u}_{*,\mathrm{2}}$. To solve for both parameters in position (1), we assume a fixed Charnock parameter derived from the dualsetup analysis (see Fig. 9b):
$$\begin{array}{ll}{\displaystyle}& {\displaystyle}{G}_{\mathrm{2}}\left(t\right)={\displaystyle \frac{{u}_{*,\mathrm{2}}(tk)}{\mathit{\kappa}}}\\ \text{(7)}& {\displaystyle}& {\displaystyle}\sqrt{{\left(\text{ln}\left({\displaystyle \frac{g}{f{\mathit{\alpha}}_{c,\mathrm{2}}{u}_{*,\mathrm{2}}(tk)}}\right)A\right)}^{\mathrm{2}}+{B}^{\mathrm{2}}},\text{(8)}& {\displaystyle}& {\displaystyle}{z}_{\mathrm{0},\mathrm{2}}={\mathit{\alpha}}_{c,\mathrm{2}}{\displaystyle \frac{{u}_{*,\mathrm{2}}^{\mathrm{2}}}{g}}.\end{array}$$ 
The forecasted wind speed is
$$\begin{array}{}\text{(9)}& {\displaystyle}{\displaystyle}{U}_{\mathrm{2},z\mathrm{2}}\left(t\right)={\displaystyle \frac{{u}_{*,\mathrm{2}}(tk)}{\mathit{\kappa}}}\text{ln}\left({\displaystyle \frac{{z}_{\mathrm{2}}}{{z}_{\mathrm{0},\mathrm{2}}(tk)}}\right).\end{array}$$
4.4 Advection model with height, roughness and orography correction (AHRO)

It is assumed that the corrections due to orography at positions further away are negligible.

Steps 1–3 from the AHR model are conducted.

The orography corrections at the downstream position are applied using the measured wind direction at (1), i.e.
$$\begin{array}{}\text{(10)}& {\displaystyle}{\displaystyle}{U}_{\mathrm{2},z\mathrm{2}}^{o}\left(t\right)={U}_{\mathrm{2},z\mathrm{2}}\left(t\right)O({x}_{\mathrm{2}},{z}_{\mathrm{2}},{\text{dd}}_{\mathrm{1}}).\end{array}$$
4.5 Statistical models
To evaluate the goodness of the forecasting techniques in Sect. 4, we use the benchmarks persistence and ARIMA.

ARIMA, denoted as ARIMA$(p,d,q)$ is a statistical model widely used in very shortterm predictions of wind speeds (Kavasseri and Seetharaman, 2009). It uses recent past values (autoregressive, AR) and recent residuals of the forecast (moving average) to predict current values. This model is suitable to analyse nonstationary processes since it uses nonseasonal differences (d) to build the forecasting model. A general equation for the ARIMA model is
$$\begin{array}{ll}{\displaystyle}U(t+k)& {\displaystyle}=\sum _{i=\mathrm{1}}^{p}{\mathrm{\Phi}}_{i}U(t+ki)+a(t+k)\\ \text{(11)}& {\displaystyle}& {\displaystyle}\sum _{j=\mathrm{1}}^{q}{\mathrm{\Theta}}_{j}a(t+kj),\end{array}$$where Φ_{i} is the ith autoregressive parameter, Θ_{j} is the jth movingaverage parameter, a(t) is the error term at time t, k is the forecasting horizon and U(t) is the value of the wind speed observed at the time t. Here we build a new ARIMA model for each period. To test the stationarity of the time series, we first look at its autocorrelation function and its partial autocorrelation function. If it has positive autocorrelations out to a high number of lags, we include an order of differentiation d. To test if this order is sufficient, we look at the residuals of the differentiated time series and perform a unit root test using the Dickey–Fuller test (Dickey and Fuller, 1979). To determine the order p and q of the ARIMA model, we compute the autocorrelation function and the partial autocorrelation function of the stationarized time series, following the method explained in Cadenas and Rivera (2007). The model chosen for every period is the one which minimises the residuals. For every individual set, the previous hour of observations is used to derive the AR and moving average parameters using the method by Box and Jenkins (1976).

Persistence is a particular case of the ARIMA model in which q=0, p=1, d=0 and the AR coefficient is set to 1 since it assumes that the previous and the current values are highly correlated. Our predicted wind speed is defined as
$$\begin{array}{}\text{(12)}& {\displaystyle}{\displaystyle}U\left(t\right)=U(tk).\end{array}$$
We evaluate the accuracy of the 5 min forecast of wind speeds based on the described advection techniques and compare it with the results of the statistical methods persistence and ARIMA. To do so, three criteria are employed, namely the rootmeansquare error (RMSE), the mean bias error (MBE) and the maximum absolute error (MaxAE). Table 3 includes the RMSE, MBE and MaxAE for all periods stated in Table 2. Minimum values are indicated in bold. The improvements of the advection models over the benchmarks persistence and ARIMA are shown in Table 4. Values corresponding to best performance are indicated in bold.
For neutral conditions (periods 2, 3, 4, 5, 7 and 8), the advection model with corrections performs in general better than the statistical forecasting models. The improvement over persistence using the best calibrated advection model for each period ranges from 21 to 38 % (see Table 4). Compared to the benchmark ARIMA the improvement ranges from 4 to 28 %. As an example, the distribution of errors produced by all models for period 7 can be seen in Fig. 11. The forecasting error ϵ is defined as ${\mathit{\u03f5}}_{i}={U}_{\mathrm{p},i}{U}_{\text{ob},i}$, where U_{ob,i} is the actual observation for a time position t_{i} and U_{p,i} is the forecast for the same period. The statistical methods show a broader distribution of errors. This is because ARIMA and persistence fail to predict the phase of the events since they construct their predictions according to the previous observations.
The dependency of the forecasting errors on the mean wind speed of the downstream observation for all advection models is shown in Fig. 12. From there it can be inferred that the local topographic corrections are required since the AH model overestimates wind speeds. For the neutral periods 4, 5, 7 and 8, the forecasting accuracy of the AHR model is higher than that of any other advection model. In those periods, introducing the orography correction results in an underestimation of the wind speed, as can be seen in the MBE of those periods. For wind speeds close to 16 m s^{−1} and neutral conditions, AHRO produces smaller errors (Fig. 12a). Therefore for period 2, which has a higher mean wind speed, introducing the orography correction results in a more accurate forecast than that of any other model. This is because both the roughness change correction and the orography corrections are estimated with mean neutral profiles, whose mean wind speed at the forecasting height is also close to 16 m s^{−1}. For period 3, the one with the highest wind speed, the increasing underprediction of AHRO and AHR with wind speed results in AH predicting better than the other models.
For stable cases (periods 1, 6, 9 and 10) the performance of the advectionbased models is quite similar to the performance seen in neutral cases. During periods 6 and 10, the AHR model produces smaller errors than the statistical models. Figure 13 shows the comparison of the observed and forecasted wind speeds for all models during period 6. The figure shows that there is more scatter for the persistence and ARIMA models than for the advection models. For wind speeds above 8 m s^{−1} a high scatter between the advection models and the observations is also found.
For periods 1 and 9, all advection models show larger errors than the conventional models. This is because these periods are characterised by higher wind speeds than periods 6 and 10. The effect of the mean wind speed in the forecasting error of stable cases is shown in Fig. 12b. Above 6 m s^{−1} the forecast error tends to increase with wind speed. For higher wind speeds, the forecast originates from further upstream positions and consequently higher heights. If we now look at the differences between the PPI observations at 2950 m and the estimation of wind speeds using Eqs. (1) and (3) from the dualsetup observations (see Fig. 10), we can see that the differences are more pronounced at further distances. Thus, it is difficult to accurately predict the magnitude of the wind speed during stable conditions and high wind speeds, due to the increasing height in the observations at further positions, the differences in the dualsetup and PPI observations, and the assumption of neutral stabilities during stable conditions. Although we include the shear in our advection models we are not considering the atmospheric stability.
However, in quantifying the errors for the best fitted advection model in both stable and neutral cases, the RMSE of the stable cases is in general smaller than those of the neutral periods because during stable conditions the inflow is less turbulent. For stable cases, disregarding the periods of high wind speed (1 and 9), the best calibrated advection models give improvements over persistence of 21–26 % and over ARIMA of 24–28 %.
When looking at the forecast of wind speeds during period 7 (see Fig. 14), we can see that the advection models are able to forecast the phase of the events, but the forecast does not contain as many fluctuations as the observed wind speed at the downstream position. To analyse if this is due to the model or to the nature of the observations, the dependency of the level of fluctuations on the horizontal reconstructed wind speeds with the distance of the measurements is investigated. In Fig. 15, the ensemble average of the standard deviation of U, computed for every hour and elevation angle during periods when all measurements are available, is displayed. The standard deviation observed by the lidar is higher closer to the coast. We attribute this to a combination of two sources: sitespecific conditions and measurement artifacts. In the first source we consider the higher roughness length close to the coast, compared to positions further offshore, and the topographic effects. In the second source we include the different height in the observations for the different ranges and the different arc length used for the reconstruction of horizontal wind speeds from the lidar. Since the arc length used for the measurement increases with the distance, the reconstruction of wind speeds acts as a lowpass filter for further distances. This filtering effect deteriorates the prediction of the magnitude of the events, and consequently influences the maximum absolute error.
This paper evaluated the use of wind lidar observations for a very shortterm forecast of nearcoastal winds, using wind speed advectionbased models. From our analysis on periods with neutral atmospheric conditions, the bestfitted advectionbased model with corrections showed an improvement over the benchmarks persistence and ARIMA of 21–38 and 4–28 %, respectively. Due to the different turbulent conditions experienced at every range position, the forecasting technique was not able to predict the turbulence of the fluctuations. We attribute these differences partly to the presence of the coast increasing the turbulence level as the flow approaches and the lowpass filtering inherent in the wind speed reconstruction from the lidar measurements. During stable periods, we could only produce an accurate prediction of the magnitude of the wind speed during low wind speeds. This is a reflection of the increasing difficulty of predicting winds as (i) the observation height increases at further positions, (ii) there are differences in the dualsetup and PPI observations and (iii) neutral stability is assumed during stable conditions, due to a lack of a precise estimate of the offshore stability.
In this paper the forecasting horizon is limited to 5 min due to the maximum range of the lidar measurements (6 km) and the high wind speeds experienced during the measurement campaign. A longrange lidar system with a maximum range of 10 km could forecast wind speeds of up to 17 m s^{−1}, thus generating forecasts with a horizon of 10 min. Since commercially available ultrarange lidars can now measure up to 30 km (Kameyama et al., 2012), the forecasting horizon for this application could be extended up to 30 min ahead.
The corrections applied in our advectionbased models to forecast the magnitude of the wind speed observations are necessary due to the tilted trajectories and local effects of the coastline and the cliff. However, the corrections are not perfect. The results are based on a limited number of dualsetup measurements and it is clear that we could not find a model with a zero mean bias error. The best performing advection model depends on the wind speed and stability. Despite all these limitations, we showed that lidars (i) provide rangeresolved information to derive sitespecific effects influencing the wind speed and (ii) are promising candidates for very shortterm wind power applications since they can forecast wind speeds with more accuracy than the benchmarks persistence and ARIMA. To use an advectionbased wind speed forecasting technique, one could better benefit from horizontal trajectories that do not require height corrections. Additionally, applying this technique in pure offshore areas improves the results since no corrections due to local effects are required. An operational lidarbased forecasting system on an offshore wind farm would need no corrections at all. Thus, it is reasonable to expect that the forecasting performance of such a system would be better than the best results we have achieved since the many corrections might not have benefited the forecasting accuracy.
Very shortterm wind power forecasts typically use statistical techniques that learn from the wind speed and power data at the location of interest and surroundings. Based on our results, a longrange lidar system is likely to decrease the uncertainty in the prediction of offshore wind power, especially during ramp events, i.e. large variation in wind speed within a short period of time, for which statistical methods do not perform well.
Our analysis is a first input component to a decisionmaking model that may include spot market prices, scheduled supply and demand, and balancing costs. Thus, here it is not intended to quantify the economic impact of using a lidarbased wind speed forecast. However, as the balancing costs are proportional to the rootmeansquare error, it can be assumed that they will decrease. In particular, as in most of the periods analysed the maximum absolute error is lower than that of the benchmarks, using a lidarbased wind speed forecast might have a positive impact on integrating offshore wind power into the grid.
Further research will focus on using longrange, remotely sensed wind speed observations to predict the power produced by a single wind turbine or a wind farm.
Lidar data can be requested from Alfredo Peña at DTU Wind Energy (aldi@dtu.dk).
LV conducted the research work and wrote the paper. AP and MC extensively contributed to the modelling of the coastal effects and the use of lidar observations for forecasting wind speed, respectively. LvB and MK supervised the research work and contributed to the structure of the paper. All coauthors participated in the outline and review of the paper.
The authors declare that they have no conflict of interest.
This project has received funding from the European Union's Horizon 2020
research and innovation programme under the Marie SkłodowskaCurie grant
agreement no. 642108. Funding from the ForskEL program to the project “RUNE”
no. 12263 and from the Ministry of Science and Culture of Lower Saxony to the
project “ventus efficiens” (ZN3024, MWK Hannover) is
acknowledged.
Edited by: Julie
Lundquist
Reviewed by: two anonymous referees
Alexiadis, M., Dokopoulos, P., Sahsamanoglou, H., and Manousaridis, I.: Shortterm forecasting of wind speed and related electrical power, Sol. Energy, 63, 61–68, https://doi.org/10.1016/S0038092X(98)000322, 1998. a
Barthelmie, R., Murray, F., and Pryor, S.: The economic benefit of shortterm forecasting for wind energy in the UK electricity market, Energ. Policy, 36, 1687–1696, https://doi.org/10.1016/j.enpol.2008.01.027, 2008. a
Beck, H. and Kühn, M.: Dynamic data filtering of longrange Doppler LiDAR wind speed measurements, Remote SensingBasel, 9, 561, https://doi.org/10.3390/rs9060561, 2017. a
Box, G. E. P. and Jenkins, G. M.: Time series analysis: forecasting and control, J. Time Ser. Anal., 31, p. 303, https://doi.org/10.1111/j.14679892.2009.00643.x, 1976. a
Cadenas, E. and Rivera, W.: Wind speed forecasting in the South Coast of Oaxaca, México, Renew. Energ., 32, 2116–2128, https://doi.org/10.1016/j.renene.2006.10.005, 2007. a
Cadenas, E. and Rivera, W.: Short term wind speed forecasting in La Venta, Oaxaca, México, using artificial neural networks, Renew. Energ., 34, 274–278, https://doi.org/10.1016/j.renene.2008.03.014, 2009. a
Charnock, H.: Wind stress on a water surface, Q. J. Roy. Meteor. Soc., 81, 639–640, https://doi.org/10.1002/qj.49708135027, 1955. a
Damousis, I. G., Alexiadis, M. C., Theocharis, J. B., and Dokopoulos, P. S.: A fuzzy model for wind speed prediction and power generation in wind parks using spatial correlation, IEEE T. Energy Conver., 19, 352–361, https://doi.org/10.1109/Tec.2003.821865, 2004. a
Danish Wind Industry Association (DWIA): Wind Energy Production as a Percentage of Total Electricity Consumtion 2005–2017, available at: http://www.windpower.org/en/knowledge/statistics/the_danish_market.html (last access: 22 January 2018), 2018. a
Dickey, D. A. and Fuller, W. A.: Distribution of the estimators for autoregressive time series with a unit root, J. Am. Stat. Assoc., 74, 427–431, https://doi.org/10.2307/2286348, 1979. a
EPEXSPOT: Intraday Lead Times, available at: https://www.epexspot.com/en/productinfo/intradaycontinuous/intraday_lead_time (last access: 25 September 2017), 2017. a
Floors, R., Lea, G., Pena Diaz, A., Karagali, I., and Ahsbahs, T.: Report on RUNE's Coastal Experiment and First InterComparisons Between Measurements Systems, Tech. rep., DTU Wind Energy E0115(EN), DTU Wind Energy, Roskilde, Denmark, 2016a. a
Floors, R., Peña, A., Lea, G., Vasiljević, N., Simon, E., and Courtney, M.: The RUNE Experiment – a database of remotesensing observations of nearshore winds, Remote SensingBasel, 8, 884, https://doi.org/10.3390/rs8110884, 2016b. a, b, c, d, e
Frehlich, R.: Scanning Doppler lidar for input into shortterm wind power forecasts, J. Atmos. Ocean. Tech., 30, 230–244, https://doi.org/10.1175/JTECHD1100117.1, 2013. a
Geostyrelsen: available at: http://download.kortforsyningen.dk/content/dhm2007overflade16mgrid (last access: 5 April 2017), 2016. a
Giebel, G., Brownsword, R., Kariniotakis, G., Denhard, M., and Draxl, C.: The StateoftheArt in ShortTerm Prediction of Wind Power – A Literature Overview, Technical Report, EU Project ANEMOS, 109 pp., available at: http://orbit.dtu.dk/getResource?recordId=274635&objectId=1&versionId=1 (last access: 28 May 2018), 2011. a
GonzalezAparicio, I. and Zucker, A.: Impact of wind power uncertainty forecasting on the market integration of wind energy in Spain, Appl. Energ., 159, 334–349, https://doi.org/10.1016/j.apenergy.2015.08.104, 2015. a
Hill, D. C., McMillan, D., Bell, K. R. W., and Infield, D.: Application of autoregressive models to U.K. wind speed data for power system impact studies, IEEE T. Sustain. Energ., 3, 134–141, https://doi.org/10.1109/TSTE.2011.2163324, 2012. a
Holttinen, H., Meibom, P., Orths, A., Hulle, F. V., Lange, B., O'Malley, M., Pierik, J., Ummels, B., Tande, J. O., Estanqueiro, A., Gomez, E., Söder, L., Strbac, G., Shakoor, A., Ricardo, J., Smith, C., Milligan, M., and Ela, E.: Design and Operation of Power Systems with Large Amounts of Wind Power, Final Summary Report, IEA Wind Task 25, available at: https://community.ieawind.org/task25/viewdocument/designandoperationofpowersyste2?CommunityKey=4aa822101b2e43c5b37b1cdf11020dc8 (last access: 28 May 2018), 2016. a
Høyer, J. L. and She, J.: Optimal interpolation of sea surface temperature for the North Sea and Baltic Sea, J. Marine Syst., 65, 176–189, https://doi.org/10.1016/j.jmarsys.2005.03.008, 2007. a
Ibanez, E. and Milligan, M.: Impact of transmission on resource adequacy in systems with wind and solar power, in: IEEE Power and Energy Society General Meeting, 1–5, https://doi.org/10.1109/PESGM.2012.6343955, 2012. a
Kameyama, S., Sakimura, T., Watanabe, Y., Ando, T., Asaka, K., Tanaka, H., Yanagisawa, T., Hirano, Y., and Inokuchi, H.: Wind sensing demonstration of more than 30 km measurable range with a 1.5 µm coherent Doppler lidar which has the laser amplifier using Er,Yb:glass planar waveguide, Proc. SPIE, 8526, https://doi.org/10.1117/12.977330, 2012. a
Kavasseri, R. G. and Seetharaman, K.: Dayahead wind speed forecasting using fARIMA models, Renew. Energ., 34, 1388–1393, https://doi.org/10.1016/j.renene.2008.09.006, 2009. a, b
Kraus, E. B.: AtmosphereOcean Interaction, Oxford University Press, New York, p. 275, 1972. a
Larson, K. A. and Westrick, K.: Shortterm wind forecasting using offsite observations, Wind Energy, 9, 55–62, https://doi.org/10.1002/we.179, 2006. a
Li, G. and Shi, J.: On comparing three artificial neural networks for wind speed forecasting, Appl. Energ., 87, 2313–2320, https://doi.org/10.1016/j.apenergy.2009.12.013, 2010. a
Mikkelsen, T., Angelou, N., Hansen, K., Sjöholm, M., Harris, M., Slinger, C., Hadley, P., Scullion, R., Ellis, G., and Vives, G.: A spinnerintegrated wind lidar for enhanced wind turbine control, Wind Energy, 16, 625–643, https://doi.org/10.1002/we.1564, 2013. a
Monfared, M., Rastegar, H., and Kojabadi, H. M.: A new strategy for wind speed forecasting using artificial intelligent methods, Renew. Energ., 34, 845–848, https://doi.org/10.1016/j.renene.2008.04.017, 2009. a
Peña, A., Floors, R., Sathe, A., Gryning, S. E., Wagner, R., Courtney, M. S., Larsén, X. G., Hahmann, A. N., and Hasager, C. B.: Ten years of boundarylayer and windpower meteorology at Høvsøre, Denmark, Bound.Lay. Meteorol., 158, 1–26, https://doi.org/10.1007/s1054601500798, 2016. a
Peña, A., Mann, J., and Dimitrov, N.: Turbulence characterization from a forwardlooking nacelle lidar, Wind Energ. Sci., 2, 133–152, https://doi.org/10.5194/wes21332017, 2017. a
Sanchez, R. and Rørbæk, K.: Metocean Buoy Deployment, Technical Report, DHI, 2016. a
Schlueter, R. A., Sigari, G., and Costi, A.: Wind array power prediction for improved operating economics and reliability, IEEE T. Power Syst., 1, 137–142, https://doi.org/10.1109/TPWRS.1986.4334859, 1986. a
Simon, E. and Courtney, M.: A Comparison of SectorScan and Dual Doppler Wind Measurements at Høvsøre Test Station – One Lidar or Two?, Tech. rep., DTU Wind Energy E0112(EN), DTU Wind Energy, Roskilde, Denmark, 2016. a, b
Smith, S. D.: Wind stress and heat flux over the ocean in gale force winds, J. Phys. Oceanogr., 10, 709–726, https://doi.org/10.1175/15200485(1980)010<0709:WSAHFO>2.0.CO;2, 1980. a
Stawiarski, C., Träumner, K., Knigge, C., and Calhoun, R.: Scopes and challenges of dualDoppler lidar wind measurements – An Error Analysis, J. Atmos. Ocean. Tech., 30, 2044–2062, https://doi.org/10.1175/JTECHD1200244.1, 2013. a
Taylor, G. I.: The spectrum of turbulence, P. Roy Soc. AMath. Phy., 164, 476–490, https://doi.org/10.1098/rspa.1938.0032, 1938. a
Torres, J. L., Garcia, A., De Blas, M., and De Francisco, A.: Forecast of hourly average wind speed with ARMA models in Navarre (Spain), Sol. Energy, 79, 65–77, https://doi.org/10.1016/j.solener.2004.09.013, 2005. a, b
Troen, I. and Lundtang Petersen, E.: European Wind Atlas, Risø National Laboratory, Roskilde, p. 656, 1989. a
Wang, Q., MartinezAnido, C. B., Wu, H., Florita, A. R., and Hodge, B. M.: Quantifying the economic and grid reliability impacts of improved wind power forecasting, IEEE T. Sustain. Energ., 7, 1525–1537, https://doi.org/10.1109/TSTE.2016.2560628, 2016. a
Wharton, S., Newman, J., Qualley, G., and Miller, W.: Measuring turbine inflow with verticallyprofiling lidar in complex terrain, J. Wind Eng. Ind. Aerod., 142, 217–231, https://doi.org/10.1016/j.jweia.2015.03.023, 2015. a