Measurement and Prediction of Wind Fields at an Offshore Site by Scanning Doppler LiDAR and WRF

LiDAR-based wind speed measurements have seen a significant increase in interest in wind energy. However, reconstruction of wind speed vector from a LiDAR-measured radial wind speed is still a challenge. Furthermore, for extensive application of LiDAR technology, it can be used as a means to validate simulation and analytical models. To that end, this study employed scanning Doppler LiDAR for assessment of wind fields at an offshore site and compared Weather Research and Forecasting (WRF)-based mesoscale simulations and several wake models with the measurements. Firstly, the effect of carrier-to-noise-ratio (CNR) and data availability on the quality of scanning LiDAR measurements was evaluated. Analysis of vertical profiles show that the average wind speed is higher for wind blowing from the sea than that blowing from the land. Furthermore, profiles obtained from the WRF simulation also show a similar tendency in the LiDAR measurements in general, though it overestimates the wind speeds at higher altitudes. A method for reconstruction of wind fields from plan-position indicator (PPI) and range height indicator (RHI) scans of LiDAR-measured line of sight velocities was then proposed and first used to investigate the effect of coastal terrain. An internal boundary layer with strong shear could be observed to develop from the coastline. Finally, the flow field around wind turbine was measured using PPI scan and used to validate wake models.


Introduction
Measurement and collection of accurate wind data is important for a range of wind energy applications, including wind resource evaluation for prospective wind farm sites, optimization of farm layouts, and control of turbines. The general practice in the wind energy industry is to use instruments mounted at a single height or at multiple heights on meteorological towers in order to collect wind data. However, due to structural and cost constraints, such towers are usually about 50 m to 100 m in height; thus, they cannot measure the wind field across the rotor of recent utility-scale turbines with typical hub heights of 100 to 130 m and rotor diameters of 150 to 200 m. As a result, remote sensing techniques and in particular LiDAR technologies are getting increasingly popular in wind energy research due to their ability to measure wind speeds over large regions and higher altitudes. Another crucial advantage of LiDARs over conventional meteorological masts is that the former provides flexibility regarding transportation and installation. This can be particularly important for offshore sites, where installation of a meteorological tower with either a fixed or floating platform can be prohibitively expensive.
Doppler LiDAR measures wind speed by emitting a laser beam that gets Doppler-shifted after being backscattered from aerosol particles, which are assumed to be moving at the same speed as the wind. This Doppler shift is used to compute the wind speed along the beam direction, which is known as the line of sight (LOS) velocity or radial velocity (u r ). Various methods have been proposed

Measurement and Prediction Methods
This section first describes the test facility and the measurement devices. This is followed by discussions about three scan modes and the methods for the retrieval of velocity vectors from LOS velocity data. Finally, the details about the WRF (Weather Research and Forecasting) for the wind field prediction of the site is presented.

Test Site for Wind Field Measurement and Prediction
The measurements were collected from the Choshi offshore wind energy test facility located about 3.5 Km offshore of Choshi city of Chiba prefecture, Japan. As shown in Figure 1, the facility consists of a 2.4 MW wind turbine (MWT92/2.4) with a rotor diameter D = 92 m and a hub height h = 80 m, and a meteorological tower located 285 m east from the turbine. The meteorological tower has a height of 95 m, with 3D sonic anemometers (SAT-600A of Sonic Corporation) located at 40 m, 60 m, and 80 m. In addition to this, cup anemometers and wind vanes are installed between heights 20 m and 90 m at an interval of 10 m. The tower has a platform at a height of 15 m above the mean sea level. Two LiDARs, WindCube V1 (Leosphere, Orsay, France) and a scanning Doppler LiDAR (WindCube100S), are mounted on this platform. V1 is a profiling LiDAR; therefore, it can only measure vertical profiles at the installation point. WindCube 100S is capable of performing 3D volumetric scans [17]. Except for the validation purpose, the current study primarily used WindCube100S. Note that measurements from the V1 LiDAR were already validated against measurements from a sonic anemometer [18].

Wind Field Measurement by Scanning Doppler LiDAR
In order to retrieve vertical profiles of wind speeds and directions, line-of-site (LOS) wind speeds ( ) were measured in the Doppler Beam Swinging (DBS) scanning configuration, shown schematically in Figure 2. In this configuration, the LiDAR beam is swung from north to east to south to west and to vertical directions. The wind speed can then be calculated at each measurement height using: where , , are wind speed components in west to east, north to south, and vertical directions, respectively. , , , , are radial wind speeds along the LOS measured in east, west,

Wind Field Measurement by Scanning Doppler LiDAR
In order to retrieve vertical profiles of wind speeds and directions, line-of-site (LOS) wind speeds (u r ) were measured in the Doppler Beam Swinging (DBS) scanning configuration, shown schematically in Figure 2. In this configuration, the LiDAR beam is swung from north to east to south to west and to vertical directions. The wind speed can then be calculated at each measurement height using: where u, v, w are wind speed components in west to east, north to south, and vertical directions, respectively. u rE , u rW , u rN , u rS , u rV are radial wind speeds along the LOS measured in east, west, north, south, and vertical directions, respectively. Finally, φ is the elevation angle and, except for the Atmosphere 2020, 11,442 4 of 20 vertical scan, it is set to 62 • . Note that some LiDARs in DBS mode do not perform vertical scans and instead compute the vertical velocity from the following relation: w = u rE + u rN + u rW + u rS 4 sin φ . (2) conditions (aerosol concentration, presence of cloud, or mist layers) strongly influence the availability of the LiDAR data. Therefore, it is necessary to evaluate how availability influences the quality of measured wind speed. To this end, a 10-min availability of wind speed data is defined as [19]: where is availability, is the maximum number of data that the LiDAR can measure in 10 min, and is the actual number of data measured by the LiDAR in a 10-min time period. For the scanning LiDAR, is 140. When is small, the number of samples for computing the 10-minaveraged wind speed is also small, and this may result in a higher error in the mean wind speed. In addition to its application as a profiling LiDAR with a very high measurement range, scanning LiDAR can also perform conical and vertical slice scans [17]. Figure 3 shows two commonly employed scan modes, i.e., range height indicator (RHI) and plan-position indicator (PPI), of scanning Doppler LiDAR. In the range height indicator (RHI) scan, the azimuth angle is fixed, and the elevation angle is varied to sweep the beam through the vertical slice. In the plan-position indicator (PPI) scan, the LiDAR beam sweeps over a range of azimuth angles, while maintaining the elevation angle at a constant value. Usually, in an RHI scan, high resolution samples are acquired in the direction of the elevation angle, whereas, in a PPI scan, high resolution samples are acquired in the azimuthal direction [20]. Multiples of such scans by varying azimuth angle for RHI mode (or elevation angle for PPI mode) are necessary to construct a three-dimensional flow field. However, care should be taken in that multiple RHI or PPI scans can result in a significantly low time resolution, thus making it difficult to construct reliable turbulence fields from scanning LiDAR measurements. Using the DBS mode, the vertical velocity profiles were collected for a period of four months between March and July of 2016. It should be noted that the scanning speed of the scanning LiDAR is 20 s per cycle, whereas the V1 LiDAR completes one scan cycle every five seconds. Consequently, the number of data collected by the scanning LiDAR is 1/4th of the data collected by the V1 LiDAR. This is because the range of the scanning LiDAR is 3000 m, while the maximum range of the V1 LiDAR is 186 m, thus requiring a longer processing time for the former. Furthermore, atmospheric conditions (aerosol concentration, presence of cloud, or mist layers) strongly influence the availability of the LiDAR data. Therefore, it is necessary to evaluate how availability influences the quality of measured wind speed. To this end, a 10-min availability of wind speed data is defined as [19]: where η is availability, N max is the maximum number of data that the LiDAR can measure in 10 min, and N mes is the actual number of data measured by the LiDAR in a 10-min time period. For the scanning LiDAR, N max is 140. When η is small, the number of samples for computing the 10-min-averaged wind speed is also small, and this may result in a higher error in the mean wind speed.
In addition to its application as a profiling LiDAR with a very high measurement range, scanning LiDAR can also perform conical and vertical slice scans [17]. Figure 3 shows two commonly employed scan modes, i.e., range height indicator (RHI) and plan-position indicator (PPI), of scanning Doppler LiDAR. In the range height indicator (RHI) scan, the azimuth angle is fixed, and the elevation angle is varied to sweep the beam through the vertical slice. In the plan-position indicator (PPI) scan, the LiDAR beam sweeps over a range of azimuth angles, while maintaining the elevation angle at a constant value. Usually, in an RHI scan, high resolution samples are acquired in the direction of the elevation angle, whereas, in a PPI scan, high resolution samples are acquired in the azimuthal direction [20]. Multiples of such scans by varying azimuth angle for RHI mode (or elevation angle for PPI mode) are necessary to construct a three-dimensional flow field. However, care should be taken in that multiple RHI or PPI scans can result in a significantly low time resolution, thus making it difficult to construct reliable turbulence fields from scanning LiDAR measurements.  As stated earlier, in order to retrieve the velocity vector (V) from the LiDAR-measured line of site velocity ( ), additional processing of the data is required. The method for the estimation of V from the PPI and RHI scan data is explained next. Assuming that the vertical component of wind speed w is very small, it was neglected in this study. and components can then be expressed as: with being the wind direction. The LOS velocity r measured by the scanning LiDAR is the projection of V in the laser beam direction; therefore, V can be expressed as: Here, is the azimuthal angle and is the elevation angle. The wind direction at every point in a scan is assumed to be the same as that measured by the wind vane on the meteorological tower or the nacelle. With measured by the LiDAR, , , and known, computation of V is then straightforward from Equation (5). Once V is known, it can be substituted into Equation (4) to compute and v. An alternative approach to retrieve the velocity vector (V) from can be to sweep the laser beam in a conical scan for a complete 360°. The azimuth angle for which the magnitude of is maximum can be considered as the wind direction. Figure 4 shows an example of as a function of collected using a PPI scan. Peaks can be observed around azimuth angles 30° and 210°. Of the two peaks, the azimuth angle for which the magnitude of is higher is used to estimate the wind direction. Note that, for the current LiDAR, is positive when the wind blows away from the LiDAR and is negative otherwise. In order to avoid interference due to the meteorological tower and other installations on the platform, the region 90°≤ ≤ 180° is not scanned in the current study. With the wind direction known, V can be computed from Equation (5). This essentially is a form of VAD approach (see Section 1). Wind speeds reconstructed using this approach are also presented in the next section. As stated earlier, in order to retrieve the velocity vector (V) from the LiDAR-measured line of site velocity (u r ), additional processing of the data is required. The method for the estimation of V from the PPI and RHI scan data is explained next. Assuming that the vertical component of wind speed w is very small, it was neglected in this study. u and v components can then be expressed as: with α being the wind direction. The LOS velocity u r measured by the scanning LiDAR is the projection of V in the laser beam direction; therefore, V can be expressed as: Here, θ is the azimuthal angle and φ is the elevation angle. The wind direction at every point in a scan is assumed to be the same as that measured by the wind vane on the meteorological tower or the nacelle. With u r measured by the LiDAR, φ, α, and θ known, computation of V is then straightforward from Equation (5). Once V is known, it can be substituted into Equation (4) to compute u and v.
An alternative approach to retrieve the velocity vector (V) from u r can be to sweep the laser beam in a conical scan for a complete 360 • . The azimuth angle for which the magnitude of u r is maximum can be considered as the wind direction. Figure 4 shows an example of u r as a function of θ collected using a PPI scan. Peaks can be observed around azimuth angles 30 • and 210 • . Of the two peaks, the azimuth angle for which the magnitude of u r is higher is used to estimate the wind direction. Note that, for the current LiDAR, u r is positive when the wind blows away from the LiDAR and is negative otherwise. In order to avoid interference due to the meteorological tower and other installations on the platform, the region 90 • ≤ θ ≤ 180 • is not scanned in the current study. With the wind direction known, V can be computed from Equation (5). This essentially is a form of VAD approach (see Section 1). Wind speeds reconstructed using this approach are also presented in the next section.

Wind Field Prediction by WRF
Considering that mesoscale simulations are often used to conduct wind resource assessment of wind farm sites, accuracy of mesoscale model is crucial for accurate power prediction. In the current paper, LiDAR measurements were used to validate results of mesoscale simulations. For the mesoscale simulation, the Weather Research and Forecasting (WRF Ver. 3.4.1) model developed by the National Center of Atmospheric research (NCAR) was employed [21]. The computation domain has five levels of horizontal nesting, i.e., 18 km-6 km-2 km-667 m-222m, with the higher resolution in the region where LiDAR measurements were conducted (cf. Figure 1). Other than wind resource assessment for wind energy [18], WRF-based simulations are extensively used in wind engineering applications [22,23]. The results from the three grid resolutions of 2 km, 667 m, and 222 m were compared with the measurements. As the initial and boundary conditions of the simulations, National Center for Environmental Prediction (NCEP) final operational global analysis data (NCEP-FNL 1° × 1°, 6-hourly) data was used [24]. However, for the sea surface temperature, Operational Sea Surface Temperature and Ice Analysis (OSTIA) data was used [25]. OSTIA provides near-real time global sea surface temperature at the grid resolution of 1/20° (~ 6 km).

Results and Discussion
The LiDAR measurements are presented and discussed in this section. Section 3.1 presents vertical wind profiles measured using DBS mode and compares those by mesoscale simulations. Section 3.2 validates the proposed method for velocity retrieval from volumetric scan. In Section 3.3, the effect of coastal terrain on the near-shore marine ABL is characterized. The flow field around the wind turbine is measured and compared with wake models in Section 3.4.

Measurement and Prediction of Vertical Wind Profiles
Before presenting the comparison and analysis of vertical wind profiles, the quality of LiDAR measurements is assessed. To that end, the effect of carrier-to-noise-ratio (CNR) on the quality of measurement, particularly for different availability limits ( ), is discussed. Here, is the lower limit of the availability (η) set for the scanning LiDAR data (cf. Equation (3)). Figure 5 shows the distribution of CNR for the measurements at the height of 140 m. As is clear from the figure, the minimum threshold value of CNR is −25 dB. However, in order to investigate the effect of CNR, the minimum threshold values are increased to −23 dB and −20 dB. Figure 6 presents comparisons of 10min-averaged horizontal velocities measured by scanning LiDAR and V1 LiDAR, for these three CNR limits and for 80%. A comparison of statistical parameters for an 20% and 80% and for all three CNR limits is summarized in Table 1. It can be appreciated from the comparison that the

Wind Field Prediction by WRF
Considering that mesoscale simulations are often used to conduct wind resource assessment of wind farm sites, accuracy of mesoscale model is crucial for accurate power prediction. In the current paper, LiDAR measurements were used to validate results of mesoscale simulations. For the mesoscale simulation, the Weather Research and Forecasting (WRF Ver. 3.4.1) model developed by the National Center of Atmospheric research (NCAR) was employed [21]. The computation domain has five levels of horizontal nesting, i.e., 18 km-6 km-2 km-667 m-222 m, with the higher resolution in the region where LiDAR measurements were conducted (cf. Figure 1). Other than wind resource assessment for wind energy [18], WRF-based simulations are extensively used in wind engineering applications [22,23]. The results from the three grid resolutions of 2 km, 667 m, and 222 m were compared with the measurements. As the initial and boundary conditions of the simulations, National Center for Environmental Prediction (NCEP) final operational global analysis data (NCEP-FNL 1 • × 1 • , 6-hourly) data was used [24]. However, for the sea surface temperature, Operational Sea Surface Temperature and Ice Analysis (OSTIA) data was used [25]. OSTIA provides near-real time global sea surface temperature at the grid resolution of 1/20 • (~6 km).

Results and Discussion
The LiDAR measurements are presented and discussed in this section. Section 3.1 presents vertical wind profiles measured using DBS mode and compares those by mesoscale simulations. Section 3.2 validates the proposed method for velocity retrieval from volumetric scan. In Section 3.3, the effect of coastal terrain on the near-shore marine ABL is characterized. The flow field around the wind turbine is measured and compared with wake models in Section 3.4.

Measurement and Prediction of Vertical Wind Profiles
Before presenting the comparison and analysis of vertical wind profiles, the quality of LiDAR measurements is assessed. To that end, the effect of carrier-to-noise-ratio (CNR) on the quality of measurement, particularly for different availability limits (η lim ), is discussed. Here, η lim is the lower limit of the availability (η) set for the scanning LiDAR data (cf. Equation (3)). Figure 5 shows the distribution of CNR for the measurements at the height of 140 m. As is clear from the figure, the minimum threshold value of CNR is −25 dB. However, in order to investigate the effect of CNR, the minimum threshold values are increased to −23 dB and −20 dB. Figure 6 presents comparisons of 10-min-averaged horizontal velocities measured by scanning LiDAR and V1 LiDAR, for these three CNR limits and for η lim = 80%. A comparison of statistical parameters for an η lim = 20% and 80% and for all three CNR limits is summarized in Table 1. It can be appreciated from the comparison that the quality of measurements degrades when the minimum threshold value of the CNR is increased.
In particular, the RMSE increases significantly when the minimum threshold of the CNR is increased for the same η lim . This is due to the distribution of CNR, as shown in Figure 4. Setting a higher value for the minimum threshold of the CNR means discarding data below that CNR.
Atmosphere 2020, 11, x 7 of 20 quality of measurements degrades when the minimum threshold value of the CNR is increased. In particular, the RMSE increases significantly when the minimum threshold of the CNR is increased for the same lim . This is due to the distribution of CNR, as shown in Figure 4. Setting a higher value for the minimum threshold of the CNR means discarding data below that CNR.  Atmosphere 2020, 11, x 7 of 20 quality of measurements degrades when the minimum threshold value of the CNR is increased. In particular, the RMSE increases significantly when the minimum threshold of the CNR is increased for the same . This is due to the distribution of CNR, as shown in Figure 4. Setting a higher value for the minimum threshold of the CNR means discarding data below that CNR.   Next, Table 2 summarizes a comparison of 10-min-averaged horizontal velocities measured by scanning LiDAR and V1 LiDAR for four different availability limits (η lim ) for the minimum threshold value of CNR of −25 dB. Note that, for η lim = 0%, even if one data is measured in a 10-min time slot, then that time slot is considered valid. It can be appreciated from the table that, for lower η lim , the variation is larger and the coefficient of determination (R 2 ), as a result, is also smaller. However, the difference between R 2 for η lim = 20% and that for η lim = 80% does not change too much. A further comparison between the 100S and V1 LiDARs is presented quantitatively using root mean square error (RMSE) in Figure 7. The RMSE increases with decreasing η lim .  Next, Table 2 summarizes a comparison of 10-min-averaged horizontal velocities measured by scanning LiDAR and V1 LiDAR for four different availability limits ( ) for the minimum threshold value of CNR of −25 dB. Note that, for 0%, even if one data is measured in a 10-min time slot, then that time slot is considered valid. It can be appreciated from the table that, for lower , the variation is larger and the coefficient of determination ( ), as a result, is also smaller. However, the difference between for 20% and that for 80% does not change too much. A further comparison between the 100S and V1 LiDARs is presented quantitatively using root mean square error (RMSE) in Figure 7. The RMSE increases with decreasing .  The availability of scanning Doppler LiDAR data is further analyzed as a function of height. Figure 8a shows data availability in the form of number of data measured every 10 min, i.e., as the function of height. It is obvious from the figure that the is higher at the lower height and reduces with height. Figure 8b shows the vertical profiles of the availability ( ) of 10-min-averaged wind speeds. is given by: The availability of scanning Doppler LiDAR data is further analyzed as a function of height. Figure 8a shows data availability in the form of number of data measured every 10 min, i.e., N mes as the function of height. It is obvious from the figure that the N mes is higher at the lower height and reduces with height. Figure 8b shows the vertical profiles of the availability (η 10 ) of 10-min-averaged wind speeds. η 10 is given by: possible number of 10-min-averaged data for the given time period. As shown in the figure, if lim is set high, 10 becomes smaller, and vice versa. For example, when lim is set to 80%, 10 is around 40% at a height of 500 m, whereas, when lim is set to 20%, 10 is around 70% at the same height. Therefore, although setting a higher lim improves the quality of the average wind speed data (see Table 2), the number of data decreases significantly. This effect is more pronounced for higher altitudes.  Figure 9 presents one-month-averaged boundary layer profiles for the wind blowing from the land and that from the sea obtained from the LiDAR measurements, as well as from the mesoscale simulations. The definition of the wind blowing from the land and from the sea is based on the Here, N 10,η lim is the number of 10-min-averaged data for η ≥ η lim , and N 10, max is the maximum possible number of 10-min-averaged data for the given time period. As shown in the figure, if η lim is set high, η 10 becomes smaller, and vice versa. For example, when η lim is set to 80%, η 10 is around 40% at a height of 500 m, whereas, when η lim is set to 20%, η 10 is around 70% at the same height. Therefore, although setting a higher η lim improves the quality of the average wind speed data (see Table 2), the number of data decreases significantly. This effect is more pronounced for higher altitudes. Figure 9 presents one-month-averaged boundary layer profiles for the wind blowing from the land and that from the sea obtained from the LiDAR measurements, as well as from the mesoscale simulations. The definition of the wind blowing from the land and from the sea is based on the coastline feature of the test site (cf. Figure 1). The simulation results show a good agreement with the measurements for the wind blowing from the sea, with the maximum error of 2.3%. Regarding the wind blowing from the land, the velocity profile obtained from the simulation reproduces the measured profile up to a height of 300 m. However, for higher altitudes, the simulation overestimates the measurement by up to 10%. The largest difference between the simulation results (due to grid resolution) was observed at 490 m, and it was about 0.05%. However, the differences between simulations with the three grid resolutions are negligible.
In Figure 9, the wind speed is higher when the wind is blowing from the sea than when it is blowing from the land. For example, around the hub height level (≈ 80 m), the wind speed when the wind is blowing from the land is 6.8 m/s, whereas that of the wind blowing from the sea is 9.2 m/s, with the latter being 35% higher. This trend is also observed in the diurnal cycle of vertical profiles in Figure 10. The figure shows hourly averaged vertical profiles, which are further averaged over the period of one month. With very few exceptions, profiles for the wind blowing from sea have higher wind speeds. Thus, the effect of wind direction on power production can be significant for near-shore wind farms and must be taken into account while deciding the layout of such farms. Note that this work does not consider boundary layer stratification. Dividing the measured data into stable, neutral, and convective boundary layers and then analyzing the wind blowing from sea and land may present a further detailed picture of the ABL at the site, as shown by Kikuchi et al. [18].    coastline feature of the test site (cf. Figure 1). The simulation results show a good agreement with the measurements for the wind blowing from the sea, with the maximum error of 2.3%. Regarding the wind blowing from the land, the velocity profile obtained from the simulation reproduces the measured profile up to a height of 300 m. However, for higher altitudes, the simulation overestimates the measurement by up to 10%. The largest difference between the simulation results (due to grid resolution) was observed at 490 m, and it was about 0.05%. However, the differences between simulations with the three grid resolutions are negligible. In Figure 9, the wind speed is higher when the wind is blowing from the sea than when it is blowing from the land. For example, around the hub height level ( 80 m), the wind speed when the wind is blowing from the land is 6.8 m/s, whereas that of the wind blowing from the sea is 9.2 m/s, with the latter being 35% higher. This trend is also observed in the diurnal cycle of vertical profiles in Figure 10. The figure shows hourly averaged vertical profiles, which are further averaged over the period of one month. With very few exceptions, profiles for the wind blowing from sea have higher wind speeds. Thus, the effect of wind direction on power production can be significant for near-shore wind farms and must be taken into account while deciding the layout of such farms. Note that this work does not consider boundary layer stratification. Dividing the measured data into stable, neutral, and convective boundary layers and then analyzing the wind blowing from sea and land may present a further detailed picture of the ABL at the site, as shown by Kikuchi et al. [18].  Figure 11 shows the comparison of simulated and measured hourly averaged boundary layer profiles, which are further averaged over the period of one month. The comparisons show that during the morning hours the difference between the simulations and the measurements is significant, with  Figure 11 shows the comparison of simulated and measured hourly averaged boundary layer profiles, which are further averaged over the period of one month. The comparisons show that during the morning hours the difference between the simulations and the measurements is significant, with simulations underestimating the wind speed in the lower region of the boundary layer (below 300 m) and overestimating it in the higher region (above 300 m). The largest error of almost 20% can be observed at the height of 140 m. On the other hand, during the day time, the agreement of the simulations with the measurements is significantly high. Therefore, although the long-term average data from mesoscale simulations are reliable, short-term statistics-like hourly averaged data-should be used with caution.
Atmosphere 2020, 11, x 11 of 20 simulations underestimating the wind speed in the lower region of the boundary layer (below 300 m) and overestimating it in the higher region (above 300 m). The largest error of almost 20% can be observed at the height of 140 m. On the other hand, during the day time, the agreement of the simulations with the measurements is significantly high. Therefore, although the long-term average data from mesoscale simulations are reliable, short-term statistics-like hourly averaged datashould be used with caution.

Validation of velocity vector retrieved from PPI and RHI scan data
This section presents validation of the two approaches -proposed in Section 2.2-for retrieving velocity vector from LOS velocity fields obtained in the PPI and RHI scans. The first approach, which requires wind direction information (cf. Equations (4) and (5)), is validated against measurements from the anemometer mounted on the nacelle of the neighboring wind turbine (see Figure 1). A PPI scan is employed to collect scanning LiDAR measurements, as shown in Figure 12. The LiDAR beam is swept at ∅ 13.2° and for an azimuth angle range of 40°, such that the beam passes through the nacelle anemometer, and the wind turbine is in the middle of the azimuth range. In order to obviate the effect of the turbine wake on the nacelle measurement, the comparison was made when the turbine was stopped.

Validation of Velocity Vector Retrieved from PPI and RHI Scan Data
This section presents validation of the two approaches -proposed in Section 2.2-for retrieving velocity vector from LOS velocity fields obtained in the PPI and RHI scans. The first approach, which requires wind direction information (cf. Equations (4) and (5)), is validated against measurements from the anemometer mounted on the nacelle of the neighboring wind turbine (see Figure 1). A PPI scan is employed to collect scanning LiDAR measurements, as shown in Figure 12. The LiDAR beam is swept at ∅ = 13.2 • and for an azimuth angle range of 40 • , such that the beam passes through the nacelle anemometer, and the wind turbine is in the middle of the azimuth range. In order to obviate the effect of the turbine wake on the nacelle measurement, the comparison was made when the turbine was stopped.
Atmosphere 2020, 11, x 11 of 20 simulations underestimating the wind speed in the lower region of the boundary layer (below 300 m) and overestimating it in the higher region (above 300 m). The largest error of almost 20% can be observed at the height of 140 m. On the other hand, during the day time, the agreement of the simulations with the measurements is significantly high. Therefore, although the long-term average data from mesoscale simulations are reliable, short-term statistics-like hourly averaged datashould be used with caution.

Validation of velocity vector retrieved from PPI and RHI scan data
This section presents validation of the two approaches -proposed in Section 2.2-for retrieving velocity vector from LOS velocity fields obtained in the PPI and RHI scans. The first approach, which requires wind direction information (cf. Equations (4) and (5)), is validated against measurements from the anemometer mounted on the nacelle of the neighboring wind turbine (see Figure 1). A PPI scan is employed to collect scanning LiDAR measurements, as shown in Figure 12. The LiDAR beam is swept at ∅ 13.2° and for an azimuth angle range of 40°, such that the beam passes through the nacelle anemometer, and the wind turbine is in the middle of the azimuth range. In order to obviate the effect of the turbine wake on the nacelle measurement, the comparison was made when the turbine was stopped.  Since the nacelle anemometer-measured wind speed is influenced by the rotor or the shape of the nacelle, the measurement of the nacelle anemometer is first compared against the wind speed measured by cup anemometers installed at a height of 80 m on the meteorological tower. The comparison showed that on average, the wind speed measured by the nacelle anemometer is 5.7% higher than that obtained from the meteorological tower measurement [19]. Therefore, the measurement of the nacelle anemometer is corrected using the velocity ratio 1.057. Figure 13 shows a comparison of 10-min-averaged wind speeds retrieved from the scanning LiDAR measurement against the measurements of the nacelle-mounted anemometer. The time series and scatter plots show a good agreement between the two measurements with R 2 = 0.979, RMSE = 0.562, and slope being 0.976.
Atmosphere 2020, 11, x 12 of 20 Since the nacelle anemometer-measured wind speed is influenced by the rotor or the shape of the nacelle, the measurement of the nacelle anemometer is first compared against the wind speed measured by cup anemometers installed at a height of 80 m on the meteorological tower. The comparison showed that on average, the wind speed measured by the nacelle anemometer is 5.7% higher than that obtained from the meteorological tower measurement [19]. Therefore, the measurement of the nacelle anemometer is corrected using the velocity ratio 1.057. Figure 13 shows a comparison of 10-min-averaged wind speeds retrieved from the scanning LiDAR measurement against the measurements of the nacelle-mounted anemometer. The time series and scatter plots show a good agreement between the two measurements with 0.979, RMSE 0.562, and slope being 0.976.   Figures 14 and 15 show comparisons of the wind speed and wind direction obtained by scanning Doppler LiDAR (retrieved using the VAD technique) against the wind speed and wind direction measured by the met mast tower. The comparison height is 90 m. Both the time series and scatter plot show a good agreement between the two measurements. The coefficient of determination is 0.99 for both wind speed and wind direction, and the RMSE of wind speed is 0.31 m/s, whereas that of wind direction is 5.6°. However, it should be mentioned that the complete conical scan required by this approach can significantly reduce the data sample. Furthermore, due to uncertainty in the LiDAR measurement, it may be necessary to perform multiple scans for the same elevation angle in order to obtain accurate wind directions. This makes the approach less attractive for constructing the  Figures 14 and 15 show comparisons of the wind speed and wind direction obtained by scanning Doppler LiDAR (retrieved using the VAD technique) against the wind speed and wind direction measured by the met mast tower. The comparison height is 90 m. Both the time series and scatter plot show a good agreement between the two measurements. The coefficient of determination is R 2 = 0.99 for both wind speed and wind direction, and the RMSE of wind speed is 0.31 m/s, whereas that of wind direction is 5.6 • . However, it should be mentioned that the complete conical scan required by this approach can significantly reduce the data sample. Furthermore, due to uncertainty in the LiDAR measurement, it may be necessary to perform multiple scans for the same elevation angle in order to obtain accurate wind directions. This makes the approach less attractive for constructing the velocity field. The current study employs the former approach validated in this section for velocity field reconstruction.
Atmosphere 2020, 11, x 12 of 20 Since the nacelle anemometer-measured wind speed is influenced by the rotor or the shape of the nacelle, the measurement of the nacelle anemometer is first compared against the wind speed measured by cup anemometers installed at a height of 80 m on the meteorological tower. The comparison showed that on average, the wind speed measured by the nacelle anemometer is 5.7% higher than that obtained from the meteorological tower measurement [19]. Therefore, the measurement of the nacelle anemometer is corrected using the velocity ratio 1.057. Figure 13 shows a comparison of 10-min-averaged wind speeds retrieved from the scanning LiDAR measurement against the measurements of the nacelle-mounted anemometer. The time series and scatter plots show a good agreement between the two measurements with 2 = 0.979, RMSE = 0.562, and slope being 0.976.

Measurement and Prediction of Wind Field in Near-Shore Boundary Layer
In order to evaluate the effect of coastal terrain on the near-shore marine ABL, the wind field between the meteorological tower and the coastline was measured. RHI scans were conducted towards the north and for an elevation angle range of 0° to 60°. The wind direction information required by Equation (5) to construct wind speeds from LiDAR-measured LOS velocities ( ) was obtained from a wind vane mounted on the meteorological tower. Furthermore, for the spatial averaging, the measurement data was divided into vertical grids along the north south direction. The size of each block in the grid was set to 50 50 m . Average wind speeds were computed from all the measurement points in that block and for a fairly long time span of three hours.

Measurement and Prediction of Wind Field in Near-Shore Boundary Layer
In order to evaluate the effect of coastal terrain on the near-shore marine ABL, the wind field between the meteorological tower and the coastline was measured. RHI scans were conducted towards the north and for an elevation angle range of 0 • to 60 • . The wind direction information required by Equation (5) to construct wind speeds from LiDAR-measured LOS velocities (u r ) was obtained from a wind vane mounted on the meteorological tower. Furthermore, for the spatial averaging, the measurement data was divided into vertical grids along the north south direction. The size of each block in the grid was set to 50 × 50 m 2 . Average wind speeds were computed from all the measurement points in that block and for a fairly long time span of three hours. Figure 16 shows the mean velocity field in the vertical plane between the observation tower and the northern coast. The left side of the color map is 1000 m from the coast, and the LiDAR is 3500 m from the coast. Strong shear can be observed in the region close to the coast. The coastline has a vertical cliff of about 50 m above the sea level. This enforces the development of an internal boundary layer (IBL) that persists approximately up to 1500 m from the LiDAR position (almost 2 km from the shore).

Measurement and Prediction of Wind Field in Near-Shore Boundary Layer
In order to evaluate the effect of coastal terrain on the near-shore marine ABL, the wind field between the meteorological tower and the coastline was measured. RHI scans were conducted towards the north and for an elevation angle range of 0° to 60°. The wind direction information required by Equation (5) to construct wind speeds from LiDAR-measured LOS velocities ( ) was obtained from a wind vane mounted on the meteorological tower. Furthermore, for the spatial averaging, the measurement data was divided into vertical grids along the north south direction. The size of each block in the grid was set to 50 50 m . Average wind speeds were computed from all the measurement points in that block and for a fairly long time span of three hours.   simulations and the measurements is around 10% up to 500 m. Note that the height where the error line is broken in Figure 17a corresponds to the location where no measurement data was available. For Figure 17b,c, which are further from the coastline, LiDAR measurements are well reproduced by the profiles predicted by the WRF simulations up to a height of 250 m. However, as was the case for the wind blowing from the land (cf. Figure 9 c), for higher altitudes, the simulations overestimate the measurement, for the current measurement period. Until now, results from mesoscale simulations were mainly validated against tower-mounted anemometers or profiling LiDAR with a lower range. Therefore, it was only possible to confirm the accuracy of the simulation results at 200 m or below. By using scanning Doppler LiDAR, together with the proposed method for the reconstruction of the velocity field from LOS velocities, it is now possible to validate the mesoscale simulations against measurements for higher altitudes and different locations in the IBL.

Measurement and Prediction of Wind Turbine Wake
In this section, the flow field around the wind turbine is measured in order to characterize the wind turbine wake. Measuring the wind turbine wake using scanning Doppler LiDAR poses several challenges. In addition to a low sampling rate, PPI scans conducted using LiDAR can only provide measurements in conical planes. An alternative can be to conduct multiple RHI scans, but this method can only be applied when LiDAR is located in the wake of the turbine, and the wind direction does not vary too much. In this study, multiple PPI scans are performed for a sequence of elevation angles. Figure 18 shows a schematic of this scan approach. From the measured LOS velocities, only those at the hub height level are extracted (the grey plane in Figure 18). Wind speed is then constructed from the LOS velocity using Equation (5).  • from the north were averaged to produce these profiles. It can be appreciated that the simulations are able to produce more or less similar velocity magnitudes as the measurements. However, for the profiles at 900 m from the coast, i.e., closest to the coast, the difference between the measurement and simulations is significant. Difference between the simulations and the measurements is around 10% up to 500 m. Note that the height where the error line is broken in Figure 17a corresponds to the location where no measurement data was available. For Figure 17b,c, which are further from the coastline, LiDAR measurements are well reproduced by the profiles predicted by the WRF simulations up to a height of 250 m. However, as was the case for the wind blowing from the land (cf. Figure 9c), for higher altitudes, the simulations overestimate the measurement, for the current measurement period. Until now, results from mesoscale simulations were mainly validated against tower-mounted anemometers or profiling LiDAR with a lower range. Therefore, it was only possible to confirm the accuracy of the simulation results at 200 m or below. By using scanning Doppler LiDAR, together with the proposed method for the reconstruction of the velocity field from LOS velocities, it is now possible to validate the mesoscale simulations against measurements for higher altitudes and different locations in the IBL.

Measurement and Prediction of Wind Turbine Wake
In this section, the flow field around the wind turbine is measured in order to characterize the wind turbine wake. Measuring the wind turbine wake using scanning Doppler LiDAR poses several challenges. In addition to a low sampling rate, PPI scans conducted using LiDAR can only provide measurements in conical planes. An alternative can be to conduct multiple RHI scans, but this method can only be applied when LiDAR is located in the wake of the turbine, and the wind direction does not vary too much. In this study, multiple PPI scans are performed for a sequence of elevation angles. Figure 18 shows a schematic of this scan approach. From the measured LOS velocities, only those at the hub height level are extracted (the grey plane in Figure 18). Wind speed is then constructed from the LOS velocity using Equation (5).   Figure 19 visualizes the mean flow field around the wind turbine reconstructed from the PPI scans. A velocity deficit, which is a common characteristic of a turbine wake, is clearly visible in the figure. Furthermore, a dual peak is also observed in the near-wake region up to ⁄ = 2. It can be appreciated that the wake region persists approximately 7D downstream from the wind turbine. Note that, for ⁄ ≈ 2, the azimuth angles are approximately perpendicular to the incoming flow and, therefore, the measurement error is high in this region.
Several other studies have performed scanning Doppler LiDAR-based measurements of wind turbine wakes. However, reconstruction of velocity fields from LiDAR-measured LOS velocities poses a great challenge. Therefore, most of the studies simply tried to analyze and interpret LOS velocities, and did not perform any further processing of the measured data [26][27][28]. As is also mentioned in the Introduction, the assumption of constant wind direction assumption during the measurement period and the application of multiple LiDARs are also receiving a lot of interest in the wind energy community [29,30]. The method of wind turbine wake measurement with a single Scanning Doppler LiDAR for arbitrary directions proposed in this study is therefore different from those in the literature.
Finally, the measurement data is compared against four analytical wake models proposed by Jensen [31], Frandsen et al. [32], Bastankhah and Porte-Agel [33], and Qian and Ishihara [34]. As with Figure 18. Schematic of the multiple PPI scans for the measurement and reconstruction of flow field around wind turbine. Wind speed data is extracted in the grey horizontal plane at the hub height level. Figure 19 visualizes the mean flow field around the wind turbine reconstructed from the PPI scans. A velocity deficit, which is a common characteristic of a turbine wake, is clearly visible in the figure. Furthermore, a dual peak is also observed in the near-wake region up to x/D = 2. It can be appreciated that the wake region persists approximately 7D downstream from the wind turbine. Note that, for x/D ≈ 2, the azimuth angles are approximately perpendicular to the incoming flow and, therefore, the measurement error is high in this region.   Figure 19 visualizes the mean flow field around the wind turbine reconstructed from the PPI scans. A velocity deficit, which is a common characteristic of a turbine wake, is clearly visible in the figure. Furthermore, a dual peak is also observed in the near-wake region up to ⁄ = 2. It can be appreciated that the wake region persists approximately 7D downstream from the wind turbine. Note that, for ⁄ ≈ 2, the azimuth angles are approximately perpendicular to the incoming flow and, therefore, the measurement error is high in this region.
Several other studies have performed scanning Doppler LiDAR-based measurements of wind turbine wakes. However, reconstruction of velocity fields from LiDAR-measured LOS velocities poses a great challenge. Therefore, most of the studies simply tried to analyze and interpret LOS velocities, and did not perform any further processing of the measured data [26][27][28]. As is also mentioned in the Introduction, the assumption of constant wind direction assumption during the measurement period and the application of multiple LiDARs are also receiving a lot of interest in the wind energy community [29,30]. The method of wind turbine wake measurement with a single Scanning Doppler LiDAR for arbitrary directions proposed in this study is therefore different from those in the literature.
Finally, the measurement data is compared against four analytical wake models proposed by Jensen [31], Frandsen et al. [32], Bastankhah and Porte-Agel [33], and Qian and Ishihara [34]. As with Several other studies have performed scanning Doppler LiDAR-based measurements of wind turbine wakes. However, reconstruction of velocity fields from LiDAR-measured LOS velocities poses a great challenge. Therefore, most of the studies simply tried to analyze and interpret LOS velocities, and did not perform any further processing of the measured data [26][27][28]. As is also mentioned in the Introduction, the assumption of constant wind direction assumption during the measurement period and the application of multiple LiDARs are also receiving a lot of interest in the wind energy community [29,30]. The method of wind turbine wake measurement with a single Scanning Doppler LiDAR for arbitrary directions proposed in this study is therefore different from those in the literature.
Finally, the measurement data is compared against four analytical wake models proposed by Jensen [31], Frandsen et al. [32], Bastankhah and Porte-Agel [33], and Qian and Ishihara [34]. As with most other wake models, the wake models validated in this study require one or more of the following three input parameters: mean inflow velocity, inflow turbulence intensity, and thrust coefficient (C T ) of the turbine. The inflow wind speed and turbulence intensity measured by the met mast at the hub height levels are used. Figure 20a,b show the time series of the 10-min-averaged wind speed and turbulence intensity. The average inflow velocity during the comparison period was 8.2 m/s, and the turbulence intensity was 3.5%. From the C T − u characteristic of the turbine, a thrust coefficient value of 0.8 was used [35]. Figure 20c shows one-minute-averaged wind directions and yaw directions of the turbine during the period when the wind turbine wake data was collected. Figure 20d shows the yaw misalignment. The distribution of the yaw misalignment during the evaluation period was roughly symmetrical on the two sides of 0 • and the maximum yaw error was ±10 • . The wake model proposed by Qian and Ishihara [34] is able to account for the yaw misalignments. Therefore, for an accurate comparison, the yaw error distribution is used to construct weighted average velocity profiles from this wake model.
Atmosphere 2020, 11, x 16 of 20 most other wake models, the wake models validated in this study require one or more of the following three input parameters: mean inflow velocity, inflow turbulence intensity, and thrust coefficient ( ) of the turbine. The inflow wind speed and turbulence intensity measured by the met mast at the hub height levels are used. Figure 20a,b show the time series of the 10-min-averaged wind speed and turbulence intensity. The average inflow velocity during the comparison period was 8.2 m/s, and the turbulence intensity was 3.5%. From the characteristic of the turbine, a thrust coefficient value of 0.8 was used [35]. Figure 20c shows one-minute-averaged wind directions and yaw directions of the turbine during the period when the wind turbine wake data was collected. Figure 20d shows the yaw misalignment. The distribution of the yaw misalignment during the evaluation period was roughly symmetrical on the two sides of 0° and the maximum yaw error was ±10°. The wake model proposed by Qian and Ishihara [34] is able to account for the yaw misalignments. Therefore, for an accurate comparison, the yaw error distribution is used to construct weighted average velocity profiles from this wake model. Note that Jensen's model, which assumes top-hat distribution of velocity deficit, is simply a function of and downstream distance from the turbine. The model proposed by Frandsen et al. [32] is based on the conservation of mass and momentum but is also a function of and assumes top-hat velocity deficit distribution. In the model by Bastankhah and Porte-Agel [33], they employed Gaussian distribution to represent velocity deficit profiles. Although the model of Qian and Ishihara also employed Gaussian distribution for velocity deficit, the model and the wake growth rates are the function of both and inflow turbulence intensity. Note that Jensen's model, which assumes top-hat distribution of velocity deficit, is simply a function of C T and downstream distance from the turbine. The model proposed by Frandsen et al. [32] is based on the conservation of mass and momentum but is also a function of C T and assumes top-hat velocity deficit distribution. In the model by Bastankhah and Porte-Agel [33], they employed Gaussian distribution to represent velocity deficit profiles. Although the model of Qian and Ishihara also employed Gaussian distribution for velocity deficit, the model and the wake growth rates are the function of both C T and inflow turbulence intensity. Figure 21 shows the horizontal profiles of normalized velocity deficit at the hub height in the wake of the wind turbine. Here, the normalized velocity deficit is defined as: where U 0,h is the average inflow velocity which is equal to 8.2 m/s during the current comparison period, U is the streamwise velocity, and ∆U is the velocity deficit. Profiles obtained from scanning LiDAR are compared with those predicted by the wake models. The models of Jensen and Frandsen et al. underestimate the velocity deficit-compared to the measurement-around the wake center. Similar tendency has been observed in other studies on comparison of these two wake models with the large-eddy simulations (LES) (see, e.g., Ref. [33]). The performance of these models can be further improved by tuning the wake growth rate parameter, as done in some other studies (see Ref. [36] for the review), but the top-hat assumption of velocity deficit will never allow them to reproduce the Gaussian nature of the measured profiles. It can be appreciated that the velocity-deficit profiles predicted by the models of Bastankhah and Porte-Agel and Qian and Ishihara show similar trend as the measured profiles. The model of Bastankhah and Porte-Agel, however, overestimates the velocity deficit at x/D =1 and significantly underestimates the profiles further downstream. The discrepancy can be attributed to the fact that, in this study, the wake growth rate parameters were not obtained from fitting them to the current measurements, instead the parameters from the closest LES case in Ref. [33] were used. The model of Qian and Ishihara shows the most favorable agreement with the measured profiles, though even this model visibly overestimate or underestimate the velocity deficit compared to the measured profiles. The velocity deficit is maximum at x/D = 3, though a deficit can be observed even at x/D = 6. At x/D = 1.0, close to the rotor, a horizontal profile with a dual peak can be observed in the measurement. This dual peak is the result of an uneven distribution of force on the rotor in the radial direction. The near-wake flow is significantly affected by the specific blade aerodynamics, especially when the ambient turbulence is low. This has also been observed in experimental studies as shown by Ishihara and Qian [37]. However, these differences between the wake model and measurements in the near wake region are not very relevant when considering engineering applications of the model. The reason being that the turbine spacing in actual wind farms is always larger than 3D, and, in most cases, the streamwise spacing is larger than 7D [38].  Figure 21 shows the horizontal profiles of normalized velocity deficit at the hub height in the wake of the wind turbine. Here, the normalized velocity deficit is defined as:

Conclusions
The current work evaluated the performance of scanning Doppler LiDAR for wind resource measurement and analysis for wind energy applications. The measurement results were also compared with WRF simulations and wake models.
Firstly, 10-min-averaged wind speeds measured by scanning LiDAR in DBS mode were validated against the measurements of an existing V1 LiDAR. It was observed that the quality of measurements degraded when the minimum threshold value of the CNR was increased. As a result, the agreement between the scanning LiDAR and the V1 LiDAR measurements was better when the availability of the scanning LiDAR was higher. However, even with a lower availability of 20% (30 data per 10 min), fairly good 10-min-averaged wind speed data could be produced. Vertical wind profiles were measured in DBS-mode and then investigated for wind blowing from the land and that from the sea. It was found that the wind speed was higher when the wind was blowing from the sea than when it was blowing from the land. The simulated wind profiles by WRF showed a similar tendency to the LiDAR measurements, but they overestimated the measurements at higher altitudes when the wind was blowing from the land.
A method for the reconstruction of wind speeds from conical and slice scans (also known as PPI and RHI) was proposed in this study. The method assumes the wind direction at every point in a scan to be the same as that measured by another wind vane. The proposed velocity retrieval method with RHI and wind directions obtained from the wind vane mounted on the meteorological tower was then used to investigate the effect of coastal terrain on the ABL by measuring the flow field between the LiDAR and the coastline. The coastal terrain was responsible for the development of an IBL resulting in strong shear in the region close to the coast. It was also found that this IBL persisted for approximately 2000 m offshore.
Finally, the flow field around the wind turbine was measured using multiple PPI scans together with wind directions obtained from the wind vane mounted on the meteorological tower and processed to characterize the wind turbine wake. A velocity deficit in the wind turbine wake was clearly visible along with a dual peak in the near-wake region in the flow field visualization. Horizontal profiles of velocity deficit in the wake of the wind turbine showed common characteristics observed in earlier numerical and experimental studies. The LiDAR measurement was also used to validate the wake models. Two of the Gaussian wake models showed better agreement with the measurements than the models with top-hat velocity distribution.