Assimilating high-resolution winds from a Doppler lidar using an ensemble Kalman filter with lateral boundary adjustment

Monitoring severe weather, including wind shear and clear air turbulence, is important for aviation safety. To provide accurate information for nowcasts and very short-range forecasts up to an hour, a rapid-update prediction system has been developed, with a particular focus on lateral boundary adjustment (LBA) using the local ensemble transform Kalman filter (LETKF). Due to the small forecast domain, limited-area forecasts are dominated by the lateral boundary conditions from coarse-resolution global forecasts. To effectively extend the forecast lead time for the small domain, a new LBA scheme using the LETKF has been developed and assessed with three sea-breeze front cases. Observing system simulation experiments for high-resolution winds from a simulated Doppler lidar were performed with the Japan Meteorological Agency Nonhydrostatic Mesoscale Model at a horizontal resolution of 400 m and 15-minute update cycle. The results indicate that the LBA improved the forecast significantly. In particular, the 1-hour wind-speed forecast with the LBA is as accurate as the 15-minute forecast without the LBA. The assimilation of Doppler lidar high-resolution wind data with the LBA is a promising approach for very short-range forecasts up to an hour with a small domain, such as for aviation weather.


Introduction
For aviation safety, reliable weather forecasts with a high spatio-temporal resolution are crucial (e.g. Frech et al., 2007). It is particularly important to monitor real-time weather conditions, such as the clear air turbulence, wind shear and aircraft wake vortex, during takeoff and landing (Proctor et al., 2004). For real-time monitoring of aviation weather, Doppler lidars are used at the Haneda and Narita international airports in Japan (Yamamoto, 2009), Hong Kong International Airport (Shun and Chan, 2008) and other international airports. Doppler lidar is a powerful tool for observing three-dimensional wind fields at a high spatiotemporal resolution in clear air conditions and is suitable for aviation purposes (Ishii et al., 2007;Iwai et al., 2008Iwai et al., , 2011. Assimilating the high-resolution observations into a numerical model enables short-range numerical forecast guidance for 1Á12 hours, as, for example, in Rapid Refresh (RAP, http://rapidrefresh.noaa.gov/), an operational regional numerical weather prediction (NWP) system at the National Centers for Environmental Prediction.
Dynamical downscaling using a limited-area model (LAM) is an efficient and common approach to highresolution NWP. Generally, this approach requires an initial condition and lateral boundary conditions (LBCs). The LBCs are typically derived from a coarse-resolution global NWP model; the accuracy of the downscaling forecast depends on the initial and boundary data (Laprise et al. 2000). The domain size and the forecast length determine the degree to which the LBCs affect the model solution in the domain interior (Warner et al., 1997). When the forecast *Corresponding author. email: msawada@mri-jma.go.jp Tellus A 2015. # 2015 M. Sawada et al. This is an Open Access article distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), allowing third parties to copy and redistribute the material in any medium or format and to remix, transform, and build upon the material for any purpose, even commercially, provided the original work is properly cited and states its license. domain is small, the LAM forecast is rapidly dominated by the LBCs (Errico and Baumhefner, 1987;Gustafsson et al., 1998). The slight change in the coarse outer-domain forecasts may cause critical differences in the inner-domain forecasts. Moreover, the downscaling forecast using the coarser resolution model may have a difficulty in representing the small-scale phenomena because the LBCs do not include the mesoscale evolution and the domain size is not sufficiently large to self-generate within a limited domain. In this case, modifying the LBCs is a way to extend the predictability of the small forecast domain.
Data assimilation has been widely used to improve the initial conditions and the LBCs at the initial time. Several studies have included LBCs as control variables in fourdimensional variational data assimilation (4DVAR) to optimise the initial conditions and, simultaneously, the LBCs (Zou and Kuo, 1996;Kawabata et al., 2007;Gustafsson, 2012). Gustafsson (2012) showed that including LBCs as control variables in 4DVAR improved the forecasts, especially in the case of strong inflow from the lateral boundary. However, these previous studies did not consider changing the forecasted LBCs beyond the assimilation window. By modifying the LBCs for the entire forecast period, we may extend the predictability of the downscaled prediction, particularly when the forecast domain is very small. Hereafter, the scheme to modify the LBCs in the forecast period is referred to as lateral boundary adjustment (LBA). Fig. 1 illustrates how the LBA works. When the LBCs gradually deviate from the natural state (e.g. lagging of the front passage and phase error of the low/high pressure systems), the forecast errors grow because the LAM forecast approaches the inaccurate LBCs (Fig. 1a). If we apply the LBA approach and optimise the LBC in the forecast, the forecast errors are closer to the modified LBCs (Fig. 1b). If the modified LBCs have smaller forecast errors than those in the original LBCs, the analysisÁforecast cycle with the LBA is expected to suppress the error growth. This idea is similar to the 'quasi-outer-loop' algorithm (Yang et al., 2012), which adjusts the ensemble mean and recentres the ensemble perturbations around the adjusted more accurate mean state. We extend the recentreing of ensemble perturbations for the entire forecast period without the corresponding (future) observations (Section 2). This approach will not be a good choice if the background error covariance suddenly changes.
The purpose of this study is to develop a rapidly updated NWP system combined with LBA for a very short-range forecast ( Â1 hour) and to assess its impact on forecasts in an observing system simulation experiment (OSSE) with Doppler lidar. We use the local ensemble transform Kalman filter (LETKF, Miyoshi and Aranami, 2006;Hunt et al., 2007;Miyoshi and Kunii, 2012) to develop the LBA scheme. The paper is organised as follows. Section 2 describes the newly developed LBA scheme, followed by the experimental design in Section 3. Section 4 presents the results. Finally, Section 5 provides concluding remarks and discussions.

Lateral boundary adjustment
The LBA scheme corrects the ensemble mean of the LBCs for the entire forecast period using the LETKF. First, we introduce a basic equation of LETKF, which is used for providing an initial condition of forecastÁanalysis cycle. Then a specific equation for LBA proposed in this study is introduced.
In the LETKF, the analysis equation for the ensemble mean state is (e.g. Hunt et al., 2007;Miyoshi, 2010) as follows: where x is the model state variable and dX is the ensemble perturbation matrix (each column corresponds to each member). w i is the weight matrix or transformation matrix. e P a is the analysis error covariance matrix in the space spanned by analysis ensemble perturbations and R is the observation error covariance matrix. H is the observation operator, H is its tangent linear and y8 is the observation. The overbar denotes the ensemble mean. Superscripts a and f represent analysis and forecast, respectively, and the subscript i is the time index. The second term on the righthand side of eq. (1) is known as the analysis increment. Practically, the analysis ensembles for the initial condition are constructed by adding analysis ensemble perturbations dX f i w i to the forecast ensemble mean " x f i (Bishop et al., 2001;Hunt et al., 2007;Miyoshi, 2010). The x-axis represents the forecast time, and the y-axis is the forecast error. The stars and filled circles indicate analysis (initial conditions) and forecast (first guess), respectively. Solid arrows represent the forecast ensemble mean trajectory, and the dotted lines correspond to the LBC ensemble mean trajectory. Grey ellipses denote the ensemble spread for each analysis and forecast state.
Following this, the LBA equation is derived based on eq. (1) by assuming that the transformation matrix w i remains constant for the entire forecast period. This assumption can be made, because the forecast time is very short (several hours at most) and the transformation matrix cannot be obtained satisfactorily for the entire forecast period without the corresponding (future) observations y n (n!i). Additionally the computational cost can be saved by reducing the matrix operation. Based on this assumption, we obtain the analysis ensemble mean for the entire forecast period: " x a n ¼ " Here, the subscript n is the index for the forecast time step from i'1 to N (end of forecast). Moreover, eq. (3) is identical to eq. (1) at n 0i (analysis time). The ensemble mean of the LBCs for the forecast period (" x f n ) is modified using eq. (3). The lateral boundary perturbations (LBPs) are derived from the ensemble downscaling forecasts of coarse-resolution global NWP. Thus the LBCs (x a n ) are obtained: x a n ¼ " x a n þ dX f n .  . Previous studies have applied the LETKF to the JMA-NHM with promising results (Miyoshi and Aranami, 2006;Seko et al., 2011;Saito et al., 2012). In the present study, the nature run was created by downscaling the JMA operational mesoscale analysis at a 10-km resolution to a 0.4-km resolution, at 3-hourly intervals. A horizontal resolution of 0.4 km was necessary to resolve sea-breeze front sharply in this case. The simulation was run for 168 hours, starting at 09 LT on June 16, 2007. The nature run data consisted of meteorological fields every 15 minutes. The computational domain contained 72 )72 horizontal grid points and 40 vertical layers. For the lateral boundaries, a sponge layer with eight grids on each boundary was employed. The eastern and western boundaries covered the ocean and the land regions, respectively (Fig. 2b).

Experimental design of OSSE
To create the LBPs for the analysisÁforecast cycles with LETKF, ensemble downscaling forecast experiments were performed using the JMA operational 1-week ensemble prediction system (EPS) with a horizontal resolution of 1.258 and 6-hourly intervals. The unperturbed and the first 10 positive perturbed members were used for the initial and LBCs, downscaled to horizontal resolutions of 12 km and 2 km, with one-way nesting. By subtracting the first 10 positively perturbed members from the unperturbed member, 10 negatively perturbed members were obtained; the entire forecast ensemble, with a 0.4-km mesh, included 21 members. The forecast length and output interval, in the experiment with a 2-km mesh, were the same as those in the TRUE.
In the OSSE, the lidar observations were simulated from the TRUE fields through interpolation at the observation locations, assuming a lidar specification similar to that used at the Sendai airport (Iwai et al., 2008). The simulated lidar performed plane position indicator (PPI) scans at 20 distinct elevation angles (i.e. 18, 28, 48, 68, 88, 108, 128, 148,178, 208, 258, 308, 358, 408, 458, 508, 558, 608, 708 and 808) with an azimuth resolution of 188 and a radial resolution of 180 m from 360 m to 10 080 m (55 points) in the radial direction (a horizontal range of observations is shown in Fig. 2b). The radial resolution used by Iwai et al. (2008) was 90 m but the number of elevation angles increased by decreasing radial resolution to produce spatially uniform observation data. The observed height is limited to below an altitude of 2 km. The total number of lidar observations was 12 860. A simulated lidar observation was generated every 15 minutes, for 168 hours. Gaussian random errors, with standard deviations of 1 m s (1 , were added to the radial velocity. The simulated lidar data were generated with the same grid spacing as the forecast experiment such that the representativeness error, unlike the real situation, was not included. The LETKF data assimilation scheme  was adjusted to the JMA-NHM through several modifications, including the implementation of vertical hybrid coordinates and the introduction of LBPs in the forecastÁanalysis cycles. The state vector incorporated air density, u, v, w, temperature, pressure, water vapour mixing ratio, surface pressure, ground temperature, underground temperature and surface precipitation. The assimilation cycle for the initial conditions was conducted with a horizontal resolution of 0.4 km and an assimilation window of 15 minutes. The 15-minute forecast with a horizontal resolution of 0.4 km was used as a first guess. Simulated observation data were assimilated during the analysis. The ensemble size was set to 21, and a multiplicative inflation factor of 10% was used. The localisation parameters of LETKF were set to 2 km in the horizontally and 0.5 km in the vertically. The LBA was conducted using only the ensemble perturbations with the 2-km mesh (LBPs), because the ensemble perturbations with the 0.4-km mesh for the entire forecast period could not be obtained before the analysisÁ forecast cycle was completed. The horizontal localisation parameter of the LBA was set to 50 km, modifying the LBCs, even though the lidar data were far from the lateral boundary; the horizontal localisation parameter of 2 km is too short to modify the LBCs. The localisation radius of influence is 2 ffiffiffiffiffiffiffiffiffiffi 10=3 p times the localisation parameter (Gaspari and Cohn, 1999;Hamill et al., 2001). The vertical localisation parameter was fixed at 0.5 km. A superobservation method was applied in order to reduce the computational time for the assimilation cycle and the random instrumental error. Super-observations are generally applied to higher-resolution observation data (180 m, in the present study) than the model grid spacing (2 km). The simulated radial velocity was interpolated to the 2-km horizontal grid using the Cressman scheme (Seko et al., 2004). Table 2 presents an outline of the three experiments (and additional sensitivity experiments) performed in order to examine the effect of LBA. The first experiment was a pure ensemble forecast from the JMA 1-week EPS without assimilation (NOAS), aimed at obtaining the impact through comparison with the two subsequent experiments. The second experiment was an analysisÁforecast cycle using the simulated lidar data to correct only the initial conditions (AS). The third experiment was an analysisÁforecast cycle with LBA, aimed at correcting both the initial conditions and the LBCs (ASLB). This final experiment required twice the assimilation processing for each analysis cycle, one for the initial conditions and another for the LBCs (LBA). The model setting (resolution and physical parameterisations) in the three experiments were the same as those in TRUE. Three study cases were selected (Table 3)

A case study with and without the LBA
As noted above, previous studies have suggested that the forecast error in the LAM is restricted by the LBCs. To overcome the limitation of the LBCs, this study attempted to optimise the LBC by applying the LETKF to forecasted variables of the LBC in the OSSE framework. First, one case study was described in detail. Then, a composite of the three cases was analysed. Fig. 3 shows simulated wind fields at a 125-m height for the case 14:00Á18:00 LT on June 16, 2007, for TRUE, NOAS, AS and ASLB. Vectors and shaded regions represent the horizontal wind and eastÁ west wind speed at a 125-m height, respectively. Ensemble means in NOAS and analysis means in AS and ASLB are displayed. In this case, a low-pressure system moved from the southwest to the east coast of Japan, which caused the prevailing northwesterly winds (Fig. 2a). The prevailing northwesterlies gradually weakened in the afternoon; the sea-breeze front passed the Sendai airport at 17:00 LT in the surface weather observation. In the TRUE wind field, the sea-breeze front appeared at the eastern boundary at 15:30 LT, intruded over land at 17:00 LT and reached the Sendai airport at 17:15 LT (not shown). In the NOAS experiment, the sea-breeze had already penetrated the inland area at 15:00 LT because the synoptic westerlies were weaker than those of TRUE due to the weaker low-pressure system. In the AS experiment, the synoptic northwesterlies, which were not formed in NOAS, were formed successfully, and the sea-breeze was not observed at 15:00 LT. At 16:00 LT, the sea-breeze penetrated the boundary, arising from the excessive easterlies of the LBCs, as shown for NOAS in Fig. 3. The position of the sea-breeze front was obscure at 17:00 LT. The passage of the sea-breeze front at the Sendai airport was successfully captured. However, the front overpenetrated near the northern and southern lateral boundaries due to excessive easterlies in the LBC at 18:00 LT. The ASLB experiment also captured the right timing and position of the sea-breeze intrusion similar to AS. At 15:00 LT, the excessive easterlies at the eastern lateral boundary became weaker than in AS. At 16:00 LT, the excessive  easterlies along the northeastern lateral boundary were reduced through the assimilation cycle. However, the seabreeze near the southeastern lateral boundary was not captured because it was located outside the observation range of the lidar. One hour later, the sea-breeze front was within the observation range; the front structure was then well captured. The sea-breeze front was clearly formed at 18:00 LT and, unlike AS, the excessive penetration of the sea-breeze near the northern and southern lateral boundaries did not occur. The easterlies at the eastern lateral boundary were weaker than in the TRUE fields, possibly because the LBA uses a long horizontal localisation parameter so that both the westerlies and easterlies near the centre were assimilated. The vertical structures of the simulated sea-breezes are examined from eastÁwest cross sections at y 00 km for 1-hour intervals in Fig. 4. Thick lines represent contour lines of eastÁwest wind speeds of 0 m s (1 , indicating the boundary between the sea-breeze and synoptic flows for the TRUE fields. The sea-breeze appeared at the eastern lateral boundary approximately 16:00 LT, and its front, with a thickness of less than 400 m, travelled inland at a speed of 2Á3 m s (1 , passing the Sendai airport (x00 km) at 18:00 LT in the TRUE data. The sea-breeze (easterly) with a thickness of 600Á800 m penetrated inland in the NOAS experiment at 15:00 LT and remained until 18:00 LT. The NOAS experiment failed to reproduce the sea-breeze thickness and the timing of its intrusion. The sea-breeze position and thickness were well captured between 17:00 and 18:00 LT in both AS and ASLB because of the assimilation of lidar data. These cases reproduced the synoptic westerlies, which were not formed in NOAS. Compared with AS, a stronger seabreeze near the surface was found in ASLB, similar to the TRUE fields.
To investigate the extent to which the LBA reduces the forecast and analysis errors, the horizontal map of root- mean-square errors (RMSEs) for both wind speed and direction are shown for the NOAS, AS and ASLB experiments in Fig. 5. The RMSE was calculated as follows: where U and V are ensemble mean eastÁwest and northÁ south wind speeds, respectively, in NOAS, AS or ASLB and Ut and Vt are TRUE fields. z and t are indexes for the vertical layer and time step, respectively. Here, NZ is 20, which corresponds to a 1-km height; NT is 16 for the period 14:15Á18:00 LT, with an interval of 15 minutes (16 assimilation cycles). The NOAS experiment had the largest RMSE over the entire domain, especially along the southern lateral boundary, where the excessive progression of the sea-breeze front was observed. In both first guess and analysis fields of AS, the RMSEs became 2Á3 m s (1 smaller inside the lidar observation range than near the lateral boundaries. The RMSEs of the first guess near the inflow lateral boundaries were close to those in NOAS because both experiments employed the same LBCs. The RMSE near the lateral boundaries was remarkably smaller in the first guess and analysis of ASLB than in those of AS. In the southwest quadrant, the RMSEs of both the first guess and analysis of ASLB decreased by more than half of the AS experiment. Another important difference is that the RMSE in the domain interior decreased in the first guess of ASLB. For the first guess, the LBA had a positive impact on the forecast not only near the lateral boundaries but also in the inner area.
To evaluate the temporal changes of forecast error in the three experiments for 14:00Á18:00 LT on June 16, 2007, domain-averaged RMSEs of wind speed below 1 km are shown in Fig. 6. The lines represent both the first guess and analysis RMSEs. In the NOAS experiment, the RMSE reached approximately 8.5 m s (1 at 15:15 LT, gradually decreasing thereafter to 6 m s (1 . The RMSE for the analysis in AS was smaller than that in the NOAS. However, the error rapidly grew from the lateral boundaries (Fig. 5b), approaching the NOAS error. This pattern is consistent with the findings of previous studies (c.f. Errico et al., 1993). The assimilation reduced the temporally averaged RMSEs of the first guess and analysis for the period 14:15Á18:00 LT in AS by 6.19 m s (1 and 4.03 m s (1 , respectively. When the LBA was applied, the RMSEs for the analysis and first guess of ASLB became 2.39 m s (1 (39%) and 0.91 m s (1 (23%) smaller, respectively, than those in the AS. In ASLB, the RMSE for the first guess increased by less than 1 m s (1 from the analysis for the entire period. It is evident that the RMSE reduction for the first guess arises from the LBA, indicating that this has a positive effect on the extended forecast.

Extended forecast experiment
The effectiveness of the LBA in extending the forecast lead time was quantitatively investigated with additional experiments. The forecasting processes were extended from 15 minutes to 1 hour at each analysis time (15-minute intervals) for the three cases with and without the LBA. Hereafter, these experiments are called F1h_AS and F1h_ASLB. The model setting (resolution and physical parameterisations) in the extended forecast experiments were the same as  Domain-averaged RMSE of wind speed from surface to 1-km height for the analysisÁforecast cycle in the NOAS (black), AS (blue) and ASLB (red). The x-axis denotes the forecast date and the y-axis is the averaged RMSE. Circles and triangles represent analysis and first guess. those in the TRUE. Figure 7 shows the time-series of the domain-averaged wind-speed RMSEs for each forecast at different initial times (14:00 LT, 06:00 LT and 05:00 LT on June 16, 17 and 19, 2007, respectively). The circles represent the initial time, and blue and red lines correspond to F1h_AS and F1h_ASLB, respectively. The values are normalised by the NOAS RMSE at each forecast time. These panels demonstrate that the error growth rate depends significantly on the specific cases, especially in the F1h_AS scenario. In case-1 of the F1h_AS simulation, the normalised RMSE was 0.6 at the initial stage (analysis), exceeding unity (the NOAS RMSE) in the 30-minute forecast. In case 2 of F1h_AS, the RMSE rapidly increased and exceeded the NOAS RMSE after 15 minutes. This result may be caused by a dynamical inconsistency between the assimilated interior flows and the LBC. The RMSE grew relatively slower in case 3 of F1h_AS, approaching the NOAS RMSE. A particularly interesting result is that the RMSE became significantly smaller in F1h_ASLB than in F1h_AS for all initial times and cases. In almost all cases, the RMSE at 1 hour was smaller than the NOAS RMSE except for a few initial points in case-1. This difference is because the wind direction changed from a synoptic westerly to a seabreeze (easterly) at the southeastern boundary, outside of the lidar observation range, during the period 16:15Á17:00 LT (Fig. 3). Even though the lidar data without the information on the wind direction change were assimilated, the LBA could not reproduce suitable LBCs, resulting in a large RMSE in that region. Figure 8 illustrates the normalised error growth by the NOAS RMSE. The two lines with one standard deviation bar are the averaged RMSE for the 36 forecast experiments, which consist of three cases, with each case with 12 1-hour forecast experiments from different initial times. In F1h_AS, the normalised RMSE was less than 70% of the NOAS RMSE at the initial time and increased to unity after 15 minutes, which means the RMSE is comparable with the NOAS RMSE. The normalised RMSE in F1h_ AS was approximately constant for the 30-to 60-minute forecasts because the LBC controls the domain interior during this period. The normalised RMSE in F1h_ASLB at the initial time was less than 40% of the NOAS RMSE and grew quite slowly with time, reaching a normalised RMSE of 0.4890.23. The 1-hour forecast error in F1h_ASLB was smaller than at the initial time in F1h_AS. It is evident that the LBA contributes to extending the forecast lead time in this study.
The horizontal distributions of normalised RMSE in F1h_AS and F1h_ASLB are presented in Fig. 9. Shaded regions represent the normalised RMSE, and contours represent the RMSE for the 15-minute (Fig. 9a, c) and 60minute (Fig. 9b, d) forecast lengths. The 15-minute forecast of F1h_AS shows that the normalised RMSE became less than 0.8 within half of the lidar observation range and was comparable with that of the NOAS RMSE along the lateral boundaries. For the 1-hour forecast, there was no improvement in nearly the entire domain, indicating that the effects of assimilated initial conditions quickly disappear. The RMSE of the 15-minute forecast in F1h_ASLB was less than 0.4 of the NOAS RMSE near the lateral boundaries and within the lidar observation range, representing a significant improvement compared with F1h_AS. The RMSE distribution for the 1-hour forecast was similar to that of the 15-minute forecast, revealing that the forecast error did not grow quickly through the lateral boundaries in F1h_ASLB. In this study, the LBA was effective for the 1-hour forecast length and suppressed the error growth throughout the entire forecast domain.

Additional experiments on the observation range and number of instruments
The impacts of the number of instruments and observation range on the forecast error were investigated by performing several sensitivity experiments. Two and four lidar observation sites were deployed (Fig. 2b), referred to as AS2p and AS4p, respectively. The model setting (resolution and physical parameterisations) in the sensitivity experiments were the same as those in the TRUE. Domain-averaged RMSEs are presented for the observation range of approximately 10 km and 5 km from the lidar in Fig. 10. RMSEs are averaged for the three cases. As the number of lidars increased, the RMSEs of analysis decreased in all experiment cases. In the analysis of AS, the RMSE reduction in response to the increase in lidars was larger than in ASLB. However, the RMSEs for the first guess of AS did not decrease as suggested by the analysis when the number of lidars increased. Despite the AS4p covering the entire extent of the lateral boundaries, the forecast error grew quickly and became comparable to the NOAS RMSE after 15 minutes, which suggests that the initial condition correction was not sufficient to reduce the forecast error. The RMSE for the first guess was smaller in ASLB1p than in AS4p, implying that the LBA effectively uses past observation information because the error reduction was larger in ASLB than in AS even though there were fewer observations in ASLB.
To investigate the impacts of the observation range radius, we performed sensitivity experiments with half the observation range of AS and ASLB, referred to as AS_h and ASLB_h, respectively. The RMSE for the analysis was 0.5Á1.5 m s (1 larger in AS_h and ASLB_h than in AS and ASLB because the observation coverage was small in the AS_h and ASLB_h experiments. The RMSE for the first guess of AS_h was similar to that of AS, which was caused by saturation of the forecast error produced by the LBC. In contrast, the RMSE for the first guess of ASLB_h became 1.2Á1.6 m s (1 larger than in ASLB. The effects of the number of instruments in AS_h and ASLB_h were nearly qualitatively the same as in AS and ASLB. In these cases, the results imply that the LBA, rather than increasing the number of lidars and observation range, is required to extend the forecast lead time because the RMSE was dominated by the lateral boundaries rather than the initial conditions.

Conclusion and discussion
We have proposed an LBA scheme using the LETKF and evaluated its effectiveness at extending the forecast lead time for three cases in which a sea-breeze was formed by performing OSSE using Doppler lidar. The RMSEs for the first guess and analysis of AS decreased by 1.0 and 3.2 m s (1 (14 and 44%), respectively, compared with the RMSE of NOAS. By applying the LBA, the RMSEs were 3.4 and 4.1 m s (1 (47 and 57%) smaller than that in NOAS and were 2.4 and 0.9 m s (1 (39 and 23%) smaller than in AS, respectively. Compared with the AS, the RMSE was clearly reduced in the first guess because the LBA reduced the RMSE both in the domain interior and near the lateral boundaries. Forecast experiments for 1 hour were performed at each analysis time (12 initial times for one case) for the three cases to investigate the extent to which the LBA extends the forecast lead time. The domain-averaged RMSE of ASLB for the 1-hour forecast for all cases decreased by 3.5 m s (1 (50%) of the AS RMSE, which is less than or comparable to the RMSE of the 15-minute forecast for AS.
In previous studies and most operational centres, data assimilation is primarily used for the improvement of the initial conditions. This adjustment is not always beneficial in a LAM forecast, especially with a narrow domain, because the observation information is quickly lost from the lateral boundaries (Warner et al., 1997). We accomplished an extended forecast for a narrow domain under the OSSE framework based on the assumption of linearity of the analysis increment during a forecast period ( Â1 hour in this study). To the best of the authors' knowledge, this method is the first LETKF attempt to mitigate the error growth from LBCs, and the findings indicate that the proposed LBA is a promising approach for prediction for narrow domains. The combination of prediction, Doppler lidar and data assimilation with LBA will be useful for providing not only aviation hazard guidance but also severe weather (tornado and torrential rain) and situational awareness information.
In the operational sense, a caution may be needed with regard to the computational cost. The ASLB experiment requires double assimilation for each analysis, increasing the computational demand compared with AS (one assimilation for each analysis). However, the assimilation uses only less than 5% of the total CPU time of 20-member ensemble runs in this configuration. Therefore, the additional assimilation computation is not essential.
Another concern is that the remarkable improvement after applying the LBA may be related to the non-normalised LBP. Here, the LBP was produced from the dynamical downscaling from the JMA global ensemble forecast, in which the spread is too large because the forecast time is different from the small-domain forecast. The large LBP causes the ensemble spread to increase, resulting in an excessive correction by the observations. We may need to consider generating an appropriate LBP (c.f. Torn et al., 2006;Saito et al., 2012). This will remain an important future issue.
The assumption of constant background error covariance matrix during the whole forecast period is also shortcoming because it may lose the flow-dependent characteristics, leading to filter divergence. Throughout this work, the assumption of constant dX f i is valid at least for the selected dry cases. Additionally, the lidar observation is only assimilated in the limited-area domain, possibly leading to the significant improvement. It is necessary to evaluate a combinational effect between the lidar observation and other instruments, such as surface observation, which could reduce the impact of assimilation with lidar observation. The LBA is effective only in some rare cases with a large forecast displacement error such as in this study.
Similarly to the LBA approach, parameter estimation techniques based on data assimilation could be applied to produce the LBCs at appropriate forecast time. Parameter estimation generally provides optimal parameters (LBCs in the present case) at the analysis time using the present observations. If the idea of the four-dimensional ensemble Kalman filter (Hunt et al., 2004) is applied to parameter estimation, parameter values at the later forecast times can be estimated with the present observations; this is analogous to the LBA approach. Similar assumptions are made by the running-in-place method (Kalnay and Yang, 2010), and quasi-outer-loop method (Yang et al., 2012). These methods improve state estimates using the 'future' observations within the assimilation window. The parameter estimation and the LBA can be considered as a similar approach.
Although the preliminary results from this case study provide a promising approach to exploring how the LBCs improve LAM predictions, future work is required to practically implement the proposed LBA. Because only three study cases with very specific condition (only assimilated lidar data without other observations and large forecast displacement error) and the OSSE framework were used in this study, a more comprehensive evaluation of the LBA effects on the real forecast are needed to show statistical confidence. Moreover, nature includes moist processes, which have stronger nonlinearity compared with dry processes, such as sea-breezes. This method should be verified for cases including precipitation to understand the case dependency.