## Classi-Fly: Inferring Aircraft Categories from Open Data using Machine Learning

## Martin Strohmeier Vincent Lenders

firstname.lastname@armasuisse.ch armasuisse Science and Technology Thun, Switzerland

## Matthew Smith Ivan Martinovic

firstname.lastname@cs.ox.ac.uk Department of Computer Science, University of Oxford Oxford, United Kingdom

## ABSTRACT

In recent years, air traffic communication data has become easy to access, enabling novel research in many fields. Exploiting this new data source, a wide range of applications have emerged, from weather forecasting to stock market prediction, or the collection of information about military and government movements. Typically these applications require knowledge about the metadata of the aircraft, specifically its operator and the aircraft category.

armasuisse Science + Technology , the R&amp;D agency for the Swiss Armed Forces, has been developing Classi-Fly, a novel approach to obtain metadata about aircraft based on their movement patterns. We validate Classi-Fly using several hundred thousand flights collected through open source means, in conjunction with ground truth from publicly available aircraft registries containing more than two million aircraft. Classi-Fly obtains the correct aircraft category with an accuracy of over 88%, demonstrating that it can improve the meta data necessary for applications working with air traffic communication. Finally, we show that it is feasible to automatically detect specific flights such as police and surveillance missions.

## CCS CONCEPTS

· Information systems ! Clustering and classification; Business intelligence; · Applied computing ! Aerospace;

## KEYWORDS

aviation, aircraft classification, open datasets, air traffic control, air traffic management, open source intelligence

## 1 INTRODUCTION

Crowdsourced aircraft trajectory data has gained significant importance in recent years enhancing data collection for several scientific fields and opening up new research opportunities. With the largescale and open availability of globally distributed ground sensors, hundreds of research publications across different subjects have been exploiting this novel data source. 1 Most, if not all, of these investigations require knowledge about the broad category of an aircraft, e.g., whether it is a commercial airliner, a fighter jet, or a small private airplane; and its use case (e.g., government, business, surveillance, or military).

1 A regularly maintained list of examples can be found at https://opensky-network. org/community/publications.

This development has been enabled by the rise of softwaredefined radios (SDRs), which have become readily available and affordable over the past decade and greatly reduced the barriers to entry, meaning that users can now take part in crowdsourced sensor networks with little cost. As such, it is now straightforward to collect all air traffic surveillance communication directly from aircraft, providing information about flight trajectories, including its position, velocity and unique identifiers in order to tell them apart and track them over time. Based on this principle, many commercial ventures have been created, for example Flightradar24 and FlightAware, who in turn provide sophisticated data-based services to aviation stakeholders.

The first area, epidemiological modeling of the possible spread of COVID-19, was of crucial interest early in the stages of the pandemic. The utility of flight data for this purpose was illustrated for example in widely circulated studies such as [6] but has been knowntobeuseful in the context of pandemics for much longer (e.g., [22]). Knowing aircraft category and operator can help improve the accuracy of these models, i.e. in estimating the number of passengers on each tracked aircraft.

In the wake of COVID-19, the interest in trajectory-based research has accelerated even further. Relevant research can be separated into two different areas of modeling, pandemic and economic .

The second main area, economic modeling, uses flights either as an indicator of economic activity (at a given airport, region, or globally) as illustrated in [23] or as a direct measure of the impact on the aviation sector (in particular cargo and passenger transport). Here, the speed of these indicators is the crucial advantage, as they allow to 'nowcast' the economy faster than traditional methods. Examples of such use of data provided by OpenSky can be found in the Bank of England's quarterly Monetary Policy Report [4] or the National Statistics Office of Denmark [37]. Again, obtaining accurate metadata about the operators and purpose of the tracked aircraft is crucial for improving the accuracy of these models.

Crucially, contrary to the trajectories, there is no helpful information about the aircraft types, categories or operators broadcast by the aircraft themselves. To solve this issue, researchers rely broadly on external databases, maintained through a mix of crowdsourcing by aviation enthusiasts and official databases provided by a few countries' aviation agencies such as the US Federal Aviation Administration (FAA) [11]. Unfortunately, these sources are of only limited use as they are incomplete and outdated for many of the world's aircraft, which poses a major challenge to applications in this area.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01061v2-figure-c0170da0-2816-4520-9cff-f9cfab21c882.jpg)

ATC Ground Stations

Figure 1: Representation of ADS-B and SSR systems.

In this paper, we present a machine-learning approach to solve this problem called Classi-Fly. We first deal with the challenges of incomplete and unreliable ground truth by verifying the categories of our dataset manually. Using this dataset, we show that it is feasible to automatically classify aircraft into different operator categories based purely on their learned flight movement patterns. Contrary to other classification approaches, Classi-Fly works on features directly derived from the trajectories, which cannot be altered by the aircraft operator, a crucial advantage over classical methods.

In this paper, we make the following contributions:

- GLYPH&lt;15&gt; On a dataset of 182,325 flights, we show that it is feasible to automatically learn the category of a given aircraft with over 88% accuracy based solely on its flight behaviour. We build a model based on features derived from this behaviour and compare the accuracy of four different classifiers.
- GLYPH&lt;15&gt; We discuss the implications of our method, including potential countermeasures, and analyze a case study of previously unidentified aircraft with sensitive mission profiles.
- GLYPH&lt;15&gt; Using this approach, we classify previously unknown aircraft into different categories, effectively deriving metadata information for these aircraft, which can be used for popular research applications from open-source intelligence to epidemiological modeling.

The remainder of this paper is structured as follows: Section 2 describes the necessary background on air traffic control and tracking. Section 3 describes the fundamentals of our approach. 4 introduces our data collection process. Section 5 describes our experimental design before Sections 6 and 7 present the performance of the analyzed classifiers on the ground truth and new data, respectively. Finally, Section 8 and 9 cover the discussion and related work before Section 10 concludes.

## 2 BACKGROUND

This section provides the necessary background to how aircraft tracking works. Fig. 1 shows the wireless communication links of two considered technologies, which are explained in the following.

## 2.1 Surveillance Technologies in Aviation

There are two main surveillance technologies used for cooperative tracking of civil aircraft. Secondary Surveillance Radar (SSR) uses the so-called transponder Modes A, C, and S, which provide digital target information (altitude, squawk identification) compared to traditional primary radar (PSR). Aircraft transponders are interrogated on the 1030 MHz frequency and reply with the desired information on the 1090 MHz channel (see Fig. 1, right.) With the newer Automatic Dependent Surveillance-Broadcast (ADS-B) protocol (see Fig. 1, left), aircraft regularly broadcast their own identity, position, velocity and other information such as intent or emergency codes. These broadcasts do not require interrogation; position and velocity are automatically transmitted at 2 Hz [27].

## 2.2 Aircraft Identifiers in Air Traffic Communication

A 24-bit address assigned by the International Civil Aviation Organization (ICAO) to every aircraft is transmitted via both ADS-B and SSR. Crucially, this identifier is different to an aircraft squawk or call sign . Squawks, of which only 4096 exist, are allocated locally and not effective for continuous tracking. The call sign can be set separately through the flight deck for every flight. Call signs of private aircraft typically consist of the aircraft registration number, commercial airliners use the flight number, and military and government operators often use special call signs depending on their mission. In contrast, the ICAO identifier is globally unique and provides an address space of 16 million; while the transponder can be re-programmed by engineers, the identifier is not easily (or legally) changed by a pilot. These characteristics make it ideal for continuous tracking over a prolonged period of time.

## 2.3 Required Data Mining Capabilities

Aircraft tracking is the act of obtaining live or delayed positional information on aircraft by passive actors. Their motivations range from traditional hobbyist planespotting enthusiasm over military and business interests to environmental science. Where traditionally most spotters have conducted their trade purely using visual means, i.e., seeing and recognizing the aircraft near an airport, modern software-defined radio (SDR) technology has made accurate, fast and scalable tracking of aircraft feasible for anyone.

There are two options to exploit SDRs: install their own personal receivers or use the SDR data aggregated by web tracking services. While a single receiver with a radius of up to 600 km can already provide interesting results, the insights are increased considerably with a larger network. Both live tracking data and the required metadata are easily accessible on-line as discussed in Section 4.

## 3 CLASSI-FLY

Classi-Fly is a machine-learning approach to categorize an aircraft purely on its behaviour, i.e. the way it moves over time rather than relying on self-reported information such as call signs and / or fallible databases.

More concretely, Classi-Fly analyses collected trajectories of aircraft and breaks them down into 12 principal features, based on flight duration, position, velocity and acceleration. Using these features, it learns the behaviour of 8 different aircraft categories,

Table 1: Description of the ground truth dataset, comprising 9880 randomly selected aircraft with a maximum of 25 flights.

| Flights /   | Aircraft   | States / Flight   | Duration / Flight [s]   |
|-------------|------------|-------------------|-------------------------|
| Mean        | 20.3148    | 152.62            | 4669                    |
| Median      | 25         | 79                | 1897                    |
| Total       | 200,710    | 30,633,219        | 937,223,660             |

from commercial to military and government aircraft. This in turn enables the user to automatically classify aircraft improve a knowledge base that is often incomplete and erroneous.

## 3.1 Applications

During the COVID-19 pandemic, aircraft trajectory-based research has accelerated significantly. Scientific modeling has used aircraft trajectory data to measure the impact of the implemented mitigation measures on the economy [12, 16], nature [18] and air traffic itself [36]. The growing interest is reflected in the popularity of new datasets, which have seen significant uptake in a short time [25].

To support these use cases, Classi-Fly can also contribute directly towards open data initiatives such as the OpenSky Network aircraft metadata database, which in turn is used for a wide variety of research applications (e.g., [28, 31, 33]).

Beyond this, journalists and scientific researchers alike have previously examined aircraft trajectories to better understand the security and privacy of the global air traffic system. For example, reporters have used such data to uncover federal surveillance aircraft [2]. More scientific use cases in the literature include the analysis of mergers and acquisitions data of public companies [33], or the in-depth analysis of privacy leaks [31].

## 3.2 Traditional Acquisition of Aircraft Categories

The standard way to do obtain the required aircraft category uses the owner (or operator) and the aircraft model (e.g., corporate jets such as a Gulfstream), from which the use case and category of the aircraft may be inferred with good certainty. For instance, the Swiss Air Force operates aircraft for military operations, whereas business jets are likely used for corporate purposes.

However, in a large percentage of cases there is no meta information available for observed aircraft. This makes it much more difficult to identify the category of an aircraft. A recent study found that around 15% of all transponder-equipped aircraft could not be found using publicly available data [28]. Typically, these aircraft are from countries that do not provide an open aircraft registry. Furthermore, they may be outdated entries or registered recently.

## 3.3 Advantage of Behavioural Features

In order to supplement and verify unreliable and incomplete metadata collected through external sources, Classi-Fly uses exclusively behavioural features based on an aircraft's trajectories. The main advantage of this approach is that these features cannot be trivially altered or spoofed by the aircraft operator. This is contrary to any classification based on the content of their communication.

There is additional information about the capabilities of an aircraft provided by the Mode S Enhanced Surveillance (EHS) protocol features used by some aircraft. While interesting in theory, we have decided to not use these for our classification task for the following reasons: First, they, too, can easily be changed and manipulated by the aircraft operator at will. More crucially however, these communication options are not consistently used, over 50% of aircraft do not broadcast any information besides position and velocity.

Non-behavioural features, which are based on content broadcast by an aircraft, are primarily its transponder code, the call sign and the squawk code, i.e., the aircraft identifiers described in Section 2.2. The latter two are set by the pilot (often in accordance with local customs and air traffic controllers) and can thus be adapted practically at will. The transponder code is also not reliable in many situations, e.g. when the US American Air Force One presented as a different, non-existent aircraft [14].

Requirements for Classi-Fly include robustness against malicious actors and intentional manipulation of the communication data by the aircraft operator. Consequently, all features that could be manipulated were excluded in our design. This leaves only behavioural features, which cannot be altered by an aircraft at all (e.g. increasing maximum possible velocity) or at least not without significant cost in terms of time and resources (e.g., diversions, distraction flights or other large changes to the mission pattern that make an aircraft look like a different class). Please note that while the content of ADS-B messages is fundamentally not secured, the behavioural features based on location and speed can be verified reliably using physical-layer methods such as multilateration [34].

## 4 DATA COLLECTION

We now describe the processes for the collection of fine-grained tracking data and for obtaining aircraft ground truth from public sources. All data used in this work has been openly available and is thus already accessible to researchers on an ever growing scale.

## 4.1 The OpenSky Network

OpenSky is a crowdsourced network [27], which is used as the backbone of our data collection. As of January 2020, the OpenSky Network consists of more than 2500 registered sensors streaming data to its servers. The network has currently received and stored over 16 trillion ATC messages, adding over 20 billion messages by more than 50,000 different aircraft every day. As a non-profit, research-oriented network, OpenSky offers open access to its data to academic researchers and has been used for a large number of publications spanning many different domains from aviation security to climate change research. Detailed information about the history, infrastructure and use cases of OpenSky are found in [27].

Data Acquisition and Pre-Processing. Aircraft trajectories can be retrieved from the OpenSky Network for free for universities, flight authorities, and other non-profit research institutions. 2 The available data goes back several years, for which it offers dense coverage

2 https://opensky-network.org/data/impala

of Europe and the USA. More recently, it has spread to other continents, although coverage in Africa in particular is still lacking as it is based on volunteers to provide the locally broadcast aircraft communication. We obtained about 200,000 such aircraft trajectories for our ground truth and another 180,000 for the different classification categories.

The raw data is obtained from OpenSky via an Impala shell and consists of so-called state vectors , which describe the state of every observed aircraft, i.e., its position, altitude, and velocity in increments of one second. All state vectors were then separated into flights , by dividing the positional data messages received by all aircraft as follows: Each positional state which is more than 10 minutes older than the next and is at an altitude of less than 2500 m is considered an arrival state, and hence a finished flight. Note that not all flights seen by OpenSky are necessarily complete, if a flight begins or finishes outside the coverage area, the first/last message will constitute the end point of the flight. We did not differentiate between complete or incomplete flights in order to maximize the robustness of our approach. OpenSky conducts some additional processing to filter out erroneous messages and transmission-induced noise as well as potentially maliciously altered data [29].

## 4.2 Aircraft Behavioural Ground Truth

To facilitate the feature selection in the next section, we required ground truth on the average flight and movement behaviour of aircraft. We first retrieved the positional data of 9880 randomly selected aircraft seen by OpenSky in the year 2017 to be able to obtain the average values as boundaries for our features. This data was capped at maximum of 25 flights per aircraft, which resulted in more than 200,000 collected flights, with an average duration of 4669 seconds and a total number of more than 30 million analyzed state vectors. Table 1 provides the details of the ground truth dataset.

We then used these randomly selected aircraft to learn the average aircraft behaviour with regards to its flight features, which are discussed in Section 5.2. For each feature, we quantized the data set into q quantiles and learned these quantiles' specific bounds. These are then used to model the relative behaviour of different aircraft categories for our classification task.

## 4.3 Aircraft Metadata Ground Truth

There are several public sources which provide meta-information on aircraft based on their identifiers: the aircraft registration or a unique 24-bit address provided by ICAO. This typically includes the aircraft model (e.g., Airbus A320) and the owner/operator (e.g., British Airways), which we exploited to label our aircraft category ground truth.

We have used the following openly available sources to collect and verify the ground truth for our work:

- GLYPH&lt;15&gt; The OpenSky Network has recently released an aircraft database complimenting its tracking efforts with crowdsourced metadata on over 495,000 aircraft. Available here: https://opensky-network.org/aircraft-database
- GLYPH&lt;15&gt; Another non-profit project, Airframes.org, is a valuable source, offering comprehensive metadata about 609,000 aircraft identifiers. This includes background knowledge such

Table 2: Description of the evaluation data set.

| A/C Category   |   Aircraft |   Ratio [%] | Flights   | States [x1000]   |
|----------------|------------|-------------|-----------|------------------|
| Business       |       1000 |        16.6 | 36,119    | 5196             |
| Commercial     |       1000 |        16.6 | 48,590    | 12,465           |
| Fighter        |        921 |        15.3 | 6918      | 751              |
| Small Utility  |        440 |         7.3 | 16,071    | 3360             |
| Surveillance   |        403 |         6.7 | 15,384    | 4571             |
| Tanker         |        402 |         6.7 | 7657      | 1125             |
| Trainer        |       1080 |        18   | 23,778    | 5602             |
| Transport      |        768 |        12.8 | 27,808    | 5067             |
| Sum            |       6014 |       100   | 182,325   | 38,142           |

as pictures and historical ownership information (available at http://airframes.org).

- GLYPH&lt;15&gt; Furthermore, the plane spotting community actively maintains many separate databases with spotted aircraft. They usually operate SSR receivers and enrich the received data with information such as operator, model, or registration manually. The database structure of Kinetic Avionic's BaseStation software has become the de facto standard format and is also used to exchange and share their databases in forums and discussion boards. Our database version used stems from November 2017, containing 455,457 rows of aircraft data.
- GLYPH&lt;15&gt; For aircraft registered in the USA, the FAA provides a daily updated database of all owner records, online and for download. These naturally exclude any sensitive owner information but overall contain over 320,000 clean and well-organised records as of January 2018 (available at https://registry.faa.gov/aircraftinquiry/).
- GLYPH&lt;15&gt; Lastly, web services such as FlightAware and FlightRadar24 provide online access to more than a million aircraft IDs (available at http://www.flightaware.com and http://www. flightradar24.com).

When considering all these databases together, we had access to metadata for 2,180,803 unique aircraft identifiers; this snapshot for our work was taken in January 2018.

Note that these sources are naturally noisy, since they rely on compiling many separate smaller databases, are often (partly) crowdsourced and change over time; aircraft are frequently registered, de-registered and transferred globally. Due to the number of aircraft involved in the experiments in this paper we could not verify the model and operator of every aircraft by hand (i.e., by following their behaviour on web trackers and ensure consistency with the existing database). Nonetheless, this is a realistic situation for anyone looking to accurately categorize aircraft and requires an approach which is robust to such noise fluctuations.

## 4.4 Aircraft Category Extraction

Based on the trajectory data provided by OpenSky and the collected metadata, we obtained flight behaviour data for eight different aircraft categories described here in brief:

- GLYPH&lt;15&gt; Business jets: Business stakeholders typically fly jets capable of 4-20 passengers. Gulfstream's G-range, Cessna's Citation jets and Bombardier's Learjet and Challenger aircraft are amongst the most popular choices. However, this category also comprises smaller and larger aircraft as long as they are operated for business use.
- GLYPH&lt;15&gt; Small utility aircraft ('general aviation'): This aircraft group comprises a large variety of aircraft used privately and in commercial operations, which we class as so-called general aviation aircraft. Typical examples are the Cessna 172 and 182, the most sold aircraft models in the world.
- GLYPH&lt;15&gt; Commercial airliners: A large group that makes up a vast majority of passenger miles in the air. It is defined by the operator, i.e. a commercial airline that conducts scheduled transport, typically with large aircraft seating 50 or more passengers (e.g., Airbus 320 or Boeing 737).
- GLYPH&lt;15&gt; Military fighter aircraft: Fighters are designed primarily for air-to-air combat. Relatively few of these are equipped with ADS-B transponders; our group consists mainly of Eurofighters, Tornados and F15/16 aircraft.
- GLYPH&lt;15&gt; Military trainer aircraft: This category includes smaller jet and turboprop aircraft used as training vehicle for military pilots by air forces and navies around the world. Representative examples of such trainer aircraft are the Northrop T-38 Talon or the Pilatus PC-21.
- GLYPH&lt;15&gt; Military tanker aircraft: These aircraft are capable of refuelling other aircraft in the air and provide essential operational capabilities. By far the most representative example in our dataset is the Boeing KC-135 Stratotanker.
- GLYPH&lt;15&gt; Military transport aircraft: These are large aircraft used by the military to transport troops or equipment. Generally slower than aircraft intended for air fighting, they share some similarities with tanker aircraft. In our dataset these are represented mainly through the McDonnell Douglas/Boeing C-17 Globemaster III.
- GLYPH&lt;15&gt; Civil surveillance Aircraft: These aircraft are used by police agencies for surveillance purposes. They are typically small utility aircraft with special equipment and exhibit particular behaviour during their missions.

We note that these categories are not determined solely on aircraft model but instead on their use cases as defined by the operator (i.e., military or not). Indeed, there is also overlap in some military aircraft models, for example Multi Role Tanker Transport (MRTT) aircraft fulfil several roles.

Knowledge of these categories can help with a number of use cases. With the exception of the commercial airliners and small utility aircraft, all are directly potentially sensitive aircraft categories. Commercial airliners and business jets are required as input for research on economic activity (for example [23]), while the latter are also particularly interesting for investment banking studies [33]. Civil surveillance aircraft as a category have played a role in uncovering clandestine operations by state and non-state-actors [2], with small utility aircraft being the category that most of such surveillance aircraft masquerade as. In the military context, telling apart unidentified commercial aircraft from potential threats can make a difference in highly volatile situations such as the accidental

Table 3: Top origin countries of the main dataset.

Table 4: Top origin countries of unknown aircraft.

| Country     |   Aircraft |   [%] |
|-------------|------------|-------|
| USA         |       2916 |  48.5 |
| Germany     |        816 |  13.6 |
| China       |        287 |   4.8 |
| UK          |        239 |   4   |
| Australia   |        212 |   3.5 |
| Netherlands |        160 |   2.7 |
| Belgium     |        119 |   2   |
| Canada      |        110 |   1.8 |

| Country    |   Aircraft |   [%] |
|------------|------------|-------|
| UK         |        121 |  11.4 |
| Austria    |         96 |   9   |
| Germany    |         71 |   6.6 |
| China      |         67 |   6.3 |
| Czech Rep. |         59 |   5.5 |
| Ireland    |         53 |   5   |
| Australia  |         43 |   4   |
| Brazil     |         40 |   3.8 |

shooting down of Ukraine International Airlines Flight 752 in Iran in January 2020 [1]. Differentiating between the categories of fighters, tankers, trainers, and transport aircraft serves as an additional piece of information and can provide tactical advantages.

## 5 EXPERIMENTAL DESIGN

We describe the features used to determine aircraft behaviour and explain the evaluation data set used to predict aircraft categories.

## 5.1 Evaluation Data Sets

To select our evaluation data set, we first queried the full sample of aircraft seen by OpenSky in January 2018, which spanned 87,000 aircraft in total. This sample was then classified into eight different categories based on operator and model metadata (see Section 4.4).

For small utility aircraft, the available pool was larger, however, due to the fact that many surveillance aircraft share the same aircraft model (in particular Cessna 182's [2]), manual inspection of all aircraft and their tracks was required to accurately label the ground truth. For the abundant business and commercial categories, we picked random 1000 aircraft to represent their category.

We aimed to obtain 1000 aircraft per category, however, for five of the subcategories (in particular those comprising military and surveillance aircraft) there are fewer aircraft with reliable identification and the necessary transponder equipment required to obtain the detailed flight behaviour data. Thus, we picked all available aircraft for fighters, surveillance aircraft, tankers, trainer and transport aircraft.

Thus, the main data set used for our classification experiments consists of 6014 aircraft overall, each with a maximum of 50 flights. Table 2 provides the breakdown of all aircraft categories as well as the number of flights and individual state vectors used to obtain the classification features. The lowest number of flights (6918) and messages (751,000) could be obtained for the 921 fighter aircraft, presumably due to their comparatively rare use. At the upper end, the 1000 commercial aircraft were seen on 48,590 flights with over 12 million messages, illustrating the high utilization of commercial airliners. Overall, more than 185,000 flights and almost 40 million messages were processed to obtain the behavioural features. Finally, Table 3 shows the main countries of origin of our dataset, with the US making up just under half of all aircraft, followed by several European countries, China, Australia and Canada.

Table 5: Description of features, based on quantization of each behavioural feature into q parts.

| ID                       | Name               | Feature Description                                                                           | Avg. RMI           |
|--------------------------|--------------------|-----------------------------------------------------------------------------------------------|--------------------|
| Flight Level             | Flight Level       | Flight Level                                                                                  | Flight Level       |
| f 1 ; :::; f q           | Duration           | Proportion of an aircraft's flight durations falling into q quantiles.                        | 10.42%             |
| f q + 1 ; :::; f 2 q     | Bounding Box       | Proportion of a aircraft's flight areas as bounded by a box falling into q quantiles.         | 11.88%             |
| State Vector Level       | State Vector Level | State Vector Level                                                                            | State Vector Level |
| f 2 q + 1 ; :::; f 3 q   | Altitude           | Proportion of altitude values recorded for the aircraft falling into q quantiles.             | 11.84%             |
| f 3 q + 1 ; :::; f 4 q   | Heading            | Proportion of heading values recorded for the aircraft falling into q quantiles.              | 8.66%              |
| f 4 q + 1 ; :::; f 5 q   | X-Velocity         | Proportion of X-velocity values derived for the aircraft falling into q quantiles.            | 16.63%             |
| f 5 q + 1 ; :::; f 6 q   | Y-Velocity         | Proportion of Y-velocity values derived for the aircraft falling into q quantiles.            | 13.67%             |
| f 6 q + 1 ; :::; f 7 q   | Vertical Rate      | Proportion of vertical rate values recorded for the aircraft falling into q quantiles.        | 15.81%             |
| f 7 q + 1 ; :::; f 8 q   | Heading Speed      | Proportion of heading speed values derived for the aircraft falling into q quantiles.         | 13.57%             |
| f 8 q + 1 ; :::; f 9 q   | X-Acceleration     | Proportion of X-acceleration values derived for the aircraft falling into q quantiles.        | 14.53%             |
| f 9 q + 1 ; :::; f 10 q  | Y-Acceleration     | Proportion of Y-acceleration values derived for the aircraft falling into q quantiles.        | 14.67%             |
| f 10 q + 1 ; :::; f 11 q | Vertical Acc.      | Proportion of vertical acceleration values derived for the aircraft falling into q quantiles. | 15.93%             |
| f 11 q + 1 ; :::; f 12 q | Heading Acc.       | Proportion of heading acceleration values derived for the aircraft falling into q quantiles.  | 16.07%             |

Unknown Aircraft. We further obtained all features described in Section 5.2 from 1066 unknown aircraft, i.e., aircraft sending messages with identifiers where no metadata was available from any of the structured sources. We use the communication received from these identifiers to gain insights on the category of their aircraft. Naturally, we consider that there will be some noise in this dataset, which we will not be able to fully solve due to the lack of ground truth. Thanks to OpenSky's sanity checks, wronglyreceived identifiers caused e.g. by transmission or decoding errors have already been filtered out.

We have several hypotheses and explanations for the absence of these unknown aircraft from available public sources:

Based on the 24-bit identifier, if truthful, it is possible to obtain the country the aircraft is nominally registered in, by comparing it with the official ranges defined by the ICAO [17]. Table 4 shows the main countries of origin, ranging from several European countries to China, Brazil and Australia. We find that the distribution is different to the main dataset (albeit with a small sample size), in particular the lack of US aircraft is noteworthy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01061v2-figure-fc30f915-acd4-4483-b49a-f503fa4667e1.jpg)

- (1) Sensitivity: Highly sensitive military or state aircraft are excluded from public records in most countries. Depending on their missions, their country, and their use cases, hobbyist plane spotters may not be able to fill these gaps with information gleaned through traditional planespotting.
- (3) No records available: Many countries' aviation authorities do not maintain a consistent and well-kept database in the first place. In others, such as Germany, privacy regulations are extremely strict, preventing aircraft records from finding their way into the public domain.
- (2) Novel aircraft: Depending on the quality of the public or private records, aircraft in many countries take several weeks or months until they turn up in public databases.

(4) Wrongtransponder ID: Finally, there are occurrences, where the transponder ID setting of an aircraft does not match the public records, creating discrepancies in the metadata.

Figure 2: Feature correlation matrix. 0 indicates no correlation, 1 and -1 positive and negative correlation, respectively. Wecan observe strong clusters of (anti-)correlations around both acceleration and velocity features.

## 5.2 Feature Extraction

We selected 12 different features, divided into two categories: flight level and state vector features. We explain these categories in the following; a full list of the chosen features is presented in Table 5.

Flight Level Features . These features contain information about the aircraft behaviour at the highest level, namely the distribution of the durations of all its flights as well as the distribution of the area covered by the obtained flights of the aircraft. The distribution is represented using the percentages of all flights falling into the

0.25

Figure 3: Relative mutual information. Colors indicate the different physical feature groups.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01061v2-figure-7ae696d6-2dd0-4f64-ae63-0554e159f33a.jpg)

chosen number of quantiles q based on the average bounds obtained from the random sample in Section 4.2.

State Vector Features . These features contain information at the level of the collective state vectors, i.e., the distributions of all of the aircraft's message content containing the heading, velocity, vertical rate and altitude states. The distribution is based on the training dataset as described in Section 4.2 and represented as percentage of states falling into the chosen number of quantiles q .

feature. Similarly, very high acceleration and deceleration are positively correlated in all three axes (X, Y, vertical), which reflects the capabilities and actual behavior of military fighter jets.

There are three different types of state vector features based on their physical function: positional features, velocity features, and acceleration features (or the first and second derivative of the position with respect to time). Positional features include the altitude and heading values of their aircraft. The actual position in longitude and latitude values itself is not relevant, as it does not generalize to be a distinguishing feature across aircraft models and continents. Velocity features comprise the horizontal velocity in all three spatial dimensions as well as the speed with with the heading values of the aircraft change. Finally, acceleration features are derived with respect to time from all four of the velocity features.

## 5.3 Feature Analysis

Feature Correlation . Fig. 2 shows the correlation between the features calculated on the evaluation dataset. An illustrative example is given by the heading quantiles. Here the first quantile (i.e., the ratio of no to few changes in aircraft direction) is strongly negatively correlated with all other heading quantiles. This suggests that many aircraft only ever have either few changes in directions such as commercial airlines, which stay in a straight line for most of their flight duration. Aircraft that have more or many directional changes in their flights can be clearly differentiated on this

Beyond this, we can see strong relationships mainly between the horizontal velocity and acceleration features, aircraft with many values in high X-velocity and acceleration bins also exert this behaviour in the Y-direction. On the other hand, many aircraft either fall into long flights with constant middling speeds (e.g., commercial aircraft), or instead exert many very low and very high speed and acceleration values over the course of their flights, typical for fighter jets or trainer aircraft. On the other hand, few relationships can be observed on the flight durations.

Feature Quality . To obtain a clearer view on how the classification works and to identify potentially detracting features, we estimated their quality. There is a given amount of uncertainty associated with the aircraft category-its entropy. This amount depends both on the number of classes (i.e., aircraft categories) and the distribution of the samples between them. As each feature reveals a certain amount of information about the aircraft category, this amount can be measured through the mutual information (MI). In order to measure the mutual information relative to the entire amount of uncertainty, the relative mutual information (RMI) is used. RMI measures the percentage of entropy removed from the aircraft category ( cat ) when a feature ( F ) is known [5].

The RMI is defined as

$$R M I ( c a t, F ) = \frac { H ( c a t ) - H ( c a t | F ) } { H ( c a t ) } \, \quad \quad ( 1 )$$

where H A ' ' is the entropy of A and H A B ' j ' denotes the entropy of A conditional on B . In order to calculate the entropy of a feature

Figure 4: Accuracy of the classification based on the number of flights f min (left) and feature quantiles q (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01061v2-figure-f3e4febd-2896-4fa5-87d6-acbf29760fca.jpg)

it has to be discrete. As most features are continuous we perform discretization using an Equal Width Discretization (EWD) algorithm with 20 bins [9]. This algorithm typically produces good results without requiring supervision. As outliers may have a drastic effect on the RMI computation, we use the 1 st and 99 th percentile instead of the minimal and maximum values to compute the bin boundaries in order to prevent large distortions. A high RMI indicates that the feature is distinctive on its own, but it is important to consider the correlation between features as well when choosing a feature set. Additionally, features may be more distinctive when combined, even when they are not particularly useful on their own.

Table 6: Summary of model optimization and training results across four types of classifier.

Figure 3 shows the RMI for each of our selected behavioural features, the colors indicating their physical feature group (positional, velocity, acceleration, or flight level). Overall, the velocity and acceleration features (red and blue, respectively) share the most information with the aircraft category, with many of these having an RMI of 15% or more. The positional and flight level features are relatively less distinctive, which suggests that for example the distribution of heading values or the overall flight durations are more common to any aircraft mission than a consistent behavioural feature of a category. However, we choose to keep all features for our classification to produce the best results.

## 5.4 Effects of Number of Flights and Feature Quantiles

Finally, we more closely examined the effects of two feature parameters on the accuracy of the classification: the number of flights f min collected for each aircraft's feature creation, and the number of quantiles q , into which the state vector features were divided. Fig. 4 illustrates these relationships, by training a Support Vector Machine with varying values of f min and q .

The minimum number of flights required to be create a feature vector has a significant effect on classification accuracy. With no lower bound, the overall classification accuracy is fairly poor at 61%. Such a result is likely due to the classifier accounting for lots of edge cases, making it less generalized. Performance quickly increases to over 80% with 5 collected flights; increasing the number of flights per aircraft further, the accuracy increases to over 85% at 30 flights and 88.1% at 50 flights. However, by raising f min to 50, the training set size decreases substantially-we found an f min of 30 to be a reasonable balance between data set size and accuracy.

| Model Type               | Parameters                                                                          | Avg. Accuracy   | Avg. TP Rate   |
|--------------------------|-------------------------------------------------------------------------------------|-----------------|----------------|
| Support Vector Machine   | Kernel Function: Cubic C=4.795 Multiclass Method: One vs. All Standardized Input    | 86.0%           | 84.4%          |
| K-Nearest Neighbors      | Distance Metric: City Block Distance Weight: Inverse Weight # Neighbors: 4          | 84.6%           | 83.9%          |
| Decision Tree            | Criterion: Gini Index Max Splits: 1297                                              | 71.4%           | 67.9%          |
| Ensemble (Random Forest) | Method: AdaBoost Decision Tree # Learners: 402 Max Splits: 125 Learning Rate: 0.792 | 85.9%           | 84.0%          |

Table 7: Summary of metrics for evaluation runs of each classifier. Metrics are averaged across scores from each class.

|               | Accuracy   | Precision   | TP Rate   | TN Rate   |
|---------------|------------|-------------|-----------|-----------|
| SVM           | 85.3%      | 82.6%       | 84.4%     | 97.5%     |
| KNN           | 84.2%      | 81.1%       | 84.3%     | 97.2%     |
| Decision Tree | 71.4%      | 67.9%       | 65.1%     | 94.4%     |
| Random Forest | 86.7%      | 82.5%       | 86.1%     | 97.7%     |

All results were obtained with q = 10 and represent the mean of 100 classifications.

The number of feature quantiles is also related to classification accuracy. Intuitively, as the number of quantiles increases, the risk of overfitting may increase. With the minimum of q = 5 the accuracy was 84%, increasing to 87% at q = 10, and only increasing marginally thereafter until leveling off at q = 40 and 87.8%. Further increases to q = 50 show no positive effect. As such, we found 10 quantiles to be a good balance of accuracy and generality. For the analysis, f min of 30 was used, with scores averaged over 100 repetitions.

(a) SVM Evaluation Confusion Matrix

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01061v2-figure-6f2f5de6-9385-40ac-8c1d-c3b17ac10bf5.jpg)

(b) KNN Evaluation Confusion Matrix

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01061v2-figure-7afda6d7-ea01-42a6-95c9-22270102d039.jpg)

(c) RF Evaluation Confusion Matrix

Figure 5: Confusion Matrices of observations on the evaluation set for the Support Vector Machine, Random Forest and KNearest Neighbors classifiers.

## 6 MODEL SELECTION

We now compare the results of four classifiers for this classification task, implemented in Matlab: Decision Tree, Random Forests, Support Vector Machines and K-Nearest Neighbors. With the settings obtained in the previous section for f min , we retained 3519 aircraft. We trained each classifier on 80% of our training set (2859 instances), retaining 20% (713 instances) for evaluation. We used 5-fold cross-validation on our training set to reduce overfitting.

## 6.1 Classifier Training and Optimization

The results of the classification show whether aircraft categories can be distinguished purely on their movement behaviour. Based on our analysis above, we used a minimum number of flights f min = 30 and number of feature quantiles q = 10. We trained a range of model types to allow for comparison: Support Vector Machines (SVM), K-Nearest Neighbors (KNN), Decision Tree and Ensemble methods (here Random Forests, RF). To find reasonable parameters for the models, we carried out a randomized search for 30 iterations.

Decision tree is the worst performing model of the four; it performed very well in identifying commercial and training aircraft but much worse in all other categories. This could be due to the limited ability of the model to handle labels with fewer training instances, with the tree instead mainly able to classify the more frequent labels.

As shown in Table 6, SVM, KNN and RF offer similar performance during training across both accuracy and true positive rate. SVM performs slightly better across both metrics.

## 6.2 Classifier Evaluation

After training, we evaluated our classifier performance on the unseen test portion of our dataset. As shown in Tab. 7, performance is in line with performance during training in Tab. 6. As before modelling with a decision tree offered poor performance relative to the other approaches, having a high TPR for the three most populous classes and relatively low TPR elsewhere. As such we do not consider this classifier further.

The remaining three classifiers have similar performance across the metrics in Tab. 7. In Fig. 5 we can see the number of observations

predicted correctly by the SVM, KNN and RF models. As indicated by the performance metrics, each classifier performs similarly with slight biases towards certain classes. More specifically, KNN made more errors when classifying commercial aircraft than the SVM and RF models but performed better on the trainer and transport classes than the RF.

Of the three classifiers, RF has the lowest 'spread' of misclassification, i.e. misclassified business aircraft fell into four categories rather than across all the other seven. This suggests the model is better generalized than the others, which have quite a few cases of single instance misclassifications in unusual places, e.g. SVM and KNN classifying a surveillance aircraft as commercial.

All classifiers made similar classification errors. Transport aircraft were particularly susceptible to this, with transport often being predicted when the aircraft was actually business, or actual transport aircraft being misclassified. This is likely due to the low number of transport aircraft instances as well as their similarity to other aircraft movements or multi-purpose nature, i.e flying one-off or irregularly timed routes between special-purpose airports [13].

These errors also highlight the distinctiveness of some aircraft compared to others. We can see that trainer aircraft are rarely misclassified and when they are, the predicted class falls into one of four other labels. Business aircraft, however, have misclassifications across a range of labels. This could be a result of business aircraft being used for a wide range of purposes which in turn might result in flights similar to other categories.

If a single classifier is needed, the RF or SVM models provide equally good across-the-board performance in comparison to KNN. However, any of the three models would perform quite well and could be used to construct a meta-classifier. Further training examples in the lower population classes would help to explore whether certain classifiers perform better for some classes than others, helping to better assess the benefit of a meta-classifier. However, this would need to be done with care as individual classifiers good at identifying certain types of aircraft might be outvoted.

## 7 ANALYSIS OF UNKNOWN AIRCRAFT

Table 8 shows the classification of approximately 1000 aircraft, about which there was no data available in any publicly accessible database at the time of our snapshot. All selected aircraft had at least 10 flights and 500 state vector data points available for their feature creation, to reduce the amount of noise to a minimum and ensure that these are consistently used aircraft identifiers. To obtain categories for these aircraft, we used the random forest classifier trained on the known aircraft data as described above. As an ensemble classifier it provides confidence scores, i.e., the percentage of times a sample has been classified as a particular category. We used these scores as a cut off threshold, i.e., any sample classified with a score of less than 0.5 in any of the eight classes was judged as too low to provide useful insights. Taking this into account, 52.3% of all aircraft were classified confidently into one group. Table 8 shows the full results.

The commercial aircraft could overwhelmingly be verified manually using the most current online source, FlightRadar24, as having been put into service after the time our metadata snapshot was

Table 8: Classification of unknown aircraft.

| Aircraft Category   | Aircraft   | Percentage   |
|---------------------|------------|--------------|
| Business            | 116        | 10.9%        |
| Commercial          | 316        | 29.6%        |
| Fighter             | -          | -            |
| Small Utility       | 49         | 4.6%         |
| Surveillance        | 74         | 6.9%         |
| Tanker              | -          | -            |
| Trainer             | 2          | 0.2%         |
| Transport           | -          | -            |
| Other               | 509        | 47.7%        |
| Sum                 | 1066       | 100 %        |

Figure 6: Example of seven flight trajectories from a previously unknown surveillance aircraft detected in Croatia. Each colour is an individual flight and all flights clearly exhibit ''circling' features exhibited by surveillance aircraft.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.01061v2-figure-11607b08-1931-4e07-ab21-286748eeb03f.jpg)

taken in January 2018. Indeed, of the 316 aircraft, 305 were classified correctly, with the 11 misclassifications being larger business jets. The new airliners in this set included, for example, 9 Boeing Dreamliners delivered to Norwegian in the first half of 2018 [7] or new aircraft in China, the biggest growth market for commercial aviation.

We further find that a large number of aircraft are seen by the classifier as business and small utility aircraft (10.9% and 4.6% respectively). This is plausible, as information on such private aircraft is not necessarily well-publicized, potentially even sensitive, and many countries other than the English-speaking world either do not require such aircraft to be on a public register or even do not publish any aircraft register at all. While we can naturally still not verify the accuracy of the classification, many such classified aircraft are regulars at typical business airports (e.g., Farnborough, UK or Teterboro, US), improving our confidence. The final large group was made up of surveillance aircraft (6.9%), whose sensitivity provides a clear motivation for not publishing their meta information. We discuss a detailed case study on such an aircraft in the next section. There was a small minority of aircraft classified as trainer aircraft (0.2%). Finally, no military fighter, transport, or tanker aircraft were found in this dataset.

Detection of Surveillance Aircraft. We now take an example case study of 'interesting' aircraft categories being detected by ClassiFly. As first laid out in [2] for the case of the United States, there are many federal surveillance operations conducted with undercover aircraft. Such aircraft register with inconspicuous call signs, registrations and transponder identifiers, thus evading all databasebased detection and making behavioural classification necessary. In the case of surveillance aircraft, some of the defining features can be visualized well just with the basic trajectory data. Most notably, such aircraft do not fly from point A to point B in a relatively straightforward fashion. Instead, they fly to a target area (say, a border area or an ongoing high-profile criminal incident) where they circle steadily in search of persons or objects, shifting altitude and center position occasionally.

This case study shows that our approach generalizes across different countries and their surveillance institutions and is able to detect surveillance aircraft around the globe.

Fig. 6 visualizes this precise behavior. It shows seven flights from an example classified with very high confidence (RF score of 0.91) as surveillance aircraft. While no information about this aircraft is available, as it does not appear in any database, it clearly exhibits the patterns of an aircraft used for surveillance of a narrow area, which are picked up by the classifier. The number of flights with this nature clearly increases confidence in the correctness of this prediction. On further analysis, we found that these flights were conducted in Croatia and we speculate that they are related to military and anti-terrorist missions.

## 8 DISCUSSION

We now discuss the limitations of Classi-Fly, possible refinements, and potential countermeasures to our approach.

## 8.1 Limitations

The greatest limitation of Classi-Fly is the inherent non-specificity of some categories. For example, it is difficult to identify the precise use case of a business aircraft; besides business travel, the same Gulfstream G550 could be used for transport of goods for the military or people for leisure, which can pose a potential threat to the validity of the results. Likewise, there may be differences between countries and institutions not captured in these categories and the ground truth. As long these do not pertain to behaviour, our methods work well but it is conceivable that different global regions or military institutions may also exhibit differences in behaviour with their aviation hardware.

However, with further research into potential subcategories and howtodefine them based on metadata such as the operator or owner or the airports frequented, this could be mitigated and their different behavior learned. This applies also to currently neglected aircraft categories such as unmanned aerial vehicles (UAV, or drones) and ultralight aircraft (ULAC), which will become transponderequipped in larger numbers in the future and are a major interest for researchers and the aviation industry.

## 8.2 Refinement

Besides improving the ground truth metadata on the examined aircraft categories, other, non-behavioral features can be integrated into Classi-Fly. As many wireless standards (not only in aviation) give manufacturers a large amount of freedom over the actual softand hardware implementations, differences emerge that can be used as classification features.

Besides these approaches, it is possible to add a host of features derived from the actual message content sent out by the aircraft. In a non-adversarial setting where the aircraft operators do not actively seek to obfuscate their identity (beyond excluding it from public databases), this would greatly improve classification accuracy.

On the physical layer, [19] proves that it is feasible to distinguish aircraft transponders based on anomalies in the frequency stability of their messages. On the data link layer, research has exploited differences in the transponders' random backoff algorithms [32].

Overall, we assume that certain uncommon aircraft may be individually identifiable through the combination of features. Future work will thus consider the possible granularity that several approaches can provide if they are combined and further quantify the privacy impact for aircraft owners and operators.

## 8.3 Countermeasures

As our classification approach is agnostic to any non-behavioral features, it is more difficult to apply effective countermeasures against it. Aircraft could deliberately change their behavior to avoid detection and classification. However, this forces the aircraft into not being able to fulfil its intended function freely, for example surveillance aircraft not circling their target, or military fighter jets deliberately flying slowly. This limits the potential benefit of such an option.

The most effective countermeasure as concluded by the literature consists of the randomisation of the aircraft's ICAO identifier, making it difficult to continuously track the same aircraft over time. If done globally for all aircraft, and in conjunction with other pseudonymisation measures regarding the registration, it could effectively thwart consistent aircraft tracking and by extension also Classi-Fly. However, the cat may largely be out of the bag already; with the current widespread availability of comprehensive aviation data there is sufficient input available for training.

Related work [33] has looked at countermeasures to the basic enabling mechanisms of aircraft tracking, which is generally based on the ICAO identifier or other directly identifying information broadcast voluntarily by the aircraft (such as its registration). There are two popular privacy-preserving approaches to aircraft tracking found in the aviation industry: the first consists of not displaying aircraft on popular web feeds (such as FlightRadar24 or FlightAware), the second comprises the use of shell companies to hide the real owners of an aircraft and thus undermine the collection of accurate metadata. Both ideas, while certainly popular, are ultimately not effective against a moderate threat model that allows for the independent localization of aircraft [33].

## 9 RELATED WORK

The classification of objects or subjects based on wireless communication has been a popular field of research, in particular with a focus on security and privacy aspects. Exemplary studies outside the aviation domain range from the mobility states of humans [24]

to the classification of intruders (people, soldiers, vehicles) in a military setting [3].

Regarding such trajectories, the authors in [21] used GPS-based tracks of cab drivers to study their behavior and classify them into high-earning and average-earning drivers through the use of angularity and travel time features. Using taxi tracks with a different focus, further work attempted to uncover anomalous trajectories in a dataset by comparing and isolating tracks which are few and different from the majority [39].

The closest related academic research is the classification of different types of ground vehicles. Vehicle type classification is an important signal processing task with widespread military and civilian applications in intelligent transportation systems [10]. Several data types have been used for vehicle classifications, collected for example from acoustic or seismic [30, 38] sensor sources. While these may be applicable in the aviation domain, too, our work focuses on the trajectories for classification.

Most closely related in the vehicle domain, the authors in [35] distinguish two classes of vehicles (trucks and passenger cars) using GPS data extracted from mobile traffic sensors with a misclassification rate of 4.6%. The main features are based on the vehicles' acceleration and deceleration behavior. Our work transfers this idea into three dimensions and applies it to the very different speeds and vehicle types found in aviation.

Finally, the closest non-academic work related to our approach is the successful attempt of investigative journalists to uncover unknown surveillance aircraft in the USA, which was presented at DEFCON 25 [15]. The authors report on the background of socalled spy aircraft, which are identified using a machine learning approach on aircraft flight data pre-processed by a large commercial tracking website. While we follow a similar basic approach concerning such surveillance aircraft in this work, we systematically analyze the effectiveness and validity of applying machine learning to aircraft behavior. In order to do this, we process a large open data set, and discuss requirements on features and number of flights. Furthermore, we generalize this approach to many aircraft categories.

In the aircraft domain, wireless classification has focused on traditional non-cooperative PSR communication as the medium. Such work exists for both military [20] and commercial aircraft [40] and exploits for example Doppler signatures [8] and high resolution range profiles [40] to identify the type of aircraft seen by the radar. However, primary radars are prohibitively expensive and thus widely inaccessible for research. As they are being replaced globally with the more accurate and cost-efficient ADS-B, we choose to focus on this cheap and openly available source of aircraft trajectories.

## 10 CONCLUSION

In this work, we presented Classi-Fly, a machine learning method to infer the categories of aircraft, both anonymous and known, based purely on their movement behavior. We validate our approach using publicly available flight data, comprising several hundred thousand flights with tens of millions of states in conjunction with meta information obtained from publicly available aircraft registries. Our results show that we can obtain the category of an aircraft with a likelihood of almost 90%, based on features obtained from 30 flights or fewer. In cases where no metadata is publicly available for an aircraft, we show that our approach can be used to create this data, which is necessary for many research projects based on air traffic communication. Finally, we have examined a case study showing that it is possible to automatically discover sensitive aircraft in a large data set using Classi-Fly, including police, surveillance and military aircraft.

Future work in this area can focus on defining even more aircraft categories relevant for new and existing research applications. Specifically, we plan to extend the existing categories to include UAV and other non-standard aircraft such as gliders or ULAC. This requires these aircraft categories to have sufficiently broad equipage with ADS-B transponders or alternatives such as FLARM but can be a worthwhile expansion of coverage in many countries where FLARM is popular. 3 Differences in behaviour across institutions and global regions could also be of much further interest.

## REFERENCES

- [1] Ahmad Ali Abin, Shahabedin Nabavi, and Mohsen Ebrahimi Moghaddam. 2020. It is time for AI to prevent unintentionally disastrous human errors. Authorea (June 2020). https://doi.org/10.22541/au.159162125.57899841
- [3] Anish Arora, Prabal Dutta, Sandip Bapat, Vinod Kulathumani, Hongwei Zhang, Vinayak Naik, Vineet Mittal, Hui Cao, Murat Demirbas, Mohamed Gouda, et al. 2004. A line in the sand: a wireless sensor network for target detection, classification, and tracking. Computer Networks 46, 5 (2004), 605-634.
- [2] Peter Aldhous. 2017. BuzzFeed News Trained A Computer To Search For Hidden Spy Planes. This Is What We Found. Buzzfeed News (Aug. 2017). https://www. buzzfeed.com/peteraldhous/hidden-spy-planes
- [4] Bank of England, Monetary Policy Committee. 2020. Monetary Policy Report . Technical Report. https://www.bankofengland.co.uk/-/media/boe/files/ monetary-policy-report/2020/may/monetary-policy-report-may-2020.pdf
- [6] Isaac I Bogoch, Alexander Watts, Andrea Thomas-Bachli, Carmen Huber, Moritz UG Kraemer, and Kamran Khan. 2020. Potential for global spread of a novel coronavirus from China. Journal of travel medicine 27, 2 (2020), taaa011.
- [5] Roberto Battiti. 1994. Using mutual information for selecting features in supervised neural net learning. IEEE Trans. on Neural Networks 5, 4 (1994), 537-550.
- [7] Breaking Travel News. 2018. Norwegian sees increase in passenger numbers for April. http://www.breakingtravelnews.com/news/article/ norwegian-sees-increase-in-passenger-numbers-for-april/
- [9] James Dougherty, Ron Kohavi, and Mehran Sahami. 1995. Supervised and unsupervised discretization of continuous features. In Machine Learning Proceedings 1995 . Elsevier, 194-202.
- [8] Barry D Bullard and Patrick C Dowdy. 1991. Pulse Doppler Signature of a Rotary-wing Aircraft. IEEE Aerospace and Elec. Systems Mag. 6, 5 (1991), 28-30.
- [10] Marco F Duarte and Yu Hen Hu. 2004. Vehicle classification in distributed sensor networks. J. Parallel and Distrib. Comput. 64, 7 (2004), 826-838.
- [12] Dario Garcia-Gasulla, Sergio Alvarez Napagao, Irene Li, Hiroshi Maruyama, Hiroki Kanezashi, Raquel Perez-Arnal, Kunihiko Miyoshi, Euma Ishii, Keita Suzuki, Sayaka Shiba, et al. 2020. Global Data Science Project for COVID-19 Summary Report. arXiv preprint arXiv:2006.05573 (June 2020). https://arxiv.org/ abs/2006.05573
- [11] Federal Aviation Administration. 2020. FAA Registry - Aircraft Inquiry. https: //registry.faa.gov/aircraftinquiry/
- [13] Thomas L Gibson. 2002. The Death of 'Superman': The Case Against Specialized Tanker Aircraft in the USAF . Technical Report. Air University, Maxwell, Alabama.
- [15] Jason Hernandez, Sam Richards, and Jerod MacDonald-Evoy. 2017. Tracking Spies in the Skies. Presented at DEFCON 25.
- [14] Dareh Gregorian. 2018. Trump's first trip to a combat zone was secret fi/exclam\_question except on Twitter. NBC News (Dec. 2018). https://www.nbcnews.com/politics/donald-trump/ trump-s-first-trip-combat-zone-was-secret-except-twitter-n952136
- [16] Stefano Maria Iacus, Fabrizio Natale, Carlos Satamaria, Spyridon Spyratos, and Michele Vespe. 2020. Estimating and Projecting Air Passenger Traffic during the COVID-19 Coronavirus Outbreak and its Socio-Economic Impact. arXiv preprint arXiv:2004.08460 (April 2020). https://arxiv.org/abs/2004.08460

3 FLARM is a cooperative low-cost collision avoidance system developed for the gliding community [26]; its signals are collected by some web trackers.

|      | Registration of Aircraft Addresses with Mode S Transponders . Technical Report NACC/DCA/3, WP/05.                                                                                                                                                                                                                                                                                                                          |
|------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------|
| [18] | Punta Cana, Dominican Republic. Thomas Lecocq, Stephen P Hicks, Koen Van Noten, Kasper van Wijk, Paula Koelemeijer, Raphael SM De Plaen, Frederick Massin, Gregor Hillers, Robert E Anthony, Maria-Theresia Apoloner, et al. 2020. Global quieting of high-frequency seismic noise due to COVID-19 pandemic lockdown measures. Science (Jul 2020). https://science.sciencemag.org/content/early/2020/07/22/science.abd2438 |
| [19] | Mauro Leonardi, Luca Di Gregorio, and Davide Di Fausto. 2017. Air Traffic Security: Aircraft Classification Using ADS-B Message Phase-Pattern. Aerospace 4, 4 (2017), 51.                                                                                                                                                                                                                                                  |
| [20] | H Lin and AA Ksienski. 1981. Optimum frequencies for aircraft classification. IEEE Trans. Aerospace Electron. Systems 5 (1981), 656-665.                                                                                                                                                                                                                                                                                   |
| [21] | Liang Liu, Clio Andris, and Carlo Ratti. 2010. Uncovering cabdrivers behavior patterns from their digital traces. Computers, Environment and Urban Systems 34, 6 (2010), 541-548.                                                                                                                                                                                                                                          |
| [22] | Liang Mao, Xiao Wu, Zhuojie Huang, and Andrew J Tatem. 2015. Modeling monthly flows of global air travel passengers: An open-access data resource. Journal of Transport Geography 48 (2015), 52-60.                                                                                                                                                                                                                        |
| [23] | Sam Miller, Helen Susannah Moat, and Tobias Preis. 2020. Using aircraft location data to estimate current economic activity. Scientific reports 10, 1 (2020), 1-7.                                                                                                                                                                                                                                                         |
| [24] | M Mun, Deborah Estrin, Jeff Burke, and Mark Hansen. 2008. Parsimonious Mobility Classification Using GSM and WiFi traces. In Proceedings of the Fifth Workshop on Embedded Networked Sensors (HotEmNets) .                                                                                                                                                                                                                 |
| [25] | Xavier Olive, Martin Strohmeier, and Jannis Lbbe. 2020. Crowdsourced air traffic data from The OpenSky Network 2020. Zenodo (may 2020). https://doi.org/10. 5281/zenodo.3782598                                                                                                                                                                                                                                            |
| [26] | Christoph G Santel, Paul Gerber, Simon Mehringskoetter, Verena Schochlow, Joachim Vogt, and Uwe Klingauf. 2014. How Glider Pilots Misread the FLARM Collision Alerting Display: A Laboratory Study. Aviation Psychology and Applied Human Factors 4, 2 (2014), 86.                                                                                                                                                         |
| [27] | Matthias Schaefer, Martin Strohmeier, Vincent Lenders, Ivan Martinovic, and Matthias Wilhelm. 2014. Bringing Up OpenSky: ALarge-scale ADS-B Sensor Net- work for Research. In Proceedings of The 13th ACM/IEEE International Conference on Information Processing in Sensor Networks (IPSN) . 83-94.                                                                                                                       |
| [28] | Matthias Schaefer, Martin Strohmeier, Matthew Smith, Markus Fuchs, Vincent Lenders, Marc Liechti, and Ivan Martinovic. 2017. OpenSky Report 2017: Mode S and ADS-B Usage of Military and Other State Aircraft. In Digital Avionics Systems Conference (DASC), 2017 IEEE/AIAA 36th . IEEE, 1-10.                                                                                                                            |
| [29] | Matthias Sch¨fer, a Martin Strohmeier, Matthew Smith, Markus Fuchs, Vincent Lenders, and Ivan Martinovic. 2018. OpenSky report 2018: assessing the integrity of crowdsourced mode S and ADS-B data. In 2018 IEEE/AIAA 37th Digital Avionics Systems Conference (DASC) . IEEE, 1-9.                                                                                                                                         |
| [30] | James F Scholl, Loren P Clare, and Jonathan R Agre. 1999. Seismic attenuation characterization using tracked vehicles . Technical Report. Rockwell International Corp Thousand Oaks Ca Science Center.                                                                                                                                                                                                                     |
| [31] | Matthew Smith, Daniel Moser, Martin Strohmeier, Vincent Lenders, and Ivan Mar- tinovic. 2018. Undermining privacy in the aircraft communications addressing and reporting system (ACARS). Proceedings on Privacy Enhancing Technologies 2018, 3 (2018), 105-122.                                                                                                                                                           |
| [32] | Martin Strohmeier and Ivan Martinovic. 2015. On Passive Data Link Layer Fingerprinting of Aircraft Transponders. In Proceedings of the First ACMWorkshop on Cyber-Physical Systems-Security and/or PrivaCy (CPS-SPC) . ACM, 1-9.                                                                                                                                                                                           |
| [33] | Martin Strohmeier, Matthew Smith, Vincent Lenders, and Ivan Martinovic. 2018. The Real First Class? Inferring Confidential Corporate Mergers and Government Relations from Air Traffic Communication. In IEEE European Symposium on .                                                                                                                                                                                      |
| [34] | Security and Privacy (EuroS&P) Martin Strohmeier, Matt Smith, Matthias Sch¨fer, a Vincent Lenders, and Ivan Martinovic. 2017. Crowdsourcing security for wireless air traffic communications. In Cyber Conflict (CyCon), 2017 9th International Conference on . IEEE, 1-18.                                                                                                                                                |
| [35] | Zhanbo Sun and Xuegang Jeff Ban. 2013. Vehicle Classification Using GPS Data. Transportation Research Part C: Emerging Technologies 37 (2013), 102-117.                                                                                                                                                                                                                                                                    |
| [36] | Toyotaro Suzumura, Hiroki Kanezashi, Mishal Dholakia, Euma Ishii, Sergio Al- varez Napagao, Raquel Prez-Arnal, and Dario Garcia-Gasulla. 2020. The Impact of COVID-19 on Flight Networks. arXiv preprint arXiv:2006.02950 (June 2020). arXiv:2006.02950 [cs.SI] https://arxiv.org/abs/2006.02950                                                                                                                           |
| [37] | United Nations Department of Economic and Social Affairs. 2020. Us- ing experimental statistics to monitor of the impact of COVID-19 in Denmark. https://covid-19-response.unstatshub.org/data-solutions/ using-experimental-to-monitor-the-impact-of-covid19-in-denmark/                                                                                                                                                  |
| [38] | Huadong Wu, Mel Siegel, and Pradeep Khosla. 1998. Vehicle sound signature recognition by frequency vector principal component analysis. In Instrumentation and Measurement Technology Conference, 1998. IMTC/98. Conference Proceedings. , Vol. 1. IEEE, 429-434.                                                                                                                                                          |
| [39] | IEEE Daqing Zhang, Nan Li, Zhi-Hua Zhou, Chao Chen, Lin Sun, and Shijian Li. 2011. iBAT: Detecting Anomalous Taxi Trajectories from GPS Traces. In Proceedings of the 13th International Conference on Ubiquitous Computing . ACM, 99-108.                                                                                                                                                                                 |

- [40] Anthony Zyweck and Robert E Bogner. 1996. Radar Target Classification of Commercial Aircraft. IEEE Trans. on Aerospace and Elec. Sys. 32, 2 (1996), 598-606.
Martin Strohmeier Matthew Smith Vincent Lenders Ivan Martinovic
Classi-Fly: Inferring Aircraft Categories from Open Data using Machine Learning
In recent years, air traffic communication data has become easy to access,
enabling novel research in many fields. Exploiting this new data source, a wide
range of applications have emerged, from weather forecasting to stock market
prediction, or the collection of information about military and government
movements. Typically these applications require knowledge about the metadata of
the aircraft, specifically its operator and the aircraft category.
  armasuisse Science + Technology, the R\&D agency for the Swiss Armed Forces,
has been developing Classi-Fly, a novel approach to obtain metadata about
aircraft based on their movement patterns. We validate Classi-Fly using several
hundred thousand flights collected through open source means, in conjunction
with ground truth from publicly available aircraft registries containing more
than two million aircraft. Classi-Fly obtains the correct aircraft category
with an accuracy of over 88%, demonstrating that it can improve the meta data
necessary for applications working with air traffic communication. Finally, we
show that it is feasible to automatically detect specific flights such as
police and surveillance missions.