Analysis of Mountain Wave Effects on a Hard Landing Incident in Pico Aerodrome Using the AROME Model and Airborne Observations

A hard landing incident in Pico Aerodrome (LPPI) involving an Airbus A320-200 aircraft is investigated using airborne observations and forecasts of the AROME (Applications of Research to Operations at Mesoscale) model. A second flight is also analyzed. The severity of the wind shear during both flights is quantified using the intensity factor “I” that is based on aerial data and recommended by ICAO (International Civil Aviation Organization). During Flight 1, 36% of the landing phase (below 2100 ft) occurred under “severe” wind shear conditions and 16% occurred under “strong” conditions. Upstream characteristics included southwest winds, stable stratification and a Froude number close to 1. According to the AROME model, these circumstances triggered the development of vertically propagating mountain waves, with maximum vertical velocities above 400 ft/min and exceeding 200 ft/min in the flight path. These conditions, together with the severe wind shear, may have caused the incident. During the second flight, a wake with lee vortices and reversed flow developed in the region of the flight path, which is consistent with a low upstream Froude number and/or with the flow regime diagram of previous studies. During the approach phase of this flight, “severe” wind shear conditions were absent, with “strong” ones occurring 4% of the time. It predominantly displayed “light” conditions during 68% of this phase. As a result of the comparison between “I” and the AROME turbulence indicators, preliminary thresholds are proposed for these indexes. Lastly, this study provides an objective verification of AROME wind forecasts, showing a good agreement with airborne observations for wind speeds above 10 kt, but a poor skill for weaker winds.


Introduction
In a stably stratified atmosphere, the airflow towards a topographic obstacle may trigger gravity (or buoyancy) waves downstream of the obstacle, known as mountain waves [1]. The flow response to the obstacle depends on several factors like the mountain shape and the Froude number [2]. At a low Froude number, several phenomena have been observed: the flow may be blocked upwind [3] or it may also pass around the obstacle rather than over it (flow splitting phenomenon) [2,4]. In association with the flow splitting phenomenon, a low-level wake may develop on the leeside of the obstacle [4][5][6][7]. The wake can be expressed as a quasi-steady pair of counter-rotating vortices circulating about vertical Internacional-Azores Airlines, S.A.) Airbus A320-200 aircraft. Two flights were considered, where the first one represents the hard landing incident (henceforth referred to as Flight 1). Upon touchdown, the aircraft bounced and registered a vertical acceleration above 2 G. No fatalities or serious injuries were reported. The second one, labeled Flight 2, will serve as a control parameter that is representative of a flight in which no incidents were reported.
This study has three goals. The first goal is to characterize the nature of the orographic phenomenon involved during the two flights. The second is to provide an assessment of the agreement between AROME wind forecasts and airborne observations. The third objective is to establish a correspondence between AROME turbulence indicators and the wind shear severity registered by the aircraft. For this purpose, five turbulence indicators (Brown, Ellrod TI1, Ellrod TI2, CAT1 and EDR or Eddy Dissipation Rate) are computed from AROME forecasts. The wind shear perceived by the aircraft is quantified using the wind shear intensity factor "I", as originally suggested by Woodfield and Woods [23]. This assessment is also of great importance because the Portuguese Institute for Sea and Atmosphere (IPMA) routinely provides AROME forecasts of the wind and EDR for Pico and Faial Islands to SATA on an experimental basis. This test phase started in November 2017.

Data and Methodology
This section characterizes the wind from aerial data gathered by SATA A320-200 aircraft during both flights. The method of calculation of the wind shear intensity factor "I", which is a function of these airborne measurements, is established. A brief description of the AROME model is also provided. The Froude number and turbulence indicators, which are based on forecasts, are then defined. Lastly, definitions for the RMSE (root mean square error) and MAE (mean absolute error) concerning wind magnitude and direction are included.

SATA Airborne Measurements
For each flight, wind intensity and direction were recorded approximately every 4 seconds along with the corresponding altitude and geographical coordinates (latitude and longitude). Only the approach and landing phases (approximately final 10 minutes of flight or altitudes below 5000 ft) were perused since these are the most pertinent to understanding the incident. The average trajectory for both flights along with Pico's orography in the AROME model are mapped in Figure 2. Forecasts of mountain wave induced turbulence rely on algorithms based on outputs of numerical weather prediction (NWP) models. Many of these algorithms rely on wind, temperature and turbulent kinetic energy [15]. Therefore, the accuracy of these fields is a key issue. However, an accurate simulation of such atmospheric variables is particularly challenging in the atmospheric boundary layer (ABL) due to its unsteady and turbulent nature [16]. Several factors may influence the model accuracy, namely, unresolved topographic features [17], turbulence closure schemes [18,19] as well as the representation of turbulence mixing [20,21]. The dynamical core also plays a key role. Ricard et al. [22] for instance, in analyzing the kinetic energy spectra of AROME and Meso-NH models (which have the same physical parameterizations), concluded that a model using a Semi-Implicit Semi-Lagrangian (SISL) scheme has a coarser effective resolution than a model using an Eulerian explicit scheme (less efficient).
The current article focuses on the orographic flow associated to Pico Mountain in two different synoptic backgrounds, using forecasts from the operational AROME (Applications of Research to Operations at Mesoscale) model and airborne measurements collected by SATA (SATA Internacional-Azores Airlines, S.A.) Airbus A320-200 aircraft. Two flights were considered, where the first one represents the hard landing incident (henceforth referred to as Flight 1). Upon touchdown, the aircraft bounced and registered a vertical acceleration above 2 G. No fatalities or serious injuries were reported. The second one, labeled Flight 2, will serve as a control parameter that is representative of a flight in which no incidents were reported.
This study has three goals. The first goal is to characterize the nature of the orographic phenomenon involved during the two flights. The second is to provide an assessment of the agreement between AROME wind forecasts and airborne observations. The third objective is to establish a correspondence between AROME turbulence indicators and the wind shear severity registered by the aircraft. For this purpose, five turbulence indicators (Brown, Ellrod TI1, Ellrod TI2, CAT1 and EDR or Eddy Dissipation Rate) are computed from AROME forecasts. The wind shear perceived by the aircraft is quantified using the wind shear intensity factor "I", as originally suggested by Woodfield and Woods [23]. This assessment is also of great importance because the Portuguese Institute for Sea and Atmosphere (IPMA) routinely provides AROME forecasts of the wind and EDR for Pico and Faial Islands to SATA on an experimental basis. This test phase started in November 2017.

Data and Methodology
This section characterizes the wind from aerial data gathered by SATA A320-200 aircraft during both flights. The method of calculation of the wind shear intensity factor "I", which is a function of these airborne measurements, is established. A brief description of the AROME model is also provided. The Froude number and turbulence indicators, which are based on forecasts, are then defined. Lastly, definitions for the RMSE (root mean square error) and MAE (mean absolute error) concerning wind magnitude and direction are included.

SATA Airborne Measurements
For each flight, wind intensity and direction were recorded approximately every 4 s along with the corresponding altitude and geographical coordinates (latitude and longitude). Only the approach and landing phases (approximately final 10 min of flight or altitudes below 5000 ft) were perused since these are the most pertinent to understanding the incident. The average trajectory for both flights along with Pico's orography in the AROME model are mapped in Figure 2.  ) and AROME orography. Shading shows elevation in meters. Pink marker shows the aerodrome location. "AB" and "CD" are endpoints of cross-sections analyzed in subsequent sections.

Wind Shear Intensity Factor "I"
The severity of wind shear intensity, as perceived by the aircraft, will be measured according to the wind shear intensity factor "I", detailed in ICAO's (International Civil Aviation Organization) Manual on Low-level Wind Shear [24]. Following Guan and Yong's study [25], wind shear is classified as "light", "moderate", "strong" or "severe" based on the variation of the airspeed with time and on the airspeed proportion . The intensity factor "I" itself may be expressed as Figure 2. Approximate flight path (blue dashed line) and AROME orography. Shading shows elevation in meters. Pink marker shows the aerodrome location. "AB" and "CD" are endpoints of cross-sections analyzed in subsequent sections.

Wind Shear Intensity Factor "I"
The severity of wind shear intensity, as perceived by the aircraft, will be measured according to the wind shear intensity factor "I", detailed in ICAO's (International Civil Aviation Organization) Manual on Low-level Wind Shear [24]. Following Guan and Yong's study [25], wind shear is classified as "light", "moderate", "strong" or "severe" based on the variation of the airspeed with time dV TAS dt and on the airspeed proportion . The intensity factor "I" itself may be expressed as follows: The computation of the true airspeed, V TAS , is based on the following relation: where GS represents the ground speed of the aircraft and was computed by employing the three-dimensional distance formula between any two points described by their altitude and geographic coordinates (latitude and longitude). Note that ∆V TAS is simply the difference between true airspeeds at the n th and (n + 1) th instants. The head or tailwind component (see Equation S2.1 and Figure S2, from Supplementary Materials) is then subtracted from the ground speed to obtain the true airspeed. Once dV TAS dt and ∆V TAS V TAS are known, the wind shear intensity may be classified by consulting the proper boundaries specified in Figure S3 from the Supplementary Materials. In both flights, this was performed for altitudes approximately below 2000 ft (610 m), which according to Guan and Yong [25], represents the appropriate range for low-level wind shear considerations.

Description of the AROME Model
The AROME model that runs operationally at IPMA uses a horizontal grid spacing of 2.5 km. It uses a SISL scheme and employs a three-class ice scheme with ECMWF (European Center for Medium-Range Weather Forecasts) radiation parameterizations, as was described by Seity et al. [26]. The representation of turbulence in the ABL is based on a prognostic TKE (Turbulent Kinetic Energy) equation [27], combined with a diagnostic mixing length. AROME initial and boundary conditions emanate from the ARPEGE (Action de Recherche Petite Échelle Grande Échelle) model. At the time of Flight 1, the operational AROME model was integrated with 46 levels, while at the time of Flight 2, additional vertical levels were incorporated, leading to a total of 60 levels.

Froude Number Definitions
The Froude number is a dimensionless parameter that may be utilized to describe the interaction mechanism between kinetic and potential energies as air flows past orography, according to Sutherland [28]. The Froude number can be interpreted as the ratio of the inertial to gravity forces in the flow (Lynch and Cassano [29]). This ratio may also be interpreted physically as the ratio between the mean flow velocity and the shallow water gravity wave speed (Holton [30]). The Froude number has also been applied to infer the type of flow phenomenon that may occur in the lee of a mountain when paired with a flow regime diagram, as was done by Sheridan and Vosper [31].
Three definitions of the Froude number were considered in this study. The first is the classical form, the second is based on Sheppard's dividing-streamline concept [32] and the third is applicable when a temperature inversion is present in the vicinity of the upstream region of the mountain. The exact upstream location at which Froude numbers were computed in both flights is indicated by point C in Figure 2. The classical Froude number, Fr, may be defined as: where U is upstream flow speed, N is the Brunt-Väisälä frequency and h m is the mountain height [6,7,33]. Vosper et al. [34] provide an additional definition, F L , that is based on a horizontal length scale L, such that F L = U/NL. In this study, h m is 1200 m, which is the maximum height of Pico Mountain according to AROME topography and L = 21 km. Some authors [2,4,35], alternatively, refer to the non-dimensional mountain height M as: where N would now represent an average value for the Brunt-Väisälä frequency below h m and U nm a cross-mountain wind speed average also under h m , as is done similarly by Jiang et al. [35]. Several authors summarize their results from laboratory, numerical and theoretical studies in a schematic regime diagram showing that for Fr > 1 mostly small-amplitude gravity waves are expected, while for Fr < 1, flow splitting, lee vortices and/or wave breaking may occur depending also on the horizontal aspect ratio of the obstacle. This is seen in Figure 2 from [4] and in Figure 13 from [2], the latter of which is shown here as Figure A1a.
For the parameter F L , whenever F L ≤ 1 vertically propagating waves are expected, but are evanescent if F L > 1 [34].

Froude Number (Fr h ) in Terms of the Dividing-Streamline Height h c
Under stable conditions and with varying wind speed and stratification, according to Sheppard [32] and Heinze et al. [6], it is possible to determine the dividing-streamline height h c (also called the critical height) as follows: where Heinze et al. [6] define V h as the mean horizontal velocity as a function of the horizontal velocity components u and v: V h = √ u 2 + v 2 . It is assumed in Equation (5) that all the kinetic energy of the flow is converted into potential energy [6]. As was done in [6], Equation (5) will be solved iteratively for h c via the trapezoidal rule to numerically compute the integral.
As noted by Heinze et al. [6], if Fr > 1, all fluid parcels, even near the surface, can flow over the obstacle. On the other hand, when Fr < 1, the fluid may stagnate in the upstream region below a certain critical height h c , passing primarily around the orographic obstacle rather than over it (flow splitting phenomenon). Above h c , the flow would be able to ascend the obstacle [6]. An illustration of flow splitting and of the critical height h c is provided in Figure 3.
So, when flow splitting occurs, an expression for the Froude number in terms of the dividing-streamline height may then be obtained, as is noted by Trombetti and Tampieri [36]: According to Wooldridge et al. [37], Equation (6) may also be used to find the critical height h c if Fr h is replaced by Fr (defined by Equation (3)), which is applicable when Fr < 1.
According to Wooldridge et al. [37], Equation (6) may also be used to find the critical height ℎ if Fr is replaced by Fr (defined by Equation (3)), which is applicable when Fr < 1.

Inversion Froude Number (Fr )
The inversion Froude number may be computed in the presence of an inversion by adopting Vosper's [38] definition: where is the inversion height, ∆ is the magnitude of the potential temperature difference across the inversion, is the mean potential temperature below the inversion base and is the cross-mountain wind speed (averaged below ). See Figure 3 for a visual representation of some of these parameters.
In the presence of a temperature inversion, several authors [5,38,39] provide a regime diagram for a two-layer flow over mountains as a function of the inversion Froude number (Equation (7)) and of the ratio , as is shown in Figure A1b from Appendix A.

Inversion Froude Number (Fr i )
The inversion Froude number may be computed in the presence of an inversion by adopting Vosper's [38] definition: where z i is the inversion height, ∆θ is the magnitude of the potential temperature difference across the inversion, θ m is the mean potential temperature below the inversion base and U nm is the cross-mountain wind speed (averaged below z i ). See Figure 3 for a visual representation of some of these parameters.
In the presence of a temperature inversion, several authors [5,38,39] provide a regime diagram for a two-layer flow over mountains as a function of the inversion Froude number (Equation (7)) and of the ratio h m z i , as is shown in Figure A1b from Appendix A.

Turbulence Indicators
Five turbulence indicators commonly used in aviation meteorology applications [40,41] are used in this study: Brown, Ellrod TI1, Ellrod TI2, CAT1 and EDR. These parameters were calculated using AROME forecasts.

Brown Index Φ
According to Gill and Buchanan [40] and Sharman et al. [41], the Brown Index Φ may be defined as: with u and v denoting the zonal and meridional components of the wind respectively. The vertical component of the absolute vorticity is represented by ζ a , which is the sum of the vertical component of the relative vorticity and the Coriolis frequency f . The form applied in this work, Φ ε , was based on Equation (8) and is expressed as energy dissipation in the following manner: where S V is the vertical wind shear.

Ellrod TI Indexes
As is mentioned in [40], the Ellrod indicators can be used to quantify turbulence. The first version, Ellrod TI1, is defined as: where S v is the vertical wind shear. The Ellrod TI2 index is similar to the Ellrod TI1, except that it incorporates a convergence term in the form of − ∂u ∂x + ∂v ∂y , as is shown in Equation (11):

CAT1 Indicator
The CAT1 turbulence indicator, as is denoted in this study, is referred to as "MOS (Model Output Statistics) probability forecast" by Sharman et al. [41] and is defined as: where V h is the horizontal wind speed.

EDR Indicator
The Eddy Dissipation Rate, ε, may also be used as a quantifier for turbulence. The method of calculation for this parameter in the current study followed the empirical formula delineated in Frech et al. [42]: where L e = 311 m, which is the chosen length scale [42]. TKE represents the turbulent kinetic energy, which is a prognostic variable from the AROME model. The EDR indicator itself, however, was defined as the cube root of ε according to ICAO [43]: This study presents EDR because it is the standard measure to quantify turbulence in aviation applications. Nevertheless, for other scientific communities, it is relevant to show TKE forecasts, which are presented in the supplementary materials.

RMSE and MAE for Wind Speed and Direction
The quantitative comparison between observed (from aircraft) and forecast (from AROME) wind data is established using the root mean square error (RMSE) and the mean absolute error (MAE). The calculation of the wind speed error, RMSE WSPD , is performed following the method employed by Grubišić et al. [5] for a dataset with n points, as follows: Atmosphere 2019, 10, 350 MAE WSPD is similarly defined as follows: The root mean square error and the mean absolute error for the wind direction, RMSE WDIR and MAE WDIR respectively, are determined using the method described in Jiménez et al. [44]: where ∆ WDIR (i) in degrees is given by: The comparison between model and observation was done according to the following two-step approach in order to find corresponding values for a given airborne observation: The absolute difference between all AROME levels for each recorded flight altitude, expressed as Altitude Flight − Altitude AROME , is computed and the level with the minimum difference is selected as the best match. This process is repeated for each available flight altitude.

2.
Once every flight altitude is paired with an AROME level, a latitude and longitude pair must be chosen at each one of these correspondences. Following a similar logic, the best coordinate pair from the model was deemed to minimize the distance relative to the flight coordinates at each altitude.

Results and Discussion
This section begins by illustrating results based on aerial observations. The wind shear intensity factor "I" is then analyzed. Next, outcomes based primarily on AROME forecasts like Froude numbers and turbulence indicators are discussed. Finally, the accuracy of AROME wind predictions is assessed by comparing them with data that was collected by SATA aircraft, using objective scores.

Aerodrome Wind Data
At the LPPI aerodrome, wind was predominantly from the southwest at two instants close to the time of both flights (Table 1). A headwind (W HD ) of 12/13 kt (6/7 ms −1 ) and a crosswind (W CH ) of 14/16 kt (7/8 ms −1 ) were registered on the day of Flight 1. In regard to gusts, the aircraft experienced a headwind (GW HD ) of 18/20 kt (9/10 ms −1 ) and a crosswind (GW CH ) of 21/24 kt (11/12 ms −1 ) incoming from the left.

Characterization of Wind Profiles
The wind speed profile based on airborne observations along its path (see Figure 2 and Figure S1) is presented in Figure 4, for Flights 1 and 2. During Flight 1, the aircraft experienced strong winds oscillating between 10 kt (5 ms −1 ) and nearly 40 kt (21 ms −1 ). At 440 ft (107 m), for instance, the aircraft experienced a wind of 38 kt (19.5 ms −1 ). During Flight 2, in the approach phase, the wind speed reaches its maximum of nearly 22 kt (11 ms −1 ) at 600 ft (183 m), at a distance of 1.67 NM of the aerodrome. Closer to the aerodrome, the wind decreases to 5 kt (3 ms −1 ) near the surface. Above 1000 ft (305 m) at a distance of nearly 4 NM of the aerodrome, the wind is weak (below 8 kt or 4 ms), due to the wake in the lee of Pico Mountain, as will be shown in the next sections. During both flights, the wind was predominantly from the southwest below 1500 ft or 457 m ( Figure S4). 3.1.3. Wind Shear Intensity Factor "I" Figure 5a,b detail the relative frequency distributions in both flights for altitudes only below 2100 ft (640 m) in order to focus on the effects of low-level wind shear as is suggested in [25]. During Flight 1, about 36% of the approach phase was carried out under "severe" wind shear conditions and over 15% under "strong". In other words, more than 50% of it experienced "strong" to "severe" wind shear, which is consistent with the occurrence of the hard landing incident. During Flight 2, "severe" wind shear conditions were absent, with only 4% being classified as "strong" and nearly 29% as "moderate" (Figure 5b). Table 2 lists wind shear levels of the first three lowest levels as functions of and ∆ for both flights. By comparing them, it is clear that within the three closest levels to the ground, the hard landing incident (Flight 1) exclusively experiences "severe" wind shear, whereas in Flight 2, a "light" intensity is prevalent. 3.1.3. Wind Shear Intensity Factor "I" Figure 5a,b detail the relative frequency distributions in both flights for altitudes only below 2100 ft (640 m) in order to focus on the effects of low-level wind shear as is suggested in [25]. During Flight 1, about 36% of the approach phase was carried out under "severe" wind shear conditions and over 15% under "strong". In other words, more than 50% of it experienced "strong" to "severe" wind shear, which is consistent with the occurrence of the hard landing incident. During Flight 2, "severe" wind shear conditions were absent, with only 4% being classified as "strong" and nearly 29% as "moderate" (Figure 5b). Table 2 lists wind shear levels of the first three lowest levels as functions of dV dt and ∆V V for both flights. By comparing them, it is clear that within the three closest levels to the ground, the hard landing incident (Flight 1) exclusively experiences "severe" wind shear, whereas in Flight 2, a "light" intensity is prevalent.
3.1.3. Wind Shear Intensity Factor "I" Figure 5a,b detail the relative frequency distributions in both flights for altitudes only below 2100 ft (640 m) in order to focus on the effects of low-level wind shear as is suggested in [25]. During Flight 1, about 36% of the approach phase was carried out under "severe" wind shear conditions and over 15% under "strong". In other words, more than 50% of it experienced "strong" to "severe" wind shear, which is consistent with the occurrence of the hard landing incident. During Flight 2, "severe" wind shear conditions were absent, with only 4% being classified as "strong" and nearly 29% as "moderate" (Figure 5b). Table 2 lists wind shear levels of the first three lowest levels as functions of and ∆ for both flights. By comparing them, it is clear that within the three closest levels to the ground, the hard landing incident (Flight 1) exclusively experiences "severe" wind shear, whereas in Flight 2, a "light" intensity is prevalent.

Synoptic Analysis
During Flight 1, the flow in Pico Island was influenced by the presence of a depression located northwest of the Azores and by a subtropical anticyclone to its west/southeast. Consequently, southwesterly winds of around 20 kt (10 ms −1 ) were visible south of the island, according to an ECMWF analysis (Figure 6a). For Flight 2, an anticyclone was present south of the Azores, causing southwesterly winds of approximately 10 kt (5 ms −1 ) south of Pico Island (Figure 6b).

Mesoscale Analysis
During Flight 1, according to AROME forecasts, the upstream conditions are characterized by a stable stratification, with N varying from 0.012 to 0.015 s −1 and by the absence of an inversion ( Figure  S5a). The Froude number varies between 0.7 and 0.9 below 400 m, reaching a maximum value of 1.2 at about 1600 m (Figure 7), which coincides with a wind speed of nearly 18 ms −1 ( Figure S5c). Above this level, Fr varies between 1 and 1.2. Consistently, Fr and Fr h are 0.95 and 0.9, respectively (Table 3). Therefore, according to previous studies [2,4], no wave breaking, upstream blocking or lee vortex formation are expected (see Figure A1a). The critical height defined in Section 2 is 107 m or 60 m, depending on the method. Consequently, it is expected that air will flow nearly horizontally around the mountain below this level, but will overcome the obstacle above it. Moreover, the parameter F L is always below 0.08 ( Figure S5c) and so, according to the analysis of F L and the Froude number, vertically propagating waves are expected [34].

Synoptic Analysis
During Flight 1, the flow in Pico Island was influenced by the presence of a depression located northwest of the Azores and by a subtropical anticyclone to its west/southeast. Consequently, southwesterly winds of around 20 kt (10 ms -1 ) were visible south of the island, according to an ECMWF analysis (Figure 6a). For Flight 2, an anticyclone was present south of the Azores, causing southwesterly winds of approximately 10 kt (5 ms -1 ) south of Pico Island (Figure 6

Mesoscale Analysis
During Flight 1, according to AROME forecasts, the upstream conditions are characterized by a stable stratification, with varying from 0.012 to 0.015 s -1 and by the absence of an inversion  Table 3. Froude number and related results by flight for Point C in Figure 2 with h c1 , h c2 , z i in meters.  (3)) 0.95 0.42 (Equation (6)) 0.90 0.31 (Equation (7)) ∞ 0.92 (Equation (5) (Table 3). According to the flow regime diagram presented in Figure A1b, flow splitting and lee vortex formation with reversed flow are expected.

Characterization of Orographic Flow during Flight 1
The horizontal wind and vertical component of the absolute vorticity predicted by the AROME model are depicted in Figure 8 for Flight 1. At 30 m, the flow is partially deflected laterally (mostly on the west side), which is denounced by the curvature of the wind vectors. This deflection is still present at 200 m but is barely visible at 300 m (not shown), while at 500 m ( Figure 8b) it is no longer discernible at all. This implies that the critical altitude (h c ) should lie somewhere between 200 and 300 m. From Table 3, it may be seen that h c1 = 107 m and h c2 = 60 m, revealing that the estimate based on Equation (5) is more accurate. However, it underestimates h c , possibly due to an insufficient vertical resolution of the AROME model.
Until 2800 ft (853 m), a pair of positive and negative vertical vorticities is visible over the mountain and on the leeward region, illustrated in Figure 8a for a height of 30 m. However, no signs of reversed flow or vortex formation are present, which is consistent with the results of Bauer et al. [2] and Epifanio [4].
The horizontal wind speed at 680 ft (207 m) shows a stagnation area on the upstream region of Pico Mountain and a maximum wind speed (>40 kt or 21 ms −1 ) over it with weak winds on the leeward side (Figure 9a). This pattern is coherent with previous studies like in Bauer et al. [2] with their Figure 2.a2. Near Pico Aerodrome, AROME forecast winds with speeds above 35 kt (18 ms −1 ). Ergo, during the approach to the aerodrome, the aircraft passes a region of elevated wind shear. Figure 9b presents the maximum and minimum vertical velocities in the layer from the surface up to 6000 ft (1829 m), illustrating a wave pattern over Pico and Faial Islands. Figure 10a  The horizontal wind and vertical component of the absolute vorticity predicted by the AROME model are depicted in Figure 8 for Flight 1. At 30 m, the flow is partially deflected laterally (mostly on the west side), which is denounced by the curvature of the wind vectors. This deflection is still present at 200 m but is barely visible at 300 m (not shown), while at 500 m (Figure 8b) it is no longer discernible at all. This implies that the critical altitude (ℎ ) should lie somewhere between 200 and 300 m. From Table 3, it may be seen that ℎ = 107 m and ℎ = 60 m, revealing that the estimate based on Equation (5) is more accurate. However, it underestimates ℎ , possibly due to an insufficient vertical resolution of the AROME model.
Until 2800 ft (853 m), a pair of positive and negative vertical vorticities is visible over the mountain and on the leeward region, illustrated in Figure 8a for a height of 30 m. However, no signs of reversed flow or vortex formation are present, which is consistent with the results of Bauer et al. [2] and Epifanio [4]. During Flight 1, the mean wind speed at point C (from Figure 2) is 12.76 ms-1 ( ) and the average value of the buoyancy frequency is 0.0129 s -1 ( ). Since the horizontal scale ( ) for Pico Mountain is of the order of 21 km, the condition for vertically propagating gravity waves ( > ) applies [1,29]. See also Figure 13.3 in Cassano and Lynch [29], which illustrates the alternating pattern between upward and downward motion associated with vertically propagating gravity waves.  During Flight 1, the mean wind speed at point C (from Figure 2) is 12.76 ms −1 (U) and the average value of the buoyancy frequency is 0.0129 s −1 (N). Since the horizontal scale (L) for Pico Mountain is of the order of 21 km, the condition for vertically propagating gravity waves (N 2 > U 2 k 2 ) applies [1,29]. See also Figure 13.3 in Cassano and Lynch [29], which illustrates the alternating pattern between upward and downward motion associated with vertically propagating gravity waves. 6000 ft (1829 m) for Flight 1. Both fields are from AROME forecasts. Shading in (a) represents the wind speed in kt and the pink dot identifies the approximate aircraft position at this level. Pink dashed line identifies the approximate aircraft trajectory. Shading in (b) represents intensity in ft/min, positive values indicate upward vertical velocity and negative values represent downward vertical velocity. The pink dots are the points C and D shown in Figure 2, which are the endpoints of the cross-section (dashed pink line) used in the next figures.

Characterization of Orographic Flow during Flight 2
The horizontal wind and vertical component of the absolute vorticity predicted by the AROME model at several heights are presented in Figure 11 for Flight 2. Figure 11a depicts a stagnation area in the windward direction of Pico Mountain and the occurrence of a flow splitting phenomenon in addition to a wake with reversed flow in the lee of the hill. This wake (at 40 m) consists of a pair of vortices, one with a clockwise rotation and negative absolute vorticity and another with a counter-clockwise rotation and positive vorticity. The wake structure exemplified in Figure 11a, which is visible below 600 m (not shown), resembles the results obtained by Smith and Grubišić [45] (p.3747) in their study of Hawaii's wake, which is described as comprising "two large, nearly steady vortices". A similar structure with a pair of counter-rotating vortices downstream of the obstacle was also documented by Smolarkiewicz and Rotunno [7] (see their Figure 1c) and by Epifanio [4]

Characterization of Orographic Flow during Flight 2
The horizontal wind and vertical component of the absolute vorticity predicted by the AROME model at several heights are presented in Figure 11 for Flight 2. Figure 11a depicts a stagnation area in the windward direction of Pico Mountain and the occurrence of a flow splitting phenomenon in addition to a wake with reversed flow in the lee of the hill. This wake (at 40 m) consists of a pair of vortices, one with a clockwise rotation and negative absolute vorticity and another with a counter-clockwise rotation and positive vorticity. The wake structure exemplified in Figure 11a, which is visible below 600 m (not shown), resembles the results obtained by Smith and Grubišić [45] (p. 3747) in their study of Hawaii's wake, which is described as comprising "two large, nearly steady vortices". A similar structure with a pair of counter-rotating vortices downstream of the obstacle was also documented by Smolarkiewicz and Rotunno [7] (see their Figure 1c) and by Epifanio [4] (their Figure  3c). Moreover, the regime flow diagram of Figure A1b in Appendix A, which was used by Schär and Smith [39] and Grubišić et al. [5], accurately estimates this situation in which a wake displays reversed flow.
At higher levels, between 600 m and 850 m (Figure 11b,c), no lee vortices are present and the dipole of positive and negative vorticities is weaker and shifted to the east relative to lower levels, mainly at 850 m. At 1100 m, the incoming flow is mainly from the west and overcomes the mountain (Figure 11d). Finally, the analysis of Figure 11 suggests that the critical height in Flight 2 should be approximately 850 m, which is close to the estimate of the dividing streamline concept (see Table 3).
The horizontal wind speed at 40 m is depicted in Figure 12a, which shows that near the aircraft position, the wind speed varies between 15 and 20 kt (8-10 ms −1 ). It is also visible from this figure that the aircraft during the approach phase (pink dashed line in Figure 12a) has passed through the wake. Figures 12a and 13 also show signs of a vertically propagating mountain wave, however, less energetic than the one seen on the day of Flight 1 (compared to Figure 10) and propagating to lower levels (up to 25000 ft or 7620 m), as is illustrated in Figure S6b. During Flight 2, the maximum vertical velocity on the leeward side of Pico Mountain (Figure 13b) is less than 150 ft/min (0.8 ms −1 ), while during Flight 1 the maximum vertical velocity was above 400 ft/min (2 ms −1 ). Upon propagation of this wave to the aerodrome (Figure 13a), it is seen that only a weak downdraft (less than 100 ft/min or 0.5 ms −1 ) lies atop it with no hints of updrafts along the flight path. In Figure 13b, the reversed flow associated with the lee vortices is also visible until about 2000 ft (610 m). displays reversed flow.
At higher levels, between 600 m and 850 m (Figure 11b,c), no lee vortices are present and the dipole of positive and negative vorticities is weaker and shifted to the east relative to lower levels, mainly at 850 m. At 1100 m, the incoming flow is mainly from the west and overcomes the mountain (Figure 11d). Finally, the analysis of Figure 11 suggests that the critical height in Flight 2 should be approximately 850 m, which is close to the estimate of the dividing streamline concept (see Table 3). The horizontal wind speed at 40 m is depicted in Figure 12a, which shows that near the aircraft position, the wind speed varies between 15 and 20 kt (8-10 ms -1 ). It is also visible from this figure that the aircraft during the approach phase (pink dashed line in Figure 12a) has passed through the wake. Figures 12a and 13 also show signs of a vertically propagating mountain wave, however, less energetic than the one seen on the day of Flight 1 (compared to Figure 10) and propagating to lower levels (up to 25000 ft or 7620 m), as is illustrated in Figure S6b. During Flight 2, the maximum vertical velocity on the leeward side of Pico Mountain (Figure 13b) is less than 150 ft/min (0.8 ms -1 ), while during Flight 1 the maximum vertical velocity was above 400 ft/min (2 ms -1 ). Upon propagation of this wave to the aerodrome (Figure 13a), it is seen that only a weak downdraft (less than 100 ft/min or 0.5 ms -1 ) lies atop it with no hints of updrafts along the flight path. In Figure 13b, the reversed flow associated with the lee vortices is also visible until about 2000 ft (610 m).

Turbulence Indicators
Turbulence indicators Ellrod TI2 and CAT1 at 66 m are presented in Figure 14a,b, for Flight 1. Both indicators show maximum values leeward of Pico mountain, with TI2 exceeding 18 × 10 s -2 and CAT1 exceeding 8 × 10 ms -2 in the proximity of the aerodrome where the incident occurred. The cross-section along the flight path shows that for a distance less than 5 NM (9.3 km) of the runway, between 300 and 1600 ft, TI2 and CAT1 reach values above 18 × 10 s -2 and 8 × 10 ms -2 , respectively (Figures 14 c,d). Comparing the indicators with the wind shear Intensity Factor "I" (Figure 15a,b), the aircraft experienced severe wind shear at these levels, suggesting that these values for TI2 and CAT1 could be tested as thresholds indicative of severe wind shear conditions. Ellrod TI1 (not shown in Figure 14) has a similar spatial pattern to TI2, but with lower values (Figure

Turbulence Indicators
Turbulence indicators Ellrod TI2 and CAT1 at 66 m are presented in Figure 14a,b, for Flight 1. Both indicators show maximum values leeward of Pico mountain, with TI2 exceeding 18 × 10 −6 s −2 and CAT1 exceeding 8 × 10 −3 ms −2 in the proximity of the aerodrome where the incident occurred. The cross-section along the flight path shows that for a distance less than 5 NM (9.3 km) of the runway, between 300 and 1600 ft, TI2 and CAT1 reach values above 18 × 10 −6 s −2 and 8 × 10 −3 ms −2 , respectively (Figure 14c,d). Comparing the indicators with the wind shear Intensity Factor "I" (Figure 15a,b), the aircraft experienced severe wind shear at these levels, suggesting that these values for TI2 and CAT1 could be tested as thresholds indicative of severe wind shear conditions. Ellrod TI1 (not shown in Figure 14) has a similar spatial pattern to TI2, but with lower values (Figure 15a). The Brown index is maximum over land near the LPPI aerodrome (not shown in Figure 14), exceeding 20 × 10 −9 s −3 near the surface and decreasing drastically with height, reaching a value of 13 × 10 −9 s −3 at 600 ft ( Figure 15c).
The EDR indicator shows a different pattern when compared to others. EDR shows higher values over land than over sea, depicting a maximum value above 0.4 m 2/3 s −1 near the top of Pico Mountain (Figure 16a). The higher values over land are consistent with higher roughness lengths (than over sea), which contribute to a higher mechanical production of turbulence [16]. At 66 m, EDR reaches a value between 0.35 and 0.4 m 2/3 s −1 near the aerodrome (corresponding to moderate turbulence). Figure 16b shows that along the flight path, EDR decreases with height, varying from 0.35-0.4 m 2/3 s −1 near the surface to 0.15-0.2 m 2/3 s −1 at 2000 ft (610 m). It is interesting to mention that the thresholds applied in Figure 16 are lower than those presented in ICAO [43] but are those proposed recently [46], after a study by Sharman et al. [47] that showed that the thresholds in reference [43] were overestimated. According to reference [46], turbulence is severe when EDR ≥ 0.45 m 2/3 s −1 and is moderate for EDR between 0.2 and 0.45 m 2/3 s −1 . Comparing EDR to the wind shear Intensity Factor "I" (Figure 15d), it is clear that in the presence of severe wind shear conditions, EDR has values between 0.2 and 0.25 m 2/3 s −1 above 700 ft. For heights under 200 ft (61 m), EDR is close to 0.3, which suggests an underestimation of the severe wind shear conditions. the thresholds applied in Figure 16 are lower than those presented in ICAO [43] but are those proposed recently [46], after a study by Sharman et al. [47] that showed that the thresholds in reference [43] were overestimated. According to reference [46], turbulence is severe when EDR ≥ 0.45 m 2/3 s -1 and is moderate for EDR between 0.2 and 0.45 m 2/3 s -1 . Comparing EDR to the wind shear Intensity Factor "I" (Figure 15d), it is clear that in the presence of severe wind shear conditions, EDR has values between 0.2 and 0.25 m 2/3 s -1 above 700 ft. For heights under 200 ft (61 m), EDR is close to 0.3, which suggests an underestimation of the severe wind shear conditions.    During Flight 2, the index Ellrod TI2 shows maximum values in the leeward region of the mountain and exhibits a value of 9 × 10 s -2 near the aircraft position (Figure 17a). The index CAT1 shows maximum values above 5 × 10ms -2 downstream of Pico Aerodrome, close to the flight path (Figure 17b). Along the flight path itself, Ellrod TI2 reaches values above 9 × 10 s -2 below 500 ft (152 m), while CAT1 attains values between 4 × 10ms -2 and 8 × 10 ms -2 below 900 ft or 274 m (Figures 17 (c), (d) and S7), which are predominantly "moderate" wind shear conditions. The TI1 indicator is very similar to TI2 and the Brown index decreases again drastically with height, exhibiting values below 4 ×10 -9 s -3 above 300 ft ( Figure S7).  Figure S7), which are predominantly "moderate" wind shear conditions. The TI1 indicator is very similar to TI2 and the Brown index decreases again drastically with height, exhibiting values below 4 ×10 −9 s −3 above 300 ft ( Figure S7). Analysis of EDR depicts the occurrence of "light" turbulence near the aerodrome for Flight 2 ( Figure 18). The "I" factor is mostly consistent with this categorization as wind shear conditions reach intensities that are predominantly "light" for heights between 200-500 ft (61-152 m). Nevertheless, in some instants the wind shear was classified as "moderate", while the EDR values above 400 ft are below 0.1 m 2/3 s -1 suggesting the absence of turbulence or wind shear ( Figure S7d).
This preliminary scrutiny of turbulence indicators suggests that CAT1, Ellrod TI2 and EDR could be applied together to arrive at a classification for the hazard level of turbulence. Moreover, preliminary thresholds for these indicators are proposed (Table A1), which should be subject to further validation, using more flight data.
As expected from their similarity in mathematical construction, both TKE and EDR evolve analogously during both flights (Figures S8a,b). In Flight 1, when the EDR indicator attains values of around 0.3 m 2/3 s -1 (coinciding with "severe" wind shear conditions according to the "I" factor), TKE has values close to 4 m 2 s -2 . During Flight 2, an EDR hovering near 0.1 m 2/3 s -1 translated to a TKE value of about 0.5 m 2 s -2 .
In order to certify that the turbulence was indeed predominantly greater for Flight 1 than for Flight 2, independently of the discrepant mean flow speeds in each case, the dimensionless turbulence intensity parameter ( = √ / ), defined by Stull [16], was also plotted ( Figure S9). It Analysis of EDR depicts the occurrence of "light" turbulence near the aerodrome for Flight 2 ( Figure 18). The "I" factor is mostly consistent with this categorization as wind shear conditions reach intensities that are predominantly "light" for heights between 200-500 ft (61-152 m). Nevertheless, in some instants the wind shear was classified as "moderate", while the EDR values above 400 ft are below 0.1 m 2/3 s −1 suggesting the absence of turbulence or wind shear ( Figure S7d).

Objective Verification of Wind Forecasts from AROME
The accuracy of AROME wind predictions is assessed by comparing them with data that was collected by SATA aircraft during Flights 1 and 2. Absolute errors of wind speed and direction are also mapped for each observation and forecast pair in Figure 19. The color coding for wind magnitude and its direction is based on ICAO guidelines [43] for operationally acceptable levels of forecast errors. For wind speed, an absolute value error of up to 5 kt (2.6 ms -1 ) was deemed operationally desirable. For the direction, this value is 20°. For Flight 1, there is a good agreement between observations and forecasts.
Values for the root mean square error (RMSE) and mean absolute error (MAE) regarding wind speed and direction are displayed in Table 4, showing that wind speed RMSE for Flight 1 is only slightly above operational recommendations (about 30% larger) and the RMSE for wind direction is 11.4°, well within the acceptable limits. The Flight 1 MAE for wind speed shows that the 2.6 ms -1 recommendation is exceeded by only 7% and the wind direction MAE is again within limits with a value of 8.1°.
During Flight 2, the RMSE and MAE for wind speed are 4.43 ms -1 and 3.77 ms -1 respectively, slightly over 30% larger than Flight 1 errors. For Flight 2 wind direction errors, the RMSE and MAE are around 83.3° and 65.4°. It is important to mention that when the wind speed is above 10 kt, the wind direction errors are inferior to 45° (see Figure S10). However, during Flight 2 at altitudes above 1200 ft, the wind is below 10 kt ( Figure S4b) and at these levels, large errors (>120°) are found in the wind direction (Figures 4b,S10). The limitations of NWP models in accurately forecasting the wind direction for weak winds was previously documented by Jiménez et al. [44], who showed that RMSE increases with decreasing wind speeds, reaching values over 70° for speeds under 2 kt (1 ms -1 ). It is also important to note that during Flight 2, a wake developed in the lee of Pico mountain in the vicinity of the flight path. In the wake region, the wind varies strongly over short distances (see Figure 12), implying that a small phase error may cause large errors in wind speed and/or wind direction. According to Ricard et al. [22], the effective resolution of the AROME model is about 9 times the grid spacing (i.e. 22 km), which may be insufficient to correctly resolve the wind variations within the lee wake. This preliminary scrutiny of turbulence indicators suggests that CAT1, Ellrod TI2 and EDR could be applied together to arrive at a classification for the hazard level of turbulence. Moreover, preliminary thresholds for these indicators are proposed (Table A1), which should be subject to further validation, using more flight data.
As expected from their similarity in mathematical construction, both TKE and EDR evolve analogously during both flights ( Figure S8a,b). In Flight 1, when the EDR indicator attains values of around 0.3 m 2/3 s −1 (coinciding with "severe" wind shear conditions according to the "I" factor), TKE has values close to 4 m 2 s −2 . During Flight 2, an EDR hovering near 0.1 m 2/3 s −1 translated to a TKE value of about 0.5 m 2 s −2 .
In order to certify that the turbulence was indeed predominantly greater for Flight 1 than for Flight 2, independently of the discrepant mean flow speeds in each case, the dimensionless turbulence intensity parameter (TI = √ TKE/V h ), defined by Stull [16], was also plotted ( Figure S9). It is clear from these results that the TI of Flight 1, with two maxima slightly above 0.25 (near 200 ft or 61 m and 1200 ft or 366 m), mostly exceeds that of Flight 2, with a maximum slightly above 0.15 (near 300 ft or 91 m).

Objective Verification of Wind Forecasts from AROME
The accuracy of AROME wind predictions is assessed by comparing them with data that was collected by SATA aircraft during Flights 1 and 2. Absolute errors of wind speed and direction are also mapped for each observation and forecast pair in Figure 19. The color coding for wind magnitude and its direction is based on ICAO guidelines [43] for operationally acceptable levels of forecast errors. For wind speed, an absolute value error of up to 5 kt (2.6 ms −1 ) was deemed operationally desirable.

Summary and Conclusions
The present article focuses on the orographic flow associated with Pico Mountain in the Azores in two different synoptic environments, using SATA airborne measurements and AROME forecasts. Aerial measurements were used to compute the wind shear intensity factor "I" according to Guan and Yong [25] and to categorize its severity according to four classes: "light", "moderate", "strong" and "severe".
During the first scenario, a hard landing incident (Flight 1) occurred in Pico's aerodrome. On this day, according to the AROME model, the upstream conditions were characterized by southwest winds with speeds of up to 40 kt (20 ms -1 ) in the ABL and stable stratification, which translated to a Froude number slightly above 1 near the mountain top. These conditions triggered the development of vertically propagating mountain waves, which created an undulation pattern of upward and downward motion downstream of Pico Mountain, leading to maximum vertical velocities above 400 ft/min (2 ms -1 ) and above 200 ft/min in the flight path. Based on ICAO Annex 3 guidelines [43], a mountain wave is deemed "severe" whenever downdrafts exceed 600 ft/min (3 ms -1 ) and "moderate" whenever they vary between 350 to 600 ft/min (1.75-3.0 ms -1 ). However, it is important to note that the AROME model may have underestimated the intensity of the vertical velocities, so this pattern of oscillating updraft/downdraft may have strongly contributed to the hard landing event. Furthermore, large wind speeds (38 kt or 20 ms -1 ) were registered close to 440 ft (107 m) and the wind shear intensity factor "I" based on airborne data showed that 36% of the approach phase (below 2100 ft or 640 m) occurred under "severe" wind shear conditions and 16% under "strong". At the aerodrome, METAR observations showed southwest winds (18-21 kt or 9-11 ms -1 ) with gusts between 28 and 31 kt (14 and 16 ms -1 ).
During the second flight, the upstream conditions were characterized by southwest winds with speeds below 20 kt (10 ms -1 ) in the PBL and a temperature inversion in the layer between 860 to 1100 m. In this range, the Froude number is less than 0.5. According to previous studies, these conditions promote upstream stagnation and flow splitting with the formation of lee vortices in low levels with reversed flow, which was predicted by the AROME model. AROME forecasts for Flight 2 also showed a pattern of vertically propagating waves, but with smaller vertical velocities and shorter Values for the root mean square error (RMSE) and mean absolute error (MAE) regarding wind speed and direction are displayed in Table 4, showing that wind speed RMSE for Flight 1 is only slightly above operational recommendations (about 30% larger) and the RMSE for wind direction is 11.4 • , well within the acceptable limits. The Flight 1 MAE for wind speed shows that the 2.6 ms −1 recommendation is exceeded by only 7% and the wind direction MAE is again within limits with a value of 8.1 • . It is important to mention that when the wind speed is above 10 kt, the wind direction errors are inferior to 45 • (see Figure S10). However, during Flight 2 at altitudes above 1200 ft, the wind is below 10 kt ( Figure S4b) and at these levels, large errors (>120 • ) are found in the wind direction (Figure 4b, Figure S10). The limitations of NWP models in accurately forecasting the wind direction for weak winds was previously documented by Jiménez et al. [44], who showed that RMSE WDIR increases with decreasing wind speeds, reaching values over 70 • for speeds under 2 kt (1 ms −1 ). It is also important to note that during Flight 2, a wake developed in the lee of Pico mountain in the vicinity of the flight path. In the wake region, the wind varies strongly over short distances (see Figure 12), implying that a small phase error may cause large errors in wind speed and/or wind direction. According to Ricard et al. [22], the effective resolution of the AROME model is about 9 times the grid spacing (i.e., 22 km), which may be insufficient to correctly resolve the wind variations within the lee wake.

Summary and Conclusions
The present article focuses on the orographic flow associated with Pico Mountain in the Azores in two different synoptic environments, using SATA airborne measurements and AROME forecasts. Aerial measurements were used to compute the wind shear intensity factor "I" according to Guan and Yong [25] and to categorize its severity according to four classes: "light", "moderate", "strong" and "severe".
During the first scenario, a hard landing incident (Flight 1) occurred in Pico's aerodrome. On this day, according to the AROME model, the upstream conditions were characterized by southwest winds with speeds of up to 40 kt (20 ms −1 ) in the ABL and stable stratification, which translated to a Froude number slightly above 1 near the mountain top. These conditions triggered the development of vertically propagating mountain waves, which created an undulation pattern of upward and downward motion downstream of Pico Mountain, leading to maximum vertical velocities above 400 ft/min (2 ms −1 ) and above 200 ft/min in the flight path. Based on ICAO Annex 3 guidelines [43], a mountain wave is deemed "severe" whenever downdrafts exceed 600 ft/min (3 ms −1 ) and "moderate" whenever they vary between 350 to 600 ft/min (1.75-3.0 ms −1 ). However, it is important to note that the AROME model may have underestimated the intensity of the vertical velocities, so this pattern of oscillating updraft/downdraft may have strongly contributed to the hard landing event. Furthermore, large wind speeds (38 kt or 20 ms −1 ) were registered close to 440 ft (107 m) and the wind shear intensity factor "I" based on airborne data showed that 36% of the approach phase (below 2100 ft or 640 m) occurred under "severe" wind shear conditions and 16% under "strong". At the aerodrome, METAR observations showed southwest winds (18-21 kt or 9-11 ms −1 ) with gusts between 28 and 31 kt (14 and 16 ms −1 ).
During the second flight, the upstream conditions were characterized by southwest winds with speeds below 20 kt (10 ms −1 ) in the PBL and a temperature inversion in the layer between 860 to 1100 m. In this range, the Froude number is less than 0.5. According to previous studies, these conditions promote upstream stagnation and flow splitting with the formation of lee vortices in low levels with reversed flow, which was predicted by the AROME model. AROME forecasts for Flight 2 also showed a pattern of vertically propagating waves, but with smaller vertical velocities and shorter vertical propagation compared to Flight 1. According to airborne observations, Flight 2 occurred under wind shear conditions predominantly rated as "light" (68%) or "moderate" (29%), with its most intense class being "strong" (only during less than 4% of the landing phase).
Moreover, the pattern of several turbulence indicators derived from AROME forecasts was analyzed for both flights. All these indicators forecast higher values for Flight 1 than for Flight 2, which is consistent with the wind shear intensity factor "I" derived from airborne observations. As a result of the comparison between "I" and the AROME turbulence indicators, preliminary thresholds have been proposed for three indexes. In the future, however, more flights should be perused in order to verify the validity of such thresholds. The use of other wind shear and turbulence indicators derived from airborne data should also be tested.
Lastly, this study provides an objective verification of AROME wind forecasts, showing a good agreement with airborne observations for wind speeds above 10 kt, but a poor skill for weaker winds. For Flight 1, characterized by wind speeds over 30 kt (15 ms −1 ) during several instants, AROME forecasts show a good agreement with observations in terms of wind speed and direction, with an RMSE of 3.4 ms −1 and 11.4 • , respectively. The MAE score is 2.8 ms −1 and 8.1 • for wind speed and direction, respectively. During Flight 2, in the approach phase, the maximum wind speed is 22 kt (11 ms −1 ). In this case, AROME has an RMSE of 4.4 ms −1 and 83.3 • for wind speed and wind direction. MAE revealed values of 3.8 ms −1 and 65.4 • for the wind speed and direction. In particular, when the wind speed is below 10 kt (5 ms −1 ), large errors (>120 • ) in the wind direction can be found. This is a limitation of NWP models that was previously documented by Jiménez et al. [44] for the WRF model.
Strong ageostrophic winds often develop in coastal gaps or channels when an along-gap pressure gradient is created [48]. A gap flow may also contribute to reinforce the winds downstream of the strait between Faial and Pico islands. This hypothesis should be subject to future research.
Supplementary Materials: The following are available online at http://www.mdpi.com/2073-4433/10/7/350/s1, Figure S1: (a) Enlarged view of Pico Mountain and its Aerodrome (LPPI), (b) Location of Pico Aerodrome (LPPI) relative to Pico Mountain, Figure S2: Cross and headwind components, Figure S3: Boundaries for the classification of wind shear intensity, Figure S4: Wind profile along flight path according to aerial observations and AROME for (a) Flight 1, and (b) Flight 2, Figure S5: Absolute temperature and N profiles at point C for (a) Flight 1 and (b) Flight 2, Figure S6: (a) Vertically propagating mountain waves predicted by AROME during Flight 1, and (b) Flight 2, Figure S7: AROME turbulence indicators vs. the wind shear intensity factor I along the flight path during Flight 2. Figure S8: Variation of Turbulent Kinetic Energy (TKE) and Eddy Dissipation Rate (EDR) with Height for (a): Flight 1 and (b) Flight 2. Figure S9: Turbulence Intensity for both Flights. Figure