Aviation Turbulence: Dynamics, Forecasting, and Response to Climate Change

Atmospheric turbulence is a major hazard in the aviation industry and can cause injuries to passengers and crew. Understanding the physical and dynamical generation mechanisms of turbulence aids with the development of new forecasting algorithms and, therefore, reduces the impact that it has on the aviation industry. The scope of this paper is to review the dynamics of aviation turbulence, its response to climate change, and current forecasting methods at the cruising altitude of aircraft. Aviation-affecting turbulence comes from three main sources: vertical wind shear instabilities, convection, and mountain waves. Understanding these features helps researchers to develop better turbulence diagnostics. Recent research suggests that turbulence will increase in frequency and strength with climate change, and therefore, turbulence forecasting may become more important in the future. The current methods of forecasting are unable to predict every turbulence event, and research is ongoing to find the best solution to this problem by combining turbulence predictors and using ensemble forecasts to increase skill. The skill of operational turbulence forecasts has increased steadily over recent decades, mirroring improvements in our understanding. However, more work is needed—ideally in collaboration with the aviation industry—to improve observations and increase forecast skill, to help maintain and enhance aviation safety standards in the future.


Introduction
Atmospheric turbulence is a major aviation hazard, causing damage to aircraft and injury to passengers and crew. Turbulence is part of the chaotic atmosphere, and the chaotic nature poses a great challenge in understanding and forecasting turbulence. Aviation turbulence is incompletely understood and difficult to forecast, making it a significant hazard. The USA National Transportation Safety Board (NTSB) records the average number of air-carrier-related injuries as 58 per year (FAA 2017b). However, Sharman et al. (2006) suggest that this number is an underestimate, as not all injuries are reported. They state that in the period 1980-2008, there were 234 turbulence accidents, resulting in 298 serious injuries and three fatalities on United States operated air carriers. 184 of the serious injuries involved flight attendants and 114 involved passengers (FAA 2017b). These turbulence injuries will come at a cost to the airlines in two ways: (1) through compensation being paid, which amounts to over $10 million per year, and (2) through lost working days by injured cabin crew, which is over 7000 days per year (Sharman and Lane 2016). Although most of the flight is spent at the cruise phase (around 33,000-39,000 ft), this is also the most vulnerable part of the flight as passengers and crew are unbuckled, so encountering any turbulence would be more likely to result in injury. As a result, most incidents occur above 10,000 ft (Sharman et al. 2006). Kauffmann (2002) analysed, in detail, the impact that turbulence has on the aviation industry. He presents figures from the USA Federal Aviation Administration (FAA), who state that there were 342 reports of turbulence affecting flights on major air carriers over the period 1981-1997, with three fatalities, 80 serious injuries, and 769 minor injuries. These figures are similar to those of Sharman et al. (2006) and, again, may underestimate the real statistics of turbulence injuries. These values could also be higher in the future, as climate change is likely to increase the frequency of turbulence around the world, particularly in the mid-latitudes (Williams and Joshi 2013;Williams 2017;Storer et al. 2017). The cost of turbulence to the aviation industry is significant and comes from many sources, one of which is preventing aircraft from flying on the optimum route. This cost is as much as $16 million a year, as it is estimated that 5% of flights are forced to fly non-optimal routes (Search Technology 2000). More importantly, it is also suggested that 15% of the diversions could have been avoided with improved turbulence detection. One form of turbulence detection could be using LIght Detection And Ranging (LIDAR) technology on the front of an aircraft. The current RADAR technology used by aircraft is unable to detect turbulence unless hydrometeors are present, whereas for clear-air turbulence, they are absent. However, LIDAR technology is able to sense clearair turbulence using non-hydrometeor particles. This technology could warn pilots of turbulence ahead, enabling them to divert or put on the seatbelt sign, so passengers and crew are aware of the upcoming danger. At the moment, however, Kauffmann (2002) calculates that LIDAR technology costs more to install than the saving would be, and therefore, rolling it out is not worth the investment. Satellites can be used to avoid some kinds of turbulence, particularly turbulence associated with convection (Mecikalski et al. 2007;Francis and Batstone 2013). There is also the potential for the new generation of satellites (e.g., GOES-16) with higher spatial and temporal resolutions to improve the avoidance of turbulence. This improvement could come from better identification of deep convection and the ability to resolve gravity waves that would otherwise be invisible to the older generation GOES satellites (Wimmers et al. 2018;Nunez 2018).
Forecasting turbulence is another possibility to mitigate injuries and damage. Turbulent eddies in the atmosphere occur on scales ranging from the planetary scale down to millimeters, but only eddies of approximately 100 m in size impact aviation (Sharman et al. 2006). This is a scale that is not explicitly simulated in forecasts, because the numerical models have resolutions too coarse to resolve the individual eddies. Therefore, turbulence diagnostics are used. These diagnostics generally work on the principle that the energy associated with turbulence on aviation-affecting scales cascades down from the larger scales that can be explicitly resolved by numerical models.
The different kinds of turbulence that impact aviation are outlined in Fig. 1. Vertical wind shear instabilities, mountain waves, and convection are the three main sources that we consider in this study. Turbulence in and near clouds can also cause injuries to passengers and crew, but can easily be detected visually by pilots and using on-board RADAR. Boundary-layer turbulence is not considered here, because it influences only a small portion of the flight after take-off and before landing.
This review paper examines the current understanding of the different sources of aviation turbulence (Sect. 2), how turbulence might change in the future because of climate change (Sect. 3), and the current methods of turbulence forecasting (Sect. 4). The paper concludes with a summary and discussion (Sect. 5).

Clear-Air Turbulence
Clear-air turbulence (CAT) is defined as highaltitude aircraft bumpiness in regions devoid of significant cloudiness and away from thunderstorm activity (Chambers 1955). Far from mountains, CAT is generally accepted to result from shear instabilities. Wind shear is, therefore, a major source of CAT and is one of the best understood sources. Figure 1 indicates this type of turbulence and its association with the jet stream. To understand why shear causes turbulence, we must define the Richardson number: where N 2 is the Brunt-Väisälä frequency squared, U is horizontal wind speed, z is altitude, g is gravitational acceleration, and h is potential temperature. The Richardson number is a nondimensional number with the numerator representing the stratification and the denominator representing the vertical wind shear. It follows from theoretical considerations that instability occurs when Ri \ 0.25, so instability is favoured by large vertical wind shear (denominator) and weak stratification (numerator). In computational calculations using gridded data, numerical models rarely reach Ri = 0.25 due to the coarse resolutions, and therefore, thresholds of turbulence are model specific. To overcome this, Williams (2017) and Storer et al. (2017) chose thresholds based on the distribution of turbulence in the atmosphere. For example, they assume that severe turbulence is found in 0.01% of the atmosphere, and therefore, they take the top 0.01% (99.9-100%) of the probability distribution to be severe turbulence. Therefore, each threshold is specific to each model and resolution. It is possible for turbulence to be produced when the environmental Richardson number is much larger than the theoretical critical value, if a local effect reduces the Richardson number locally. For example, gravity waves can cause CAT by reducing the Richardson number (Ri) locally in an environment that would not normally produce turbulence, initiating the Kelvin-Helmholtz shear instability and leading ultimately to overturning and breaking billows. The various sources of gravity waves are discussed by Williams et al. (2003Williams et al. ( , 2005Williams et al. ( , 2008. In particular, gravity waves can be produced by convection (as discussed in Sect. 2.2) and spontaneous loss of geostrophic balance as the flow evolves, as described by the Lighthill-Ford theory (Lighthill 1952;Ford 1994;Knox et al. 2008;McCann et al. 2012). There is a direct connection in certain circumstances that links deformation to the generation of inertia-gravity waves (Knox et al. 2008). In some cases, this could explain the success of deformation-based CAT diagnostics such as Ellrod and Knapp (1992) TI1 and TI2 (Knox 1997).
It is also possible for the environmental Richardson number to be reduced on a much larger scale, in regions of the atmosphere with particularly strong wind shear such as the jet streams. Strong vertical wind shears around the jet stream increase the denominator in Eq. (1) which, therefore, decreases the Richardson number until it reaches a critical value and turbulence is produced. Therefore, understanding the behaviour of the jet stream will help researchers and forecasters to understand how CAT may change. This is one of the reasons for the success in forecasting turbulence of CAT diagnostics containing vertical wind shear [e.g., Colson-Panofsky index (Colson and Panofsky 1965), TI1, TI2 (Ellrod and Knapp 1992)]. Kim et al. (2016) studied the impact that the North Atlantic Oscillation (NAO) has on aviation turbulence. The NAO is a measure of the relative strength of the Icelandic low and the Azores high. The positive NAO phase implies a stronger than normal Icelandic low and Azores high, so there is a strong pressure difference between the two. In the negative NAO phase, the opposite is true and the pressure difference is weaker. In the positive NAO phase, we see a stronger jet stream and a more northerly path, whereas the negative NAO phase is associated with a weaker jet stream and a more southerly path. This change to the jet stream, therefore, has an impact on turbulence for trans-Atlantic flights. Kim et al. (2016) used wind-optimal routes to find the fastest possible flight path between London (LHR) and New York (JFK). The study found that eastbound flights fly more frequently through regions of CAT than westbound flights, and, therefore, experience more turbulence in both the positive NAO phase and the negative NAO phase. The reason for this is that eastbound flights utilise the jet stream more to benefit from strong tailwinds, so the wind-optimal routes would fly in stronger vertical shear regions more for eastbound flights than westbound flights. Westbound wind-optimal routes avoid the strong headwinds of the jet stream, and would also avoid the stronger vertical shear associated with the jet stream and, therefore, encounter less turbulence.
According to Kim et al. (2016), in the positive NAO phase, westbound flights experience more moderate-or-greater (MOG) CAT than in the negative NAO phase, because some of the westbound flights detour northward to be on the cyclonic shear side of the northerly shifted jet stream, which is more susceptible to MOG turbulence. In contrast, eastbound flights in the negative NOA phase fly through the cyclonic shear side of the southerly shifted jet. Therefore, the study suggests that westbound flights are more prone to MOG CAT in the positive NAO phase, whereas eastbound flights are more prone to MOG CAT in the negative NAO phase. This information is important for flight planners, as for example in a positive NAO phase, an eastbound flight would fly on the anticyclonic shear side of the jet stream, still using the tailwind to reduce flight time. However, a westbound flight that would normally fly north to avoid the strong headwinds should fly south and avoid the cyclonic shear side of the jet stream, reducing flight time but also the chances of strong turbulence. Information such as this can be used to avoid the strongest turbulence, while still attempting to reduce flight times and, therefore, fuel consumption.
Another example of jet stream behaviour was studied by Trier et al. (2012), who discuss how moist convection influences the upper level jet. This topic has been studied by Trier and Sharman (2009) for warm Mesoscale Convective Systems (MCSs). Latent heat release perturbs an anticyclone, and this mechanism accounts for nearly all of the magnitude of the upper level jet. This is a similar mechanism found in a cold weather outbreak, so a mid-latitude cyclone (like the MCS) enhances the downstream anticyclone. The study suggests that, although not as dominant as MCSs, mid-latitude cyclones account for 30-50% of the strength of the southerly jet stream. Without the moist convection and the perturbation of the anticyclone, the wind shear, and therefore, resultant CAT would not be as strong. This is different to CIT, because it is the indirect effect of convection on the jet stream and the increased shear that causes turbulence. The other mechanism is associated with strong convection in the mid-latitude cyclone. With cloud tops below the flight level, the cloud tops generate gravity waves induced by the convection. This mechanism is discussed further in Sect. 2.2.

Convective Turbulence
Understanding the relationship between buoyancy and shear is very important in understanding where and why turbulence forms. Lane et al. (2012) explore our current understanding of near-cloud turbulence or Convectively Induced Turbulence (CIT). They explain that the FAA regulations at the time were not sufficient for avoiding severe turbulence. For example, guideline 5 states: ''Do avoid by at least 20 miles (laterally) any thunderstorm identified as severe or giving an intense radar echo, especially under the anvil of a large cumulonimbus''. Guideline 6 states: ''Do clear the top of a known or suspected severe thunderstorm by at least 1000 ft altitude for each 10 knots of wind speed at the cloud top''.
A possible source of CIT is unstable upper tropospheric thunderstorm outflow, similar to that described in Trier and Sharman (2009). Those authors proposed a mechanism for turbulence formation after completing simulations using a convectionpermitting model. Their proposed mechanism is the formation of strong vertical wind shear in the outflow regions of Mesoscale Convective Systems (MCS), which we know reduces Ri and leads to Kelvin-Helmholtz instability, which is a well-known source of turbulence. In the simulations, they also found that strong vertical wind shear created regions of strong static instability. This was caused by differentially advecting equivalent potential temperature gradients, which were influenced by the adiabatic cooling in the convective updrafts.
The above static instabilities are different from the traditional CIT mechanisms that generally result from high-frequency gravity wave breaking (Lane et al. 2003;Lane and Sharman 2008) or from reductions in Ri as a result of propagating gravity waves . Propagating gravity waves induced by convection are indicated as a source of turbulence in Fig. 1. Lane et al. (2003) studied the traditional mechanisms forming turbulence above deep convection, and found that gravity waves formed when the overshooting top returns back down. As these waves propagate up, the changes in vertical wind shear and buoyancy with height can change the wavelength of gravity waves and the tilt. The phase speed remains the same above the jet as within the jet based on linear theory, and the decreasing speed above the jet could reach the critical level ). This is found when UðzÞ À c ¼ 0, where U(z) is vertical wind shear and c ¼ x=k, with x being the wave frequency and k being the horizontal wave number. Therefore, the gravity wave can break if wave frequency is low, or wind shear is high.
Since wind shear is an important part of the gravity wave breaking, the jet stream plays an important role. The smaller the jet width, the larger the wind shear, and therefore, the closer to the cloud top the gravity wave breaking would be. Lane et al. (2012) also showed turbulence can be found in cirrus bands (banding caused by thermal instabilities), and also ducted gravity waves that propagate far from the updraft, reducing the Richardson number and causing turbulence much farther away than the 20 km outlined in the guidelines. The new guidelines from FAA (2017a) no longer reference the acceptable height above the thunderstorm pilots can fly. Guideline 2 now states: ''Don't attempt to fly under a thunderstorm even if you can see through to the other side. Turbulence and wind shear under the storm could be hazardous''. Guideline 3 states ''Don't attempt to fly under the anvil of a thunderstorm. There is a potential for severe and extreme clear-air turbulence''. Guideline 14 states ''Do avoid by at least 20 miles any thunderstorm identified as severe or giving an intense radar echo. This is especially true under the anvil of a large cumulonimbus''.

Mountain-Wave Turbulence
Mountain-wave turbulence is similar to Convectively Induced Turbulence, as gravity waves produced by the terrain (instead of convection) propagate and then break when a critical vertical wind shear value is reached (Nastrom and Fritts 1992;Wurtele et al. 1996). This is also indicated in Fig. 1 and shows its similarity to CIT. A study by Epifanio and Qian (2008) used a high-resolution ensemble to study turbulence produced by breaking mountain waves and found that their results were consistent with the previous work by Fritts et al. (1994Fritts et al. ( , 1996, suggesting that shear is an important factor in the breaking of gravity waves. Fritts et al. (1996) ran simulations to understand the turbulence mechanism in mountain waves, and found that the dominant source of instability, and therefore turbulence, was the wind shear. This came from the mean wind field, and the differential vertical advection of the mean shear by the wave field. Similar to convective turbulence, the vertical wind shear and its interaction with propagating gravity waves is an important Vol. 176, (2019) Aviation Turbulence: Dynamics, Forecasting, and Response to Climate Change 2085 mechanism in turbulence production. This understanding will help forecasters in the future. Wolff and Sharman (2008) show that MOG turbulence is typically found where topographic heights are above 1.5 km. They also analysed the low-level wind direction when mountain-wave turbulence is produced, and concluded that the wind direction must be within 30 of the perpendicular to the mountain range orientation. For north America, this makes the Rocky Mountains a prime location for MOG MWT, as the mountain range is oriented northsouth and frequently experiences westerly low-level flow with mountain heights above 1.5 km. Greenland is also an area with high mountain ranges with the potential to cause turbulence for trans-Atlantic flights.

Climatology and Response to Climate Change
Jaeger and Sprenger (2007) used reanalysis data to understand upper tropospheric clear-air turbulence in the Northern Hemisphere, because of its role in stratosphere-troposphere exchange (Shapiro 1976(Shapiro , 1978(Shapiro , 1980. The authors create a climatology of the Ellrod and Knapp (1992) index (TI1), Brunt-Väisälä frequency squared (N 2 ), Richardson number (Ri), and potential vorticity (PV). Negative potential vorticity is a turbulence predictor as it is related to isentropic inertial instability. The study found that TI1 is largest north of the jet stream, and Kelvin-Helmholtz instability indicated by Ri is near the jets, which is what we would expect with larger wind shear in that region. Symmetric instability is most frequent south of the jets, and is particularly associated with anticyclonic jets. Hydrostatic instability is only slightly dependent on the jet position, and is most common over land where convection and gravity wave activity are most prevalent (mountain wave and deep convection gravity waves). The study also showed that winter (DJF) has the highest turbulence frequency, which follows the understanding that the jet stream is stronger in winter. The study also found long-term trends in the frequency of diagnosed turbulence over the reanalysis period. Figure 2 shows that over the reanalysis period, there is an increasing trend of turbulence in the northern Hemisphere. They were also able to relate inter-annual variability to the North Atlantic Oscillation (NAO) and Pacific/North American flow, which we know from Sect. 2.1 influences CAT.
We know that the climate system is changing due to anthropogenic forcing, and these changes may have an impact on turbulence in the future. Collins et al. (2013) showed that the changes in temperature are not uniform around the world. Importantly for turbulence, the upper troposphere and lower stratosphere respond differently to anthropogenic forcing. The tropical upper troposphere is predicted to warm faster than the tropical surface, due to an increase in latent heat release. Latent heat is released during convection, and in a warmer climate the atmosphere can hold more moisture, and as a result more convection and subsequent latent heat release will warm the troposphere. The lower stratosphere, however, will cool with the increasing greenhouse gases (Fels et al. 1980). This cooling is related to changes in infrared radiation.
The upper tropospheric and lower stratospheric changes discussed above lead to an increase in the equator-to-pole temperature gradient at flight cruising levels, as shown in Fig. 3. Particularly in the RCP8.5 simulation (right), where greenhouse gas emissions are highest, we see the warming in the tropics, and the cooling at the poles, increasing the equator-topole temperature gradient. This increase in meridional temperature gradient will induce a thermal wind response, resulting in an increase in vertical wind shear and, therefore, turbulence in the mid-latitudes. These changes to the jet stream are shown in Fig. 4, which is taken from Delcambre et al. (2013). The changes to the jet stream are predominant in the midlatitudes and at airline cruise altitudes, making their impact the largest in the busiest flight regions around the world. In addition to modifying turbulence, the increased wind speeds are also expected to modify flight times (Williams 2016). Another impact of anthropogenic forcing is the 20th-century release of chlorofluorocarbons (CFCs), which destroy atmospheric ozone (O 3 ). The loss of ozone reduces the lower stratospheric temperature, further increasing the equator-to-pole temperature gradient. This effect enhances the predicted change arising from carbon dioxide. Several studies have reported that climate change will act to increase clear-air turbulence in the future, according to climate model simulations (Williams and Joshi 2013;Williams 2017;Storer et al. 2017). The first study to look at this (Williams and Joshi 2013) focused on north Atlantic moderate-or-greater turbulence and showed that it would increase in frequency with climate change by around 40-170%. Williams (2017) then furthered the study to see how climate change might influence turbulence in five strength categories from light to severe, finding that all would increase in frequency with climate change. Storer et al. (2017) then extended the work further to see how climate change will impact CAT globally, through all four seasons, and also at two flight levels in eight geographic regions. They found that CAT will increase globally in all four seasons. This is shown in Fig. 5, with the largest and most statistically significant increases confined to the mid-latitudes around the jet stream. Some regions locally will see increases of several hundred per cent for example at 200 hPa (39,000 ft.) the north Atlantic will see an increase of 143.3% in moderate turbulence and 181.4% in severe turbulence. The study also showed that at 200 hPa (1) in winter, severe CAT by 2050-2080 will be as common as moderate CAT in the control period, and (2) for a range of turbulence strengths from light to moderate-to-severe, summertime CAT by 2050-2080 will be as common as wintertime CAT in the control period. Aircraft manufacturers and forecasters need to be prepared to prevent an increase in damage to aircraft and injuries to passengers and crew. Currently, there is a gap in the literature looking at the response of Convectively Induced Turbulence and Mountain-Wave Turbulence to climate change. However, we know that deep convection can produce turbulence, and therefore, it is likely that CIT would increase if convection was increased. Price and Rind (1994) showed an increase of 5-6% of global lightning activity with every 1 C of global warming. Lightning activity has been used to characterise turbulence events (Gill and Stirling 2013;Meneguz et al. 2016), and therefore, if we see an increase in lightning, we would expect to see an increase in convection and possibly an increase in CIT as a result. Price and Rind (1994) show that lightning activity and convection will increase particularly around the tropics, which is the area where CAT is not expected to increase (Fig. 5). Therefore, combining CAT and CIT, we would see a global increase in all aviation turbulence. This trend of increasing lightning with climate change is supported by Reeve and Toumi (1999) who show an increase in lightning activity of 40% for every 1 K of average land wet-bulb temperature. It could be possible to study the change in CIT to climate change. Convective precipitation accumulation was used by Gill and Stirling (2013) as a convective indicator for turbulence forecasts, and therefore, assessing the changes of this quantity in climate models could indicate how climate change might impact CIT.

Forecasting
Currently, the World Area Forecast Centres (WAFCs) in London (Met Office) and Washington (NOAA) produce operational turbulence forecasts for aviation. This includes a t þ 24 h significant weather (SIGWX) chart four times a day (0000, 0006, 0012, 0018 UTC) which displays multiple aviation hazards including icing, CAT, as well as the location of convective cloud. The WAFCs also produce a gridded turbulence forecast four times a day at five pressure levels with a lead time t þ 6 h to t þ 36 h (Gill 2014).

Clear-Air Turbulence
Using diagnostics to forecast clear-air turbulence allows flight planners to avoid turbulent regions. Fahey (1993) showed that strategic planning can lead to not only a reduction in injuries, but also costs. The use of diagnostics is also the only way to operationally forecast turbulence, as the turbulent eddies that affect aviation are smaller than the resolution of global atmospheric models. Turbulence predictors that have a deformation term are particularly good at forecasting CAT. Knox et al. (2008) tried to understand the previous work that linked deformation to Kelvin-Helmholtz instabilities. This included frontogenesis that could initiate Kelvin-Helmholtz instabilities, or mesoscale waves that would then break down and form turbulence (Mancuso and Endlich 1966). They found that neither could fully explain the relationship between deformation and turbulence, so instead found a link between deformation and inertia-gravity wave generation via the Lighthill-Ford theory. In certain atmospheric environments close to shear instability (low Richardson number), gravity waves emitted by imbalance destabilize the atmosphere, locally reducing Ri to below the critical value of 0.25 and leading to Kelvin-Helmholtz instabilities, generating waves which break down and form turbulence (Miles and Howard 1964;Dutton and Panofsky 1970). This explains the success of empirical diagnostics such as Turbulence Index 1 (TI1) from Ellrod and Knapp (1992): where u is the eastward wind speed and v is the westward wind speed. TI1 found in equation (2) is one of the best shear turbulence forecast diagnostics, and has been found to forecast up to 75% of all CAT cases. This is why a deterministic TI1 forecast is routinely used by the WAFCs. McCann et al. (2012) furthered the work by Knox et al. (2008) to try and improve the forecast based on spontaneous imbalance. They made an algorithm that could be used operationally by improving the forecast below FL200 (which was a problem with the current diagnostics). The apparent high bias in turbulent kinetic energy dissipation above the tropopause is reduced by including the turbulent kinetic energy production from the environment and the locally produced turbulent kinetic energy. McCann et al. (2012) also converted the turbulent kinetic energy into an eddy dissipation rate (EDR), which is now a standard measure of turbulence in the atmosphere as it does not depend on aircraft size.
There are many turbulence diagnostics, each with its own its strengths and weaknesses. All forecast CAT and are unable to forecast convective turbulence or mountain-wave turbulence. Since each individual diagnostic is unable to forecast every CAT event, Sharman et al. (2006) generated the Graphical Turbulence Guidance (GTG) system for forecasting turbulence. They selected multiple turbulence diagnostics and then used PIlot REPorts (PIREPS) to understand the diagnostics' performances. A weighting system can then be added depending on the skill of each diagnostic to produce the best overall forecast possible. The study showed that this was an improvement over the single turbulence diagnostic, and is, therefore, worth doing. Over the years, this forecast has been improved, and more recently, a mountainwave turbulence predictor has been added to the GTG system. At the moment CIT is not forecast by the GTG system, but it has been shown that it is a useful addition (Gill and Stirling 2013).

Convective Turbulence
Convection is one of the main turbulence sources. It is not yet forecast using numerical models operationally, but, instead, forecasters manually put the location of convective clouds on the significant weather (SIGWX) charts used by pilots and flight planners. Gill and Stirling (2013) showed that using a convective diagnostic from numerical weather prediction output can forecast many convective turbulence events. Therefore, combining the shear turbulence predictors and the convective predictors (similar to GTG) offers the promise to greatly improve the forecast skill. The convective indicators used in the study included convective precipitation rate, convective precipitation accumulation, and Convective Available Potential Energy (CAPE). Gill and Stirling (2013) concluded that it is possible to improve turbulence forecasts by not only combining different predictors for the same turbulence source (e.g., TI1, Brown index, etc.) but also combining predictors for different turbulence sources as well (e.g., clear-air turbulence and convective turbulence). They suggest that more work is needed in this area, but there is the potential to greatly increase the forecast skill by adding in these extra turbulence diagnostics and forecasting a greater proportion of turbulence events.

Mountain-Wave Turbulence
An algorithm for forecasting mountain-wave turbulence was described by Turner (1999). Mountain waves are a special case of gravity waves and they are smaller than the numerical resolution of the models, so the wave stress has to be parameterised. To do this, the surface winds and stability are used to work out a wind stress, and then, the wind stress is used to work out the vertical profile. The wind stress is computed for every grid point, taking into account the abruptness of orography (as that is the location of the mountains). The surface stress is passed upwards from the surface (unless a hydraulic jump, critical level, or saturation stress is diagnosed), and therefore, less stress is transferred to the next level. Above the maximum point that gravity wave stress can be sustained, the waves break and thus turbulence forms, so the models can use the wave stress at each model level to calculate the location and height at which the gravity wave stress is too high. This is, therefore, the turbulence location in the model. For mountain-wave formation, the higher the gravity wave propagates, the more likely it is to break. This is because the maximum gravity wave stress which the atmosphere can hold decreases with increased wind shear and lower air density. Mountain-wave turbulence is found to be associated most strongly with the mid-latitude westerlies and mountains perpendicular to the flow. The new GTG3 system includes mountain-wave turbulence diagnostics alongside the typical CAT diagnostics (Sharman and Pearson 2017). These authors have also used a different approach to forecasting MWT and use a combination of lowlevel wind combined with CAT diagnostics. The authors produced 14 new MWT diagnostics, which showed an increased forecast skill in MWT prone areas, supporting their use in the new GTG system.

Ensemble Forecasting
One of the main problems with the turbulence forecasting methods outlined in Sects. 4.1,4.2,and 4.3 is that they do not convey uncertainty. To resolve this issue, Gill and Buchanan (2014) and Buchanan (2016) trialled the use of ensemble forecasting for aviation turbulence. An ensemble is a collection of forecast runs, each of which is considered to be equally likely. Using an ensemble, the uncertainty in the forecast can be found. In these studies, the Met Office Global and Regional Ensemble Prediction System (MOGREPS) was used, which was made operational in 2008 (Bowler et al. 2008). Multiple Vol. 176, (2019) Aviation Turbulence: Dynamics, Forecasting, and Response to Climate Change turbulence predictors were considered, including the Dutton (1980) index, the Ellrod and Knapp (1992) indices, the Brown (1973) index, the Lunnon index (Roach and Dixon 1970), the Richardson number, and two convective predictors. Using the Derived Equivalent Vertical Gust (DEVG) as a truth in the verification, the studies found that using the ensemble forecast, they increased the forecast skill for most of the turbulence diagnostics. The studies by Gill and Buchanan (2014) and Buchanan (2016) then combined the predictors in a manner similar to Sharman et al. (2006), using an iterative scheme to maximise the forecast skill. Again, the studies showed that the ensemble forecast was more skillful than a single model deterministic forecast. A probabilistic forecast would also be more useful for pilots and flight planners, as understanding the uncertainty of the forecast can help them to make the best decision possible to maintain safety of passengers and crew, while still flying the optimum routes and reducing flight times and fuel consumption. Further work is needed in ensemble forecasting before it can be used operationally, such as studying whether multimodel ensembles can improve the forecast. In addition, more research is needed to include a mountain-wave turbulence predictor alongside the convective and shear turbulence predictors. However, it is clear that ensembles can be beneficial and should eventually be used operationally in the future.

Forecast Verification
Turbulence forecast verification is difficult, because objective verification data sets are limited. Some previous studies resolved this issue using PIREPs (Kim and Chun 2011;Schwartz 1996) to identify regions of turbulence based on a semi-quantitative scale from light to extreme, but these can be unreliable (Kane et al. 1998). PIREPs are subjective, in the sense that a more experienced pilot may catagorise an event as moderate, but an inexperienced pilot may record it as severe. In addition, PIREPs are aircraft-dependent, so a smaller aircraft will experience stronger turbulence than a larger aircraft in the same volume of air. PIREPs also have poor spatial reliability as they tend to be located in turbulence, so null turbulence events are rarely recorded (Kane et al. 1998).
To avoid the above problems, there are two main aircraft-independent measures, which are the Eddy Dissipation Rate (EDR) and the Derived Equivalent Vertical Gust (DEVG). The measures are calculated using high-resolution automated aircraft data, available at 4 s intervals. The EDR is calculated following Haverdings and Chan (2010) using: where r w is the running mean of the vertical wind shear over a 10 second time window, V a is the true airspeed, and x 1 and x 2 are cut-off frequencies set at 0.1 and 2 Hz, respectively. The DEVG (Tenenbaum 1991;Gill 2014Gill , 2016Kim et al. 2017) is calculated following Truscott (2000) using: where jDnj is the peak modulus value of the deviation of the aircraft acceleration from 1g in units of g, m is the total mass of the aircraft (metric tonnes), V is the calibrated airspeed at the time of the observation (knots), and A is an aircraft-specific parameter that varies with flight conditions and can be calculated using the following: and: where H is the altitude (thousands of feet), m is the reference mass of the aircraft (metric tonnes), and parameters c 1 -c 5 depend on the aircraft's flight profile as outlined by Truscott (2000). Unlike PIREPs, EDR and DEVG are aircraftindependent, and therefore, all observations from multiple different aircraft can be combined to create a consistent observational database. Also unlike PIR-EPs, EDR and DEVG report all results including null reports. Because it is automated, the entire flight is recorded, and therefore, the only limit in coverage is the location the aircraft fly. Therefore, if every commercial aircraft had this capability, then all turbulence events could be recorded, and forecast verification improved. There are limitations to using the DEVG data set, as aircraft manoeuvres and active control techniques can enhance or dampen vertical accelerations of aircraft, leading to over-or underrepresentation of the vertical gusts (WMO 2003). However, despite these limitations, using an aircraftindependent measure (such as EDR or DEVG) is preferable to pilot reports (PIREPs).

Summary
Atmospheric turbulence is a major aviation hazard and we do not yet fully understand it. As a result, it causes many injuries to passengers and crew all around the world, with millions of dollars lost due to aircraft damage, planes being in the wrong place from a diverted flight, loss of hours from cabin crew, and compensation. This study has reviewed how turbulence that affects aviation can come from three main sources. The first is shear turbulence and is mainly found around the jet stream. The second is Convectively Induced Turbulence (CIT), and this forms around deep convection. However, there are multiple mechanisms that produce turbulence from convection, one of which is similar to the mechanism for the third source of turbulence known as Mountain-Wave Turbulence, which occurs when orography perturbs the atmosphere and produces gravity waves that can lead to turbulence.
This study has also reviewed how climate change is expected to have a large impact on clear-air turbulence in the future, with several studies finding that CAT will increase with climate change, particularly around the mid-latitudes in all seasons.
This makes turbulence forecasting very important in the future, to maintain the safety of passengers and crew. Currently, only one turbulence index known as the Ellrod and Knapp (1992) Turbulence Index 1 (TI1) is used by the two World Area Forecast Centres (WAFCs), but others have been included in the Graphical Turbulence Guidance (GTG) forecast. GTG combines many indices but is not yet used fully operationally in the WAFCs. Ensemble forecasting has also been developed to improve the forecasts, as currently only a deterministic forecast is created and that does not show the uncertainty in the forecast. The use of next-generation satellites (e.g., GOES-16) could help to improve nowcasting of turbulence, because with higher spatial and temporal resolutions, convection and gravity waves are more clearly resolved, aiding avoidance on the short time scales.
We know that the current methods of aviation turbulence forecasting are imperfect, and not all turbulence events are predicted. Therefore, continuing research to improve the forecast systems should be a high priority in the future. In addition, the aviation industry must make sure that the aircraft of the future are up to the challenge. This is because turbulence is not the only aviation hazard, with take-off weight restrictions due to high temperatures or aircraft icing also major hazards that might change in the future. In addition, the airline industry's desire to choose the most time-effective route to reduce fuel consumption and costs may force some pilots to fly into more turbulence regions. By making sure that all these aspects are taken into account, new and better working processes can be developed to make sure that the impact aviation has on the climate, and the impact climate change may have on aviation, are both reduced. With the volume of air travel also expected to increase over time, resulting in more congested skies, the research conducted now is vital to making sure that the industry is prepared for the future.