arxiv: v2 [cs.sy] 6 Feb 2011

Transcription

1 Arcraft Proxmty Maps Based on Data-Drven Flow Modelng Erwan Salaün 1, Georga Insttute of Technology, Atlanta, GA, Maxme Garel 2, Massachusetts Insttute of Technology, Cambrdge, MA, arxv: v2 [cs.sy] 6 Feb 2011 Adan E. Vela 3, and Erc Feron 4 Georga Insttute of Technology, Atlanta, GA., Wth the forecast ncrease n ar traffc demand over the next decades, t s mperatve to develop tools to provde traffc flow managers wth the nformaton requred to support decson makng. In partcular, decson-support tools for traffc flow management should ad n lmtng controller workload and complexty, whle supportng ncreases n ar traffc throughput. Indeed, the growth of ar transportaton s condtoned on the ablty to mantan acceptable safety. Whle many decson-support tools exst for short-term traffc plannng, few have addressed the strategc needs for medum- and long-term plannng for tme horzons greater than 30 mnutes. Ths paper seeks to address ths gap through the ntroducton of 3D arcraft proxmty maps that evaluate the future probablty of presence of at least one or two arcraft at any gven pont of the arspace. Three types of proxmty maps are presented: () presence maps that ndcate the local densty of traffc, () conflct maps that determne locatons of potental conflcts and ther correspondng probabltes, and () outler proxmty maps that show the probablty of conflct due to arcraft not belongng to 1 Postdoctoral Fellow, School of Aerospace Engneerng, 270 Ferst Drve, erwan.salaun@gatech.edu. 2 Postdoctoral Assocate, Laboratory for Informaton and Decson Systems, 77 Massachusetts Avenue, mgarel@mt.edu. The work was realzed when Maxme was a PhD canddate n the School of Aerospace Engneerng of Georga Insttute of Technology. 3 Ph.D. Canddate, School of Mechancal Engneerng, 801 Ferst Drve, aevela@gatech.edu. 4 Professor of Aeronautcs and Astronautcs, School of Aerospace Engneerng, 270 Ferst Drve, feron@gatech.edu 1

2 domnant traffc patterns. These maps provde traffc flow managers wth nformaton relatng to the complexty and dffculty of managng an arspace. The ntended purpose of the maps s to antcpate how arcraft flows nteract, and how outlers mpact the domnant traffc flow for a gven tme perod. Ths formulaton s able to predct whch crtcal regons may be subject to conflcts between arcraft, thereby requrng careful montorng and addtonal effort to manage the arspace. These probabltes are computed usng a generatve arcraft flow model. Tme-varyng flow characterstcs, such as geometrcal confguraton, speed, and probablty densty functon of arcraft spatal dstrbuton wthn the flow, are determned from archved Enhanced Traffc Management System data, usng a talored clusterng algorthm. Arcraft not belongng to flows are dentfed as outlers. I. Introducton ACCORDING to the predcted ar traffc growth, current arspace capacty lmts,.e. the maxmum number of arcraft allowed n an arspace, wll be ncreasngly exceeded throughout the Natonal Arspace System (NAS) n the comng decades [1]. Currently, the capacty of an arspace s measured accordng to the peak number of arcraft allowed n the consdered sector wthn any one mnute perod [2]. When the predetermned thresholds are exceeded, congeston alerts are trggered. In practce, capacty thresholds are set by the ablty of controllers to handle traffc and to mantan manageable workloads for worst-case plannng over a wde range of traffc patterns [3]. Hence, current arspace capacty thresholds are ndependent of any changes n the traffc confguratons that may occur due to weather or other extraneous traffc perturbatons. Such a ll-defned method for establshng and measurng arspace capacty often leads to under-utlzaton of the arspace under nomnal condtons. Or conversely, when convectve weather s present, over-utlzaton s not uncommon. To accommodate hgh levels of ar traffc throughput and to effectvely manage arspace utlzaton whle mantanng safety, ntrnsc descrptors of traffc and the arspace are requred to support arspace management n the redesgn or adjustment of flow confguratons. These measures 2

3 are often referred to as complexty measures, snce they am at provdng an mage of the dffculty or complexty of managng the consdered arspace and traffc confguraton. There has been a sgnfcant volume of research related to the estmaton of ar traffc control complexty, and ts applcaton to determne arspace capacty. One method for defnng and establshng ar traffc control complexty s dynamc densty [4, 5]. Ths metrc s a weghted sum of the local densty of arcraft; the number of headng, alttude, and speed changes; and a varety of other factors. The weghtng factors are determned and valdated through human-n-the-loop experments. Dynamc densty makes use of nstantaneous arcraft confguratons, and therefore requres accurate predcton of traffc to hold any sgnfcant forecastng value. The mpact of uncertantes on arcraft trajectores on predctng dynamc densty has been studed n [5, 6], and was demonstrated to cause short-comngs n forecastng complexty. Another ntrnsc measure of ar traffc control complexty, proposed n [7], analyzes the stablty of a nonlnear dynamcal system correspondng to the consdered traffc pattern. Ths analyss method s attractve for ts ablty to seamlessly account for arcraft poston uncertantes. An alternatve approach toward managng complexty s presented n [8]. The authors propose a dynamc sectorzaton of the arspace based on ar traffc controller workload, whch leads to the generaton of another knd of complexty map. Recently, sgnfcant effort has been placed nto understandng the probablty of conflct and the assocated resoluton-commands requred to ensure separaton. For example, several tools are presented n [9 11] to determne the probablty of conflct for a gven geometrc traffc confguraton, wth and wthout arcraft poston uncertanty. In [12], an nput-output approach s proposed that generates complexty maps based on the conflct-resoluton commands requred to accommodate dsturbances such as the entrance of an exogenous arcraft nto the arspace. In [13], the authors use arcraft poston and ntent nformaton to compute the probablty of presence of the arcraft at any pont n the arspace wthn a gven tme frame. In [14, 15], the authors propose a complexty metrc that s a functon of the local average velocty vector, the densty of arcraft and the proxmty of an arcraft to the consdered pont of the arspace, n order to evaluate the solutons computed by ther traffc flow management (TFM) algorthm. Nevertheless, a key pece of nformaton formng the foundaton of the prevously cted works s 3

4 knowledge of the ntal arcraft postons and trajectores ( known wth or wthout uncertanty). Such nformaton, whle useful for short-term ar traffc management, e.g. plannng a safe trajectory for an ntruder arcraft through the arspace (e.g. [13]), s not well-suted for md-term or long-term ar traffc management, that s, for tme horzons greater than 30 mnutes. Indeed, from a Traffc Flow Management perspectve, strategc plannng requres knowledge of ar traffc patterns and flows characterstcs (e.g. geometrcal confguraton, flght plans, dstances between two consecutve arcraft), and not the current poston and ntent nformaton of ndvdual arcraft. As explaned n [16], consderng arcraft flows enables the computaton of more predctable and perturbatonresstant estmates than consderng ndvdual arcraft. Representng ar traffc as flows provdes a more meanngful and hgher-level vew of ar traffc to traffc flow managers. In fact, such a representaton s also consstent wth an ar traffc controllers understandng of a sector. As explaned n [17, 18], one such way ar traffc controllers abstract sectors s accordng to domnant traffc flow patterns and traffc flow nteractons. Few works have developed TFM tools wth ths vewpont n mnd. The authors of [16] provde a lst of flow features to descrbe flow pattern, such as the number of flows, the major flow and ts sze, the number of crossng flows, etc. In [19], these crtera allow the authors to estmate and predct sector capacty as a functon of the traffc flow pattern, and to study the mpact of severe weather. Whle the lst of flow features establshed n [16] provdes a noton of the traffc complexty (gven as a scalar number), more nformaton may be needed by traffc flow managers to dentfy regons of hgh complexty. In [20], a smplfed approach to generate complexty maps, takng nto account the nfluence of automated an conflct-resoluton algorthm and the flows characterstcs, s proposed. However, only two flows are consdered, and the extenson of the results of [20] to a realstc arspace (.e. multple ntersectons wth multple flows) s not obvous and may be computatonally expensve. Buldng upon the prelmnary work n [21], ths paper proposes a new knd of ar traffc complexty maps amed at supportng traffc flow management. These maps, referred to as (arcraft) proxmty maps, provde an estmaton of probablty of proxmty of one or more arcraft to a regon of arspace. For smplcty, the expresson probablty of proxmty s used throughout the paper nstead of the correct (and more rgorous) expresson probablty an arcraft s wthn a gven 4

5 neghborhood of a pont P n the arspace. Three knds of proxmty maps are developed, each wth a dfferent probablty metrc. They are as follows: Presence maps depct the probablty of presence of at least one arcraft from a flow n a gven volume of the arspace. These maps enable the dentfcaton of the man flows of arcraft and the regons of hgh arcraft densty. Conflct maps provde the probablty of presence of at least two arcraft from two dfferent flows n a gven volume of the arspace. These maps hghlght regons of arspace wth hgh probablty of conflct or close proxmty between arcraft, thereby ndcatng regons of the arspace that requre addtonal attenton from the ar traffc controllers. Outler proxmty maps gve the probablty of presence of an arcraft from the domnant traffc and an outler arcraft (defned n Secton II E) n a consdered volume of the arspace. These maps ad traffc flow managers n evaluatng the mpact of arcraft outlers on the domnant ar traffc flows. Regons of the maps wth hgh values for the probablty of presence, conflct, or outler proxmty are nterpreted to be crtcal regons that requre montorng. The proposed proxmty maps take nto account flow characterstcs n a manner consstent wth TFM (.e. aggregate flow descrptons such as average mles-n-tral), thereby makng them an ntutve and effectve strategc decson-support tool. To enable the creaton of the presence and proxmty maps, a consstence model of arspace s requred. However, because a descrbng arspace s qute complex, we ntroduce a dmensonally-reduced model. The dmensonally-reduced model serves as a generatve model that s consstent wth the tools we are proposng. That s, the ntal modelng provdes a framework n whch a parametrc model can be establshed. For ths purpose, the clusterng algorthm ntroduced n [22] s used to process archved Enhanced Traffc Management System (ETMS) data to generate a realstc ar traffc flow model,.e. the generatve model. The resultng model accounts for the arcraft spatal dstrbutons and arcraft arrval dstrbutons for each flow. Ultmately, the methodology for generatng the traffc model provdes a data-drven framework on whch proxmty maps are based. The process, startng from archved data to the 5

6 Fg. 1 Work-process flow from hstorcal data to the generaton of proxmty maps. generaton of proxmty maps, s llustrated n Fg. 1, wth an example of clustered trajectores and presence maps. The usage of proxmty maps extends beyond dsplayng crtcal nformaton. They may be used for decson makng n traffc flow management. Examples nclude dynamcally adjustng flow characterstcs; redefnng or creatng new flght routes; or vsualzng the nfluence of rogue arcraft (smlar to [12]) or severe weather (as n [15]) located n a specfc arspace regon. The tools presented n ths artcle may also be utlzed as decson support tools n conjuncton wth the Flow-Based Route Plannng method used n [14, 15], or other such tools. The remander of the paper s organzed as follows. Secton II presents a methodology for developng a traffc model through the extracton of ar traffc flows and outlers from archved flght data. The resultng model serves as a parametrc generatve model for computaton of the proxmty maps. In the next secton, an emprcal analyss of the clusterng s completed. The analyss s completed to support assumptons and smplfyng calculatons used n the proxmty maps. Secton IV ntroduces the noton of the resdual dstance of an arcraft and descrbes how presence, conflct and outler proxmty maps are generated from the flow and outler models. In Secton V, example proxmty maps are generated for a representatve regon of arspace. Fnally, n Secton VI, our conclusons are presented. 6

7 II. Arspace Modelng: Transformaton of flght data to flow and outler models Ths secton ntroduces a methodology to buld mathematcal models of arcraft operatons wthn an arspace derved from ETMS data. The resultng model serves as a generatve parametrc model to create the proposed proxmty maps. To precede, useful defntons are frst ntroduced. Next, the ETMS data set s descrbed, followed by a bref dscusson of the clusterng algorthm used to generate the traffc models. The clusterng algorthm s then appled to an typcal arspace, and analyss of the resultng traffc flow and outler model s provded. A. Defntons Consder a 3D arspace A R 3, whose arbtrarly shaped boundares are fxed and ndependent of the alttude. The coordnates of a pont P = (x, y, z) of the arspace are expressed n an East- North-Up frame R 0, where the x, y axes are expressed n nautcal mles (NM) and the alttude z s expressed n feet (ft); a smple lattude-longtude to flat earth projecton s used, wth the orgn of the arspace s located at the center of the arspace. For the scope and scale of the arspace (approxmately 400 NM n the largest dmenson), the error ntroduced by ths approxmaton s nsgnfcant n relaton to the topologcal descrptors need to understand generate the proxmty maps. As customary, Flght Levels (FL) are used for the z-coordnate. Flght levels are dscretzed at 100 ft ncrement, such that FL 320 corresponds to an alttude of 32,000 ft. The followng terms, used throughout the artcle, some of whch are depcted n Fg. 2, are defned as follows: A trajectory s a sequence of l ponts defned by a set of 3D coordnates. All trajectores are resampled to have same number of ponts, l = 8. Trajectores are spatally sampled, nstead of temporally sampled, n order to allow for dentfcaton of spatally smlar trajectores. A cluster, C, s a labelng of m trajectores that are dentfed as smlar by the clusterng algorthm. A track T s the centrod of the cluster C,.e the trajectory center of mass of all the trajectores n the cluster. The track T can be nterpreted as the representatve trajectory for cluster C. 7

8 The wndow W k of cluster C s located at the k th pont along the assocated track T. The wndow W k s vertcally orentated wth ts face algned accordng to the headng of T at the k th pont towards the (k + 1) th pont. The dmensons of each wndow are determned by the extrema of the k th ponts correspondng to each trajectory wthn the cluster labelng. At each wndow W k, the local coordnate frame ( a k, b k, c k ) located at pont kth pont on track T k. The vertcal axs, c k, s algned wth the global z coordnate axs, whle a k track T, such that t ponts towards W k+1. algn wth the A box B k, k = 1... l 1, s the polyhedron defned by the convex hull of two consecutve wndows, W k and W k+1. A tube s the successon of all the boxes of a gven cluster. For every wndow W k, the sampled lateral dstrbuton of arcraft, f a, s computed from the k orthogonal projectons of the k th pont of the trajectores n C onto the a k axs of wndow W k. For every wndow W k, the sampled vertcal dstrbuton of arcraft, f c k s computed from the alttudes of the k th pont of the trajectores n C along the c k axs. The arcraft arrval rate parameter λ(t), s the expected number of arcraft enterng the arspace regon durng a gven tme perod. The arrval rate for each cluster C, λ (t), s defned smlarly. For the remander of the paper, the consdered tme perod s 15 mnutes long. A flow F s a 3-dmensonal tube that contans the trajectores of cluster C. The tube s defned by a successon of l wndows W k, k = 1... l. Assocated wth a flow are ts arcraft speed dstrbuton, arcraft arrval rate, and ts arcraft lateral and vertcal dstrbutons. B. ETMS Data Set The data used to construct the arspace traffc model s taken from the Enhanced Traffc Management System (ETMS) data set. Cleveland center s selected for the study because of ts sgnfcance 8

9 Fg. 2 A 3D representaton of an ascendng flow wth ts correspondng tube. to the NAS. Indeed, a large share of the traffc to and from the North-East of the Unted States fles through the center. Furthermore, as noted n [23], Cleveland sector s one of the most congested and delay-prone Centers n the contnental Unted States, and centrally to located to a number of hgh traffc regons. The ETMS data ncludes arcraft radar tracks (longtude, lattude, alttude) taken at approxmately 1 mnute ntervals. Durng the 123 days (May to August 2005) covered by the data, all 526,840 arcraft trajectores wth at least one pont over FL 250 are consdered; ths subset of the data covers the majorty of en route arcraft. The dates 08/04/2005, 08/24/2005, and 08/26/2005 are not ncluded n the clusterng analyss, as the correspondng data was not avalable. ETMS data s notorous for ts nconsstences n alttudes, therefore a smple flterng algorthm on reported alttudes was used. After ths ntal flterng process, arcraft wth nconsstences n reported postons and alttudes are excluded. Flght data s deemed to be nconsstent f t exhbts physcally unrealzable trajectores, e.g. rapd changes n alttude or poston. After flterng out nconsstent trajectores, a clean data-set of 338,060 trajectores remans. Statstcs not requrng the alttude or the poston nformaton, such as the arcraft arrval rate nto the arspace are computed on the entre en route data set, whle the clusterng of trajectores s run on the clean subset. Fgure 3 9

10 presents a randomly sampled densty plot of 110,000 clean trajectores over Cleveland center (only 110,000 flghts are plotted because of memory lmtatons). Major ar traffc flows are clearly dentfable as the darkest regons of the arspace; many correspondng to defned jetways. The appearance of the gray shades ndcates that some arcraft fly outsde of the major flows. The clusterng algorthm presented n the next subsecton s used to dentfy the major flows, thereby establshng to domnant traffc pattern. Fg. 3 Densty plot of a subset of the clean trajectores. C. Trajectory Clusterng The 338,060 trajectores are clustered usng a slghtly modfed verson of the algorthm presented n [22]. The algorthm s defned accordng to the followng steps: 1. Clean and format the trajectores. Inconsstent trajectores are excluded. The remanng trajectores are spatally resampled accordng to a lnear nterpolaton wth a fxed number of equally space ponts. The lnear approxmaton s suffcent as most arcraft fly along straght trajectores n ths center. 2. Augment the dmensonalty of the trajectory data by ncludng addtonal features such as headng, polar coordnates, etc. 10

11 3. Apply herarchcal clusterng. Organze and separate the trajectores by alttude and atttude to create separate data sets. Level trajectores are splt by flght level; clmbng trajectores are sorted by destnaton (or arspace regon ext) flght level; and descendng trajectores are splt by ntal flght level. 4. Normalze each feature and concatenate the data nto a sngle row vector for each flght. Each herarchcal groupng - defned by alttude and atttude - has a correspondng matrx of data, ncludng nterpolated trajectory ponts and augmented data. Each row of the matrces corresponds to a flght, and each column corresponds to a feature. 5. Apply prncpal component analyss [24] on the prevously descrbed matrces and keep the projectons onto the frst 5 prncpal components, thereby reducng the dmensonalty of the data. 6. Cluster the values of the projectons usng a densty-based clusterng algorthm (DB- SCAN [25]). 7. Obtan clusters of trajectores and outlers for each alttude and atttude category. Outlers are trajectores that are not assgned to a cluster. Here, a key extenson to the work presented n [22] s to splt the dataset nto subsets, by alttude and atttude (clmb, descent, level flght), to greater fdelty n modelng. Ths ntal categorzaton results n 136 subsets. As noted, the trajectores for each category are clustered usng the presented algorthm. The number of clusters and ntra-cluster smlarty can be adjusted by tunng nput parameters to DBSCAN (the mnmum number of ponts requred to start a cluster, and the mnmum dstance between those ponts). Wth such a large sample set of trajectores, a two-pass approach s taken towards the clusterng. The frst round dentfes both large and small clusters. A second applcaton of DBSCAN s appled to large clusters wth more dscrmnatng parameters. Followng categorzaton and clusterng, approxmatvely 80% of the trajectores are grouped nto 685 clusters, or flows. Fgure 4 presents a 2D and 3D vew of the centrods of the clusters. Blue lnes represent westbound level traffc, yellow lnes represent eastbound level traffc, green lnes 11

12 200 Centrods for all clusters 500 Centrods for all clusters [NM] Alttude [x100 ft] [NM] 0 [NM] [NM] (a) Top vew (b) Oblque vew Fg. 4 Centrods for all traffc flow clusters. represent descendng traffc, and red lnes represent ascendng traffc. Eastbound traffc (yellow lnes) appear to be solely passng through the center of the center. The reason s that Eastbound traffc s wdespread over the entre center and no clear major flow s dentfable. Therefore, for each flght level, most of the Eastbound trajectores belong to the same cluster whch s very wde. In Fgure 4(b), major arports of Cleveland center Cleveland Hopkns arport (CLE), Detrot Metropoltan Wayne County Arport (DTW) and Pttsburgh Internatonal Arport (PIT) are clearly dentfable by the clusters correspondng to descendng ar traffc. Gven the greater varablty n takeoff trajectores, n comparson to landng trajectores, few ascendng clusters from arports are dentfed. Prmary reasons for varablty nclude the dverse range n arcraft performance and plot decsons characterstcs. D. Clusterng Analyss Applyng clusterng to the flght data dentfes numerous flows over a broad range of trajectory types. To better understand the structure of ar traffc n relaton to the clusterng algorthm, results are analyzed over a broad range of parameters to determne when flow models of the arspace may lose valdty. Of partcular nterest s the proporton of outler arcraft, that s, arcraft trajectores 12

13 Number of trajectores Clustered Outlers Leveled at FL Ascendng to FL Descendng from FL Fg. 5 Dstrbuton of clustered nad outler flghts by flght level. that are not assgned a cluster labelng. Fgure 5 presents the dstrbuton of clustered and outler trajectores by flght level. Leveled at FL ndcates that the alttude s constant for the correspondng arcraft trajectory; Ascendng to FL ndcates the destnaton flght level for the arcraft trajectory wthn the center; and Descendng from FL ndcates the ntal flght level of the arcraft before startng ts descent. For the consdered flght levels over 80% of trajectores belong to a cluster. Ths result supports the hypothess that flows exst, and can be meanngful towards descrbng ar traffc. Further analyss s completed to determne the fracton of trajectores clustered accordng to tme of day, and date. Fgure 6 llustrates the dstrbuton of clustered and outler trajectores by tme of the day. Accordng to the fgure, the fracton of outlers s relatvely constant. Meanwhle, the number of trajectores clustered vares throughout the day, wth the maxmum number of clusters beng establshed durng peak hours,.e. between 11 a.m. and 2 a.m. GMT (5 a.m. and 8 p.m. local tme). Fgure 7 presents the number of trajectores dentfed as part of a cluster, or dentfed as an outler for each of the 123 days (1 = May 1, 2005). The results presented here demonstrate that the clusterng remans consstent over a broad range of parameterzatons (tme of day, date, alttude/atttude). Therefore, clusterng yelds a model that can be utlzed for subsequent complexty analyss. 13

14 Number of trajectores Clustered 4000 Outlers Tme of the day (GMT) Fg. 6 Dstrbuton of clustered and outler flghts by tme of the day Number of trajectores Clustered Outlers Day n the data set Fg. 7 Analyss of clustered and outlers flghts by day of traffc. E. Outlers Whle the vast majorty of the traffc (approxmatvely 80%) are clustered nto flows, 20% are desgnated to be outlers. Outlers are trajectores that are deemed too dfferent from trajectores assgned flow labels by the clusterng algorthm. Too dfferent s a relatve measure that s reflected by parameters of the clusterng algorthm. In order to approach a complete model of the arspace, t s vtal that non-clustered trajectores are accounted for. Outler arcraft are partcularly mportant to assessng ar traffc control complexty because they act as perturbatons to an otherwse statc system. As such, complete traffc models requre ncluson of outler arcraft that are not clustered accordng to standard traffc flows. The outler model s gven n the next secton, wth the complete parameterzed ar traffc model. 14

15 III. Creatng a Parameterzed Generatve Model The prevous secton provded a methodology for clusterng arcraft trajectores nto flows and for determnng whch arcraft can be consdered as outler (non-clustered arcraft). Usng the clusters as a bass, the model s further analyzed to create a parameterzed generatve model. Conceptually, the parameterzed generatve model provdes an nput to generate proxmty complexty maps. Parameterzng the arspace traffc model n ths fashon allows for tools to test what-f scenaros or to create proxmty maps wthout re-clusterng and analyzng the ar traffc data. To parameterze traffc, the arcraft arrval process and the speed dstrbuton for each flow are analyzed and modeled, as well as the outlers trajectores. A. Spatal Flow Dstrbuton The trajectores labeled to a partcular flow by the clusterng algorthm do not follow a sngular trajectory, hence the ntroducton of trajectory wndows and flow tubes. In fact, some flows have a large wdth, up to several hundred mles for Eastbound flght levels wth a large throughput. Fgure 8(a) presents an example cluster, wth correspondng trajectores, wndows, and tube. The vertcal and lateral arcraft dstrbutons for a specfc wndow (resp. for each wndow W k) are llustrated lateral Fgure 8(c). The vertcal and horzontal dstrbutons appear to be decorrelated. For the generatve model, the decorrelaton of vertcal and lateral dstrbuton are assumed. To support ths hypothess, Fgure 9 shows a hstogram of correlaton values between the vertcal and lateral dstrbutons, for each of the 8 wndows for the 685 flows. The dstrbuton of the correlaton values between the vertcal and lateral dstrbutons s narrow, and centered n 0. For 90% of all cluster wndows, the absolute value of the correlaton s less than Ths allows the parameterzed generatve model to use the approxmaton that the vertcal and lateral dstrbuton of arcraft trajectores are ndependent. B. Speed dstrbuton Wthn a gven flow, arcraft fly at dfferent speeds. Factors for such varaton nclude dfferences n arcraft and engne type, arcraft weght, weather, wnd, and plot preferences. Analyss of the flows was conducted to determne a dstrbuton to best model arcraft speeds. Ultmately, the t- 15

16 (a) Example flow Vertcal axs, FL Projecton of trajectores on a wndow Wndow lateral axs, NM (c) Dstrbuton analyss for a flow (b) Dstrbuton analyss at a wndow Fg. 8 Lateral and vertcal dstrbutons of arcraft postons at a wndow wthn a flow. locaton scale probablty dstrbuton was selected. The probablty dstrbuton, f V, for the arcraft speeds n flow F s f V (x) = Γ( ν+1 2 ) σ νπγ( ν 2 ) [ ] ν + ( x µ ( ν+1 σ )2 2 ), (1) ν where µ s the locaton parameter, σ > 0 s the scale parameter, ν > 0 s the shape parameter, and Γ s the standard gamma probablty dstrbuton. If x has a t-scale dstrbuton wth parameters µ, σ and ν, then x ν σ has a Student s t-dstrbuton wth ν degrees of freedom. Fgure 10 presents an example of an emprcal probablty densty functon for arcraft speeds wthn a major flow, as well as the t-locaton scale ft for the emprcal dstrbuton. 16

17 0.2 Correlaton between lateral and vertcal dstrbutons 0.15 Frequency Correlaton Fg. 9 Correlaton between lateral and vertcal dstrbutons of arcraft postons at the wndows. Fg. 10 Dstrbuton of arcraft speeds for the busest flow. C. Arcraft arrval process The clusterng method provdes a spatal representaton of the flows. To obtan a complete traffc model t s necessary that the traffc rate nto each flow s approprately descrbed from hstorcal data. Due to the sparsty of arcraft arrvals nto a specfc flow, t s dffcult to obtan meanngful and verfable dstrbutons for arcraft arrval tmes for ndvdual flows. Therefore, the arcraft arrval process s analyzed over the entre arspace regon (.e. the center level) and proportoned amongst the flows. Because the number of arcraft enterng the arspace regon s a functon of the tme of the day, the objectve of ths secton s to apply the χ 2 test to determne f the arcraft arrval process can be modeled as a non-homogeneous exponental process. Addtonally, to strengthen the modelng process, smlar hypothess testng s appled to the data to check for a Posson countng 17

18 process. Whle the Posson countng process mples an exponental nter-arrval dstrbuton, and vs versa, testng of arrval countng process yelds slghtly dfferent results. However, over the vast majorty of the day, hypothess testng reveals that both arspace models are approprate. For an arbtrary homogeneous Posson countng process, N, defned by parameter λ [events/mnute], the probablty dstrbuton s gven by [26]: P (N(t) N(t s) = k) = (λ)k e λs. (2) k! The assocated cumulatve dstrbuton for the random varable I representng event nter-arrval tmes, s gven by: P (I < a) = 1 e λa. (3) To drve the modelng process, each day s dvded nto 15 mnute ncrements, durng whch tme, the arcraft arrval process s assumed to be statc and follow exponental nter-arrval tmes or a Posson process. The χ 2 goodness-of-ft test s appled to all 15 mnute tme ntervals for both the exponental nter-arrval tme dstrbutons, as well as the Posson countng process. The null hypothess declares that the nter-arrval tmes nto the sector are exponental, wth rate λ j ; or n the case of the countng process, that the dstrbuton follows a Posson process wth parameter λ j τ. The parameter, λ j τ, s the average number of arcraft arrvals durng the j th tme nterval of length τ, over all smlar days of the week. In ths case, τ = 15mn. For a χ 2 -test result of 1, the null hypothess s rejected, mplyng that t s unlkely that the nter-arrval tmes follows an exponental dstrbuton (or Posson countng process) wth rate λ j (λ j τ); otherwse, the hypothess cannot be rejected. The α sgnfcance level was set to 0.05, whch s a standard value for ths test. For the exponental dstrbuton, the bn-sze for χ 2 -test s selected to be 2 seconds. Fgure 11(a) llustrates the traffc varaton over the day accordng to the best-ft 15 mnute parameter λ j τ for the days Wednesday and Sunday. Over the vast majorty of the 24 hour tme span consdered, for Wednesdays (although the same results hold for all other days), the null hypothess s not rejected for ether the exponental dstrbuton nor the Posson process, as llustrated n Fgure 11(b). That s, the arrval process for each 15 mnute nterval follows a best-ft exponental dstrbuton. The dfferences between the values of λ j τ for Wednesdays and Sundays ndcate daly 18

19 changes n ar traffc demand. As an example, for the tme perod between 15:15 and 15:30 Zulu tme, over all Wednesday of the week, comparson the dstrbuton of arcraft nter-arrval tmes nto the sector, wth the best-ft exponental dstrbuton, shows a close ft, as llustrated n Fg. 11(c). λτ: Average number of arcraft arrvals (τ=15 mnutes) λτ estmated, Wed λτ estmated, Sun Arcraft arrval process Tme of the day (hours, ZULU tme) (a) Evoluton of the rate parameter (Zulu). H 0 =1 reject χ 2 results for arcraft counts χ 2 results for arcraft nter arrval tmes Tme of the day (hours, ZULU tme) (b) χ 2 test results for Wednesdays. 1 Dstrbtuon of Arcraft nter arrval tmes ( 15:15 to 15:30 Z) P(T K+1 T k < 1/λ) Data (Wed) Exponental CDF, 1 / λ = nter arrval samples χ 2 test = 0, p = Inter arrval tme 1/λ (sec) (c) Example emprcal and modeled expoental dstrbutons. Fg. 11 Arcraft arrval process nto the center. 19

20 The arcraft nter-arrval process nto the center s therefore modeled as a non-homogeneous posson process over the complete day, wth statc rate parameters for any 15 mnute nterval. The rate parameter changes wth the tme of the day, and day of the week. Addtonally, t s lkely that seasonal changes exst. However for the 4 months consdered, any long-term varatons over same days of the week are small. Whle, the χ 2 test does not valdate that arcraft arrve wth an exponental dstrbuton, the test does verfy a specfc property, enablng modelng of arcraft nter-arrval tmes nto the arspace as exponental. Verfyng that the ar traffc arrval process nto an en route arspace can be modeled as an exponental dstrbuton supports an assumpton found extensvely throughout the lterature (e.g. see [20, 27, 28]). Smlar modelng statements concernng the arrval process nto each flow s not possble, as the number of arcraft arrvals per flow per 15 mnute perod s low, and thus fals to yeld strong results for χ 2 testng. For modelng purposes, we assume that the flows are ndependent, and so are ther arrval processes. The arrval process for the greater arspace s proportoned to accordng to traffc volume along each flow. The rate parameter for the flow F durng tme perod j s λ j = πj λj, where π j wth x j d, s gven by 30 π j = d=1 30 d=1 x j d, x j d, (4) beng the count of arcraft enterng the center durng the tme perod j, on data day d, and dentfed flyng on flow F. D. Outler model Outlers are arcraft that are not clustered nto flows, and they represent approxmately 20% of the traffc. To model outler trajectores, densty maps of arcraft postons are generated. The arspace s dscretzed nto 1 NM squares and 1000 ft hgh cubes. By spatally nterpolatng all arcraft trajectores, the total number of flghts passng through a partcular cube s determned. To generate a probablty dstrbuton, the arcraft counts n each cube are normalzed between 0 and 1 over the entre arspace. Fgure 12 presents densty plots for the locaton of outlers at FL 350 and FL 360. In the fgure, black ndcates a hgh densty of arcraft, and whte a low densty. Wth 20

21 excepton to the locaton where major flows pass, the dstrbuton of outlers over the arspace at a gven flght level s relatvely unform. The spatally probablty dstrbuton of outlers defnes an outler probablty model, f O (P ). Drawng an analogy wth the Kalman flter, the clusters can be nterpreted to be the man sgnal, and the outlers as a nose njected nto the sgnal. Fg. 12 Densty of outlers at FL 350 and FL 360. E. Flow and outler reconstructon Over the last two sectons, clusterng and outler analyss provded models for the spatal dstrbuton of arcraft through arspace. Furthermore, the vertcal and lateral dstrbuton of arcraft postons at each wndow provdes greater granularty n the spatal dstrbuton of arcraft. Coupled wth a model for arcraft speeds, and arcraft nter-arrval tmes nto the arspace, an ar traffc model s generated. The beneft of such a model s that t can be parameterzed to test what-f scenaros wthout reconstructon of the traffc model. That s, adjustments to the model propertes: flow rates, mean arcraft speeds, trajectory wndows, or flow trajectores can be made quckly, thereby defnng a parameterzed generatve model. The nformaton provded by parameterzed generatve model s then used to construct the proxmty maps proposed n ths paper. 21

22 IV. Generaton of arcraft proxmty maps Modelng arspace va arcraft flows and outlers enables the generaton of the presence, conflct, and outler proxmty maps, whose computaton s now presented. These maps use the nformaton generated by the parameterzed generatve model detaled n the prevous sectons. We focus frst on the presence and conflct maps, whch consder only arcraft from modeled flows. Then, we wll ntroduce outler proxmty maps. A. Defnton Consder the East-North-Up frame R 0 = ( a 0, b 0, c 0 ) defned n Secton II and fxed wth respect to the arspace A; and consder any pont P = (x p, y p, y p ) of the arspace. Gven a model of ar traffc parameterzed, the objectve s to compute maps based on the probablty P 1 (P ) that at least one arcraft fles n the proxmty of a pont P and the probablty P 2 (P ) that at least two arcraft from dfferent flows fly n the proxmty of a pont P. These probabltes wll help to generate the presence maps (related to P 1 (P )) and the conflct map (related to P 2 (P )). The probablty P 1 (P ) can be nterpreted as the probablty of presence of arcraft, and P 2 (P ) can be related to probablty of conflct. The shape of the proxmty volume s arbtrary but we consder frst the proxmty volume as the full cylnder P c (P ) centered on P, and defned as follows: P c (P ) = {(x, y, z) A (x x p ) 2 + (y y p ) 2 a 2 p and (z z p ) 2 b 2 p}, (5) wth a p, b p > 0, as llustrated n Fg. 13. Choosng a p = 2.5 NM, b p = 1000ft = 0.165NM reflects the noton of conflct between arcraft: f two arcraft are smultaneously located nsde the cylnder, they are n conflct. For any pont P, the probablty P 1 (P ) (resp. P 2 (P )) s then the probablty that at least one (resp. two) arcraft from dfferent flows are located nsde P c (P ), at any gven tme. However, consderng the proxmty volume P c (P ) gven n Eq. (5) s computatonally too expensve to generate proxmty maps n real-tme, snce the flows are not only a sngle track but a 3D dstrbuton. The cylnder s approxmated by ts crcumscrbng rectangular box. The resultng values of the probabltes P 1 (P ) and P 2 (P ) are then upper bounds on the correspondng values f P c (P ) had been consdered nstead. 22

23 Fg. 13 Proxmty volume P c(p ) for a gven pont P of the arspace. For a gven arspace geometry and flow confguraton, two three-dmensonal proxmty maps, M 1 and M 2, are generated accordng to the values P 1 (P ) and P 2 (P ) assocated wth all ponts P n the arspace. Whle the map M 1 s based on the global presence of arcraft from a flow n the arspace, the map M 2 hghlghts the regons of possble conflct between arcraft from dfferent flows, whch are the most senstve regons for ar traffc controllers. B. Inter-arcraft dstance and resdual dstance of an arcraft From Secton II, a flow F s characterzed by, among other parameters, four probablty densty functons : f a k and f c k for the lateral and vertcal dstrbutons of arcraft at the wndow W k ; f T for the nter-arrval tme dstrbuton between two consecutve arcraft (for clarty, the ndex of the day j s omtted n ths secton); f V for the arcraft speed dstrbuton. To create spatal proxmty maps, the nter-arcraft dstance,.e. the dstance between two consecutve arcraft flyng wthn the same flow F, must be ntroduced, nstead of nter-arrval tme. It wll result a complete knowledge of spatal dstrbuton of arcraft wthn the flow F, along the vectors ( a k, b k, c k ). Denote {AC k, k N } the sequence of arcraft flyng n flow F. The longtudnal dstance d between two consecutve arcraft of F s defned as the dstance between arcraft projecton on the track (or centrod) T (see defnton n Secton II and depcted n Fg. 14(a)) of two consecutve arcraft (AC k, ACk+1 ). The nter-arcraft dstance s a random varable such that d = v t, and has assocated probablty densty functon f D. Assumng that the two random varables v and t are ndependent, f D can be determned from the knowledge of f V and f T. Let d m be the mean nter-arcraft dstance. 23

24 (a) Inter-arrval dstance (b) Resdual dstance Fg. 14 Defntons of the nter-arcraft dstance and resdual dstance. The noton of resdual dstance of the arcraft s now ntroduced, as t s a key to the constructon of the proxmty maps. The noton of resdual dstance for arcraft s smlar to the concept of resdual lfe (or excess lfe ) n renewal theory (see [26] for further detals), and t s defned as follows: consder any pont A located on the track T of flow F (see Fg. 14(b)). Observng traffc at tme t 0 = 0 (we assume that the process starts at tme t = ), let AC k to cross A at any tme t t 0. The resdual dstance r k of the arcraft ACk be the next k th arcraft s the dstance between the projecton of the current poston of arcraft AC k onto T and A. Especally, r 1 s the dstance between the current poston of arcraft AC 1 and A. The varable r 1 s called the frst resdual dstance because t s related to the frst arcraft to cross A. From [26], the PDF of the resdual dstance r k of the arcraft ACk, f R k(rk ), can be derved from the PDF of the nter-arcraft dstance, f D ( d ), and s gven by the followng expresson: f R 1 (r 1 ) = 1 d m ( 1 F D (r 1 ) ), (6) f R k = f R k 1 f D, k 2, (7) where F D s the cumulatve dstrbuton functon of the nter-arcraft dstance assocated to f D, and represents the convoluton operator. To llustrate the noton of resdual dstance, consder a straght flow F and assume the same constant velocty for all arcraft wthn the flow, gven by the maxmum velocty estmate of the speed dstrbuton. From Secton II, t s assumed that the nter-arrval tme dstrbuton f T s an exponental dstrbuton, characterzed by the arrval rate λ. The PDF f D s then also an 24

25 exponental dstrbuton, characterzed by the mean nter-arcraft dstance d m : 1 d e 1 f D ( d ) = m d m d when d 0 0 when d < 0. (8) From Eq. (6), the probablty densty functon of the frst resdual dstance, f R 1, can then be wrtten as: 1 f R 1 (r 1 d e 1 ) = m d m r 1 when r 1 0 (9) 0 when r 1 < 0, whch s smlar to f D for the consdered example. For any k th arcraft AC k, the probablty densty functon of ts resdual dstance r k, f R, can be then derved from Eqs. (7) and (9). k C. Intersecton volume As mentoned above, consderng a cylnder as the proxmty volume requres too many computatons to generate results for a realstc regon of the arspace n real-tme (.e. n a few seconds). Therefore, the proxmty cylnder P c (P ) s approxmated by a crcumscrbng box, whose geometry depends of the box B k consdered. Indeed, gven any k th box B k of flow F and any pont P A, the correspondng proxmty box P k (P ) s defned such that (see Fg. 15): ts center s P, ts dmensons are a p a p b p (5NM 5NM 0.165NM), t s algned wth the centrod of B k. To compute the probabltes P 1 (P ) and P 2 (P ), the frst step s to determne the ntersecton volume V k(p ), defned to be the ntersecton of any box Bk wth Pk (P ),.e. Vk (P ) = Bk P k (P ) as depcted n Fg. 15. Consderng both the dmensons of the proxmty box (5NM 5NM 0.165NM), and the fact that arcraft trajectores are slanted just a few degrees (due to the physcal constrants of the arcraft flght dynamcs), t s approxmated for computatonal smplcty that the consdered box B k s horzontal around P k (P ). Consder the frame R k = ( a k, b k, c k ) correspondng to the box Bk, wth P as the geometrc center. Snce R k s algned wth the box B k, t s also algned wth Vk (P ). The bound of the 25

26 Fg. 15 Defnton of the ntersecton volume V k (P ). ntersecton volume V k (P ) s then defned by ts entry and ext coordnates along the three axes: A and A along a k, B and B along b k, and C and C along c k. These values can be easly determned from the equatons of the proxmty volume P k (P ) and of the box Bk. D. Presence maps Presence maps consder the probablty that one or more arcraft s located near a pont P,.e. nsde the consdered proxmty volume. Gven any pont P A, consder frst a sngle box B k, and the proxmty box P k(p ). Let P 1(P, B k) be the probablty that at least one arcraft AC from box B k be nsde P k (P ). Snce the three probablty dstrbutons for the presence of an arcraft along each axs are assumed to be ndependent (see Secton II), the probablty P 1 (P, B k ) can be wrtten P 1 (P, B k ) = P( AC P k (P )) = P( AC [AA]).P( AC [BB]).P( AC [CC]). (10) So far, only the lateral and vertcal PDFs at the wndows W k (.e. at the entrance of the boxes B k) have been defned and determned n Secton II: f a k and f c k. Snce the ntersecton volume V k(p ) s not necessarly near W k, the arcraft spatal dstrbutons wthn the flow F must be extended around the ponts A,C. Therefore, a lnear nterpolaton of the lateral (resp. vertcal) dstrbutons f a k and f a k+1 (resp. f c k and f c k+1) between the wndows W k and W k+1 s used to get f a k /A (resp. f c k /C),.e. the lateral (resp. vertcal) arcraft dstrbuton at the pont A (resp. C). The lnear nterpolaton process approxmates the smooth varaton of arcraft presence densty wthn the 26

27 flow, preservng the nature of a probablty densty functon. Any term n the rght-hand sde of Eq. (10) can be computed as a functon of the probablty densty functons prevously defned, f a k /A, f c k /C, and f R 1 : P( AC [AA]) = P( AC [CC]) = P( AC [BB]) = A A C C BB 0 f a k /A(x) dx, (11) f c k /C(z) dz, (12) f R 1 (y) dy, (13) where BB = BB. The probablty P 1(P, B k ) can then be determned from Eqs. (10) (13). The prevous defntons and equatons can naturally be extended to the l 1 boxes B k of the flow F. It leads to ntroduce the probablty, P 1 (P, F ), that there s at least one arcraft AC from flow F flyng at proxmty of the pont P. The probablty s computed from Eqs. (10) (13) takng the l 1 boxes B k nto account. Consder now the entre arspace A modeled wth N ndependent flows. The presence probablty of at least one arcraft from a flow about any pont P, P 1 (P ), can then be wrtten N N ( P 1 (P ) = 1 P( [1, N], AC P (P )) = 1 P( AC P (P )) = 1 1 P1 (P, F ) ). =1 =1 (14) E. Conflct maps Conflct maps determne the probablty that at least two arcraft from dfferent flows are smultaneously n the same regon of the arspace. More formally, the probablty P 2 (P ), that at least two arcraft from dfferent flows are n the proxmty volume surroundng P, s wrtten : P 2 (P ) = P ( (, j) [1, N] 2, j, s.t. (AC P (P )) & (AC j P j (P ) ). Usng the probablty P 1 (P ) defned n Secton IV D, P 2 (P ) can be wrtten P 2 (P ) = P 1 (P ) P( [1, N], (AC P (P )) & (AC j / P j (P ), j )) N N = P 1 (P ) P 1 (P, F ) (1 P 1 (P, F j )). =1 j=1,j 27

28 F. Outler proxmty maps Presence and conflct maps consder only arcraft that belong to flows. For a complete modelng of the traffc, outlers non clustered arcraft should be consdered because they account for 20% of the overall traffc and may have an mportant nfluence on the domnant traffc. The spatal dstrbuton of the outlers, f O (P ), s gven by the parameterzed generatve model presented n Secton III. Ths leads to the ntroducton of outler proxmty maps, denoted M O. Outlers proxmty maps estmate the nfluence of outlers on the modeled flows that represent the majorty of traffc nsde the arspace under consderaton. Therefore the probablty P O (P ), that one arcraft from the flows and one outler are smultaneously n the same neghborhood of any pont P of the arspace, s determned. Only one outler s consdered, snce the probablty of two outlers to be n proxmty too each other s neglgble. Usng a methodology smlar to the one used n Sectons IV D and IV E, P O (P ) s approxmated by 1 P O (P ) = P 1 (P ). a 2 P b P P O (P ) f O (P ) dτ, where P O (P ) s a parallelogram whose dmensons are a P a P 2b P, and whose center s P. The spatal dstrbuton of the outlers, f O (P ), can be desgned (.e. smulated) to evaluate the nfluence of the presence of rogue arcraft (consdered as external dsturbances) on a nomnal ar traffc consstng of the flows n a gven regon of the arspace: ths approach extends the works n [12, 13] from ndvdual arcraft to entre arcraft flows. V. Generaton of proxmty maps for Cleveland center ZOB Ths secton combnes both hstorcal flow parameters determned from the clusterng algorthm presented n Sectons II and III, and the methods to construct proxmty maps gven n Secton IV. The maps use the flows characterstcs computed for May and the tme perod between 2:00p.m. and 2:15p.m. local tme. Durng ths tme perod, the expected number of enterng arcraft s 65. For easer vsualzaton of the results, proxmty maps have been computed for the horzontal slce of arspace between FL 310 and FL 350. In ths regon, the mean nter-arcraft dstances for the 10 busest flows are gven n Table 1, and the centrods of the flows are plotted n red on Fgs. 16(a) and 17(a). Fgures 16 and 17 llustrate the presence and conflcts maps, respectvely. A 3D vew 28

29 and a 2D vew are presented for each flght level. The outlers proxmty maps for the consdered center are gven n Fg. 18, whch hghlghts the regons of the arspace that may be senstve to the presence of an outler. These maps provde an ntuton of the mportance of each flow, how they nteract wth each other, and how much nfluence an outler can have on the man traffc. Snce these maps provde a depcton of the probablty of proxmty of arcraft, they may be used n combnaton wth a forecasted weather map, n order to estmate the regons where bad weather may have strong mpact on the predcted traffc (as llustrated n [21]). The maps can also be used as a tool for traffc flow managers to redesgn routes or adjust flow rates n order to allow ncreased throughput whle keepng a low level of proxmty probablty. Flow number, d m (NM) Table 1 Mean nter-arcraft dstance, Cleveland center ZOB between 2:00pm-2:15pm VI. Concluson Ths paper s concerned wth computng arspace statstcal complexty metrcs usng a parameterzed generatve model drven by hstorcal flght data. Trajectores are clustered nto major flows that account for 80% of the traffc and the remanng 20% are dentfed as outlers. These flows and outlers consttute a structure of ar traffc, used as a parameterzed generatve model. Wthn the flows, arcraft spatal and temporal probablty dstrbutons are computed. Analyss of the arspace demonstrates that the arcraft arrval process nto the center can be modeled as tme-varyng Posson process. Usng the structure determned by the flow and outlers model, three-dmensonal maps are constructed to predct and vsualze the probablty of presence and probablty of proxmty of arcraft n the arspace. The proxmty maps are nformaton dsplays that are relevant to traffc flow management. Each map (presence, conflct, or outler proxmty map) dentfes and hghlghts crtcal regons of the arspace that requre sgnfcant montorng. Further, because the maps are based on parameterzed hstorcal traffc models, t s possble to consder what-f scenaros, n whch parameters (flow rates and dstrbutons) are adjusted n real-tme to consder daly traffc varatons. Ultmately, 29

30 (a) 3D vew (b) 2D vew at FL 310 (c) 2D vew at FL 350 (d) 2D vew at FL 380 Fg. 16 Presence maps of Cleveland center ZOB between 2:00pm-2:15pm. ths tool provdes an nformaton dsplay to ad n the strategc and operatonal decson makng requred n traffc flow management. Acknowledgements Ths work s funded by NASA under Grant NNX08AY52A and the FAA under Award No. 07-C- NE-GIT, Amendment Nos. 005, 010, and 020. References [1] NextGen Integraton and Implementaton Offce. Nextgen mplementaton plan. Techncal report, Federal Avaton Admnstraton,

31 (a) 3D vew (b) 2D vew at FL 310 (c) 2D vew at FL 350 (d) 2D vew at FL 380 Fg. 17 Conflct maps of Cleveland center ZOB between 2:00pm-2:15pm. (a) 2D vew at FL 310 (b) 2D vew at FL 350 Fg. 18 Outler proxmty maps of Cleveland center ZOB between 2:00pm-2:15pm. 31