Accuracy Comparison of Coastal Wind Speeds between WRF Simulations Using Different Input Datasets in Japan

In order to improve the accuracy of the wind speed simulated by a mesoscale model for the wind resource assessment in coastal areas, this study evaluated the effectiveness of using the Japan Meteorological Agency (JMA)’s latest and finest (2 km × 2 km) grid point value (GPV) data, produced from the local forecast model (LFM) as input data to the mesoscale model. The evaluation was performed using wind lidar measurements at two sites located on the coasts of the Sea of Japan and Pacific Ocean. The accuracy of the LFM–GPV was first compared with that of two products from the JMA Meso Scale Model (MSM) (5 km × 5 km): MSM-GPV and mesoscale analysis (MANAL). Consequently, it was shown that LFM–GPV exhibited the most accurate wind speeds against lidar measurements. Next, dynamical downscaling simulations were performed using the weather research and forecasting model (WRF) forced by the three datasets above, and their results were compared. As compared to the GPVs, it was found that the WRF dynamical downscaling simulation using them as input can improve the accuracy of the coastal wind speeds. This was attributed to the advantage of the WRF simulation to improve the negative bias from the input data, especially for the winds blowing from the sea sectors. It was also found that even if the LFM–GPV is used as an input to the WRF simulation, it does not always reproduce more accurate wind speeds, as compared to the simulations using the other two datasets. This result is partly owing to the tendency of WRF to overestimate the wind speed over land, thus obscuring the higher accuracy of the LFM–GPV. It was also shown that the overestimation tendency cannot be improved by only changing the nudging methods or the planetary boundary layer schemes in WRF. These results indicate that it may be difficult to utilize the LFM–GPV in the WRF wind simulation, unless the overestimation tendency of WRF itself is improved first.


Introduction
In March 2017, the New Energy and Industrial Technology Development Organization (NEDO) released the Japanese offshore wind resource map, NEDO Offshore Wind Information System (NeoWins) [1]. It provides the annual mean wind speed distributions and detailed wind climate statistics at five heights from 60 to 140 m, with a spatial resolution of 0.5 km × 0.5 km for coastal waters within 30 km from the coast, and with a spatial resolution of 0.100 • × 0.125 • for open oceans within the Japanese Economic Exclusive Zone. In addition, NeoWins offers wind information and any information regarding natural and social environmental conditions required for offshore wind energy development. Thus, NeoWins has become a powerful tool for planning offshore wind farms in Japan.
The wind information of NeoWins was produced using the weather research and forecasting model (WRF) [2], which is a mesoscale meteorological model frequently used for offshore wind resource assessment (e.g., [3][4][5][6][7][8]). The accuracy of wind speeds with a spatial resolution of 0.5 km × 0.5 km was examined using wind turbine hub-height measurements from three meteorological masts and a wind lidar at four coastal sites [9][10][11]. Subsequently, it was confirmed that the target accuracy of the map (the relative bias is within ±5% for the annual mean wind speed at a hub height) was achieved at all the sites. This result means that though the resolution of 0.5 km × 0.5 km is in the so-called "gray zone" [12], and the WRF simulation exhibits accurate wind speeds. However, considering the spread of the offshore wind development in the near future, it is important to obtain more accurate wind speed simulated by WRF, as WRF is used for dynamical downscaling for a specific site and will used for the next-generation wind resource map.
The choice of input data is crucial for a dynamical downscaling. Thus, the sensitivity of the input data to the mesoscale model simulation has often been discussed in previous studies (e.g., [13][14][15][16][17]). In NeoWins, mesoscale analysis (MANAL) [18], which is the objective analysis with a spatial resolution of 5 km × 5 km produced by the Japan Meteorological Agency (JMA), is used as initial and lateral boundary conditions to force WRF. Akimoto and Kusaka (2010) [19] showed that the use of MANAL for the downscaling provides more accurate weather conditions in the Kanto Plain, Japan, than that of the National Centers for Environmental Prediction (NCEP)-final operational global analysis (FNL) with a spatial resolution of 1 • × 1 • [20]. In addition to MANAL, two other grid point values (GPVs) can be utilized as input data to WRF: Local forecast model (LFM)-GPV [21] and meso scale model (MSM)-GPV [22], which are numerical weather prediction data from the JMA LFM and MSM, respectively. The LFM is the latest and finest nonhydrostatic model introduced in August 2012, and is designed to allow for the rapid production and frequent updating of analyses at a higher resolution (2 km × 2 km, hourly). Our previous study [23] evaluated the accuracy of wind speeds derived from the LFM-GPV and the MSM-GPV using in-situ onshore and offshore measurements obtained at the low altitudes less than 40 m above ground level (AGL). They indicated that the wind speeds of the LFM-GPV were more accurate than those of the MSM-GPV for both onshore and offshore, and presumed that the result was attributed to higher temporal and spatial resolutions and the assimilation of data from the automated meteorological data acquisition system (AMeDAS), operated by the JMA. Therefore, we expect that the accuracy of wind speed in the WRF simulation can be improved partially using the LFM-GPV as input data.
This study aimed to examine whether using the LFM-GPV would result in a more accurate WRF simulation for wind resource assessment in coastal areas, as compared to the cases using the MSM-GPV and MANAL. The wind speeds of the LFM-GPV, MSM-GPV, and MANAL, and the WRF simulations forced by them were compared with wind lidar measurements at two coastal sites in Japan. In-situ observations, input data, model configurations, and evaluation methods are described in Section 2. The accuracies of the GPV and WRF wind speeds are compared and evaluated in Section 3. Section 4 discusses the possible methods to improve the accuracy of the WRF simulation, and the conclusions of this study are presented in Section 5.  Figure 1 shows the maps of Japan and target areas, i.e., Niigata and Ibaraki, which are located on the coasts of the Sea of Japan and the Pacific Ocean, respectively. They exhibit different wind climatology and are thus suitable for the evaluation of the WRF simulation in the coasts of Japan. Table 1 shows the specifications of lidar observations at both sites. At Niigata, lidar measurements were performed on flat land at approximately 300 m inland from the coastline. The ZephIR 300 lidar [24] was used to measure the wind speed and direction at the heights of 40, 58, 80, 95, and 145 m above mean sea level (AMSL). At Hazaki, the WindCube WLS7-86 lidar [25] was installed at the tip of a 427-m long pier owned by the Hazaki Oceanographic Research Station [26] of the Port and Airport Research Institute (PARI). Vertical wind profiles from 47 to 207 m AMSL were measured with an interval of 20 m. The lidar measurements are accurately calibrated by using the different in-situ measurements:

In-Situ Observation Data
The meteorological mast at Niigata (Mito et al., 2016 [27]), and another lidar measurement at Ibaraki (Shimada et al., 2018 [28]). The heights of 95 m at Niigata and 107 m at Ibaraki are referred to as a "hub height" in this study. The lidar data of the 10-min average wind speed and direction at every height for a whole year (from October 2015 to September 2016) was used for the analyses.  Figure 1 shows the maps of Japan and target areas, i.e., Niigata and Ibaraki, which are located on the coasts of the Sea of Japan and the Pacific Ocean, respectively. They exhibit different wind climatology and are thus suitable for the evaluation of the WRF simulation in the coasts of Japan. Table 1 shows the specifications of lidar observations at both sites. At Niigata, lidar measurements were performed on flat land at approximately 300 m inland from the coastline. The ZephIR 300 lidar [24] was used to measure the wind speed and direction at the heights of 40, 58, 80, 95, and 145 m above mean sea level (AMSL). At Hazaki, the WindCube WLS7-86 lidar [25] was installed at the tip of a 427-m long pier owned by the Hazaki Oceanographic Research Station [26] of the Port and Airport Research Institute (PARI). Vertical wind profiles from 47 to 207 m AMSL were measured with an interval of 20 m. The lidar measurements are accurately calibrated by using the different insitu measurements: The meteorological mast at Niigata (Mito et al., 2016 [27]), and another lidar measurement at Ibaraki (Shimada et al., 2018 [28]). The heights of 95 m at Niigata and 107 m at Ibaraki are referred to as a "hub height" in this study. The lidar data of the 10-min average wind speed and direction at every height for a whole year (from October 2015 to September 2016) was used for the analyses.   Figure 2 shows the annual wind roses at Niigata (95 m) and Ibaraki (107 m). The prevailing wind directions are northwest and east for Niigata, and north-northeast and south-southwest for Ibaraki. The NW winds at Niigata have a higher speed than the E winds, while the NNE and SSW wind   Figure 2 shows the annual wind roses at Niigata (95 m) and Ibaraki (107 m). The prevailing wind directions are northwest and east for Niigata, and north-northeast and south-southwest for Ibaraki. The NW winds at Niigata have a higher speed than the E winds, while the NNE and SSW wind speeds at Ibaraki are similar. To investigate the accuracy of wind speeds from GPVs and WRF simulations, which is considered to depend on wind direction, it is convenient to define the land and sea sectors for wind direction. Considering the shape of the coastlines at Niigata and Ibaraki, the sectors are defined as indicated by the red dot lines in Figure 1b,c. At Niigata, the land and sea sectors are 0 • ≤ θ < 225 • and 225 • ≤ θ < 360 • . At Ibaraki, the land sector is 150 • ≤ θ < 330 • and the sea sector is 0 • ≤ θ < 150 • or 330 • ≤ θ < 360 • . When the wind directions at all of the measurement heights are classified into one of each sector, the wind is defined as that for the sector. Table 2 shows the number of samples of wind speed data for each sector. Niigata has more winds from the land sector (ENE to E), while Ibaraki has more winds from the sea sector (NNE to NE). "Other" in the table implies cases where the sector cannot be identified owing to non-uniform wind directions with height. speeds at Ibaraki are similar. To investigate the accuracy of wind speeds from GPVs and WRF simulations, which is considered to depend on wind direction, it is convenient to define the land and sea sectors for wind direction. Considering the shape of the coastlines at Niigata and Ibaraki, the sectors are defined as indicated by the red dot lines in Figure 1b,c. At Niigata, the land and sea sectors are 0° ≤ θ < 225° and 225° ≤ θ < 360°. At Ibaraki, the land sector is 150° ≤ θ < 330° and the sea sector is 0° ≤ θ < 150° or 330° ≤ θ < 360°. When the wind directions at all of the measurement heights are classified into one of each sector, the wind is defined as that for the sector. Table 2 shows the number of samples of wind speed data for each sector. Niigata has more winds from the land sector (ENE to E), while Ibaraki has more winds from the sea sector (NNE to NE). "Other" in the table implies cases where the sector cannot be identified owing to non-uniform wind directions with height.

JMA GPVs
The LFM and MSM are the current forecasting models for disaster prevention and aviation safety, operated by the JMA. The horizontal resolution for the surface and pressure layer is 2 km × 2 km (1531 × 1301 grids) in the LFM and 5 km × 5 km (817 × 661 grids) in the MSM. The number of vertical levels in the LFM is 58 from the surface to 20.2 km and that of the MSM is 48 from the surface to 21.8 km. The LFM uses forecasting results from the MSM as the lateral boundary condition, and is designed to provide shorter-range forecasting with a higher spatio-temporal resolution. In the forecast cycle, the first guesses from the LFM and MSM are combined with the latest multiple observations through data assimilation, producing the local analysis (LANAL) and MANAL, respectively. It is noteworthy that only LANAL assimilates surface data from AMeDAS to improve the reproducibility of local-scale phenomena. More details on the models and data assimilation systems are available in the JMA report [29].
The GPVs from the LFM and MSM, including LANAL and MANAL, are distributed by the Japan Meteorological Business Support Center (JMBSC). Because LANAL is a new dataset and available only after January 2018, this study used three GPVs from the LFM (hereinafter called LFM-GPV),

JMA GPVs
The LFM and MSM are the current forecasting models for disaster prevention and aviation safety, operated by the JMA. The horizontal resolution for the surface and pressure layer is 2 km × 2 km (1531 × 1301 grids) in the LFM and 5 km × 5 km (817 × 661 grids) in the MSM. The number of vertical levels in the LFM is 58 from the surface to 20.2 km and that of the MSM is 48 from the surface to 21.8 km. The LFM uses forecasting results from the MSM as the lateral boundary condition, and is designed to provide shorter-range forecasting with a higher spatio-temporal resolution. In the forecast cycle, the first guesses from the LFM and MSM are combined with the latest multiple observations through data assimilation, producing the local analysis (LANAL) and MANAL, respectively. It is noteworthy that only LANAL assimilates surface data from AMeDAS to improve the reproducibility of local-scale phenomena. More details on the models and data assimilation systems are available in the JMA report [29]. The GPVs from the LFM and MSM, including LANAL and MANAL, are distributed by the Japan Meteorological Business Support Center (JMBSC). Because LANAL is a new dataset and available only after January 2018, this study used three GPVs from the LFM (hereinafter called LFM-GPV), MSM (hereinafter called MSM-GPV), and MANAL, which are available for the study period from October 2015 to September 2016. Table 3 shows the specifications of those three GPVs. The GPVs include the islands of Japan and their surrounding waters, having horizontal resolutions of 0.025 • × 0.020 • (LFM-GPV), 0.0625 • × 0.0500 • (MSM-GPV), and 5 km × 5 km (MANAL) for the surface layer. The number of vertical layers of the LFM-GPV and MSM-GPV are both 17, and includes the surface layer and 16 pressure levels (1000, 975, 950, 925, 900, 850, 800, 700, 600, 500, 400, 300 250, 200, 150, and 100 hPa). MANAL has the same levels except for 975 hPa. The primary parameters included in the datasets are sea-level pressure, geopotential height, wind speed, temperature, relative humidity, etc. As for the LFM-GPV and MSM-GPV, only the initial values of each forecasting cycle (not forecast values) were used in this study. Table 3. Specifications of the LFM-GPV, MSM-GPV, and MANAL [18,21,22]. To extract the wind speed at an arbitrary point from the GPVs, bilinear interpolation was used for the horizontal plane. For height, a linear-log interpolation was applied in the U-logZ plane, using three wind speeds at 10 m AGL and at the AGL heights of the two lowest pressure levels, following our previous study [23].

WRF Configuration
In this study, the Advanced Research WRF (ARW) version 3.8.1 [2] was used as a meteorological mesoscale model. The 12 one-month WRF simulations were performed for whole a year from October 2015 to September 2016. Figure 1b,c show the model domains for Niigata and Ibaraki, and Table 4 shows the model configuration used for the simulations. The model configuration follows that used for the NeoWins project [1]. The roughness length used in the WRF simulation is based in the land-use data from the National Land Numerical Information (NLNI), which has 12 categories in land-use [30]. The value of roughness length was set so that it was similar to the original value for the United States Geological Survey (USGS) land-use category [31]. The soil temperature and moisture used for the Noah land surface model [32] were obtained from NCEP-FNL (1 • × 1 • , 6 hourly). The sea surface temperature (SST) based on a moderate resolution imaging spectroradiometer (MODIS) (0.02 • × 0.02 • , daily) (MOSST), developed by the National Institute of Advanced Industrial Science and Technology (AIST) and Kobe University [33], was used as the lower boundary condition. The number of vertical levels was 40 from the surface to 100 hPa, and the lowest levels were located at approximately 23, 73, 130, 199, and 287 m AGL. Four-dimensional data assimilation (FDDA) was conducted by applying the grid nudging to full levels in the outer domain, and to the levels above the top of the planetary boundary layer (PBL) in the inner domain, where the horizontal components of the wind speed (u, v), potential temperature (θ), and water vapor mixing ratio (q) are nudged toward the GPVs. As a PBL scheme, the Mellor-Yamada-Janjic (MYJ) scheme [34] was used because it has often been adopted in previous studies (e.g., [1, [35][36][37]).    Table 5 shows the input data and model grids used in the three WRF simulations. Considering the horizontal resolutions of the LFM-GPV (0.025 • × 0.020 • ), MSM-GPV (0.0625 • × 0.0500 • ), and MANAL (5 km × 5 km), the WRF simulation using the LFM-GPV as input (hereinafter called WRF-LFM) contains 1.5 to 0.5 km domains, while those using MSM-GPV and MANAL (WRF-MSM and WRF-MANAL, respectively) contains 2.5 to 0.5 km domains. In order to capture large-scale meteorological conditions from the GPVs, the outer domains were set to the rectangular area with the side of about 250 km. This study used the 3 hourly initial values of the GPVs as lateral boundary conditions to force WRF, and for data assimilation. For accuracy evaluation, hourly WRF wind speeds at observation points, expected from the inner domain using bilinear and linear-log interpolations, were compared with in-situ measurements. Table 5. Input data and model grids used in WRF-LFM, WRF-MSM, and WRF-MANAL simulations.

Evaluation Methods
The accuracies of the GPV and WRF wind speeds were evaluated using three statistical parameters: Bias, root mean square error (RMSE), and correlation coefficient. They are written as follows: where n is the number of data samples, W EST,i is the GPV or WRF wind speed, W OBS,i is the observed wind speed, and W EST and W OBS are temporal averages for each wind speed. In this study, the bias and RMSE are presented using relative values (%), i.e., those normalized by the observed mean wind speed, because a relative value is customarily used in the research of wind energy and its application.

Comparison between GPV Wind Speeds
First, we investigated the accuracy of wind speeds from the three GPVs. Note that the annual mean wind speed at Ibaraki (107 m) is 7.85 m/s, though that at Niigata cannot be specified in this study due to the agreement with the data providers. Figure 3 shows the annual biases, RMSEs, and correlation coefficients of the GPV wind speeds at Niigata (95 m) and Ibaraki (107 m). The annual biases of the LFM-GPV, MSM-GPV, and MANAL were −5.9%, −9.6%, and −7.5% at Niigata, respectively, and −4.8%, −8.8%, and −8.3% at Ibaraki, respectively. The biases were all negative, meaning that the GPVs underestimate the annual mean wind speed at both sites. Dividing the winds into those coming from the land and sea sectors, the biases of the LFM-GPV, MSM-GPV, and MANAL for the land-sector winds were −0.2%, −2.0%, and −1.7% at Niigata, respectively, and 5.2%, 3.6%, and 5.2% at Ibaraki, respectively, whereas those for the sea-sector winds were −9.5%, −14.6%, and −11.3% at Niigata, respectively, and −9.9%, −15.0%, and −15.2% at Ibaraki, respectively. All of the GPVs remarkably underestimated the wind speed for the sea-sector winds at both sites, whereas the biases of the land-sector winds were closer to zero. These results indicate that the negative biases of the annual mean wind speeds are primarily caused by the sea-sector winds. The results of the GPVs underestimating the sea-sector winds (winds having been blowing on the sea surface) and overestimating the land-sector winds (winds having been blowing on the land surface) are in good agreement with those shown in our previous study [23], although the study examined the surface winds. The only exception is the underestimation for the land-sector wind at Niigata. For this wind, the GPVs exhibited RMSEs of 41.2% to 42.5% and correlation coefficients of 0.65 to 0.69, i.e., a much lower accuracy than those of other winds. The low accuracy is attributed to the low wind speed of the predominant ENE-E winds (as shown in Figure 2a), which are primarily land breezes at night (shown later). It is speculated that the MSM-GPV and MANAL cannot reproduce such local winds well, owing to its low spatial resolutions. Figure 4 shows the vertical profiles of the GPV and the observed mean wind speeds and its difference (bias) at Niigata and Ibaraki. Here, the wind speeds were normalized by the observed mean wind speed at the hub height (95 m at Niigata and 107 m at Ibaraki). The biases were normalized by the observed mean wind speed at each height. From this figure, it is found that the three GPVs exhibited negative annual biases; that is, the GPVs underestimate the annual mean wind speeds at both the hub heights and other heights. In the bias profiles for the land-sector and sea-sector winds, the absolute values were found to increase toward the surface, including positive biases found for the land-sector winds at Ibaraki. This finding is consistent with the result of Shimada et al. (2011b) [36] that indicated a similar tendency of MANAL wind speeds using wind profiler measurements at 10 JMA stations.  Figure 4 shows the vertical profiles of the GPV and the observed mean wind speeds and its difference (bias) at Niigata and Ibaraki. Here, the wind speeds were normalized by the observed mean wind speed at the hub height (95 m at Niigata and 107 m at Ibaraki). The biases were normalized by the observed mean wind speed at each height. From this figure, it is found that the three GPVs exhibited negative annual biases; that is, the GPVs underestimate the annual mean wind speeds at both the hub heights and other heights. In the bias profiles for the land-sector and sea-sector winds, the absolute values were found to increase toward the surface, including positive biases found for the land-sector winds at Ibaraki. This finding is consistent with the result of Shimada et al. (2011b) [36] that indicated a similar tendency of MANAL wind speeds using wind profiler measurements at 10 JMA stations.   Figure 4 shows the vertical profiles of the GPV and the observed mean wind speeds and its difference (bias) at Niigata and Ibaraki. Here, the wind speeds were normalized by the observed mean wind speed at the hub height (95 m at Niigata and 107 m at Ibaraki). The biases were normalized by the observed mean wind speed at each height. From this figure, it is found that the three GPVs exhibited negative annual biases; that is, the GPVs underestimate the annual mean wind speeds at both the hub heights and other heights. In the bias profiles for the land-sector and sea-sector winds, the absolute values were found to increase toward the surface, including positive biases found for the land-sector winds at Ibaraki. This finding is consistent with the result of Shimada et al. (2011b) [36] that indicated a similar tendency of MANAL wind speeds using wind profiler measurements at 10 JMA stations.  Next, the accuracies of the three GPV wind speeds were compared. Figures 3 and 4 show that the LFM-GPV tends to exhibit higher wind speeds than the other GPVs, and hence exhibited the smallest (reduced negative) biases. The annual RMSEs of the LFM-GPV, MSM-GPV, and MANAL were 31.2%, 33.4%, and 32.7% at Niigata, respectively, and 24.7%, 24.6%, and 25.6% at Ibaraki, respectively; meanwhile, the annual correlation coefficients were 0.84, 0.81, and 0.82 at Niigata, respectively, and 0.85, 0.86, and 0.85 at Ibaraki, respectively. It is found that the two statistical parameters were not significantly different between the three GPVs, and the same holds for each of the land-sector and sea-sector winds. Therefore, the accuracy comparison between the GPVs results in the conclusion that the LFM-GPV provides the most accurate wind speeds at all the heights of both sites. The high accuracy of the LFM-GPV is likely owing to its higher spatio-temporal resolution and the data assimilation of AMeDAS data, as discussed in our previous study [23].

Comparison between GPV and WRF Wind Speeds
This subsection compares the wind speeds simulated by WRF using the three GPVs discussed in the previous subsection, as the input. Figure 5 shows the annual biases, RMSEs, and correlation coefficients for the three WRF wind speeds: WRF-LFM, WRF-MSM, and WRF-MANAL. The annual biases for the WRF-LFM, WRF-MSM, and WRF-MANAL were −2.8%, −0.1%, and 1.5% at Niigata, respectively, and 3.7%, 1.2%, and 2.7% at Ibaraki, respectively. Compared to the biases of GPVs (Figure 3), the absolute values of the biases, especially those for the sea-sector wind, were significantly smaller at both sites. When comparing Figure 4 (showing the GPV wind profiles) and Figure 6 (showing the WRF wind profiles), improvements are shown for all the observation heights. This appears to be the effect of the higher spatio-temporal downscaling with WRF. Figure 6 also shows that the biases for the land-sector wind were negative at Niigata and positive at Ibaraki for all the heights, and the bias tendency increased toward the surface. However, the former result in that the bias was negative is inconsistent with the results of previous studies [11,[35][36][37] that reported the tendency of a positive bias over land for the WRF wind speeds. To explain this discrepancy, Figure 7 shows the accuracy variations of the WRF wind speed in a day at Niigata. As shown, a negative bias appeared from 22 to 8 LST, primarily at night. Because the number of samples was significantly larger at night (Figure 7a) and the wind speed from the land sector was low (Figure  2a), the wind at night is primarily considered as a typical land breeze, blowing from the ENE or E. The result above indicates that WRF tends to underestimate the land breeze circulating at night. Meanwhile, during the day, positive biases prevail that are likely the same type of positive bias as those reported in the previous studies [11,[35][36][37]. In other words, the negative bias found at Niigata is caused by another tendency of WRF, which generally tends to overestimate the wind speed over land. Next, the accuracies of the three GPV wind speeds were compared. Figures 3 and 4 show that the LFM-GPV tends to exhibit higher wind speeds than the other GPVs, and hence exhibited the smallest (reduced negative) biases. The annual RMSEs of the LFM-GPV, MSM-GPV, and MANAL were 31.2%, 33.4%, and 32.7% at Niigata, respectively, and 24.7%, 24.6%, and 25.6% at Ibaraki, respectively; meanwhile, the annual correlation coefficients were 0.84, 0.81, and 0.82 at Niigata, respectively, and 0.85, 0.86, and 0.85 at Ibaraki, respectively. It is found that the two statistical parameters were not significantly different between the three GPVs, and the same holds for each of the land-sector and sea-sector winds. Therefore, the accuracy comparison between the GPVs results in the conclusion that the LFM-GPV provides the most accurate wind speeds at all the heights of both sites. The high accuracy of the LFM-GPV is likely owing to its higher spatio-temporal resolution and the data assimilation of AMeDAS data, as discussed in our previous study [23].

Comparison between GPV and WRF Wind Speeds
This subsection compares the wind speeds simulated by WRF using the three GPVs discussed in the previous subsection, as the input. Figure 5 shows the annual biases, RMSEs, and correlation coefficients for the three WRF wind speeds: WRF-LFM, WRF-MSM, and WRF-MANAL. The annual biases for the WRF-LFM, WRF-MSM, and WRF-MANAL were −2.8%, −0.1%, and 1.5% at Niigata, respectively, and 3.7%, 1.2%, and 2.7% at Ibaraki, respectively. Compared to the biases of GPVs (Figure 3), the absolute values of the biases, especially those for the sea-sector wind, were significantly smaller at both sites. When comparing Figure 4 (showing the GPV wind profiles) and Figure 6 (showing the WRF wind profiles), improvements are shown for all the observation heights. This appears to be the effect of the higher spatio-temporal downscaling with WRF.     Figure 6 also shows that the biases for the land-sector wind were negative at Niigata and positive at Ibaraki for all the heights, and the bias tendency increased toward the surface. However, the former result in that the bias was negative is inconsistent with the results of previous studies [11,[35][36][37] that reported the tendency of a positive bias over land for the WRF wind speeds. To explain this discrepancy, Figure 7 shows the accuracy variations of the WRF wind speed in a day at Niigata. As shown, a negative bias appeared from 22 to 8 LST, primarily at night. Because the number of samples was significantly larger at night (Figure 7a) and the wind speed from the land sector was low (Figure 2a), the wind at night is primarily considered as a typical land breeze, blowing from the ENE or E. The result above indicates that WRF tends to underestimate the land breeze circulating at night. Meanwhile, during the day, positive biases prevail that are likely the same type of positive bias as those reported in the previous studies [11,[35][36][37]. In other words, the negative bias found at Niigata is caused by another tendency of WRF, which generally tends to overestimate the wind speed over land.

Comparison between WRF Wind Speeds
As mentioned in the introduction, a primary purpose of this study was to examine whether using the LFM-GPV as the input for the WRF simulation can improve the accuracy of wind speed simulated by WRF. Thus, this subsection describes the sensitivity of the WRF simulations to different input datasets. When focusing on the land-sector wind at Ibaraki, where such an overestimation is typically shown (Figure 6b), it was found that the WRF-LFM exhibited slightly lower wind speeds than the WRF-MSM and WRF-MANAL; however, the differences were slight. Figure 8 shows the spatial distributions of the relative differences between the WRF simulations in the 100-m-height (AMSL) wind speed averaged for the times when the land-sector wind was observed at Ibaraki. Figure 8a,b show the values of the WRF-LFM minus those of the WRF-MSM and those of the WRF-LFM minus those of the WRF-MANAL, respectively. From both figures, it is found that the differences were primarily approximately ±2% over land except for the vicinity of the lateral boundaries; further, the negative values, indicating WRF-LFM < WRF-MSM/WRF-MANAL, encompass larger parts of the land. These results suggest that the LFM-GPV can facilitate a reduction of the tendency of WRF to overestimate the wind speed over land as compared to the other GPVs; however, the effect is limited and is insufficient for a complete improvement.

Comparison between WRF Wind Speeds
As mentioned in the introduction, a primary purpose of this study was to examine whether using the LFM-GPV as the input for the WRF simulation can improve the accuracy of wind speed simulated by WRF. Thus, this subsection describes the sensitivity of the WRF simulations to different input datasets. When focusing on the land-sector wind at Ibaraki, where such an overestimation is typically shown (Figure 6b), it was found that the WRF-LFM exhibited slightly lower wind speeds than the WRF-MSM and WRF-MANAL; however, the differences were slight. Figure 8 shows the spatial distributions of the relative differences between the WRF simulations in the 100-m-height (AMSL) wind speed averaged for the times when the land-sector wind was observed at Ibaraki. Figure 8a The WRF-LFM exhibits another problem for the sea-sector wind at Niigata, as shown in Figure  6a (Table 4), the tendency of the WRF-LFM to underestimate the sea-sector wind speed is likely to be caused by the grid nudging of the low wind speeds of the LFM-GPV in the WRF simulations.   (Table 4), the tendency of the WRF-LFM to underestimate the sea-sector wind speed is likely to be caused by the grid nudging of the low wind speeds of the LFM-GPV in the WRF simulations. The WRF-LFM exhibits another problem for the sea-sector wind at Niigata, as shown in Figure  6a That is, the LFM-GPV exhibits remarkably lower wind speeds than the MSM-GPV in the WRF simulation domains, when the wind blows from the sea sector at Niigata. Considering the model configuration (Table 4), the tendency of the WRF-LFM to underestimate the sea-sector wind speed is likely to be caused by the grid nudging of the low wind speeds of the LFM-GPV in the WRF simulations.

Discussion on Overestimation for Wind Speeds over Land
The accuracy comparisons in the previous section indicate that the LFM-GPV exhibited a higher accuracy in wind speed at the two coastal sites, compared to the MSM-GPV and MANAL. It was also shown that the WRF dynamical downscaling simulation using GPVs as input can improve the accuracy of the coastal wind speeds, as compared to the original GPVs. However, it is surprising to find that even if the LFM-GPV, exhibiting the highest accuracy of the three, was used for the WRF simulation, the WRF-LFM simulation did not clearly exhibit more accurate wind speeds than those using the other two GPVs. This result may be interpreted as follows: The effectiveness of using the LFM-GPV is obscured by the tendency of WRF to overestimate the wind speed over land (i.e., the land-sector wind speed in the previous section). In other words, it may be difficult to demonstrate the effectiveness of using the LFM-GPV in the WRF simulation, unless the WRF's overestimation tendency is improved. Actually, such an observation tendency can lead to an overestimation of the capacity factor of wind turbines, and consequently, to design an unprofitable layout of wind farms. In addition, a similar overestimation tendency was reported in previous studies [11,[35][36][37]. Shimada et al. (2011a [35]; 2011b [36]) showed the tendency of WRF to overestimate wind speeds in the PBL over land. Draxl et al. (2014) [37] indicated that the wind speeds simulated by WRF primarily exhibit positive biases at typical wind turbine heights in onshore areas. Kato et al. (2017) [11] speculated that the positive bias in coastal waters was caused by the advection of overestimated wind speeds from land, by demonstrating that the magnitude of the positive bias depends on the fetch from the coast. Thus, to obtain a method to reduce the overestimation, sensitivity simulations with different FDDA (nudging) options and PBL schemes were performed and the possible causes of overestimation are discussed.

Effect from Nudging Method
Grid nudging is an FDDA technique that bridges the gap between the model simulations and large-scale meteorological conditions provided from the input data. In the WRF dynamical downscaling simulation for wind resource assessments, the nudging method has frequently been used in previous studies (e.g., [16,[35][36][37]).
In WRF, the grid nudging of an arbitrary physical variable, which can be the horizontal component of the wind speed (u, v), potential temperature (θ), or water vapor mixing ratio (q), is implemented by adding an additional tendency term to the nudged variable's equation as follows: where p * α is the flux from of variable α, and the dry column pressure, p * = p s − p t ; p s , is the surface pressure, and p t is the constant pressure at the top of the model. F(α, x, t) is the physical forcing term, where x is the independent spatial variable and t is time. G α is a timescale controlling the nudging strength applied to variable α(G α = 1/t ) and is the nudging timescale. W represents the spatial and temporal weighting, where W = w xy w σ w t . The analysis quality factor, ε, ranging between 0 and 1, is based on the quality and distribution of the data used to produce the gridded analysis.α 0 is the observation analyzed to the grid and interpolated linearly in time.
For the simulations in this study (Table 4), the nudging coefficient, which corresponds to the nudging strength, is set to 0.0003 s −1 as the WRF default. Grid nudging is enabled for the entire outer domain, while it is excluded in the PBL in the inner domain; this method is called FDDA-DYNAMIC in this study. Specifically, the characteristic of FDDA-DYNAMIC is that the nudging space varies dynamically in response to the PBL height; that is, as the PBL becomes higher, the effect of nudging decreases accordingly in the inner domain, and vice versa. To be compared with the FDDA-DYNAMIC, another simulation is performed with an FDDA method, in which the nudging for the outer domain is the same as that for FDDA-DYNAMIC; however, for the inner domain, FDDA is excluded within the PBL defined below a specified height (set to 1500 m AGL in this study). This study refers to this method as FDDA-STATIC. These FDDA-DYNAMIC and FDDA-STATIC simulations were conducted and evaluated for the land-sector wind speeds at Ibaraki, as shown in Figure 6b-2. Figure 10 shows the vertical profiles of the bias in mean wind speeds for the land sector at Ibaraki, derived from the two WRF-MSM simulations using FDDA-DYNAMIC and FDDA-STATIC. No major differences are shown between the two simulations. A similar result was obtained in the spatial distribution of the relative difference between both nudging simulations in the 100 m height (AMSL) wind speed (not shown). It was therefore concluded that the tendency of WRF to overestimate the wind speed over land is not attributed to the FDDA-DYNAMIC used in Section 3.
outer domain is the same as that for FDDA-DYNAMIC; however, for the inner domain, FDDA is excluded within the PBL defined below a specified height (set to 1500 m AGL in this study). This study refers to this method as FDDA-STATIC. These FDDA-DYNAMIC and FDDA-STATIC simulations were conducted and evaluated for the land-sector wind speeds at Ibaraki, as shown in Figure 6b-2. Figure 10 shows the vertical profiles of the bias in mean wind speeds for the land sector at Ibaraki, derived from the two WRF-MSM simulations using FDDA-DYNAMIC and FDDA-STATIC. No major differences are shown between the two simulations. A similar result was obtained in the spatial distribution of the relative difference between both nudging simulations in the 100 m height (AMSL) wind speed (not shown). It was therefore concluded that the tendency of WRF to overestimate the wind speed over land is not attributed to the FDDA-DYNAMIC used in Section 3.

Effect from PBL Scheme
Next, to identify the sensitivity of the overestimated wind speed over land to the PBL scheme, further WRF simulations were conducted using different PBL schemes. Hu et al. (2006) [38] demonstrated that a positive bias is clearly found only when using the MYJ scheme; Floors et al. (2013) [39] demonstrated that the WRF simulations with different PBL schemes based on the firstand higher-order closure models performed differently for wind directions in the Danish coast. The results of both previous studies indicated that the choice in PBL scheme could affect the accuracy of onshore wind simulations. Thus, in this study, we conducted WRF-MSM simulations with the two different PBL schemes: The first-order Yonsei University (YSU) scheme (Hong et al., 2006 [40]) and the second-order Mellor-Yamada-Nakanishi-Niino Level-3 (MYNN3) scheme (Nakanishi et al., 2006 [41]). Figure 11 shows the vertical profiles of the bias for the land-sector wind at Ibaraki, derived from the WRF-MSM simulations with the YSU, MYJ, and MYNN3 schemes. No significant differences are indicated among the three simulations, and all of them shows that the positive bias increases in the downward direction. Figure 12 shows the spatial distribution of the relative difference among them in the 100 m height (AMSL) wind speed averaged for the times when the land-sector wind was observed at Ibaraki. Figure 12a,b shows the values of YSU minus MYJ and those of MYNN3 minus MYJ, respectively. Focusing on the land area, positive values prevail in both figures, although the latter exhibit larger values up to approximately 5%. This implies that the wind speeds derived from the YSU and MYNN3 schemes were higher than those from the MYJ scheme. Therefore, we conclude that WRF tends to overestimate the land-sector wind (i.e., wind blowing over land) regardless of the choice in PBL scheme.

Effect from PBL Scheme
Next, to identify the sensitivity of the overestimated wind speed over land to the PBL scheme, further WRF simulations were conducted using different PBL schemes. Hu et al. (2006) [38] demonstrated that a positive bias is clearly found only when using the MYJ scheme; Floors et al. (2013) [39] demonstrated that the WRF simulations with different PBL schemes based on the firstand higher-order closure models performed differently for wind directions in the Danish coast. The results of both previous studies indicated that the choice in PBL scheme could affect the accuracy of onshore wind simulations. Thus, in this study, we conducted WRF-MSM simulations with the two different PBL schemes: The first-order Yonsei University (YSU) scheme (Hong et al., 2006 [40]) and the second-order Mellor-Yamada-Nakanishi-Niino Level-3 (MYNN3) scheme (Nakanishi et al., 2006 [41]). Figure 11 shows the vertical profiles of the bias for the land-sector wind at Ibaraki, derived from the WRF-MSM simulations with the YSU, MYJ, and MYNN3 schemes. No significant differences are indicated among the three simulations, and all of them shows that the positive bias increases in the downward direction. Figure 12 shows the spatial distribution of the relative difference among them in the 100 m height (AMSL) wind speed averaged for the times when the land-sector wind was observed at Ibaraki. Figure 12a,b shows the values of YSU minus MYJ and those of MYNN3 minus MYJ, respectively. Focusing on the land area, positive values prevail in both figures, although the latter exhibit larger values up to approximately 5%. This implies that the wind speeds derived from the YSU and MYNN3 schemes were higher than those from the MYJ scheme. Therefore, we conclude that WRF tends to overestimate the land-sector wind (i.e., wind blowing over land) regardless of the choice in PBL scheme.

Other Possible Causes
The previous subsections indicated that WRF's tendency to overestimate the wind speed over land cannot be improved by only changing the nudging methods and PBL schemes. In other words, WRF's overestimation tendency appears to be attributed to other dynamic or thermodynamic causes in the simulation. It is thus necessary to verify any processes in the surface layer and PBL that can cause the overestimation of the wind speed over land. This final subsection discusses the possible causes of the overestimation, reviewing relevant studies in the past.
From a dynamical point of view, Hahmann et al. (2015) [6] reported that the displacement height and surface roughness length are not considered sufficiently over land in the WRF simulation. Floors et al. (2013) [39] and Varquez et al. (2014) [42] demonstrated that the correction of surface parameters based on observations could improve the accuracy of the wind conditions in WRF simulation. Shimada et al. (2015) [33] suggested using the topo-wind option (Jimenez and Dudhia, 2012 [43]) that considers an additional surface drag owing to unresolved subgrid-scale topography in the momentum equation, to reduce the surface wind speed over complex terrains. However, because this

Other Possible Causes
The previous subsections indicated that WRF's tendency to overestimate the wind speed over land cannot be improved by only changing the nudging methods and PBL schemes. In other words, WRF's overestimation tendency appears to be attributed to other dynamic or thermodynamic causes in the simulation. It is thus necessary to verify any processes in the surface layer and PBL that can cause the overestimation of the wind speed over land. This final subsection discusses the possible causes of the overestimation, reviewing relevant studies in the past.
From a dynamical point of view, Hahmann et al. (2015) [6] reported that the displacement height and surface roughness length are not considered sufficiently over land in the WRF simulation. Floors et al. (2013) [39] and Varquez et al. (2014) [42] demonstrated that the correction of surface parameters based on observations could improve the accuracy of the wind conditions in WRF simulation. Shimada et al. (2015) [33] suggested using the topo-wind option (Jimenez and Dudhia, 2012 [43]) that considers an additional surface drag owing to unresolved subgrid-scale topography in the momentum equation, to reduce the surface wind speed over complex terrains. However, because this

Other Possible Causes
The previous subsections indicated that WRF's tendency to overestimate the wind speed over land cannot be improved by only changing the nudging methods and PBL schemes. In other words, WRF's overestimation tendency appears to be attributed to other dynamic or thermodynamic causes in the simulation. It is thus necessary to verify any processes in the surface layer and PBL that can cause the overestimation of the wind speed over land. This final subsection discusses the possible causes of the overestimation, reviewing relevant studies in the past.
From a dynamical point of view, Hahmann et al. (2015) [6] reported that the displacement height and surface roughness length are not considered sufficiently over land in the WRF simulation. Floors et al. (2013) [39] and Varquez et al. (2014) [42] demonstrated that the correction of surface parameters based on observations could improve the accuracy of the wind conditions in WRF simulation. Shimada et al. (2015) [33] suggested using the topo-wind option (Jimenez and Dudhia, 2012 [43]) that considers an additional surface drag owing to unresolved subgrid-scale topography in the momentum equation, to reduce the surface wind speed over complex terrains. However, because this option is adjusted to Corporation Power Ltd., and those at Ibaraki were obtained under the project P07015 of the NEDO. The authors are grateful to all of the organizations above. Comments from the editor and anonymous reviewers are also acknowledged. We would also like to thank Editage for English language editing.

Conflicts of Interest:
The authors declare no conflict of interest.