Analysis and preprocessing of Czech Mode-S observations

Benedikt Strajnar, Slovenian Environment Agency (ARSO)

in cooperation with Alena Trojáková, CHMI

Stay report, CHMI Prague, 6-24 July 2015 1 Introduction

Modern surveillance techniques allows air traffic controls (ATC) to collect a large num- ber of flight parameters from aircraft in real time. Among the available communication techniques, Mode-S (selective mode communication between ATC and airplane) is re- alized with a ground requesting radar (the so-called secondary radar), transponders on board the aircraft, and Mode-S ground sensors to collect the replies. These systems are of interest for meteorology because it is possible to extract meteorological information. This can be done either by using mandatory Mode-S data, the Enhanced Surveillance (EHS) parameters which are available in any Mode-S system (de Haan, 2011) or by di- rectly requesting wind and temperature messages contained in a non-mandatory Mode-S data register Meteorological Routine Air Report (MRAR) (Hrastovec and Solina, 2013; Strajnar, 2012). Both Mode-S EHS and MRAR methods have their advantages and disadvantages. A great advantage of Mode-S EHS observations is that EHS parameters are available in any Mode-S equipped ATC system and from all Mode-S equipped aircraft. However, the computation is indirect, especially for temperature, which can only be computed through speed of sound equation from Mach number and airspeed (de Haan, 2011). The wind is computed from a vector difference between air speed (including heading) and ground speed. The magnetic heading in particular has been shown to include aircraft type specific biases. A correction procedure was proposed by (de Haan, 2011) and then refined in (de Haan, 2013). On the other hand, Mode-S MRAR includes direct wind and temperature observations and whose quality was shown to be the same as the quality of Automated Meteorological Data Relay (AMDAR) on the AMDAR-equipped aircraft (Strajnar, 2012). However, because this data are not mandatory in a standard Mode-S system, the ATC radars must be configured to interrogate an additional MRAR register containing this data (where this is possible at all with respect to the density of air traffic). The support for MRAR also depends on the type of the transponder on board the aircraft. In the Slovenian case (Strajnar, 2012), MRAR was available from 5% of all Mode-S reports. Since few years, Mode-S observations (both EHS and MRAR) have been collected by Air Navigation Services of Czech Republic (ANS-CZ). The data was recently made available for evaluation also to Czech Hydrometeorological Institute (CHMI). This work is devoted to a basic analysis and preprocessing of these Mode-S observations. We first describe the data and provide examples of flight profiles. As a next step, raw Mode-S MRAR and Mode-S EHS (as computed for the moment at ANS CZ) are compared to AMDAR in a colocation study. A further evaluation is provided by a comparison to ALADIN/CZ numerical weather prediction model. This comparison enables an evalua- tion for each Mode-S observation (also non-AMDAR equipped aircraft). Description of scripts, software and data can be found last in the Appendix A.

2 Data description

2.1 Mode-S data set The Czech Mode-S data includes 3 Czech Mode-S radars Praha/Ruzyně, Písek (Jince) and Buchtův kopec (Sněžné). Those provide Mode-S EHS every 4 s and Mode-S MRAR every 10 seconds. Mode-S EHS is also available form 2 additional radars near Auersberg

1 Figure 1: (left) Mode-S MRAR and (right) Mode-S EHS coverage in 24 hours (from Jan Kráčmar and Jiří Frei, ANS CZ). Note the much higher density and slightly larger extent of Mode-S EHS.

(D) and Bratislava (SK). The horizontal coverage is shown in Fig. 1. There are around 3.5 million Mode-S EHS and around 140 thousand Mode-S MRAR observations per day. The time period used for analysis is 12 June - 9 July 2015.

2.2 ICAO aircraft addresses Within Mode-S, ICAO address is provided as an unique aircraft identification. It is as- signed by flight authorities in each country. Consequently, no centralized public database is publicly available. It is possible to relate ICAO address with aircraft information by online resources (e.g. various live flight tracking or plane spotter websites). Another possibility is to retrieve aircraft information from flight plans available at ATC. We used a combination of both to identify as many ICAO addresses as possible (more than 8000 of around 9000 different aircraft). Figure 2 shows the distribution of Mode-S EHS by aircraft type (without detailed subtypes). It can be noted that most frequent type is Boeing 737 followed by Airbus 320,319 and 321. The most frequent aircraft reporting Mode-S MRAR is Canadair Regional Jet (CRJ).

2.3 Flight profiles The easiest way to observe properties of Mode-S data sets is by plotting profiles of mea- sured variables during the flight. Figure 3 shows an ascent and a descent of a CRJ aircraft which provides both Mode-S EHS and Mode-S MRAR observations. For wind, a good agreement between Mode-S EHS and MRAR can be seen, although there is a small sys- tematic difference (smaller wind speeds in EHS observations). MRAR wind time series seems to be slightly more smooth, but oscillations of about 1-2 m/s can be observed both types. In the case of temperature (right part of Fig. 3 ), EHS is much more noisy and oscillates around MRAR values. Another example of temperature profile demonstrates sensitivity of EHS temperatures on ground speed (which is an input in EHS calculations, Fig. 4). Figure shows MRAR

2 Number of Mode−S EHS observations per day

Boeing 737 Airbus A320 Airbus A319 Airbus A321 Boeing 777 Airbus A330 Embraer ERJ Bombardier Dash 8 Q400 Boeing 757 Airbus A380 Canadair CRJ Fokker 100 Airbus A340 Boeing 747 Boeing 767 ATR 72 Boeing 787 ATR 42 Airbus A300B4 Airbus A318 Fokker 70 McDonnell Douglas MD Avro RJ100 Cessna 510 Citation Mustang Pilatus PC

0.0 0.2 0.4 0.6 0.8 1.0

Number [million]

Figure 2: Number of Mode-S EHS observations from different aircraft types (only the 25 most frequent are shown). The aircraft which also provide MRAR are plotted in red. temperature profile in blue and EHS temperatures in green. The oscillations of EHS temperature in this case are large (5-10 K). There are also some unphysical values which are related to errors in ground speed values (red line). The ground speed is constantly decreasing over time but there are a few zero values in the time series, causing unrealistic temperature calculation. A filtering would be needed to remove such ground speeds from the data set. A time smoothing is also suggested to suppress oscillations in EHS temperature (as already done by de Haan (2011)). EHS temperature calculation may especially be problematic at low ground and air speeds.

3 Flight profile Flight profile

Altitude Aircraft heading Aircraft Speed

Mode-S EHS 15 − Mode-S MRAR 80 20 − 9000 10000 210 220 10 15 25 −

Altitude

Mode-S EHS 30 Mode-S MRAR −

50 60 70 Aircraft heading Wind (m/s) Altitude(m) Altitude(m) Temperature [°C] Temperature Aircraft heading (°) Aircraft Aircraft heading (°) Aircraft Aircraft speed (m/s) Aircraft Aircraft speed (m/s) Aircraft 35 − 40 − 180 190 200 45 − 0 5 20 30 40 50 60 70 80 20 30 40 0.0 0.2 0.4 0.6 0.8 1.0 4000 5000 6000 7000 8000 9000 10000 4000 5000 6000 7000 8000

Time Time

Figure 3: Evolution of (left) wind and (right) temperature (and other parameters) during ascent and descent. Mode-S EHS is shown in green, Mode-S MRAR in blue. Landing profile 25 140 3000 4000

Altitude 200 250 300 350

Aircraft Speed 10 15 20 Aircraft heading Mode-S MRAR

Altitude(m) Mode-S EHS 60 80 100 120 Aircraft heading (°) Aircraft Aircraft speed (m/s) Aircraft Temperature (degrees) Temperature 1000 2000 − 5 0 5 0 20 40 0 50 100 150

10:00 12:00 14:00 16:00 18:00 20:00 22:00 Time (minute)

Figure 4: Evolution of EHS and MRAR temperatures and other reported parameters during descent. Aircraft is ATR 42, Czech Airlines. 3 Colocation of AMDAR, Mode-S EHS and Mode-S MRAR

To validate Mode-S observations other meteorological observations which are close in space and time are needed. An ideal observational reference for Mode-S is AMDAR 4 Figure 5: Vertical and temporal distribution of colocated points. which is available on several types of aircraft. Although AMDAR and Mode-S obser- vations originate from the same sensors, they differ due to the preprocessing or simply the reporting frequency and temporal averaging. Such a comparison can not highlight possible systematic instrumental errors affecting both observation sources. To find Mode-S EHS/MRAR and AMDAR observation pairs (colocated observations) we allowed a time mismatch of 30 s, maximal height difference of 50 m and 10 km horizontal separation. Such a tolerance is essential because AMDAR reports are not available so frequently as Mode-S. Due to preprocessing (e.g. smoothing or averaging) the exact time and space of the AMDAR measurement does not fully agree with flight data from Mode-S. If more possible pairs are found within these criteria the closest in space is selected. Figure 6 shows time and space distribution of colocated observation pairs. There were around 73200 EHS - AMDAR pairs and around 4850 MRAR - AMDAR pairs, slightly depending on the variable.

3.1 General difference statistics A general comparison between AMDAR and Mode-S is shown in Fig. 7. The difference between Mode-S and AMDAR are mostly normally distributed. The spread of differences against EHS is much larger than the spread against MRAR, both for wind and temper- atures. The wind difference distribution even shows a two peak distribution. To verify this result, the statistics are shown against the subset of MRAR observations (around

5 Distribution by height Distribution by height difference 12000 15000 8000 10000 Frequency Frequency 4000 5000 0 0

0 2000 6000 10000 −40 −20 0 20 40

Height [m] Height difference [m]

Distribution by horizontal distance Distribution by time separation 10000 15000 6000 Frequency Frequency 5000 2000 0 0

0 5 10 15 −30 −20 −10 0 10 20 30

Distance [km] Time difference [s]

Figure 6: Distribution of Mode-S EHS/MRAR - AMDAR matches by (a) height, and (b) vertical, (c) horizontal and (d) time separation.

6 Mode−S EHS − AMDAR Mode−S MRAR − AMDAR 2000 10000 1500 6000 1000 Number Number 500 2000 0 0

−6 −4 −2 0 2 4 6 −6 −4 −2 0 2 4 6 Temperature [K] Temperature [K] raw data, N = 73211 , mean= 0.34 , sd= 1.29 K raw data, N = 4851 , mean= −0.09 , sd= 0.32 K maxrange= −300.4 −> 141.8 maxrange= −3.3 −> 4.7

Mode S EHS − AMDAR Mode S MRAR − AMDAR 7000 2000 5000 1500 1000 3000 Frequency Frequency 500 1000 0 0

−6 −4 −2 0 2 4 6 −6 −4 −2 0 2 4 6 Wind speed [m/s] Wind speed [m/s] raw data, N = 71296 , mean= 0.35 , sd= 2.13 m/s raw data, N = 4829 , mean= 0 , sd= 0.7 m/s maxrange= −24.2 −> 485.4 maxrange= −5.7 −> 5.1

Mode S EHS − AMDAR Mode S MRAR − AMDAR 2500 15000 2000 1500 10000 Frequency Frequency 1000 5000 500 0 0

−150 −100 −50 0 50 100 150 −150 −100 −50 0 50 100 150 Wind direction [Deg] Wind direction [Deg] raw data, N = 72500 , mean= 0.34 , sd= 18.11 Deg raw data, N = 4825 , mean= 0.43 , sd= 6.81 Deg maxrange= −180 −> 179.9 maxrange= −161.8 −> 100.6

Figure 7: Histograms of (left) Mode-S EHS and (right) Mode-S MRAR - AMDAR differ- ences. Shown are (top) temperature , (middle) wind speed and (bottom) wind direction. To compute histograms, only the displayed value range is shown and extreme values are indicated below the plots.

4850 cases, 8). Here, the calculation is made for the same aircraft reporting both kinds of observations. While most of the distributions do not change, wind speed for EHS seems to become somewhat biased. This may be linked to the type of the aircraft (in this case mostly CRJ as explained later) and the errors in either ground or air speed or magnetic heading.

3.2 Type dependent statistics The colocated observations stem from different aircraft types. Figure 9 presents the distribution of collocated data set between different aircraft types. It can be observed that colocated Mode-S EHS observations come mostly from Airbus A319, A320 and A321, and that Mode-S MRAR comes from CRJ aircraft. It is important to keep in mind that the colocation result will be limited only to these aircraft types. In Figs. 10,11,12 the differences are shown separately for aircraft types. Because the

7 Mode−S EHS − AMDAR Mode−S MRAR − AMDAR 2000 600 1500 400 1000 Number Number 200 500 0 0

−6 −4 −2 0 2 4 6 −6 −4 −2 0 2 4 6 Temperature [K] Temperature [K] raw data, N = 4780 , mean= 0.52 , sd= 1.36 K raw data, N = 4851 , mean= −0.09 , sd= 0.32 K maxrange= −12.5 −> 141.8 maxrange= −3.3 −> 4.7

Mode S EHS − AMDAR Mode S MRAR − AMDAR 600 2000 500 1500 400 300 1000 Frequency Frequency 200 500 100 0 0

−6 −4 −2 0 2 4 6 −6 −4 −2 0 2 4 6 Wind speed [m/s] Wind speed [m/s] raw data, N = 4783 , mean= 0.71 , sd= 1.67 m/s raw data, N = 4829 , mean= 0 , sd= 0.7 m/s maxrange= −6.9 −> 10.8 maxrange= −5.7 −> 5.1

Mode S EHS − AMDAR Mode S MRAR − AMDAR 2500 1200 2000 1500 800 600 Frequency Frequency 1000 400 500 200 0 0

−150 −100 −50 0 50 100 150 −150 −100 −50 0 50 100 150 Wind direction [Deg] Wind direction [Deg] raw data, N = 4773 , mean= 1.76 , sd= 11.89 Deg raw data, N = 4825 , mean= 0.43 , sd= 6.81 Deg maxrange= −179.9 −> 172.1 maxrange= −161.8 −> 100.6

Figure 8: Same as Fig. 7 except that Mode-S EHS comparison is only displayed for the subset where MRAR is also available.

8 Number of colocated Mode−S EHS by type Number of colocated Mode−S MRAR by type

Saab 2000 Piper PA44 Piper PA Saab 2000 McDonnell Douglas MD Hawker 800XPi Gulfstream G550 Gulfstream 550 Fokker 70 Fokker 100 Hawker Beechcraft 400XP Embraer ERJ 900EX Dassault Falcon 2000LX Cessna 560XLS Citation Excel Cessna 525 CitationJet CJ1 Cessna 525B CitationJet CJ3 Dassault Falcon 900EX Cessna 510 Citation Mustang Canadair CRJ Bombardier Global 6000 Bombardier Dash 8 Q400 Bombardier CRJ Bombardier Challenger 300 Boeing 787 Canadair CRJ Boeing 777 Boeing 767 Boeing 757 Boeing 747 Boeing 737 Boeing 717 Canadair Challenger 604 Beech C90GTX King Air Beech C90GTi King Air Avro RJ100 ATR 72 ATR 42 Airbus A380 Bombardier CRJ Airbus A340 Airbus A330 Airbus A321 Airbus A320 Airbus A319 Airbus A318 Airbus A310 Beech C90GTX King Air Airbus A300B4

0 5000 10000 15000 20000 0 1000 2000 3000 4000

Number of obs. (T or wind) Number of obs. (T or wind)

Figure 9: Number of (left) Mode-S EHS and (right) Mode-S MRAR by aircraft type.

MRAR colocated set is dominated only by CRJ, these plots are shown only for Mode-S EHS. The interpretation should be limited to the above mentioned Airbus types and CRJ because the other types are rarely colocated with AMDAR over the region. A large temperature error standard deviations of around 10 K can be seen for A319-A321 while their biases are reasonable but positive. A large warm bias appears for A318. In the case of wind speed, the errors of A319 seem to be larger that for the other Airbuses and they also include significant biases.

9 Mode−S EHS vs. AMDAR Mode−S EHS vs. AMDAR

Saab 2000 Saab 2000 Piper PA44 Piper PA44 Piper PA Piper PA McDonnell Douglas MD McDonnell Douglas MD Hawker 800XPi Hawker 800XPi Gulfstream G550 Gulfstream G550 Gulfstream 550 Gulfstream 550 Fokker 70 Fokker 70 Fokker 100 Fokker 100 Eurocopter EC135 P2 Eurocopter EC135 P2 Embraer Legacy 600 Embraer Legacy 600 Embraer ERJ Embraer ERJ Dassault Falcon 900EX Dassault Falcon 900EX Dassault Falcon 2000LX Dassault Falcon 2000LX Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel Cessna 560XLS Citation Excel Cessna 525 CitationJet CJ1 Cessna 525 CitationJet CJ1 Cessna 525B CitationJet CJ3 Cessna 525B CitationJet CJ3 Cessna 525A CitationJet CJ2 Cessna 525A CitationJet CJ2 Cessna 510 Citation Mustang Cessna 510 Citation Mustang Canadair CRJ Canadair CRJ Bombardier Global 6000 Bombardier Global 6000 Bombardier Dash 8 Q400 Bombardier Dash 8 Q400 Bombardier CRJ Bombardier CRJ Bombardier Challenger 300 Bombardier Challenger 300 Boeing 787 Boeing 787 Boeing 777 Boeing 777 Boeing 767 Boeing 767 Boeing 757 Boeing 757 Boeing 747 Boeing 747 Boeing 737 Boeing 737 Boeing 717 Boeing 717 Beech C90GTX King Air Beech C90GTX King Air Beech C90GTi King Air Beech C90GTi King Air Avro RJ100 Avro RJ100 ATR 72 ATR 72 ATR 42 ATR 42 Airbus A380 Airbus A380 Airbus A340 Airbus A340 Airbus A330 Airbus A330 Airbus A321 Airbus A321 Airbus A320 Airbus A320 Airbus A319 Airbus A319 Airbus A318 Airbus A318 Airbus A310 Airbus A310 Airbus A300B4 Airbus A300B4

0 20 40 60 80 100 120 140 −100 −50 0

Temperature standard deviation [K] Temperature mean [K]

Figure 10: (left) Standard deviation and (right) mean temperature difference between Mode-S EHS and AMDAR by aircraft type.

Mode−S EHS vs. AMDAR Mode−S EHS vs. AMDAR

Saab 2000 Saab 2000 Piper PA44 Piper PA44 Piper PA Piper PA McDonnell Douglas MD McDonnell Douglas MD Hawker 800XPi Hawker 800XPi Gulfstream G550 Gulfstream G550 Gulfstream 550 Gulfstream 550 Fokker 70 Fokker 70 Fokker 100 Fokker 100 Embraer Legacy 600 Embraer Legacy 600 Embraer ERJ Embraer ERJ Dassault Falcon 900EX Dassault Falcon 900EX Dassault Falcon 2000LX Dassault Falcon 2000LX Cessna 560XLS Citation Excel Cessna 560XLS Citation Excel Cessna 525B CitationJet CJ3 Cessna 525B CitationJet CJ3 Cessna 510 Citation Mustang Cessna 510 Citation Mustang Canadair CRJ Canadair CRJ Bombardier Global 6000 Bombardier Global 6000 Bombardier Dash 8 Q400 Bombardier Dash 8 Q400 Bombardier CRJ Bombardier CRJ Bombardier Challenger 300 Bombardier Challenger 300 Boeing 787 Boeing 787 Boeing 777 Boeing 777 Boeing 767 Boeing 767 Boeing 757 Boeing 757 Boeing 747 Boeing 747 Boeing 737 Boeing 737 Boeing 717 Boeing 717 Beech C90GTX King Air Beech C90GTX King Air Beech C90GTi King Air Beech C90GTi King Air Avro RJ100 Avro RJ100 ATR 72 ATR 72 ATR 42 ATR 42 Airbus A380 Airbus A380 Airbus A340 Airbus A340 Airbus A330 Airbus A330 Airbus A321 Airbus A321 Airbus A320 Airbus A320 Airbus A319 Airbus A319 Airbus A318 Airbus A318 Airbus A310 Airbus A310 Airbus A300B4 Airbus A300B4

0 5 10 15 20 −40 −20 0 20 40

Wind speed standard deviation [m/s] Wind speed mean [m/s]

Figure 11: (left) Standard deviation and (right) mean wind speed difference between Mode-S MRAR and AMDAR by aircraft type.

10 Mode−S EHS vs. AMDAR Mode−S EHS vs. AMDAR

Saab 2000 Saab 2000 Piper PA44 Piper PA44 Piper PA Piper PA McDonnell Douglas MD McDonnell Douglas MD Hawker 800XPi Hawker 800XPi Gulfstream G550 Gulfstream G550 Gulfstream 550 Gulfstream 550 Fokker 70 Fokker 70 Fokker 100 Fokker 100 Embraer Legacy 600 Embraer Legacy 600 Embraer ERJ Embraer ERJ Dassault Falcon 900EX Dassault Falcon 900EX Dassault Falcon 2000LX Dassault Falcon 2000LX Cessna 560XLS Citation Excel Cessna 560XLS Citation Excel Cessna 525B CitationJet CJ3 Cessna 525B CitationJet CJ3 Cessna 510 Citation Mustang Cessna 510 Citation Mustang Canadair CRJ Canadair CRJ Bombardier Global 6000 Bombardier Global 6000 Bombardier Dash 8 Q400 Bombardier Dash 8 Q400 Bombardier CRJ Bombardier CRJ Bombardier Challenger 300 Bombardier Challenger 300 Boeing 787 Boeing 787 Boeing 777 Boeing 777 Boeing 767 Boeing 767 Boeing 757 Boeing 757 Boeing 747 Boeing 747 Boeing 737 Boeing 737 Boeing 717 Boeing 717 Beech C90GTX King Air Beech C90GTX King Air Beech C90GTi King Air Beech C90GTi King Air Avro RJ100 Avro RJ100 ATR 72 ATR 72 ATR 42 ATR 42 Airbus A380 Airbus A380 Airbus A340 Airbus A340 Airbus A330 Airbus A330 Airbus A321 Airbus A321 Airbus A320 Airbus A320 Airbus A319 Airbus A319 Airbus A318 Airbus A318 Airbus A310 Airbus A310 Airbus A300B4 Airbus A300B4

0 10 20 30 40 50 60 −50 0 50

Wind direction standard deviation [deg] Wind direction mean [deg]

Figure 12: (left) Standard deviation and (right) mean wind direction difference between Mode-S MRAR and AMDAR by aircraft type.

Table 1: Thresholds used to generate white list of aircraft reporting Mode-S MRAR. variable number of observations mean std. temperature (K) 3000 < 1 K < 2 K wind speed (m/s) 3000 < 1 m/s < 5 m/s wind direction (deg) 3000 < 10◦ < 100◦

4 Validation against NWP

The comparison against AMDAR offers a good observations reference to Mode-S data but is limited to AMDAR-equipped aircraft. An evaluation of data from all aircraft is possible only by comparison against NWP model. This allows sufficiently large samples for each single aircraft and therefore a reliable error estimates. However, such a comparison is limited by forecast and model errors. In this study we used ALADIN/CZ forecasts of various lengths, depending on the hour of day. As ALADIN/CZ performs an analysis every 6 hours, we used forecasts of lengths 6-11 hours. Analysis was avoided because AMDAR observations are also used, which may in some cases prevent fair comparison. We assume that the method is robust with respect to slight decrease of the quality of forecasts with range. To find the model counterpart of each Mode-S observation we used screening configuration of the ALADIN model (e002) and saved obs-minus-guess departures stored in the observational database (ODB). The analysis was first performed for Mode-S MRAR. Partly following Strajnar et al. (2015), criteria were defined to regard measurements from an aircraft as of a good quality (Table 1). Because wind speed and direction are measured separately by measuring air speed and heading, criteria with respect to those wind variables were preferred over u and v wind components. To achieve this, the wind departures were first recalculated from observations and u and v component differences which are output of screening.

11 Figure 13 shows difference statistics computed from aircraft on the white list, dis- tributed by type. For example, the turboprop ATR aircraft (often flying in the Czech airspace and reporting MRAR) are not of sufficient quality and thus not displayed. The final white list includes 127 aircraft addresses for temperatures and 116 addresses for wind. The same procedure can also be applied to Mode-S EHS. The aggregation of the statistics takes longer because the data set is large (> 30 GB of data). For tempera- ture, 2036 out of 7278 aircraft pass the criteria defined in Table 1. For wind, the white lists was not yet computed because the model departures are so far only defined for u v components. Figure 14 shows the number of observations per observation type for wind and temperature and Fig. 15 displays the distributions of Mode-S EHS - ALADIN differences computed separately for each aircraft (wind is shown only as u components) within thresholds defined in Table 1. This suggests that for a large part of aircraft the quality of Mode-S EHS is also potentially acceptable. To create a white list of Mode-S EHS observations to be used in data assimilation experiments, the Mode-S - ALADIN differences should be recomputed in the form of wind speed and direction, like in the case of Mode-S MRAR.

12 Selected Mode−S MRAR by type Selected Mode−S MRAR by type

Saab 2000 Saab 2000 Raytheon 390 Premier Raytheon 390 Premier Learjet 60 Learjet 60 Learjet 35A Learjet 35A Hawker 900XP Hawker 900XP Hawker 800XP Hawker 800XP Gulfstream G200 Galaxy Gulfstream G200 Galaxy Dassault Falcon 900EX Dassault Falcon 900EX Dassault Falcon 900 Dassault Falcon 2000 Dassault Falcon 2000 Falcon 2000 Dassault Aviation Falcon 2000 Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel Cessna 560XLS Citation Excel Cessna 560XL Citation XLS+ Cessna 560XL Citation XLS+ Cessna 560XL Citation Excel Cessna 560XL Citation Excel Cessna 560 Citation Encore+ Cessna 560 Citation Encore+ Cessna 525 CitationJet CJ1+ Cessna 525 CitationJet CJ1+ Canadair CRJ700 Canadair CRJ700 Canadair CRJ Canadair CRJ Canadair CL605 Challenger Canadair CL605 Challenger Canadair CL604 Challenger Canadair CL604 Challenger Canadair CL601 Challenger Canadair CL601 Challenger Canadair Challenger 850 Canadair Challenger 850 Canadair Challenger 604 Canadair Challenger 604 Bombardier CRJ Bombardier CRJ Bombardier Challenger 300 Bombardier Challenger 300 Bombardier BD100 Challenger 300 Bombardier BD100 Challenger 300 Beech 390 Premier Beech 390 Premier AMD AMD

0.0 0.5 1.0 1.5 2.0 2.5 3.0 −1.0 −0.5 0.0 0.5 1.0

Std. of Mode−S − ALADIN temperature difference [K] Mean Mode−S − ALADIN temperature difference [K]

Selected Mode−S MRAR by type Selected Mode−S MRAR by type

Saab 2000 Saab 2000 Piper PA Piper PA Learjet 60 Learjet 60 Hawker 900XP Hawker 900XP Hawker 800XPi Hawker 800XPi Hawker 800XP Hawker 800XP Hawker 750 Hawker 750 Gulfstream G200 Galaxy Gulfstream G200 Galaxy Gulfstream G100 Gulfstream G100 Dassault Falcon 900EX Dassault Falcon 900EX Dassault Falcon 900 Dassault Falcon 900 Dassault Falcon 2000 Dassault Falcon 2000 Dassault Aviation Falcon 2000 Dassault Aviation Falcon 2000 Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel Cessna 560XLS Citation Excel Cessna 560XL Citation XLS+ Cessna 560XL Citation XLS+ Cessna 560XL Citation XLS Cessna 560XL Citation XLS Cessna 560XL Citation Excel Cessna 560XL Citation Excel Cessna 560 Citation Encore+ Cessna 560 Citation Encore+ Cessna 525 CitationJet CJ1+ Cessna 525 CitationJet CJ1+ Cessna 525B CitationJet CJ3 Cessna 525B CitationJet CJ3 Cessna 525A CitationJet CJ2+ Cessna 525A CitationJet CJ2+ Cessna 525A Citation CJ2+ Cessna 525A Citation CJ2+ Canadair CRJ700 Canadair CRJ700 Canadair CRJ Canadair CRJ Canadair CL605 Challenger Canadair CL605 Challenger Canadair CL604 Challenger Canadair CL604 Challenger Canadair CL601 Challenger Canadair CL601 Challenger Canadair Challenger 850 Canadair Challenger 850 Canadair Challenger 604 Canadair Challenger 604 Bombardier CRJ Bombardier CRJ Bombardier Challenger 300 Bombardier Challenger 300 Bombardier BD100 Challenger 300 Bombardier BD100 Challenger 300 Beech C90GTX King Air Beech C90GTX King Air Beech C90GTi KIng Air Beech C90GTi KIng Air Beech C90GTi King Air Beech C90GTi King Air Beech B200GT Super King Air Beech B200GT Super King Air Beech 300 Super King Air 350 Beech 300 Super King Air 350 Beech 200GT Super King Air Beech 200GT Super King Air AMD AMD

0 1 2 3 4 5 −2 −1 0 1 2

Std. of Mode−S − ALADIN wind speed difference [m/s] Mean Mode−S − ALADIN wind speed difference [m/s]

Selected Mode−S MRAR by type Selected Mode−S MRAR by type

Saab 2000 Saab 2000 Piper PA Piper PA Learjet 60 Learjet 60 Hawker 900XP Hawker 900XP Hawker 800XPi Hawker 800XPi Hawker 800XP Hawker 800XP Hawker 750 Hawker 750 Gulfstream G200 Galaxy Gulfstream G200 Galaxy Gulfstream G100 Gulfstream G100 Dassault Falcon 900EX Dassault Falcon 900EX Dassault Falcon 900 Dassault Falcon 900 Dassault Falcon 2000 Dassault Falcon 2000 Dassault Aviation Falcon 2000 Dassault Aviation Falcon 2000 Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel + Cessna 560XLS Citation Excel Cessna 560XLS Citation Excel Cessna 560XL Citation XLS+ Cessna 560XL Citation XLS+ Cessna 560XL Citation XLS Cessna 560XL Citation XLS Cessna 560XL Citation Excel Cessna 560XL Citation Excel Cessna 560 Citation Encore+ Cessna 560 Citation Encore+ Cessna 525 CitationJet CJ1+ Cessna 525 CitationJet CJ1+ Cessna 525B CitationJet CJ3 Cessna 525B CitationJet CJ3 Cessna 525A CitationJet CJ2+ Cessna 525A CitationJet CJ2+ Cessna 525A Citation CJ2+ Cessna 525A Citation CJ2+ Canadair CRJ700 Canadair CRJ700 Canadair CRJ Canadair CRJ Canadair CL605 Challenger Canadair CL605 Challenger Canadair CL604 Challenger Canadair CL604 Challenger Canadair CL601 Challenger Canadair CL601 Challenger Canadair Challenger 850 Canadair Challenger 850 Canadair Challenger 604 Canadair Challenger 604 Bombardier CRJ Bombardier CRJ Bombardier Challenger 300 Bombardier Challenger 300 Bombardier BD100 Challenger 300 Bombardier BD100 Challenger 300 Beech C90GTX King Air Beech C90GTX King Air Beech C90GTi KIng Air Beech C90GTi KIng Air Beech C90GTi King Air Beech C90GTi King Air Beech B200GT Super King Air Beech B200GT Super King Air Beech 300 Super King Air 350 Beech 300 Super King Air 350 Beech 200GT Super King Air Beech 200GT Super King Air AMD AMD

0 20 40 60 80 100 −10 −5 0 5 10

Std. of Mode−S − ALADIN wind direction difference [deg] Mean Mode−S − ALADIN wind direction difference [deg]

Figure 13: (left) Standard deviation and (right) mean difference between Mode-S MRAR 13 and ALADIN forecast by aircraft type for (top) temperature, (middle) wind speed and (bottom) wind direction. Selected Mode−S EHS by type Selected Mode−S EHS by type

Tecnam P2006T McDonnell Douglas MD Raytheon 390 Premier Fokker 100 Mooney M20TN Let L Embraer ERJ Learjet 60 Gulfstream G550 Embraer EMBB500 Phenom 100 Fokker 70 Fokker 100 Canadair CRJ700 Embraer ERJ Canadair CRJ Diamond DA Dassault Falcon 2000S Canadair CL605 Challenger Dassault Falcon 2000LX Cirrus SR22 Bombardier Challenger 300 Cirrrus SR22 Cessna T206H Turbo Stationair Boeing 787 Cessna T182T Turbo Skylane Cessna 680 Citation Sovereign Boeing 777 Cessna 560XL Citation XLS+ Boeing 767 Cessna 525A Citation CJ2+ Cessna 510 Citation Mustang Boeing 757 Cessna 208B Grand Caravan Cessna 172S Skyhawk SP Boeing 747 Canadair CRJ700 Canadair CRJ Boeing 737 Canadair CL605 Challenger Boeing 737 Canadair CL604 Challenger Bombardier Global 6000 Boeing 717 Bombardier Dash 8 Q400 Bombardier CRJ AirbusA320 Bombardier Challenger 300 Boeing 787 Airbus KC2 Voyager (A330 Boeing 777 Airbus A380 Boeing 767 Boeing 757 Airbus A350 Boeing 747 Boeing 737 Airbus A340 Beech 300 Super King Air 350 ATR 72 Airbus A330 ATR 42 Airbus KC2 Voyager (A330 Airbus A321 Airbus A380 Airbus A320 Airbus A350 Airbus A340 Airbus A319 Airbus A330 Airbus A321 Airbus A318 Airbus A320 Airbus A319 Airbus A310 Airbus A318 Airbus A300C4 Airbus A310 Airbus A300C4 Airbus A300B4 Airbus A300B4 Airbus Airbus

0 5 10 15 0 5 10 15

Number of temperature obs. [million] Number of wind obs. [million]

Figure 14: (left) Number of Mode-S EHS (left) temperature and (right) wind observations by aircraft type.

5 Conclusions and outlook

In this preliminary evaluation of Mode-S data available over the Czech airspace and surroundings, we demonstrated the observation quality with respect to AMDAR obser- vations and ALADIN NWP model, separately for each reporting aircraft or aircraft type. The difference in availability and quality between Mode-S MRAR and Mode-S EHS was also demonstrated. While Mode-S MRAR observations are of overall good quality even without preprocessing, they can be used for meteorological applications including data as- similation just after the basic data selection based on statistics of differences with respect to ALADIN model. In the case of Mode-S EHS a larger variability in the data and their errors was ob- served, and these could be further improved by smoothing the data or applying corrections to temperature wind speed and wind direction, as done at KNMI (de Haan, 2013). It appears that Mode-S EHS generally include a slight warm bias. The possible future steps include:

• checking for unrealistic air speeds resulting in unphysical temperatures in some cases (see Fig. 4)

• temporal smoothing of variables which appear in EHS temperature calculation (Mach number, air speed) and recomputing EHS temperatures, investigating the effect on warm bias

• computation of magnetic heading correction separately for each aircraft and re- computation of EHS wind observations, optional temporal smoothing (currently not performed at KNMI)

• repeated colocation with AMDAR and evaluation of impact of such corrections

14 Mode−S EHS − ALADIN distribution by aircraft Mode−S EHS − ALADIN distribution by aircraft 400 500 400 300 300 200 Frequency Frequency 200 100 100 0 0

1.0 1.2 1.4 1.6 1.8 2.0 −1.0 −0.5 0.0 0.5 1.0

Std. temperature difference [K] Mean temperature difference [K]

Mode−S EHS − ALADIN distribution by aircraft 400 300 200 Frequency 100 0

2.0 2.5 3.0 3.5 4.0 4.5 5.0

Std. u wind component difference [m/s]

Figure 15: Distributions of (left) standard deviation and (right) mean Mode-S EHS - ALADIN differences computed separately for each aircraft. Shown are (top) temperature and (bottom) u wind component.

15 • repeated comparison against ALADIN model, recalculation of wind speed and di- rection difference from wind components and data selection based on thresholds used in Table 1 to generate a white list for EHS observations. The data assimilation experiments should be carried out first with Mode-S MRAR and then Mode-S EHS. Recent studies suggested that improvements in the first few hours of forecasts can be expected (de Haan and Stoffelen, 2012; Strajnar et al., 2015). As the quality of EHS temperatures were shown to be less than the quality of wind and if this does not improve significantly with further correction it may be reasonable to assimilate EHS winds only on top of all Mode-S MRAR observations.

Acknowledgment

We thank ANS CZ for providing Mode-S observations used in this work. Jan Kračmar and Jiří Frei helped with description of data, discussion and reconstructing the aircraft details from flight plans. We are also grateful for organizing the visit to ANS CR. Many thank to colleagues from CHMI for support during the stay.

6 References

de Haan, S. (2011). High-resolution wind and temperature observations from aircraft tracked by Mode-S air traffic control radar. J. Geophys. Res., 116. de Haan, S. (2013). An improved correction method for high quality wind and tempera- ture observations derived from mode-s ehs. Technical report. de Haan, S. and Stoffelen, A. (2012). Assimilation of high-resolution Mode-S wind and temperature observations in a regional NWP model for nowcasting applications. Wea. Forecasting, 27:918–937. Hrastovec, M. and Solina, F. (2013). Obtaining meteorological data from aircraft with Mode-S radars. Aerospace and Electronic Systems Magazine, 28(12). Strajnar, B. (2012). Validation of Mode-S Meteorological Routine Air Report aircraft observations. J. Geophys. Res., 117. Strajnar, B., Žagar, N., and Berre, L. (2015). Impact of new Mode-S MRAR observations in a mesoscale model. J. Geophys. Res., 120(9).

A Scripts and data sets

The following scripts and data sets were prepared on server yaga, user mma233. The data sets are archived on directory /work/data/modes˜ (Table A). Scripts are located under /scripts/R˜ and /python˜ and possible their sub directories (Table 3).

16 Table 2: List of data sets. data set description colocated_AMDAR csv files containing colocated AMDAR, EHS and MRAR observational data screened_mrar_csv csv files containing screened MRAR observations screened_ehs_csv csv files containing screened EHS observations mrar_merged_1h csv files with merged MRAR observations, ALADIN departures and aircraft details ehs_merged_1h csv files with merged EHS observations, ALADIN departures and aircraft details ehs_merged_byICAO csv files with merged EHS observations, separate file for each aircraft address

Table 3: List of R and python scripts. script purpose profiles.R plot flight profiles compute_colocations.R compute colocated AMDAR and Mode-S points, merge both information and outprint analyse_amdar_colocation.R compute and plot statistics from colocated data analyse_amdar_colocation_by_type.R compute and plot type-dependent statistics from colocated data merge_ModeS_ehs_NWP.R prepare merged data sets for EHS merge_ModeS_mrar_NWP.R prepare merged data sets for MRAR create_whitelists_and_stats_mrar.R prepare statistics and create white lists for MRAR create_whitelists_and_stats_ehs.R prepare statistics and create white lists for EHS (not fully completed) create_stats_ehs_by_address.R create files with NWP difference statistics separately for each aircraft address superobs.R create smoothed super observations for MRAR (not used in this work) split_ehs_by_ICAOaddr.py split EHS data by aircraft registration create_screened_csv_ehs.py creation of screened EHS create_screened_csv_mrar.py create of screened MRAR ecma2modescsv.py combine variables into one row for a single observation, called by the two above scripts

17