the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Realtime estimation of airflow vector based on lidar observations for preview control
Ryota Kikuchi
Takashi Misaka
Shigeru Obayashi
Hamaki Inokuchi
As part of control techniques, gustalleviation systems using airborne Doppler lidar technology are expected to enhance aviation safety by significantly reducing the risk of turbulencerelated accidents. Accurate measurement and estimation of the vertical wind velocity are very important in the successful implementation of such systems. An estimation algorithm for the airflow vector based on data from airborne lidars is proposed and investigated for preview control to prevent turbulenceinduced aircraft accidents in flight. An existing technique – simple vector conversion – assumes that the wind field between the lidars is homogeneous, but this assumption fails when turbulence occurs due to a large windvelocity fluctuation. The proposed algorithm stores the lineofsight (LOS) wind data at every moment and uses recent and past LOS wind data to estimate the airflow vector and to extrapolate the wind field between the airborne twin lidars without the assumption of homogeneity. Two numerical experiments – using the ideal vortex model and numerical weather prediction, respectively – were conducted to evaluate the estimation performance of the proposed method. The proposed method has much better performance than simple vector conversion in both experiments, and it can estimate accurate twodimensional windfield distributions, unlike simple vector conversion. The estimation performance and the computational cost of the proposed method can satisfy the performance demand for preview control.
Atmospheric turbulence poses a potential risk to aircraft operation. Statistics reported by Boeing (2018) show that 322 nonfatal and 51 fatal accidents occurred worldwide in commercial jet flights from 2009 through 2018. Of the fatal accidents, the largest proportion (25.5 %) were due to loss of control inflight (LOCI). The International Air Transportation Association (2016) shows that LOCI frequently occurs when the aircraft speed is well below the stall speed; in conjunction with weather conditions, low speed is the most common factor in LOCI accidents. A total of 42 % of LOCI accidents occurred under degraded meteorological conditions affecting aircraft speed, in particular strong wind shear and atmospheric turbulence.
For both fatal and nonfatal aircraft accidents, the impact of atmospheric turbulence can be significant. The Japan Transport Safety Board (2020) has stated that accidents caused by turbulence accounted for 48 % of nonfatal aircraft accidents in Japan involving commercial airplanes from 2003 to 2012. An increase in the rate of accidents related to turbulence was reported by the Federal Aviation Administration (2006), Kim and Chun (2011), and Williams (2017). Accidents caused by convective systems such as cumulonimbus clouds have decreased due to advances in airborne radar (Airbus, 2020; Sermi et al., 2015). However, noncloud atmospheric turbulence, called clearair turbulence (CAT), cannot be detected by radar, as reported by Soreide et al. (2000), Barny (2012), and Inokuchi et al. (2009). Airborne CATobservation systems to minimize risks of turbulencerelated accidents are essential for aviation safety.
Numerical weather prediction (NWP), which is an essential tool for aircraft operation, can forecast weather conditions for days and even weeks in advance and output broaderarea weather information than can radar or lidar. However, NWP cannot explicitly resolve disturbances as small as most turbulence, leading to a very large predictive uncertainty (Sharman et al., 2006, Kim et al., 2011). Therefore, some researchers have developed an alternative approach that predicts turbulence potential by calculating turbulence indicators from NWP results; for example, Sharman et al. (2006) have developed an approach called graphical turbulence guidance (GTG) that combines such indicators. The turbulence potential can also be used to determine operational flight routes (Kim et al. 2015), but it has a large spatiotemporal gap on the scale of aircraft motion because it is based on NWP results such as the mesoscale model. It thus provides insufficient information to implement turbulence avoidance on aircraft in flight.
Recently, airborne Doppler lidar has been developed by Soreide et al. (2000), Barny (2012), Inokuchi et al. (2009), Machida (2017), and Inokuchi and Akiyama (2019). Emitted laser light is scattered by fine aerosol particles in the atmosphere; the backscattered light is condensed by telescopes and received by an optical transceiver. Since the wavelength of the received light varies according to the velocity of the aerosol particles due to the Doppler effect, wind speed can be calculated by comparing this wavelength with that of the received light (Inokuchi and Akiyama, 2019). However, when rain is too heavy, the backscattering signal is weakened due to strong attenuation by raindrops and a decrease in aerosols (Wei et al., 2019), making it difficult to measure the wind velocity at a distance. The Japan Aerospace Exploration Agency (JAXA) is researching and developing a coherent Doppler lidar capable of remotely detecting air turbulence in clearair conditions and has conducted a flight demonstration of a lidar system that can provide turbulence information to pilots (Inokuchi et al., 2009; Machida, 2017; Inokuchi and Akiyama, 2019). Inokuchi (2012) have shown observationally that airborne Doppler lidar can detect CAT in front of an aircraft in flight at altitudes of 3200 m; the lidar information can be detected 30 s before the turbulence affects the aircraft. The aircraft's flight speed in the test was 320 kt (160 m s^{−1}), so it detected CAT from a distance of about 4.8 km.
Based on advance airflow information, flight demonstrations have been carried out with the aim of providing pilots with the information they need to make decisions: whether to change course to avoid wind shear, and whether to turn on seatbeltsign lighting during cruise and altitude changes (Inokuchi and Akiyama, 2019). Although lidar systems are useful for providing onboard wind information to pilots, avoiding turbulence at high altitudes is difficult as the range of detection that facilitates pilots to be warned is short (Hamada, 2019). Gathering such information involves emitting a laser beam and receiving the scattered light from aerosol particles that are much smaller than precipitation droplets in the air. Therefore, when the number of aerosol particles that emit scattered light is small, it is difficult to measure wind information at a distance. Furthermore, as altitude increases, the aerosol density decreases, and the observation range tends to decrease accordingly. The maximum observation range and aerosol density measured at each altitude are shown in Inokuchi and Akiyama (2019).
Advance knowledge of turbulent atmospheric conditions would improve the performance of automatic aircraftvibration reduction systems. Automatic control to alleviate aircraft vibration is called gust alleviation and has been studied since the 1970s, mostly with only the help of feedback sensors such as inertial measurement units (Regan and Jutte, 2012). Recently, methods of reducing the vibrations due to turbulence with the help of preview controlling based on airborne lidar observation have been reported by Schmitt et al. (2007), Fezans et al. (2019), and Hamada (2019). The aim of the Aircraft Wing with Advanced Technology Operation (AWIATOR) project is the development of new directlift control devices and a lidar system for turbulence measurement (Schmitt et al., 2007). Another project – “Demonstration of Lidarbased CAT detection” (DELICAT) (Barny, 2012) – developed airborne ultraviolet lidar for gust and turbulence measurements. The test flights were carried out using an Airbus A340 aircraft equipped with ultraviolet lidar. In both the AWIATOR and the DELICAT experiments, the measurement range was short, because the lidar was developed for controlling the aircraft automatically.
In order to implement an airborne Doppler lidar gustalleviation system successfully, it is very important to measure the vertical wind velocity accurately. Both horizontal and vertical winds affect aircraft motion, but the effect of changing the vertical wind velocity is greater. This is because the effect of modifying the angle of attack is relatively larger than the effect of changing the horizontal wind velocity, which affects only the airspeed (Fezans et al., 2019). However, a fixed single Doppler lidar system can only detect the lineofsight (LOS) wind, providing a onedimensional piece of information; the vertical wind velocity in front of the aircraft cannot be measured by such a system (Hamada, 2019). It is necessary to perform the lidar measurements in two directions, upward and downward, to obtain the vertical wind velocity (Neininger, 2017). Figure 1 shows a representation of this concept. The vertical windvelocity vector is generated from the differences between the upward and downward LOS winds by using simple vector conversion. Unfortunately, this method is incapable of estimating the vertical wind velocity with high accuracy to control the aircraft automatically because the technique assumes homogeneity between the upward and downward lidars (Fezans et al., 2019). In this study, a fully turbulent field with atmospheric turbulence and gusts is considered; under these conditions, it is difficult to estimate the vertical wind velocity with high accuracy using simple vector conversion. In particular, the estimation accuracy of the vertical wind velocity rapidly worsens when the estimation position is located farther ahead from the aircraft.
In addition, actual lidar observations involve errors, noise, and loss of data, with negative effects on aircraft control, as reported by Misaka et al. (2015); these problems are worse at higher altitudes, where the aerosol density is smaller than it is at lower ones. Misaka et al. (2015) proposed a filtering algorithm based on a simple Kalman filter to remove windvelocity errors from lidar measurements. For preview control, it is essential to deal with the lidar errors, noise, and loss of data more carefully. An accurate airflowvectorestimation method and an efficient realtime filtering algorithm are required.
In this study, an estimation method and an airflowvectorfiltering algorithm are proposed for preview control to prevent turbulenceinduced aircraft accidents. The method works for both horizontally and vertically directed winds, and uses both upward and downward lidars. (In this study, “horizontal wind” means any headwind/tailwind component that does not include the crosswind component.) The lidar system in this paper is that also used by JAXA in its ongoing “lidarbased gust alleviation control” research project. The proposed algorithm stores the LOS wind data continually and uses recent and past LOS wind data to estimate the airflow vector and the wind field between lidars, whereas simple vector conversion utilizes only recent LOS wind data. The airflow vector is calculated by using wind data extrapolated from the horizontal and vertical wind components; the estimation accuracy of the airflow vector in front of the aircraft is improved by using such extrapolated wind data because the region between the lidars represents a nonhomogeneous one. A polynomial expression is used to extrapolate the wind field. In addition, the proposed method can estimate the twodimensional distribution of the wind field between the lidars, which simple vector conversion cannot.
Two test configurations – an ideal vortex flow field and a weather field – are calculated by an NWP system and utilized to evaluate the performance of the airflow vector. These experiments generate a large number of pseudolidar measurements along flight routes from the reference wind field for evaluation of the estimated performance. Comparing the prediction results with the reference wind field can confirm all the windfield values.
2.1 Airborne lidar specifications
The airborne lidar observation system currently under development by JAXA for preview control to prevent turbulenceinduced aircraft accidents is shown in this section. This system has airborne lidars that are aiming upwards and downwards; the angle between them is 20^{∘}, that is, 10^{∘} between the horizontal line and each lidar. The lidar sensor is shown in Fig. 2; its specifications are given in Table 1 (Inokuchi and Akiyama, 2019). Laser pulses generated by an optical transceiver are amplified by optical amplifiers (Sakimura et al., 2013) incorporated into an optical antenna and radiated into the atmosphere from optical telescopes. The heat generated by the optical amplifiers is dissipated by a watercooled chiller unit. The optical antenna is equipped with a 150 mm largeaperture telescope for long range observations and a 50 mm smallaperture telescope for vector conversion of shortrange observations. Each lidar measures the LOS wind velocity with an observational accuracy of ±0.09 m s^{−1}; the paired values are used to estimate the airflow vector in the region between the lidars. The observational resolution of each lidar is approximately 25 m. There are additional performance requirements for preview control: the estimation frequency and estimation accuracy of vertical wind velocity. The frequency of estimation must be more than 5 Hz, and the estimation accuracy of the vertical wind velocity must be better than 2.6 m s^{−1} in the LOS distance of 500 m. The control requirements are the conditions that are necessary for halving the peak variation in acceleration by control. This value has been specified using control simulations (Hamada, 2019), and Monte Carlo simulations have also been performed.
Next, an existing technique for estimating the airflow vector from a pair of LOS wind values is reviewed. The airflow vector in the region between the upward and downward lidars is conventionally estimated via simple vector conversion. This procedure is similar in concept to the vertical azimuth display approach used in general ground lidar systems (Newsom et al., 2017). The simple vector conversion is given by
where ${u}_{x}^{T}$ and ${u}_{z}^{T}$ are the horizontal and vertical windvelocity measurements at the observation time T; ${W}_{\mathrm{1}}^{T}$ and ${W}_{\mathrm{2}}^{T}$ are the LOS wind velocities of the upward and downwarddirected lidars at the observation time T; and θ is the angle between the horizontal line and each lidar, which is 10^{∘} in this study. The simple vector conversion assumes that the windfield region between the lidars is homogeneous (Newsom et al., 2017). The assumption of homogeneity seems natural: the regions between the lidars are 69.5 and 173.6 m at the LOS distances of 200 and 500 m ahead of the aircraft (Fig. 3). Nevertheless, the assumption would be wrong if a large fluctuation in wind velocity occurs, creating turbulence. In homogenous conditions, a simple vector conversion can estimate the airflow vector accurately; however, in nonhomogenous conditions, the estimation is expected to have poor accuracy.
2.2 Estimation algorithm based on extrapolation
Whereas simple vector conversion utilizes recent LOS wind data to estimate the airflow vector, our proposed method stores the LOS wind data continuously and uses both recent and past values to extrapolate the wind field in the region between the lidars where it has not been directly measured. The airflow vector is then calculated from Eq. (1) and the extrapolated horizontal and vertical components of the wind velocity. The airflowvectorestimation accuracy far ahead of the aircraft is improved relative to simple vector conversion by using the extrapolated wind data because the region between the upward and downward lidars is no longer assumed to be homogeneous; our algorithm uses a polynomial expression to extrapolate data points from both recent and past measurements, allowing it to be used in nonhomogenous wind fields. In addition, the proposed method can estimate the twodimensional distribution of the wind field between the lidars, again unlike simple vector conversion.
Figure 4 shows the overview of the proposed estimation method when a current data point and two past data points are used. When the aircraft speed is V and the time span of observation is dt, the airflow moves backwards at V×dt because the aircraft is advancing. Current observation times are denoted as T and past observation times as T−1 and T−2. The proposed method uses the current LOS wind values (${W}_{\mathrm{1}}^{T}$ and ${W}_{\mathrm{2}}^{T}$) and the past LOS wind values (${W}_{\mathrm{1}}^{T\mathrm{1}}$, ${W}_{\mathrm{2}}^{T\mathrm{1}}$ and ${W}_{\mathrm{1}}^{T\mathrm{2}}$, ${W}_{\mathrm{2}}^{T\mathrm{2}}$). The perpendicular distances between the horizontal line and each lidar are denoted as z^{T}, z^{T−1}, and z^{T−2}, respectively. Depending on the amount of past LOS wind data used, the order of the polynomial expression used in the extrapolation varies. The aerosol concentration in the upper sky is low, suggesting that there are considerable missing data and noise. A sufficient amount of past LOS wind data may not be available to estimate a highorder polynomial expression, and this could affect the robustness of the control. For this reason, a firstdegree polynomial expression is adopted in this study and used in the leastsquares method (LSM) to extrapolate the windfield values according at the horizontal line. The airflow vector is calculated by Eq. (1) using the extrapolated LOS wind. The equation used in the extrapolation method is
where
2.3 Filtering error and the lack of windvelocity data
In this study, two filtering algorithms are used to remove the error and the loss of data in airborne lidars. First, a filtering algorithm that is a simple representation of a Kalman filter with simplified Kalman gain is used; this filtering algorithm is described in detail in the study of Misaka et al. (2015). The algorithm assumes that infinite variance is used to exclude outliers and loss of data. This method uses the lidar spectrum data at each range bin; the algorithm defines the validity of the measurements during the lidar data peakdetection process. To identify the correct and incorrect LOS windvelocity values, two spectrum thresholds are defined. First, the largest and secondlargest spectrum values, k_{1st} and k_{2nd}, which are the fast Fourier transform points for the first and second spectrum peaks, respectively, are adjacent to each other; i.e., the magnitude of the distance between the largest and secondlargest spectrum values in the fast Fourier transform is equal to 1. Second, the distance between k_{1st} and the averaged spectrum peak k_{ave} is required to be less than a certain value k_{dif}, which represents the only hyperparameter in this algorithm as well as a parameter related to smoothness. k_{ave} is the index that conveys the location of the spectrum peak averaged in short ranges, e.g., 2–30 range bins from the lidar origin. Figure 5 shows a conceptual explanation of the variables of simplified Kalman gain in the cases of correct measurement and of an error peak. In this study, the filtering algorithm is carried out first when the observation data are obtained:
Secondly, a robust leastsquares estimation, based on Tukey's biweight methodology (Huber, 2008), is carried out to reduce the impact of the error in the LOS wind velocity. This method is based on the LOS wind data in contrast to the spectrum data from lidar observations in the first method. Although the filtering algorithm based on a simple Kalman filter can remove the error from the lidar spectrum data, error filtering via this algorithm is not perfect despite being useful. As error data can be a reason for miscontrol, it is essential to deal with the error and the loss of data of the lidars more carefully when the filtering algorithm is used for the preview control. Therefore, the robustness of the estimated airflow vector is secured by combining the simple Kalman filtering algorithm with the results of robust LSM, using Eqs. (2) and (3). In addition, the robust LSM estimation can employ the extrapolation algorithm effectively as per Eqs. (2) and (3). Therefore, a simpler and more robust algorithm is provided. Figure 6 explains the concept behind Tukey's biweight methodology as applied to lidar. The fundamental principle involves comparing the observed LOS wind values with the estimated ones from the polynomial expression used in the LSM. In the first step, the LOS wind is estimated using the general LSM (Eq. 2). In the second step, the difference ${d}_{j}^{T}$ between the observed LOS wind value and that estimated from the polynomial expression is found:
A permissible difference range L is defined and weights ${w}_{j}^{T}$ (${d}_{j}^{T}$) are calculated depending on where ${d}_{j}^{T}$ falls in the distance range:
Weights are assigned to each LOS windvelocity value. In the third step, a new firstdegree polynomial expression for the LSM with the weighted data is estimated as follows.
This process is repeated until the weight of the error value decreases and converges.
2.4 Filtering windvelocity noise
Lidar is subject not only to measuring errors and loss of LOS data values but also to random noise; this type of noise also leads to a poor estimation of the airflow vector. The random noise is caused by the reduced intensity of the received light due to the thin aerosol concentration in the sky. A general lidar signal consists of random noise superimposed on the spectral signal. If the signal intensity is low, peak search may only detect the random noise. (Additional randomness caused by environmental factors and data processing in lidar is considered here as randomness of the windspeed values.)
A simple spline algorithm generates a curve that passes through all sample points; therefore, it is not able to generate a smooth curve when the sample points have random noise, and a smoothing spline algorithm is often applied to remove the random noise in the lidar LOS wind values, as in the study by Woltring (1986). The curve generated by this algorithm does not pass through all sample points, and because of that, it can be smoother, even when there is random noise from lidar LOS wind measurements. The smoothing spline model minimizes the criterion function C_{p},
where y_{i} is a sample point value, s_{p}(x) is the value generated by a simple spline algorithm, v_{i} is a weighted factor, and p is the regularization parameter. The smoothest curve is generated when the criterion function C_{p} is minimized.
2.5 System flowchart
The airflowvectorestimation algorithm is a sequence of five different processes, which are summarized below. The system flowchart is shown in Fig. 7.

The filtering algorithm based on a simple Kalman filter is used to remove the error in lidar LOS winddata values.

The smoothing spline method is applied to reduce the negative effect of the random noise in LOS winddata values and extrapolates the values at positions for which no measurements can be read. This is identified as the firststep error.

Extrapolation, based on the polynomial expression, is carried out to estimate the windfield values by using current and past LOS wind data.

A robust LSM model is applied to obtain a more accurate polynomial expression. The calculation repeats until the parameter converges.

The airflow vector is calculated by Eq. (1) with the extrapolated LOS wind.
3.1 Ideal vortex model
We have conducted numerical experiments to evaluate the performance of actual airborne lidars. The ideal vortex model is defined and used to evaluate the estimated performance of the airflow vector. In this study, the Burnham–Hallock vortex model (Hinton and Tatnall, 1997) is used. The experiment generates a large number of pseudolidar values, from which the airflow vector is estimated. The estimation results are then compared with the reference windfield values of the ideal vortex model. Figure 8 shows the distribution of wind velocity generated using the Burnham–Hallock vortex model.
3.2 NWP model
The results predicted by a numerical weather model – the Japan Meteorological Agency NonHydrostatic Model (JMANHM) – are used to evaluate airflowvectorestimation performance (Saito et al., 2007; Kikuchi et al., 2015). To obtain highresolution weather prediction, a oneway multinesting technique (Kikuchi et al. 2015) is employed for downscaling purposes. The computational domain is nested four times to increase grid resolutions from 5.0 to 0.05 km gradually (in the sequence of 5.0, 1.5, 0.5, 0.15, and 0.05 km).
The 3 h mesoscale objective analysis data, collected using a mesoscale fourdimensional variational dataassimilation system at the Japan Meteorological Agency (Saito et al., 2007), are used for the initial condition of 5.0 km grid resolution. The experiment generates a large number of simulated twinlidar observation values along flight routes from the windfield data generated by JMANHM, which are more realistic than idealvortex model results. The airflow vector is estimated from the pseudolidar observations and compared with the JMANHM reference field. Figure 9 shows the distribution of the vertical windvelocity values generated by JMANHM.
3.3 Generation of pseudoerrors and noise
To confirm the effectiveness of the proposed filtering algorithms, errors and noise are generated artificially by using the parameter of the backscattering coefficient in the atmosphere and the statisticsbased coherent lidar equation (Kameyama et al., 2007). The backscattering coefficient is strongly related to the aerosol density in the atmosphere, and it has an impact on the lidar measurements and estimation performance. When the backscattering coefficient is very low, the measurement performance is worse, and the LOS wind data show errors and noise. Apart from this, the measurement performance is related to the focal distance, pulse width, and lidar power (Kameyama et al., 2007). The signaltonoise ratio (SNR) at the receiver, at each LOS distance, is calculated by using the coherent lidar equation and the detailed operating condition of JAXA's lidar (Inokuchi and Akiyama, 2019):
Here, R is the observation distance, η is the system efficiency, P_{t} is the lighttransmission power, ΔR is the resolution range, β is the backscattering coefficient, K is the atmospheric transmittance, D is the opening size of the optical antenna, h is Planck's constant, f is optical frequency, B is received bandwidth, F is focal distance, k is wave number, A_{c} is the vignetting factor of the optical antenna, and ${C}_{n}^{\mathrm{2}}$ is the atmospheric structure constant. In this study, the conditions are set according to the design specification for airborne lidars. Six atmospheric conditions are prepared in order to evaluate the filtering performance. The backscattering coefficients are (standard case) $\mathrm{1.8}\times {\mathrm{10}}^{\mathrm{8}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (a) $\mathrm{1.8}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (b) $\mathrm{1.35}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (c) $\mathrm{0.9}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (d) $\mathrm{0.45}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, and (e) $\mathrm{0.18}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$. Figure 10 shows the statistics for the error and noise as functions of SNR bandwidth.
4.1 Ideal vortex model without error and noise
The numerical experiments with the ideal vortex model have been carried out, and Figs. 11 and 12 show the distributions of the horizontal and vertical wind components that are estimated by the simple vector conversion and the proposed method. The flights start at the edge of the computational space. Figures 11 and 12 show the results after 10 and 15 s, respectively. Thus, they represent the instants of time before and during the aircraft's close approach to the vortex core. As shown in Figs. 11 and 12, the simple vector conversion method, which assumes that the wind field of the region between the lidars is homogeneous, cannot accurately reproduce the twodimensional distribution between the lidars. On the other hand, the figures confirm that the proposed method can estimate the twodimensional distribution of windfield values between the lidars. Figure 11 shows that the twodimensional distribution obtained with the proposed method is very similar to that of the reference field. In addition, the results show that the horizontal wind velocity with simple vector conversion is approximately −7 m s^{−1}, whereas that with the proposed method is −9.5 m s^{−1}; the horizontal wind velocity of the reference field is −9.0 m s^{−1} at a LOS distance of 450–500 m. Figure 12 shows that the results of the horizontal and vertical wind velocities with simple vector conversion are considerably lower than those of the reference field. The horizontal wind results show that the value obtained with the simple vector conversion is approximately −9.5 m s^{−1}, whereas that with the proposed method is approximately −3.5 m s^{−1}; the horizontal wind velocity of the reference field is approximately −4.5 m s^{−1} at a LOS distance of 450–500 m. The vertical wind results show that the value obtained with simple vector conversion is approximately −1.0 m s^{−1}, whereas that obtained with the proposed method is approximately 8.5 m s^{−1}; the vertical wind velocity of the reference field is approximately 7.0 m s^{−1} at a LOS distance of 450–500 m. Therefore, simple vector conversion has significantly large errors between the reference and estimated values. The errors in both the horizontal and vertical wind values estimated by the proposed method are much smaller than those estimated with simple vector conversion. Although the twodimensional distribution of the horizontal windfield values of the proposed method is larger than that of the reference field at a LOS distance of 450–500 m, the vertical windfield values can provide a good assessment of the reference field shown in Fig. 12. The 15 s timing in Fig. 12 is a more challenging case than others because the aircraft is positioned very close to the center of the vortex, and the wind direction changes abruptly. Although it is difficult to estimate the perfect windfield value at this time by using the proposed method, the proposed estimation method demonstrably has a much higher accuracy than simple vector conversion. Overall, the proposed method has much better performance than the simple vector conversion method, and it can estimate the twodimensional distribution of windfield values accurately, unlike the simple vector conversion method.
Next, the statistical estimation performance is evaluated using 100 pseudoroutes that are randomly generated 750 m above and below the center of the vortex core; Fig. 13 shows the results for the vertical wind values, along with the performance required for automatic control. The root mean square error (RMSE) between the referencefield value and the estimated windfield value is used for evaluating the estimation performance. Moreover, the effect of the number of past lidar observations used to determine the wind field, i.e., the past LOS wind, is checked. Simple vector conversion cannot satisfy the performance requirement at a LOS distance greater than 350 m. This means that achieving preview control using the simple vector conversion method may be difficult. At a LOS distance of 500 m, the RMSEs of the vertical wind values of the simple vector conversion and proposed methods are approximately 4.0 and 1.2 m s^{−1}, respectively. The proposed method can cater to the performance demand even if the number of past LOS wind values used is different; a lower number leads to better estimation performance.
4.2 Numerical weather prediction without error and noise
We also conducted numerical experiments with NWP values. Figs. 14 and 15 show the distributions of the horizontal and vertical wind components that are estimated by simple vector conversion and the proposed method. Figure 14 shows the results for the instants of time before and during the approach to a vertical wind fluctuation. The simple vector conversion method cannot accurately reproduce the twodimensional distribution of the wind field between the lidars. On the other hand, the proposed method can estimate the twodimensional distribution of the wind field between the lidars more accurately. Figure 15 shows that the wind velocities predicted by the simple vector conversion method are higher than the reference fields at 300–500 m of LOS distance in contrast to those of the proposed method.
Next, the statistical estimation performance is evaluated using 100 pseudoroutes that are randomly generated between 2 and 10 km altitude. Figure 16 shows the results, along with the performance requirement for automatic control. The effect of the number of past LOS wind values used is also checked. In this case, both simple vector conversion and the proposed method can satisfy the performance demand for preview control; however, the performance results of simple vector conversion are much worse than those of the proposed method. Moreover, the proposed method can estimate quite accurate windfield values. In this case, the use of a higher number of past LOS wind values leads to better estimation performance.
4.3 Ideal vortex model with error and noise
In this section, numerical experiments with error and noise in LOS wind values are conducted to evaluate the estimation performance of the proposed method. These numerical experiments show the error/noisefiltering performance difference between simple vector conversion and the proposed method with extrapolation from the past LOS wind. Six atmospheric conditions are prepared in order to evaluate the filtering performance. The backscattering coefficients are (standard case) $\mathrm{1.8}\times {\mathrm{10}}^{\mathrm{8}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (a) $\mathrm{1.8}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (b) $\mathrm{1.35}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (c) $\mathrm{0.9}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, (d) $\mathrm{0.45}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$, and (e) $\mathrm{0.18}\times {\mathrm{10}}^{\mathrm{11}}$ ${\mathrm{sr}}^{\mathrm{1}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{1}}$.
First, numerical experiments with the ideal vortex model are carried out. Figure 17 shows the LOS wind values, which include the measured data with error and noise, the reference wind, the smoothing spline, and the general spline model results. Figure 17 shows that the smoothing spline can filter the error and noise data of LOS wind values. When the general spline is used, the error can be filtered correctly by using a simple Kalman filter and a robust LSM; however, the noise cannot be filtered. Next, the statistical estimation performance is evaluated using 100 pseudoroutes that are randomly generated 750 m above and below the center of the vortex core. Figure 18 shows the results of the statistical estimation performance with error and noise. In addition, the difference due to the atmospheric conditions in the six cases with different backscattering coefficients is also checked. Simple vector conversion cannot satisfy the performance demand at a distance farther than 350 m LOS and cannot work correctly under atmospheric condition (e). The proposed method can always satisfy the performance demand except under atmospheric condition (e). It thus shows much better performance than simple vector conversion, even though it is difficult to estimate the windfield values by either method for atmospheric condition (e), which contains much larger noise levels than the other conditions.
In addition, the crossplots of the reference and the estimated vertical wind are shown as Fig. 19. In Fig. 19a and b, the results of the simple vector conversion are presented; Fig. 19c and d show the results of the proposed method. Figure 19a and c are the cases without error and noise, whereas Fig. 19b and d are the cases with error and noise. By comparing Fig. 19a and c, we can deduce that the proposed method provides a much better estimation than does simple vector conversion. The results in Fig. 19b and d are spread wider than those in Fig. 19a and c because of the noise data of LOS wind values. It is worth mentioning that the noise data have more negative effects on the result at 500 m LOS distance than at 100 and 300 m LOS. Nevertheless, comparison of Fig. 19b and d shows that the proposed method can provide more accurate estimations than the simple vector conversion method.
4.4 Numerical weather prediction with error and noise
We also carry out numerical experiments with NWP. The statistical estimation performance is conducted by using 100 pseudoroutes between 2 and 10 km altitude. Figure 20 shows the results of the statistical estimation performance with error and noise. Six different atmospheric conditions (standard, (a), (b), (c), (d), and (e), defined by their backscattering coefficients) are used. In this case, both simple vector conversion and the proposed method can satisfy the performance requirement for preview control; however, the simple vector conversion shows worse performance than the proposed method. The proposed method can estimate windfield values quite accurately and displays better performance than the simple vector conversion method. As in the previous experiment, it is difficult to estimate the windfield values for atmospheric condition (e) by using either simple vector conversion or the proposed method.
In this study, an airflowvectorestimation algorithm based on upward and downward airborne lidars has been proposed for preview control to prevent turbulenceinduced aircraft accidents. This estimation algorithm uses the technique of extrapolating the windfield values by using the LSM and the current and past LOS wind datasets to improve the accuracy of estimated wind values. Two test configurations for numerical experiments (ideal vortex flow and realistic NWP weather field values) have been used to evaluate the estimation of the airflow vector.
Numerical experiments on LOS wind estimation show that the proposed extrapolation method has much better performance than simple vector conversion methods, and it can estimate the twodimensional distribution of windfield values accurately, which simple vector conversion cannot. The estimation performance and the computational cost of the proposed method can satisfy the performance demand for preview control.
Numerical experiments with error and noise in the LOS wind data have been conducted to evaluate the performance of the proposed estimation method. These numerical experiments show that the smoothing spline model can filter noise correctly and reduce its negative effects. The proposed method performs much better than the simple vector conversion method, although it is difficult to estimate the windfield values for atmospheric condition (e) with either method. Atmospheric condition (e) has more noise than other conditions, and when the noise exceeds a certain level, it becomes difficult to estimate the air flow regardless of the method applied.
The proposed algorithm can satisfy the performance demands for preview control in both estimation performance and computational cost. It can estimate a twodimensional distribution that cannot be estimated by existing methods. This is valuable for improving the accuracy of the preview control: for example, the proposed method can cope with the critical case where the flight direction of the aircraft is at a steep angle with the aircraft either ascending or descending.
The findings of this study are subject to certain limitations. The target size of the atmospheric turbulence is assumed by the proposed algorithm to be comparable to or larger than the observation region between the lidars. Therefore, it is difficult to estimate a wind field with turbulence smaller than this. The effect on the aircraft vibration due to such minor turbulence, however, is minimal. An exception to this is aircraftgenerated wake turbulence, which still poses a safety risk. The radius of the actively fluctuating waketurbulence core is only a few meters, so the proposed method could lead to erroneous predictions. A second limitation is that the current results are obtained from numerical experiments and not from evaluations of actual observations. Currently, the lidar system is being modified to be smaller and lighter in order to suit small experimental aircraft. The onboard lidar system and realtime airflowvector estimation will be validated by flight experiments in 2021; the whole gustalleviation system, including preview control, will be demonstrated in 2022. The results of this research will be applied to this flight demonstration.
The data used in this paper are available from the corresponding author upon request.
RK, TM, and SO designed the experiments. RK performed the experiments, developed the model code, performed the simulations, and prepared the manuscript with contributions from all coauthors. HI contributed to the analysis and interpretation of data related to lidar and assisted in the preparation of the manuscript. All authors approve the final version of the manuscript and agree to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
The authors declare that they have no conflict of interest.
We acknowledge Naoki Matayoshi and Eiichi Yoshikawa of JAXA for their advice regarding data analysis of lidar.
This research has been supported by the Japan Society for the Promotion of Science (JSPS) KAKENHI (267391).
This paper was edited by Murray Hamilton and reviewed by two anonymous referees.
Airbus S.A.S.: Flight operations briefing notes: adverse weather operations – optimum use of the weather radar, available at: https://safetyfirst.airbus.com/optimumuseofweatherradar/, last access: 22 September 2020.
Barny, H.: DELICAT – Demonstration of Lidar Based Clear Air Turbulence Detection, in: Innovation for Sustainable Aviation in a Global Environment: Proceedings of the Sixth European Aeronautics Days, 30 March–1 April 2011, Madrid, Spain, 253, 2012.
Boeing Commercial Airplanes: Statistical summary of commercial jet airplane accidents, available at: http://www.boeing.com/resources/boeingdotcom/company/about_bca/pdf/statsum.pdf (last access: 4 December 2020), 2018.
Federal Aviation Administration: Preventing injuries caused by turbulence 2006, available at: https://www.faa.gov/regulations_policies/advisory_circulars/index.cfm/go/document.information/documentid/99831 (last access: 23 September 2020), 2006.
Fezans, N., Joos, H. D., and Deiler, C.: Gust load alleviation for a longrange aircraft with and without anticipation, Aeronaut. J., 10, 1–25, 2019.
Hamada, Y.: New LMIbased conditions for preview feedforward synthesis, Control. Eng. Pract., 90, 19–26, 2019.
Hinton, D. A. and Tatnall, C. R.: A candidate wake vortex strength definition for application to the NASA aircraft vortex spacing system (AVOSS), NASA Technical Reports, available at: https://ntrs.nasa.gov/api/citations/19970031683/downloads/19970031683.pdf (last access: 4 December 2020), 1997.
Huber, P. J.: Robust Statistics, Springer, Berlin, Heidelberg, Germany, 2008.
Inokuchi, H.: Detection of the clear air turbulence by an airborne doppler LIDAR, in: Proceedings of AsiaPacific International Symposium on Aerospace Technology, 13–15 November 2012, Jeju, Korean, 2012.
Inokuchi, H. and Akiyama, T.: True airspeed measured by an airborne coherent doppler lidar, in: APISAT 2019: Asia Pacific International Symposium on Aerospace Technology, 4–6 December 2019, Gold coast, Australia, p. 554, 2019.
Inokuchi, H., Tanaka, H., and Ando, T.: Development of an onboard doppler lidar for flight safety, J. Aircraft, 46, 1411–1415, 2009.
International Air Transportation Association: Safety Report 2015, available at: https://skybrary.aero/bookshelf/books/3408.pdf (last access: 4 December 2020), 2016.
Japan Transport Safety Board: 2003–2012: Aircraft accident reports, available at: http://www.mlit.go.jp/jtsb/, last access: 23 September 2020.
Kameyama, S., Ando, T., Asaka, K., Hirano, Y., and Wadaka, S.: Compact allfiber pulsed coherent Doppler lidar system for wind sensing, Appl. Optics, 46, 1953–1962, 2007.
Kikuchi, R., Misaka, T., and Obayashi, S.: Realtime flow prediction of lowlevel atmospheric turbulence, in: 33rd Wind Energy Symposium, 5–9 January 2015, Florida, USA, 2015.
Kim, J. H. and Chun, H. Y.: Statistics and possible sources of aviation turbulence over South Korea, J. Appl. Meteorol. Clim., 50, 311–324, 2011.
Kim, J. H., Chun, H. Y., Sharman, R. D., and Keller, T. L.: Evaluations of upperlevel turbulence diagnostics performance using the Graphical Turbulence Guidance (GTG) system and pilot reports (PIREPs) over East Asia, J. Appl. Meteorol. Clim., 50, 1936–1951, 2011.
Kim, J. H., Chan, W. N., Sridhar, B., and Sharman, R. D.: Combined winds and turbulence prediction system for automated airtraffic management applications, J. Appl. Meteorol. Clim., 54, 766–784, 2015.
Machida, S.: Project Overview of R&D for Onboard Turbulence Detection System, in: Proceedings of AsiaPacific International Symposium on Aerospace Technology 2017, Plenary Lecture 4, 16–18 October 2017, Seoul, Korea, 2017.
Misaka, T., Nakabayashi, F. K., Obayashi, S., and Inokuchi, H.: Filtering algorithm of airborne Doppler lidar measurements for improved wind estimation, T. Jpn. Soc. Aeronaut. S., 58, 149–155, 2015.
Neininger, B.: Trends in airborne atmospheric observations, European Meteorological Society Annual Meeting 2017, 4–8 September 2017, Dublin, Ireland, 14, EMS2017322, 2017.
Newsom, R. K., Brewer, W. A., Wilczak, J. M., Wolfe, D. E., Oncley, S. P., and Lundquist, J. K.: Validating precision estimates in horizontal wind measurements from a Doppler lidar, Atmos. Meas. Tech., 10, 1229–1240, https://doi.org/10.5194/amt1012292017, 2017.
Regan, C. D. and Jutte, C. V.: Survey of applications of active control technology for gust alleviation and new challenges for lighterweight aircraft, NASA Technical Memorandum, NASATM2012216008, available at: https://ntrs.nasa.gov/citations/20120013450 (last access: 4 December 2020), 2012.
Saito, K., Ishida, J. I., Aranami, K., Hara, T., Segawa, T., Narita, M., and Honda, Y.: Nonhydrostatic atmospheric models and operational development at JMA, J. Meteorol. Soc. Jpn., 85B, 271–304, 2007.
Sakimura, T., Watanabe, Y., Ando, T., Kameyama, S., Asaka, K., Tanaka, H., Yanagisawa. T., Hirano, Y., and Inokuchi, H.: 3.2 mJ, 1.5 µm laser power amplifier using an Er, Yb:glass planar waveguide for a coherent Doppler LIDAR, in: Proceedings of 17th Coherent Laser Radar Conference, 17–20 June 2013, Barcelona, Spain, 2013.
Schmitt, N. P., Rehm, W., Pistner, T., Zeller, P., Diehl, H., and Navé, P.: The AWIATOR airborne LIDAR turbulence sensor, Aerosp. Sci. Technol., 11, 546–552, 2007.
Sermi, F., Cuccoli, F., Mugnai, C., and Facheris, L.: Aircraft hazard evaluation for critical weather avoidance, in: IEEE Metrology for Aerospace (MetroAeroSpace), 4–5 June 2015, Benevento, Italy, 454–459, 2015.
Sharman, R., Tebaldi, C., Wiener, G., and Wolff, J.: An integrated approach to mid and upperlevel turbulence forecasting, Weather. Forecast., 21, 268–287, 2006.
Soreide, D. C., Bogue, R. K., Ehernberger, L. J., Hannon, S. M., and Bowdle, D. A.: Airborne coherent LIDAR for advanced inflight measurements (ACLAIM) flight testing of the LIDAR sensor, NASA Technical Reports, available at: https://ntrs.nasa.gov/citations/20000094555 (last access: 4 December 2020), 2000.
Wei, T., Xia, H., Hu, J., Wang, C., Shangguan, M., Wang, L., Jia, M., and Dou, X.: Simultaneous wind and rainfall detection by power spectrum analysis using a VAD scanning coherent Doppler lidar, Opt. Express, 27, 31235–31245, 2019.
Williams, P. D.: Increased light, moderate, and severe clearair turbulence in response to climate change, Adv. Atmos. Sci., 34, 576–586, 2017.
Woltring, H. J.: A Fortran package for generalized, crossvalidatory spline smoothing and differentiation, Adv. Eng. Softw., 8, 104–113, 1986.