InTech uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Engineering » Aerospace Engineering » "Advances in Air Navigation Services", book edited by Tone Magister, ISBN 978-953-51-0686-9, Published: August 1, 2012 under CC BY 3.0 license. © The Author(s).

# How to Manage Failures in Air Traffic Control Software Systems

By Luca Montanari, Roberto Baldoni, Fabrizio Morciano, Marco Rizzuto and Francesca Matarese
DOI: 10.5772/48685

Article top

## Overview

Figure 1. Aircraft representation

Figure 2. Collision cylinder definition

Figure 3. Conflict cylinder definition.

Figure 4. Modelling Cylinders Orientation

Figure 5. Impact Plane,Conflict and Collision cylinder.

Figure 6. Impact Plane, Collision Area, Conflict Area and Projected CPA definition.

Figure 7. Changes in the CPA coordinates due to projecting errors..

Figure 8. histogram of projected horizontal and vertical separations at the CPA (31 days of radar data)

Figure 9. Pa estimation for different CPAp. Aircraft established at a defined flight level or vz equals to zero (upper) and aircraft with vz close to zero (lower).

Figure 10. Pa estimation for different CPAp. Aircraft climbing/descending and different vz / vx ratios. vz / vx=0.1(upper), vz / vx=20 (lower)

# Probability of Potential Collision for Aircraft Encounters in High Density Airspaces

R. Arnaldo, F.J. Sáez, E. Garcia and Y. Portillo

## 1. Introduction

Collision risk estimation in airspace and mathematical modeling of mid-air collisions have been carried out for over more than 40 years [1]. During this period there has been a development of mathematical models for processes leading to possible collisions of aircraft flying nearby in order to estimate the risk of collision.

B. L. Marks [2] developed the principles in which a collision risk model could be developed in the early 1960s. Marks' work was modified and enhanced by P. Reich [3] and that model, later called the Reich model, has been the basis for many of the important developments in this field.

The Reich model uses information related to the probabilistic distributions of aircraft's lateral and vertical position, traffic flows of the routes, aircraft's relative velocities and aircraft dimensions to generate estimation of collision risk. However, this model does not cover adequately situations where ground controllers monitor the air traffic through radar surveillance and provide tactical instructions to the aircraft crews. Furthermore, the problem of collision risk modeling in the analysis of “high traffic density” ATC scenarios is different to that of “procedural scenarios”, which have been developed by Reich [4] and Brooker [5], amongst others. This is mainly due to the active role of Controllers in the first case. In this case positive control is used extensively to modify the planned aircraft route. This requires the inclusion in the model of “human factor response” behavior.

These "collision risk models" were initially applied in the 60s to determine safe separation standards between pairs of aircraft flying at the same altitude on parallel courses over the North Atlantic Ocean [6]. Since then, new models have been developed and continually refined and improved. They have been applied for different geographic regions (USA [7], European airspace [8]), for oceanic or radar [9] environments, and different flight regimes (for example, high-altitude cruise and landing on close [10,11,12] and ultra close spaced runways [13,14]), for specific flight phases [15] (focused for example on the separation between aircraft on final approach and landing, when flight risks are greater than during any other phase of flight), for different types of separation (vertical, longitudinal and lateral) and also for current and future operational concepts [16], such as free flight [17], airborne self separation [18],….

Most of these models, amongst them the formula proposed by Brooker [19] for mid-air collision risk, involve the aggregation of terms comprising different factors related to: initiating events which produce defective flight paths; the probability of safety defenses correcting these defective flight plans; and traffic and kinematic scalers. But, as he indicates: “it does no more than spell out the mechanisms by which collisions logically have to occur. The hard problem is how to populate the parameters in the formulation with sensible numbers”.

Risk models have also been developed for the estimation of conflict probability (understood as the probability that the distance between a pair of aircraft becomes smaller than some specified minimum separation value). Paielli and Erzberger’s [20,21] emphasis was on the development of algorithms to numerically evaluate approximations of conflict probabilities. Prandini et al. [22,23] emphasized the analysis of the problem and distinguished three sub-problems of evaluating conflict probability.

The main point of conflict probability is its clear relation to a well known safety criterion in civil aviation: the separation minimum, which puts a requirement on the air traffic management system; not to let aircraft come closer to each other than a certain minimum distance. In addition to minimum separation values, ICAO (International Civil Aviation Organization) has also defined limiting criteria for acceptable risk levels of fatal accident, and in particular, for the risk of mid air collision [24 ]. The allowed probability values for such events are of the order of one mid-air collision or physical crossing per 10^9 flight hour.

Furthermore, some effort has been also devoted to the problem of aircraft conflict detection. An excellent survey of the different conflict detection and resolution schemes has been carried out by Kuchar [25,26], where the conflict detection schemes are classified according to the modeling method used for projecting the aircraft position in the future.

According to Brooker [9], mid-air collisions derived from radar inaccuracies are very rare, so to estimate their frequency, it is necessary to model the factors that might lead to such events. But this extremely low value makes it difficult to obtain reliable empirical results from reasonably computational amount of data.

As collisions are very unlikely events most of the previous approaches to estimate collision risk were centered on simulations techniques applicable to rare event estimations such as Montecarlo simulations [27,28]. Nevertheless, simulations are not enough, as the components of the collision models have to be verifiable, i.e. match reality, and cautious. ‘Verifiable’ in the present context means that the model description can be demonstrated to match what happens in practice, and that most of the parameters in the model can be measured directly by analyzing air traffic patterns.

Some authors, like Dr. L. Burt [29], have formulated expressions that attempt to estimate Pa, distinguishing four different aircraft encounters geometries. The mathematical formulas are customized for these geometries so they are only applicable for circumstances that they have been developed for. They barely provide an estimate of the average conditional probability of collision Pa but they do not provide an individual value of Pa for each encounter. Therefore, this approach does not assess the severity of each individual potential encounter.

Other authors, such as Campos [30] have calculated the probability of coincidence for aircraft on arbitrary straight flight paths (either climbing, descending, or in level flight) with constant speed as an upper bound for the probability of collision. Although in this approach the time and distance of closest approach are used to calculate the position for maximum probability of coincidence. In reference [31] same authors illustrate the relationship between the aircraft RMS (Root Mean Square) position error and the minimum separation distance for achieving a certain Target Level of Safety (TLS) for low probability of collision.

Nevertheless, most of the researches on this field have just worked in the estimation of probabilities of conflict (before deliberate actions are taken to solve the conflicts) and how these probabilities depend on aircraft separation standards. Different current and future Air Traffic Management operational concepts have been studied under this perspective in an attempt to reduce aircraft separation standards [32,33] or with the aim of designing proper avoidance maneuvers in order to maintain the prescribed minimum separation standards among aircraft [34,35].

The previous considerations give an idea of the complexity of using stored aircraft tracks, within a given scenario and time frame, to infer safety level, collision risk probability and associated system weaknesses. In most high density airspace scenarios recorded tracks can be obtained for all aircraft flying in it, for example, from Radar Data Processing systems (RDP). In fact, this provides us with a robust data source, which could be used for safety analysis. This could include indirect information which is closely related to the “human factor response”. Despaite its importance not much effort have been devoted to the development of risk and collision models based upon the analysis of the stored aircraft tracks.

Furthermore, it has to be considered that the distribution of aircraft position errors over their intended tracks is one of the most important factors in determining route safety, and consequently it has been broadly studied. Reference [36], for instance, presents a modeling technique to compute the probability density function of position errors as the aircraft proceed along the route taking into account not only the time dependence, but also all the factors influencing an aircraft's position errors, e.g., surveillance and navigation errors, surveillance fix rate, and air traffic control procedures.

Following the research line initiated on [31,37,38] by the mentioned previous work, the authors are developing a more detailed mathematical model for both components of probability of collision in a radar ATC (Air Traffic Control) environment.

## 2. Fundamentals behind probability of collision estimation

Jaroslav Krystul [39] defines the risk as the probability of a particular adverse event occurring during a stated period of time. Usually, this is an event occurring when the system reaches a particular critical state. These events with a very small probability of occurrence are called rare events. Applying this definition to an ATC scenario, it is accepted that risk is closely related to those situations in which two aircraft are on conflict course and would not only pass closer than the prescribed horizontal and vertical separation minima but which would, in fact, collide.

The work presented here was originally inspired by the principle stated in [2] by B. L. Marks: “… the task of relating collision risk to a traffic configuration can be taken in two parts:

parts:

1. Determining the frequency with which aircraft are exposed to risk by passing close together; and

2. Determining what chance of collision is inherent in the passing”.

According to this idea, the probability of aircraft collision can be expressed as:

 P(collision)=FeR∗P(pot.coll/pot.conf)∗P(coll/pot.coll). (1)

where:

 FeR, Frequency of exposition to Risk, here is considered as the relative frequency that an aircraft would potentially violate the separation standards defined for the particular situation, here referred to as potential conflict. It is easily seen that this value increases with the traffic density.

 P(pot.coll/pot.conf) is the conditional probability of a potential collision (pot.coll) between two aircraft that have previously violated the separation standards (pot. conf). Its value depends on the encounter kinematics and uncertainties associated to predicted positions. It represents the intrinsic severity of the encounter and it is independent of the traffic density.

 P(coll/pot.coll ) is the conditional probability of collision among potential collisions having failed all the safety barriers (ATC, TCAS) which are in place to mitigate the risk.

A time horizon is established within which all aircraft positions are projected to explore existence of “potential conflicts”. In the following discussion 10 minutes look ahead time has been considered. Accordingly, the relative frequency of potential collisions among potential conflicts F(pot.coll/pot.conf) could be expressed as:

 F(pot.coll/pot.conf) =Num.ofpot.collisionsNum.ofpot.conflicts≈E[Pa]. (2)

where Num.pot.collisions is the number of aircraft that are about to collide (and will do if all safety barriers fail).

An initial expectation for probability of potential collision among potential conflicts, E(Pa), could be obtained as the relative frequency that two aircraft, on a conflict course, would not only pass closer than the prescribed horizontal and vertical separation minima, but would in fact collide. This expression provides an expected, or global, value and does not assess the severity of each individual potential encounter itself. This chapter proposes an approach to estimate the severity of the encounter using the conditional probability of a potential collision Pa for each particular aircraft encounter. This proposed approach aims at improving the previous works by:

• Providing an individual probability of collision for each individual encounter based on the: (1) geometry and kinematics of the encounter, (2) the minimum predicted lateral separation at the CPA, and (3) the minimum predicted vertical separation at the CPA.

• Taking into consideration the radar data errors and the segmentation errors.

### 2.1. Consideration of aircraft protection zones

As stated by Ennis [41], a protected zone represents a region around a given aircraft that no other aircraft should penetrate.

A simplification of the Bellantoni [42] approach for the definition of a collision surface can be made by modelling the aircraft as a cylinder of diameter xy and height z as indicated in figure 1.

### Figure 1.

Aircraft representation

Two aircraft are taken as colliding if their cylinders touch. With this bounded and closed airspace region representing the aircraft, a “collision cylinder” can be defined as a larger cylinder of twice the dimensions represented in figure 1, and defined by height 2z and radius 2xy (see figure 2).

On the other hand, all high density traffic ATC scenarios have established minimum separation standards defined by two values, the minimum horizontal (R) and vertical (H) separations. When two aircraft are closer than these distances the ATC system is considered to have failed. These values (R, H) allow us to use another cylinder shaped protection model for all aircraft which should be free of any other aircraft to fulfil this separation minima (see figure 3). This volume will be called the “conflict cylinder” as it is considered that two aircraft potentially violating these separations are exposed to risk.

### Figure 2.

Collision cylinder definition

During the en route phase of flight, for example, the conflict cylinder would be 5 nm in radius and 2,000 ft in height. However, these current minimum separation standards were determined many years ago and the method by which they were calculated is not well documented. Recently, Reynolds & Hansman [42] identified factors involved in defining the aircraft separation standards and discussed the importance of accurate state information for controllers in maintaining them. Ennis & Zhao [43] examined the physical compositions of the protected zone and presented a formal approach to the analysis of minimum separation standards.

### Figure 3.

Conflict cylinder definition.

A summary of the modelling cylinders defined so far is presented in the following table.

When civil aircraft are climbing or descending, it is considered that pitch angles are small and so, vertical and horizontal dimensions have small changes. Therefore, all the “modelling cylinders” will be considered as horizontal, as indicated on figure 4.

As all the cylinders are considered parallel, the longitudes and surfaces ratios among them will be constant when they are projected onto any plane.

 Cylinder Diameter Height Aircraft representation λxy λz Collision 2λxy 2λz Conflict 2R 2H

### Table 1.

Modelling cylinders definition.

### Figure 4.

Modelling Cylinders Orientation

## 3. Derivation of a general expression for probability of collision (Pa)

In order to obtain a general expression of Pa an impact plane is defined as a generic projection plane containing the centre of reference aircraft ACi (assumed as static) and perpendicular to vji(relative velocity vector between the two aircraft i and j involved in the proximity event). Additionally, the collision area is defined as the projection, over the impact plane, of the collision cylinder (2λxy,2λz). If the conflict cylinder is settled in ACi, where its centroid is the one of the cylinder as well, the conflict area could also be defined as the projection of the conflict cylinder (2R, 2H). The CPAP (Closest Point Of Approach Projection) is a point with coordinates y1p and z1p obtained by projecting intruder aircraft. Figure 5 shows that a conflict will occur if ACj encounters the stationary conflict area, that is, if the CPAp coordinates (y1p, z1p) are inside the conflict area. In the same way, a collision will occur if ACj encounters the stationary collision area, that is, if the CPAp coordinates are inside the collision area.

#### Figure 5.

Impact Plane,Conflict and Collision cylinder.

#### Figure 6.

Impact Plane, Collision Area, Conflict Area and Projected CPA definition.

Considering the changes in the CPA coordinates due to radar and radar data segmentation errors, the probability of potential collision for an intruder aircraft that has violated the separation standards and whose projection consequently hits within the conflict area can be calculated as:

 Pa(y1p,z1p)=∫SPCFdP1P2≈SPCOL·∫SPCFf1(y'1p−y1p',z'1p−z1p)·f2(−y'1p−z'1p)dy'1dz'1 (3)

This equation provides and individual probability of collision based on:

• geometry and kinematics of the encounter (SPCOL),

• the predicted minimum lateral separation at the CPA (y1p), and

• the predicted minimum vertical separation at the CPA (z1p).

This takes into consideration the two probability density functions stating segmentation lateral and vertical errors and the projection lateral and vertical errors characterization.

As a result, the bi-dimensional probability density function of the CPAs can be derived from previous equation as:

 fa(y1p,z1p)=∫SPCFf1(y'1p−y1p',z'1p−z1p)f2(−y'1p−z'1p)dy'1dz'1 (4)

Where:

• fa is the bi-dimensional probability density function of the CPAs,

• y1p is the minimum predicted lateral separation at the CPA,

• z1p is the minimum predicted vertical separation at the CPA,

• SPCF is the conflict area

• f2(y1,z1) is the probability density function, representing the distribution of y1p and z1p coordinates errors due to the errors in the segmentation process, and

• f1(y’1p,z’1p) is the statistically determined bi-dimensional probability density function (pdf) of the CPA’p coordinates (y’1p,z’1p) for each projected segment associated to an individual encounter.

#### Figure 7.

Changes in the CPA coordinates due to projecting errors..

Both expressions estimate the probability of potential collision, having a potential separation violation (potential conflict), for each aircraft encounter, provided that uncertainties in the projection of segmented trajectories and in the segmentation process have been characterised by associated pdfs, f1 and f2, respectively.

## 4. Results and discussion

The previous mathematical formulation is supported by the previously mentioned ad-hoc software, which has been developed by the authors for Eurocontrol in the framework of the 3D-CRM programme. This software is intended to measure the collision risk in high density ATC en route airspace, based on an analysis of the stored aircraft tracks that have flown in it within a given time frame.

With the purpose of evaluating the mathematical expressions to estimate the probability of collisions, the previously mentioned software tool has been applied to a radar data sample from the Maastricht Upper Area Control Centre (MUAC). EUROCONTROL’s Maastricht Upper Area Control Centre (MUAC) is a regional air traffic control centre providing seamless air navigation services in the upper airspace (above 24,500ft) for a large (approximately 700,000 square kilometres) multinational airspace in Europe. An advanced and complex ATC automated system named MADAP (Maastricht Automated Data Processing and Display System) is the technical enabler responsible for managing, processing and presenting in real time information relating to the air traffic flows in the whole area. MADAP performs centralized multi-radar tracking using the information provided by a large number of radars and computes a high quality air traffic situation. In MUAC, a unique horizontal separation standard of 5 NM is used throughout the total area of responsibility. The vertical separation minimum of 1000 ft. is used.

### 4.1. Empirical estimation for Pa

The general expression of expected Pa is calculated numerically from the relative frequency of potential collisions among all potential conflicts using the following equation:

 E[P(pot.coll/pot.conf)] =E[Pa]≈Num.ofpot.collisionsNum.ofpot.conflicts=1935166=5.4*10−4 (5)

Figure 8 illustrates the obtained bi-dimensional histogram of the projected horizontal and vertical separations at the CPA for the whole data period analysed. As it is shown, the number of potential conflicts are higher when encounters are between aircraft established at the same flight level (0ft vertical separation) and, as well, between aircraft having 2.5 and 5NM of lateral separation. It could also be noticed that the number of encounters having 1000ft separation is higher than for any other vertical separation except the 0ft. This is easily understood when taking into account that within the en-route airspace most of the time aircraft are in level flight (namaly always 1000ft apart between contiguous flight levels). If safety barriers have not been applied the number of collisions to happen would have been 19. The area used to compute the number of potential collisions is shown circled by a red dotted circle.

### Figure 8.

histogram of projected horizontal and vertical separations at the CPA (31 days of radar data)

### 4.2. Pa estimation for each aircraft encounter

Once the empirical general or expected value for Pa has been obtained, Pa was estimated for each particular encounter by the next expression.

 Pa(yp,zp)=4λxyλz⋅vxvx2+vz2[1+π4⋅λxyλz⋅vzvx]⋅f2*(−yp,−zp)==2λxy⋅f2y(−yp)⋅2λz⋅λ2z(−zp)⋅vxvx2+vz2[1+π4⋅λxyλz⋅vzvx] (6)

This equation provides and individual probability of collision based on:

• kinematics of the encounter (ratio vz,to,vx),

• the predicted minimum lateral separation at the CPA (yp), and

• the predicted minimum vertical separation at the CPA (zp).

It also takes into consideration the segmentation of lateral and vertical errors (f2y and f2z).

A result for Pa estimation for leveled flight encounter is shown in the upper part of figure 8. In this case when CPAp coordinates (yp,zp) are very close to the reference aircraft (ACi), Pa estimated value reaches 3*10-2. This value has a magnitude of two orders higher than the empirical expected result (5.4 10-4), but strongly decreases when predicted CPAp lays apart from ACi, resulting in values much lower than the empirical one. In the lower part of this figure, the graphs show when one or both aircraft are climbing/descending but having vz/vr ratio close to zero, it could be seen that regardeless the decrease of the maximum value of Pa (7*10-3), It is still greater than the empirical expected result for Pa. Furthermore, the probability of collision for CPAp for which yp coordinates close to zero but zp coordinates separated from the ACi remains significant. Pa estimation for encounters having two different aircraft climbing/descending (vz / vx) ratios is shown in figure 9.

### Figure 9.

Pa estimation for different CPAp. Aircraft established at a defined flight level or vz equals to zero (upper) and aircraft with vz close to zero (lower).

Despite the fact that the shape of both functions for Pa are similar to the one obtained in the lower part of figure 10 (aircraft climbing/descending and vz/vx close to zero), the maximum values for Pa are different in both cases (9*10-3 for vz/vx=0.1, and 2*10-2 for vz/vx=20), showing that Pa maximum values for CPAp close to reference aircraft (ACi) has a decreasing trend when vz/vr ratio increases. The following table summarises the results obtained from empirical and estimated Pa for the worst case, that is to say Pa for predicted CPAp=(0,0).

The results clearly shows that it is unrealistic to assign the same probability for potential collisions to all potential conflicts, independently of the predicted coordinates for CPA, no matter how these coordinates have been derived.

 Empirical result for expected Pa, E[Pa] 5.4·10-4 Estimated Pa for CPAp=(0,0) and level flight 3·10-2 Estimated Pa for CPAp=(0,0) and vz /vx ≈0 7·10-3 Estimated Pa for CPAp=(0,0) and vz /vx =0.1 9·10-3 Estimated Pa for CPAp=(0,0) and vz /vx =20 2·10-2

### Table 2.

Worst case Pa estimation

### Figure 10.

Pa estimation for different CPAp. Aircraft climbing/descending and different vz / vx ratios. vz / vx=0.1(upper), vz / vx=20 (lower)

### 4.3. Expected Pa estimation for a given scenario and traffic sample

When a collision risk analysis is applied to a representative aircraft population, using segmentation of their stored radar tracks, a 2D histogram of projected horizontal and vertical separations at the CPA can be obtained, as it is shown in figure 8. This histogram provides a first approach for expected Pa using equation (4), which is the way we used to obtain E[Pa]= 5.4*10-4, (this value can taken as reference value for Pa). If the histogram exhibits a close to uniform distribution, it can be understood that any “generic” potential conflict would became a potential collision with the same probability. It is also possible to propose a different approach to establish the expected value for Pa in a given scenario and for a given aircraft population, discussed below.

 E[Pa]=1N∑jiPa(yji,zji,rji)=λxyλzN∑ji[1+π4⋅λxyλz⋅rji]f2y(yji)⋅f2z(zji)==λxyλzN∑ji[1+π4⋅λxyλz⋅rji]f2y(yji)⋅f2zji(zji) (7)

Where Pa(yji,zji,rji) is the individual probability of each potential collision where:

• rji=vz/vx the between vertical and horizontal relative speeds,

• f2zji the probability density function applied to each aircraft encounter (between each pair of aircraft, i and j).

When this equation is applied to previous MUAC data sample, expected value for Pa results 8.2*10-4, which is slightly higher than the empirical results.

## 4. Conclusions

This chapter analyse in detail the inherent collision risk involved for each aircraft proximity event by assessing the conditional probability Pa of a potential collision between aircraft that are exposed to risk, that is to say, they are potentially going to violate the separation standards defined for a specific airspace if no corrective action is taken. The proposed approach allows the determination of the severity of each aircraft encounter as the probability of potential collision for each individual aircraft encounter in high density ATC en route airspace, based on an analysis of the stored aircraft tracks that have flown within a given time frame. The authors propose a mathematical formulation to characterise the severity of each aircraft proximity event using the convolution of the bi-dimensional probability density function of the predicted Closest Point of Approach between the aircraft involved and the distribution of lateral and vertical error in the projected position of the aircraft.The presented work aims to provide an individual probability of collision based on the geometry and kinematics of the encounter and the minimum lateral separation and the minimum vertical separation at the predicted Closest Point of Approach or CPA. The formula takes into consideration uncertainties introduced by the radar data error and the segmentation error. The results of this chapter shows that there is not the same severity for all the proximity events on which aircraft pass closer than the prescribed horizontal and vertical separation minima, and also that the expected severity for given a scenario and traffic sample can also vary depending on the kinematic characteristics of aircraft involved within this scenario. It is also considered that collision risk for high density of air traffic can be analysed from the estimation of three different factors:

• Relative frequency of exposition to risk (FeR). The value of this factor can be easily obtained from any radar data sample and strongly depends on the minimum applied horizontal and vertical separations standard and increases with air traffic density,

• Expected severity E(Pa). This value can be directly derived from individual probabilities of potential collision (Pa). Furthermore, having individual severities, it also permits additional assessment on safety (hot spots identification, etc.).

• Expected probability of failure of safety barriers (ATC, TCAS, etc.)

As the two first factors can be derived from the stored tracks of the traffic sample, using the software tool developed by the authors [38], further work is now devoted to develop the probability of failure of the ATM safety barriers. Once the probability of failure were stated and validated, it will be possible to estimate the collision risk for individual encounters, scenarios and air traffic samples. Results obtained for MUAC, with data sample used in previous discussion, exhibits a rounded value for frequency of exposition to risk of FeR=0.3. Probability of potential collision among encounters exposed to risk, Pa or its expected value E(Pa) for the same sample, oscillates between 8.2*10-4(expected) and 2*10-2(worst case). Previous results demand a probability of “safety barrier failure” lower than 0.4*10-5 and 1.7*10-7 respectively, to reach the ATM en route target level of safety of TLS=10-9. This last value is normally the one used as TLS. For instance, in reference (Eurocontrol, 2006) mid-air collision given as accident frequency (per flight) is 5.4*10-09, specifying that, among them, the frequency of fatal accident, directly caused by ATC (per flight), is 3.5*10-09.

## References

1 - R. E. Machol, 1995Thirty Years of Modelling Midair Collisions", Interfaces 25: 5 September-October 1995 (151-172)
2 - B. L. Marks, 1963Air traffic control separation standards and collision risk. Royal. Aircraft Establishment Technical Note 91February, 1963
3 - P. G. Reich, 1964A theory of safe separation standards for Air Traffic Control, Technical Report 64041, Royal Aircraft Establishment, UK
4 - P. G. Reich, 1966Analysis of Long-range Air Traffic Systems: Separation Standards. Journal of the Institute of Navigation, (19), 88, 169 and 331 (in three parts).
5 - Peter Brooker.Longitudinal Collision Risk for ATC Track Systems: A Hazardous Event Model. Journal of Navigation, 2006591pag. 55-70.
6 - ICAO1988Review of the General Concept of Separation Panel, 6th meeting, Doc 9536, 1ICAO, Montreal, December 1988.
7 - H. D. Sherali, C. Smith, A. A. Dr, S. Trani, Q. Sale, Chuanwen, Analysis of Aircraft Separations and Collision Risk Modeling. NEXTOR- National Center of Excellence for Aviation Operations Research 1998
8 - L. Burt, October, 200020003D Mathematical Model for ECAC Upper Airspace, Final Report
9 - Peter Brooker (Cranfield University).Radar Inaccuracies and Mid-Air Collision Risk: Part 2 En Route Radar Separation Minima The Journal of Navigation (2004).
10 - Brenda. D. Carpenter, M. I. T. Cambridge, MA , James. K. Kuchar, M. I. T. Cambridge, MA , Probability-based collision alerting logic for closely-spaced parallel approach. AIAA-1997222Aerospace Sciences Meeting and Exhibit, 35th, Reno, NV, Jan. 6-9. 1997.
11 - James. K. Kuchar, M. I. T. Cambridge, MA , Lee. F. Winder, M. I. T. Cambridge, Generalized philosophy of alerting with applications to parallel approach collision prevention. MA AIAA-20014052AIAA Guidance, Navigation, and Control Conference and Exhibit, Montreal, Canada, Aug. 6-9, 2001.
12 - F. Lee, . Winder, K. James, Kuchar, Evaluation of Collision Avoidance Maneuvers for Parallel Approach. Journal of Guidance, Control, and Dynamics 0731-5090 226doi:1999
13 - J. Powell, Stanford. David, Univ, CA , Sharon. Houck, Univ. Stanford, CA of. Assessment, possibility. the, a. of, collision. midair, an. during, closely. ultra, parallel. spaced, approach, A. I. A. , AIAA-20014205AIAA Guidance, Navigation, and Control Conference and Exhibit, Montreal, Canada, Aug. 6-9, 2001.
14 - W. Sharon, J. Houck, Powell. David, of. Probability, Collision. Midair, Ultra. During, Spaced. Closely, Approaches. Parallel, Journal of Guidance, Control, and Dynamics, 0731-5090 265do: 10.2514/2.5124. 2003
15 - Roger. Shepherd, Corp. Rannoch, V. A. Alexandria, Rick. A. Cassell, aircraft. reduced, risk. separation, model. V. A. A. I. A. assessment, AIAA-19973735AIAA Guidance, Navigation, and Control Conference, New Orleans, LA, Aug. 11-13, Collection of Technical Papers. Pt. 3 (A97-37001 10-63). 1997.
16 - H. A. P. Blom, G. J. Bakker, P. J. G. Blanker, J. Daams, M. H. C. Everdij, M. B. Klompstra, ‘‘, Accident Risk Assessment for Advanced ATM,’’ In: Air Transportation Systems Engineering, G.L. Donohue and A.G. Zellweger (Eds.), AIAA, 2001463480
17 - H. Blom, G. J. Bakker, B. , Klein Obbink and MB Klompstra. Free Flight safety risk modeling and simulation. Proceedings of 2nd International Conference on Research in Air Transportation ICRAT 2006at Beograd, Serbia, June 2428
18 - H. Blom, B. Klein, B. Obbink, National. Bakker, Laboratory. N. L. R. Aerospace, Safety Risk Simulation of an Airborne Self Separation Concept of Operation. AIAA-20077729th AIAA ATIO Conf, 2nd CEIAT Int’l Conf on Innov and Integr in Aero Sciences,17th LTA Systems Tech Conf; followed by 2nd TEOS Forum, Belfast, Northern Ireland, Sep. 18-20, 2007.
19 - Peter Brooker.Air Traffic Management accident risk. Part 1: The limits of realistic modeling. Safety Science 444194502006
20 - R. A. Paielli, H. Erzberger, Conflict. ‘‘, estimation. probability, free. for, A. I. A. A. J. flight’’, of Guidance, Control and Dynamics, 205885961997
21 - R. A. Paielli, H. Erzberger, Conflict. ‘‘, Estimation. Probability, to. Generalised, Flight’’. Non-Level, Traffic. Air, Quarterly. Control, Vol, 1952221999
22 - M. Prandini, J. Hu, J. Lygeros, S. Sastry, A. probabilistic, to. approach, conflict. aircraft, I. E. E. E. detection, Tr, onflict detection, IEEE Tr. on Intelligent Transportation Systems, 141992202000
23 - M. Prandini, J. Lygeros, A. Nilim, S. Sastry, “. A. Probabilistic, for. Framework, Conflict. Aircraft, A. I. A. Detection”, A-99-414, Proc. A. I. A. in, AIAA-99-4144, in Proc. AIAA Guidance, Navigation, and Control Conf., Portland, OR, August 9-11, 104710571999
24 - ICAO, Annex 11 - Air Traffic Services, 12th edition, incorporating amendments 1-38, July1998Green pages, attachment B, paragraph 3.2.1.
25 - J. Kuchar, L. Yang, “. A. Review, Conflict. of, Detection, Modelling. Resolution, I. E. E. E. Methods,”, on. Transactions, Transportation. Intelligent, Vol. Systems, 4179189December 2000
26 - Kuchar, L. Yang, “. Survey, Conflict. of, Detection, Modelling. Resolution, A. I. A. Methods”, A-97373 , Proc. A. I. A. in, AIAA-97-3732, in Proc. AIAA Guidance, Navigation, and Control Conf., New Orleans, LA, August 11-13, 1997
27 - A. P. Henk, Bart. Blom, Obbink. G. J. . Klein, Bakker. Bert, Safety risk simulation of an airborne self separation concept of operation, 7th AIAA Aviation Technology, Integration and Operations Conference (ATIO)<BR>2nd C 1820Belfast, Northern Ireland AIAA 2007September 2007.
28 - Yang. Lee, Hyun. Ji, James. Yang†, Eric. Kuchar‡, Massachusetts. Feron§, of. Institute, Cambridge. Technology, MA , A. A Real-Time Monte Carlo Implementation for Computing Probability of Conflict. AIAA Guidance, Navigation, and Control Conference and Exhibit 1619Providence, Rhode Island. August 2004
29 - L. Burt, October, 200020003D Mathematical Model for ECAC Upper Airspace, Final Report.
30 - L. M. B. C. . Campos, J. M. G. Marques, On the probability of collision between climbing and descending aircraft, ; Journal of aircraft 0021-8669CODEN JAIRAM / 4425505572007
31 - L. M. B. C. Campos, Probability of collision of aircraft with dissimilar position errors, ; Journal of aircraft 0021-8669CODEN JAIRAM, 3845935992001
32 - A. Leonard, The. M. I. T. R. E. Wojcik, Lean. V. A. Corporation Mc, 221, U. S. , Probabilistic Aircraft Conflict Analysis for a Vision of the Future Air Traffic Management System. AIAA 5th Aviation, Technology, Integration, and Operations Conference (ATIO) 2628Arlington, Virginia. September 2005
33 - Leonard A.Wojcik∗ The MITRE Corporation, McLean.Probabilistic Aircraft Conflict Analysis for a Future Air Traffic Management System. VA 22102 DOI:Journal of Aerospace Computing, Information, and Communication 6June 2009
34 - L. Rachelle, Ennis, J. Yiyuan, University. Zhao, Minnesota. of, M. N. Mpls, Defining Appropriate Inter-Aircraft Separation Requirements. AIAA 4th Aviation Technology, Integration and Operations (ATIO) Forum 2022Chicago, Illinois. September 2004
35 - Dingy. Jerry, Tomlinz. A. Claire, Programming. Dynamic, for. Approach, Conflict. Aircraft, A. I. A. Detection, AIAA Guidance, Navigation, and Control Conference, 1013Chicago, Illinois. August 2009
36 - D. E. Stepner, Modelling of Aircraft Position Errors with Independent Surveillance. 119AIAA Journal 1273. September 1973
37 - Eduardo José García González, INECO, Madrid, Spain, Francisco Javier Sáez Nieto, Polytechnic University of Madrid, Madrid, Spain, Maria Isabel Izquierdo, EUROCONTROL, Brussels, Bélgica.Identification and analysis of proximate events in high density en route airspaces Paper N° 63]. 7th USA- EUROPE ATM R&D Seminar July 0205Barcelona. 2007
38 - F. Saez, R. Arnaldo, E. Garcia, G. Mc Auley, M. Izquierdo, Development of a three dimensional collision risk model tool to asses safety in high density en-route airspaces. DOI:10.1243/09544100JAERO704
39 - Jaroslav Krystul.Modelling of stochastic hybrid systems with applications to accident risk assessment.September 2006
40 - L. Rachelle, Ennis† and Yiyuan J. Zhao‡ University of Minnesota, Minneapolis, Minnesota. Characterization of Aircraft Protected Zones. AIAA’s 3rd Annual Aviation Technology, Integration, and Operations (ATIO) Tech 1719Denver, Colorado. November 2003
41 - J. F. Bellantoni, The Calculation of Aircraft Collision Probabilities DOT-TSC-FAA-7127October’1971
42 - Reynolds T.G., Hansman R.J., Analysis of Aircraft Separation Minima Using a Surveillance State Vector Approach, MIT, 2001
43 - Ennis, R. L.; Zhao,Y.J., Defining Appropriate Inter-Aircraft Separation Requirements, AIAA’s 4th Annual Aviation Technology Integrations and Operations (ATIO) Forum # AIAA 2004-6203, September 20-22, 2004