Access the full text.

Sign up today, get DeepDyve free for 14 days.

Pongsakornsathien, Nichakorn; Bijjahalli, Suraj; Gardi, Alessandro; Symons, Angus; Xi, Yuting; Sabatini, Roberto; Kistan, Trevor

Aerospace
, Volume 7 (11) – Oct 28, 2020

/lp/multidisciplinary-digital-publishing-institute/a-performance-based-airspace-model-for-unmanned-aircraft-systems-zmeWpe8JuP

- Publisher
- Multidisciplinary Digital Publishing Institute
- Copyright
- © 1996-2020 MDPI (Basel, Switzerland) unless otherwise stated Disclaimer The statements, opinions and data contained in the journals are solely those of the individual authors and contributors and not of the publisher and the editor(s). MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations. Terms and Conditions Privacy Policy
- ISSN
- 2226-4310
- DOI
- 10.3390/aerospace7110154
- Publisher site
- See Article on Publisher Site

aerospace Article A Performance-Based Airspace Model for Unmanned Aircraft Systems Trac Management 1 1 1 1 Nichakorn Pongsakornsathien , Suraj Bijjahalli , Alessandro Gardi , Angus Symons , 1 1 , 1 , 2 Yuting Xi , Roberto Sabatini * and Trevor Kistan School of Engineering, RMIT University, Melbourne, VIC 3000, Australia; nicha.pongsakornsathien@rmit.edu.au (N.P.); suraj.bijjahalli@rmit.edu.au (S.B.); alessandro.gardi@rmit.edu.au (A.G.); s3601602@student.rmit.edu.au (A.S.); s3614421@student.rmit.edu.au (Y.X.); trevor.kistan@thalesgroup.com.au (T.K.) Airspace Mobility Solutions, THALES Australia, Melbourne, VIC 3000, Australia * Correspondence: roberto.sabatini@rmit.edu.au Received: 26 September 2020; Accepted: 21 October 2020; Published: 28 October 2020 Abstract: Recent evolutions of the Unmanned Aircraft Systems (UAS) Trac Management (UTM) concept are driving the introduction of new airspace structures and classiﬁcations, which must be suitable for low-altitude airspace and provide the required level of safety and ﬂexibility, particularly in dense urban and suburban areas. Therefore, airspace classiﬁcations and structures need to evolve based on appropriate performance metrics, while new models and tools are needed to address UTM operational requirements, with an increasing focus on the coexistence of manned and unmanned Urban Air Mobility (UAM) vehicles and associated Communication, Navigation and Surveillance (CNS) infrastructure. This paper presents a novel airspace model for UTM adopting Performance-Based Operation (PBO) criteria, and speciﬁcally addressing urban airspace requirements. In particular, a novel airspace discretisation methodology is introduced, which allows dynamic management of airspace resources based on navigation and surveillance performance. Additionally, an airspace sectorisation methodology is developed balancing the trade-o between communication overhead and computational complexity of trajectory planning and re-planning. Two simulation case studies are conducted: over the skyline and below the skyline in Melbourne central business district, utilising Global Navigation Satellite Systems (GNSS) and Automatic Dependent Surveillance-Broadcast (ADS-B). The results conﬁrm that the proposed airspace sectorisation methodology promotes operational safety and eciency and enhances the UTM operators’ situational awareness under dense trac conditions introducing a new eective 3D airspace visualisation scheme, which is suitable both for mission planning and pre-tactical UTM operations. Additionally, the proposed performance-based methodology can accommodate the diversity of infrastructure and vehicle performance requirements currently envisaged in the UTM context. This facilitates the adoption of this methodology for low-level airspace integration of UAS (which may dier signiﬁcantly in terms of their avionics CNS capabilities) and set foundations for future work on tactical online UTM operations. Keywords: Urban air mobility; air trac management; unmanned aircraft system; UAS trac management; airspace; sectorization; performance-based navigation; GNSS; ADS-B 1. Introduction The emerging need to integrate Unmanned Aircraft Systems (UAS) trac in the existing airspace bears notable challenges in terms of increased trac density and complexity in low-altitude airspace [1,2]. In 2012, the Federal Aviation Administration’s (FAA) Modernization and Reform Act prompted research into the ﬁeld of managing small UAS (sUAS) as civilian demand for these systems had Aerospace 2020, 7, 154; doi:10.3390/aerospace7110154 www.mdpi.com/journal/aerospace Aerospace 2020, 7, 154 2 of 25 grown signiﬁcantly [3]. This initiative was pertinent to urban operations [4–6] and set foundations for successive research eorts in this ﬁeld. In 2016, the FAA released Part 107, an addition to Title 14 of the Code of Federal Regulations (14 CFR), allowing routine sUAS operations in the National Airspace System (NAS) [7]. To date, UAS operations have largely been segregated from manned aircraft. In instances where operations spanned controlled airspace, separation from other aircraft has been ensured procedurally through an approval-seeking process with aviation authorities, and visually by the remote pilot in charge. Current regulations allow Visual Line of Sight (VLOS) operations within both segregated and unsegregated airspace (subject to a lengthy risk assessment process). However, large scale VLOS and autonomous Beyond VLOS (BVLOS) operations are not allowed within low altitude airspace. To prevent interference with manned aircraft, sUAS operations are limited to the open space below 400 ft (122 m) above ground level. However, with the rapid growth of civilian UAS use-cases and operations, it is evident that a dedicated infrastructural framework is required as trac volumes scale up. The UAS Trac Management (UTM) concept proposed and spearheaded by NASA aims to provide a framework for safely and eciently managing the demand and complexity of future operations [8]. The UTM framework proposed and demonstrated through various trials intends to support and enhance UAS operations in both VLOS and BVLOS modalities [9,10]. The UTM project is structured into four distinct phases or Technical Capability Levels (TCL). The progression of the project is characterized by increasing scenario complexity and required speciﬁc autonomous capabilities [6,11]. Each new TCL extends the capabilities of the previous TCL, with each sequential phase supporting a larger range of UAS from remotely piloted vehicles to fully autonomous UAS. Each capability is targeted to speciﬁc types of applications, geographical areas, and use cases that represent certain risk levels. TCL4 is characterized by complex operations in densely populated urban areas and Urban Air Mobility (UAM). The beneﬁts of UAS and manned/unmanned operations within these dense built-up areas are being captured by the NASA Advanced Air Mobility (AAM) working groups, which have already identiﬁed numerous public and commercial use-cases [12]. However, several challenges must be addressed to fully realise these beneﬁts. These include Beyond Line-Of-Sight (BLOS) operations, large-scale contingency management, in-ﬂight deconﬂiction and trajectory conformance monitoring. These new challenges and the safety-critical nature of urban operations will necessarily impose stringent requirements on the infrastructure supporting UTM operations. Some of the open gaps include airspace design, navigation/surveillance infrastructure, communication networks and trac control/management [13]. In particular, the co-design of coherent airspace structures and ecient sectorisation strategies remains an open research question, with no standardized criteria forming the basis for partitioning the future urban/suburban airspace into functioning operational structures. Airspace sectorisation is a key factor in managing air trac complexity, human operator workload, and demand-capacity balance [14]. In terms of capacity modelling, demand prediction and demand-capacity balance, a large body of research is available pertaining manned aircraft and conventional Air Trac Management (ATM) [14–19]. Recent developments include the dynamic morphing of airspace sectors to mitigate the disruptions associated with capacity saturation [20]. However, owing to the dierence in operational complexity, trac volumes, ﬂeet mixes, and supporting infrastructure, it is readily apparent that airspace design and sectorisation strategies for conventional air trac are inapplicable in their current form to the UAS/UAM trac management problem. While various sectorisation concepts and models have been proposed in the literature, as of today, none has been standardized or even agreed upon to cover the whole spectrum of UTM operations. The doctrine applied in this paper is that the required separation of UTM trac and subsequently, the airspace management strategy should be a direct consequence of the infrastructure supporting the operation, as originally argued in [14]. These include the hardware and software systems operating in the air (avionics systems) and on the ground (trac management systems) as well as the personnel exercising oversight and control over the operations. This leverages the opportunity to adopt some of the existing frameworks that have been applied or recommended for conventional air trac. Aerospace 2020, 7, 154 3 of 25 In particular, the concept of Performance-Based Operations (PBO), and more speciﬁcally Performance Based Navigation (PBN), proves suitable in safely accommodating the highly dynamic nature of operations as envisaged under the UTM framework. The key contribution of this paper is a novel UTM airspace model for urban areas designed to enhance the safety and eciency of trac spacing and to maximise exploitation of the airspace resources (capacity). The urban airspace model proposed here is speciﬁcally conceived to accommodate a high diversity in performance and operational characteristics of dierent UAM/UAS platforms. The airspace design and sectorization strategies are centred on the discretization of the airspace into elementary three-dimensional (3D) cells. The dimension of each cell is a function of the performance of the infrastructure supporting a given airspace region. As a consequence, the airspace structure and sector volumes are driven by the dynamically varying performance of the separation services provided. In this paper, particular emphasis is given to navigation performance, as it also aects other systems. However, the presented framework is generally applicable to the entire set of Communication, Navigation and Surveillance (CNS) systems deployed and available in a given scenario. For veriﬁcation purposes, the presented models are applied to realistic UTM/UAM scenarios, including an urban canyon given the signiﬁcant challenges this poses to navigation and cooperative surveillance performance. In particular, a sensitivity analysis of airspace sector capacity is performed considering the computational complexity of trajectory planning/re-planning and communication overheads associated with sector handovers. This paper builds upon the revisited Dynamic Airspace Management (DAM) concept introduced in [14] and presents some signiﬁcant evolutions from the initial concept introduced in [21], which are aimed to support an advanced management of low-altitude trac exploiting high levels of autonomy, for which new human-machine interaction strategies are being developed [22]. The paper is structured as follows: Section 2 presents a review of prior work in airspace design and partitioning strategies and the Concept of Operations (CONOPS) for UAS/UAM operations. The proposed airspace model is presented in Section 3. This includes the strategy to discretise the airspace and generate sectors, along with the algorithmic implementation of the model. Section 4 applies the proposed model to two veriﬁcation case studies. Section 5 provides the conclusions of this study and future research directions. 2. UTM Airspace Design and Operational Concepts Airspace structure design is an important aspect of UAM/UAS integration and to date remains an open research gap in the literature. Although lessons learned in conventional ATM are useful, many of its solutions cannot be scaled down and directly applied to the local management of UAM and UAS [12]. Table 1 summarizes the main airspace structure concepts that have been proposed in the literature so far. Four concepts are readily discernible: full mix, layers, zones and tubes. These are described in detail in [23,24]. In the full mix concept, the airspace is essentially unstructured, with vehicles using a direct route between their origin and destination, with optimised ﬂight altitudes and velocities to minimise fuel burn and other related costs. In the layered concept, the airspace is divided into sections stacked on top of each other. Each layer is typically dedicated to a type of aerial vehicle and a speciﬁc heading. The zone concept divides the airspace into dierent regions on the basis of a set of criteria. This typically includes risk exposure to people and property. Lastly, tube-based airspace structures are centred on ﬁxed corridors; vehicles must follow a common speed limit, thus ensuring safe separation and minimizing conﬂicts. The so-called “stacked-layer” concept is also gaining traction, with one proposal to split each layer into so-called sky lanes in urban centres that mimic the organisation of streets [13]. Each lane is essentially a corridor that encompasses a reference nominal trajectory with sucient buers to account for potential o-nominal situations. The separation assurance from other vehicles and static obstacles is maintained by onboard autonomous systems and/or remote pilots. The ﬂow of trac is essentially managed by oine or online design/redesign of the lanes. Situations in which lanes intersect are managed by a time-based prioritization for each lane, conceptually similar to how trac lights are applied to road trac. Another layered airspace concept which is analogous to Aerospace 2020, 7, 154 4 of 25 road trac is presented in [25,26], where roundabouts are used at the junction of multiple sky lanes, in place of time-based separations. Table 1. Summary of proposed airspace structure concepts. Airspace Concept Structure Type Layout Beneﬁt Save space Lanes Jang D.S. et al. [13] Avoid congestion at intersections Trac lights at intersections Simple Lanes Simple Sacharny S. et al. [25] Roundabouts at intersections Decrease trac delays at intersections Layer Sedov L. et al. [26] Lanes Consider vertical separation loss Simple Lanes Duchamp V. et al. [27] Flexible routes Corridors PBN-inspired density-based layers Lanes Balance safety and eciency metrics Tubes and cones for landing Free ﬂights Full mix Very ﬂexible ﬂights Landing/takeo strips Sunil E. et al. [23,24] Lanes First Come–First Serve (FCFS) principle Zone Octagon arcs at intersections Simple to implement Fix routes Tube Increase trac predictability Layers Vertical zone Amazon [28] Zone Simple Speed categorisations The zone concept has also been proposed in the literature as a means of dividing the airspace into operational volumes. In several proposals in the literature, the division is typically performed on the basis of the level of risk exposure to the general public. Alternatively, the division can also be performed on the basis of the services provided, and on the level of overall system performance required to support a given category of operation. Access to a particular zone is then contingent on the UAS meeting the level of system performance stipulated for that zone. This mirrors the implementation of PBN for manned aircraft, wherein the employed navigation systems are required to meet a certain level of performance depending on the region and phase of ﬂight. In fact, PBN has been endorsed by NASA [29], as a potentially viable concept to adopt and apply to the airspace structuring problem in the UTM context. This is also emphasised in the NASA AAM initiative [30]. Navigation performance for UAS can vary to a greater extent than manned aircraft. Most UAS navigation systems employ a Global Navigation Satellite System (GNSS) receiver as the primary source of positioning in a global frame. This is typically supplemented by fusion with other sensors including inertial sensors, visual sensors and lidar to obtain a full navigation state estimate. The achievable performance is dependent on the individual sensor characteristics as well as the employed fusion algorithm. Performance is also dependent on environmental characteristics that are dierent from those encountered in manned aircraft operations. For example, GNSS performance is highly degraded in urban environments owing to signal multipath and obscuration relative to conventional manned aircraft operations. Therefore, greater reliance on augmentation with visual sensors and intelligent fusion algorithms is necessary. A successful adaptation of PBN to the UTM context would support this widely varying performance. Further, the PBN concept can be extended more generally to include not only navigation performance but also (as a minimum), communication and surveillance performance. PBN would then be one element of the broader concept of Performance-Based Operations (PBO). However, there has been limited investigation so far into this line of research. Duchamp V. [27] presents a PBN-inspired approach using graph algorithm to separate UAS trac into dierent levels; congested city zones are divided into a larger number of layers and in contrast, low-density airspace is structured Aerospace 2020, 7, x FOR PEER REVIEW 5 of 25 navigation performance but also (as a minimum), communication and surveillance performance. PBN would then be one element of the broader concept of Performance-Based Operations (PBO). However, there has been limited investigation so far into this line of research. Duchamp V. [27] Aerospace 2020, 7, 154 5 of 25 presents a PBN-inspired approach using graph algorithm to separate UAS traffic into different levels; congested city zones are divided into a larger number of layers and in contrast, low-density airspace with a lower resolution. In [31], the airspace was partitioned into so-called ‘airboxes’ considering is structured with a lower resolution. In [31], the airspace was partitioned into so-called ‘airboxes’ communication performance. considering communication performance. 2.1. Urban UAM/UAS Operational Concept Urban UAM/UAS Operational Concept For On-demand Passenger Air Transportation (OPAT), Nneji et al. [32] present a CONOPS that For On-demand Passenger Air Transportation (OPAT), Nneji et al. [32] present a CONOPS that is conceptually similar to current aviation practice, with the pilot taking on the responsibilities of is conceptually similar to current aviation practice, with the pilot taking on the responsibilities of pre- pre-takeo checks, communication with ATC and post-landing checks. This type of CONOPS is takeoff checks, communication with ATC and post-landing checks. This type of CONOPS is currently currently gaining traction and industry interest [33]. This forms the foundation on which more gaining traction and industry interest [33]. This forms the foundation on which more advanced advanced concepts such as “Revolutionary Vehicle Autonomy” (RVA) and “Evolutionary Vehicle concepts such as “Revolutionary Vehicle Autonomy” (RVA) and “Evolutionary Vehicle Autonomy” Autonomy” (EVA) are based [32]. Two representative operational cases are described in this section: (EVA) are based [32]. Two representative operational cases are described in this section: flight over ﬂight over skyline (UAM air-taxi operations) and ﬂight below skyline (UAS delivery operations). skyline (UAM air-taxi operations) and flight below skyline (UAS delivery operations). These two These two cases are also addressed in the veriﬁcation studies presented in Section 4. For UAM cases are also addressed in the verification studies presented in Section 4. For UAM traffic, it is trac, it is recommended to adapt helicopter routes by bounded the trajectories as corridors [34]. recommended to adapt helicopter routes by bounded the trajectories as corridors [34]. However, However, simply adapting current helicopter routes will not respond to the current demands for simply adapting current helicopter routes will not respond to the current demands for Mobility-as- Mobility-as-a-Service (MAAS). Hence, in the foreseeable future, it is expected that the user will be able a-Service (MAAS). Hence, in the foreseeable future, it is expected that the user will be able to submit to submit a request by only specifying the intended departure and arrival points. The system will then a request by only specifying the intended departure and arrival points. The system will then calculate calculate the most suitable route to accomplish this service. Under these proposed concepts, the ﬂight the most suitable route to accomplish this service. Under these proposed concepts, the flight trajectory will be autonomously calculated and evaluated by its governing UTM system, allowing for trajectory will be autonomously calculated and evaluated by its governing UTM system, allowing for minimal human intervention. This is conceptually illustrated in Figure 1. The operation starts when minimal human intervention. This is conceptually illustrated in Figure 1. The operation starts when the passenger submits a request for a trip through a UTM service supplier interface. The UTM system the passenger submits a request for a trip through a UTM service supplier interface. The UTM system with a UTM operator in the loop generates an intended route considering factors such as airspace with a UTM operator in the loop generates an intended route considering factors such as airspace capacity, trac management, weather, time eciency, etc. After the intended routes are conﬁrmed, the capacity, traffic management, weather, time efficiency, etc. After the intended routes are confirmed, passenger boards the UAM platform, and the mission is executed autonomously. The ﬂight can span the passenger boards the UAM platform, and the mission is executed autonomously. The flight can multiple regions with diering safety objectives and supporting infrastructure e.g., a ﬂight can span span multiple regions with differing safety objectives and supporting infrastructure e.g., a flight can the airspace over suburban, urban and airport regions. Emergency scenarios call for the provision of span the airspace over suburban, urban and airport regions. Emergency scenarios call for the system features that allow human operators to initiate and execute contingency measures leading to a provision of system features that allow human operators to initiate and execute contingency safe and rapid recovery [32]. measures leading to a safe and rapid recovery [32]. UTM System Cruise TMA Airspace UTM Operator Hover Hover up down Decision making ATCO Approved Trip Ground Vertiport Vertiport Airport requests Transport Figure 1. Concept of Operations for a conventional “air taxi”’. Figure 1. Concept of Operations for a conventional “air taxi”. For UAS-based delivery operations, the envisaged CONOPs is as follows: For UAS-based delivery operations, the envisaged CONOPs is as follows: Phase 1—the remote pilot submits the intended mission profile to the UTM Service Supplier Phase 1—the remote pilot submits the intended mission proﬁle to the UTM Service Supplier (USS) and (USS requests ) and re clearance quests clto eara access nce to the aintended ccess the operational intended o airspace perationand al ai execute rspace the andmission; execute the mission; Phase 2—The USS checks for conﬂict with other operations, and assesses whether the UAS and its corresponding support infrastructure meets the minimum required performance requirements for the intended operational airspace; Upon clearance of these checks, the submitted request is approved; Aerospace 2020, 7, 154 6 of 25 Phase 3—The UAS takes-o and climbs to its cruise altitude; Phase 4—When the UAS is enroute to the delivery location, system performance is monitored and potential threats are assessed by both remote pilot and onboard/ground-based autonomous systems; Phase 5—Upon arrival at the delivery location, the remote pilot conﬁrms a clearance to land through UTM service, which also triggers a notiﬁcation to the customer; Phase 6—The UAS delivers the package. This can be accomplished in dierent ways. A cable and tether mechanism is typically employed; Phase 7—The customer retrieves the parcel; Phase 8—The remote pilot guides the unmanned aircraft back to base, again subject to a conﬂict check and clearance from the UTM service Human operators are required to be in the loop for assuring safety, especially during emergency procedures. However, during nominal conditions, the high level of automation allows the operator(s) to take on a strategic management role, with most low-level (and increasingly high-level) tasks being performed autonomously. 3. Urban Airspace Model The proposed performance-based airspace model is conceptually illustrated in Figure 2. The model is underpinned by a two-stage approach. During the oine planning stage, the reference grid is generated as a set of elementary cells according to the baseline CNS performance. In our approach, the elementary cells are parameterized as cuboids with a square base, but the methodology could easily Aerospace 2020, 7, x FOR PEER REVIEW 7 of 25 be extended to other 3D shapes. The length (r ), width (r ) and height (r ) of the cuboidal cells are x y z formulated as a function of the expected baseline CNS performance for the region: The total latency of position information is the Transaction delay between the time of applicability of the r , r , r = f C , N , S (1) x y z perf perf perf Latency time position measurement and the time of arrival of the ADS-B message for that position. where C , N , S 2 R are the set of performance metrics for communication, navigation, and perf perf perf NIC categorisation Navigational Integrity Category specifies the surveillance respectively. A summary of applicable metrics which quantify the CNS performance (only integrity of containment radius aligned with is provided in Lo T ca able lisatio 2.n The model is applied in two timeframes: oine planning and online RMIT Classification: Trusted ADS-B) horizontal position which then maps to RNP. airspace management. Tracking error (σ) The standard deviation of the tracking error. Offline planning Online active UAS management Specific CNS Baseline CNS Performance Performance Elementary cell Occupancy CNS Protection dimension Grid Volume Airspace sectorisation Interface Sector Figure 2. High-level overview of the performance-based airspace model. Figure 2. High-level overview of the performance-based airspace model. For communication systems in aviation applications, one measure of performance is the latency in staging an intervention due to human operator and equipment limitations. These are typically not measured in real-time but are determined apriori either analytically, or more frequently, using experimental data. For navigation systems, under the PBN framework, performance is quantified using four metrics [35]: Accuracy Integrity Continuity Availability System accuracy requirements are specified as an upper bound on position error, which is the difference between the estimated position and the actual position. Navigation state estimators typically output the standard deviation of the computed position solution as a measure of confidence. A 2σ value which bounds approximately 95% of the navigation errors is taken as an online indicator of accuracy. Integrity is a measure of trust that can be placed in the correctness of the positioning solution supplied by the navigation system. Integrity includes the ability of a system to provide timely and valid warnings to the user (alerts) when the system must not be used for the intended phase of flight. To ensure integrity in practice, the standard deviation of the computed position solution is inflated by a requisite factor to bound all but a small fraction of errors (corresponding to an ‘integrity’ risk or maximum allowable probability that the solution is out of bounds). The inflation also accounts for the maximum allowable false alarm rate of the system, which dictates the continuity of the system. Continuity is a measure of reliability that ensures that the system will perform Surveillance Aerospace 2020, 7, 154 7 of 25 Table 2. CNS Performance metrics. Uncertainty CNS Metric Description Component Human latency includes all human operator response Human latency time such as UTM operator decision-making and response time. Communication Transaction time Technical latency refers to datalink latency and control Technical latency link latency. The standard deviation of the position solution. This is Position inﬂated by multiplicative factors to meet accuracy, uncertainty () integrity, and continuity requirements. Position error Dilution of precision DOP is a ratio factor in positioning which is calculated Navigation sources (only for GNSS) from the satellites-receiver geometry. Drift rate (dead Drift stems from errors in acceleration and angular reckoning system) velocity measurement The total latency of position information is the delay between the time of applicability of the position Transaction time Latency measurement and the time of arrival of the ADS-B message for that position. Surveillance Navigational Integrity Category speciﬁes the integrity NIC categorisation of containment radius aligned with horizontal position (only ADS-B) Localisation which then maps to RNP. Tracking error () The standard deviation of the tracking error. For communication systems in aviation applications, one measure of performance is the latency in staging an intervention due to human operator and equipment limitations. These are typically not measured in real-time but are determined apriori either analytically, or more frequently, using experimental data. For navigation systems, under the PBN framework, performance is quantiﬁed using four metrics [35]: Accuracy Integrity Continuity Availability System accuracy requirements are speciﬁed as an upper bound on position error, which is the dierence between the estimated position and the actual position. Navigation state estimators typically output the standard deviation of the computed position solution as a measure of conﬁdence. A 2 value which bounds approximately 95% of the navigation errors is taken as an online indicator of accuracy. Integrity is a measure of trust that can be placed in the correctness of the positioning solution supplied by the navigation system. Integrity includes the ability of a system to provide timely and valid warnings to the user (alerts) when the system must not be used for the intended phase of ﬂight. To ensure integrity in practice, the standard deviation of the computed position solution is inﬂated by a requisite factor to bound all but a small fraction of errors (corresponding to an ‘integrity’ risk or maximum allowable probability that the solution is out of bounds). The inﬂation also accounts for the maximum allowable false alarm rate of the system, which dictates the continuity of the system. Continuity is a measure of reliability that ensures that the system will perform nominally without interruption for a given mission segment. Availability is the proportion of the overall mission time for which the navigation system meets accuracy, integrity and continuity requirements. For GNSS, the well-known Dilution of Precision (DOP) factors are also used as an easily obtained online measure of performance. Aerospace 2020, 7, x FOR PEER REVIEW 8 of 25 nominally without interruption for a given mission segment. Availability is the proportion of the o Aer ver ospace all 2020 miss , 7 io , 154 n time for which the navigation system meets accuracy, integrity and contin 8ui ofty 25 requirements. For GNSS, the well-known Dilution of Precision (DOP) factors are also used as an easily obtained online measure of performance. For surveillance systems, as in the case of navigation, state estimators used in radar such as the For surveillance systems, as in the case of navigation, state estimators used in radar such as the Extended Kalman Filter or other recursive algorithms output the standard deviation of tracking error. Extended Kalman Filter or other recursive algorithms output the standard deviation of tracking error. As in the case of navigation, this metric can be inﬂated to bind a required fraction of the error. As in the case of navigation, this metric can be inflated to bind a required fraction of the error. In the online stage, the speciﬁc CNS performance of each UAS is estimated and used to generate In the online stage, the specific CNS performance of each UAS is estimated and used to generate a virtual protection volume bounding the vehicle. This volume is the CNS protection volume, which a virtual protection volume bounding the vehicle. This volume is the CNS protection volume, which comprises the combined uncertainty of the CNS infrastructure supporting the UAS. This represents comprises the combined uncertainty of the CNS infrastructure supporting the UAS. This represents a protected volumetric bound which should not be breached by any obstacle to ensure a sucient a protected volumetric bound which should not be breached by any obstacle to ensure a sufficient margin for Separation Assurance and Collision Avoidance (SA/CA). The methodology to generate the margin for Separation Assurance and Collision Avoidance (SA/CA). The methodology to generate individual components of this volume is presented in Section 3.2. the individual components of this volume is presented in Section 3.2. An occupancy grid can successively be derived from the set of elementary cells and from the An occupancy grid can successively be derived from the set of elementary cells and from the active CNS protection volumes. Cells that are contained within and on the boundary of the protection active CNS protection volumes. Cells that are contained within and on the boundary of the protection volume correspond to occupied space and the remaining cells correspond to unoccupied space. This is volume correspond to occupied space and the remaining cells correspond to unoccupied space. This also illustrated in Figure 2. is also illustrated in Figure 2. The occupancy grid supports the demand-capacity balancing by not just considering the number The occupancy grid supports the demand-capacity balancing by not just considering the number of aircraft simultaneously active in a region, which was the traditional approach for human-centric of aircraft simultaneously active in a region, which was the traditional approach for human-centric ATM, but also their speciﬁc CNS performance, which supports more automated and autonomous ATM, but also their specific CNS performance, which supports more automated and autonomous SA/CA concepts. Airspace sectors are then generated as clusters of elementary cells to support ecient SA/CA concepts. Airspace sectors are then generated as clusters of elementary cells to support efficient management of trac across the urban region. These sectors need to consist of an optimal number management of traffic across the urban region. These sectors need to consist of an optimal number of of elementary cells to balance the trac complexity and the communication overhead due to sector elementary cells to balance the traffic complexity and the communication overhead due to sector handovers. All the elements of this performance-based airspace model i.e., elementary cells, CNS handovers. All the elements of this performance-based airspace model i.e., elementary cells, CNS protection volume, occupancy grid and airspace sectors are visualised in Figure 3. The following protection volume, occupancy grid and airspace sectors are visualised in Figure 3. The following sections introduce the underpinning theoretical framework and mathematical models for all these sections introduce the underpinning theoretical framework and mathematical models for all these entities. However, the CNS protection volumes will be addressed up front as they form the theoretical entities. However, the CNS protection volumes will be addressed up front as they form the theoretical basis upon which the overall framework and the elementary cells were deﬁned. basis upon which the overall framework and the elementary cells were defined. Figure 3. Overall urban airspace concept. Figure 3. Overall urban airspace concept. This airspace model provides the flexibility for implementation in the planning stage of three This airspace model provides the ﬂexibility for implementation in the planning stage of three different operational timeframes: offline, pre-tactical online and tactical online. The proposed dierent operational timeframes: oine, pre-tactical online and tactical online. The proposed airspace airspace model can be applied to achieve different goals. In particular, the airspace capacity can be model can be applied to achieve dierent goals. In particular, the airspace capacity can be determined determined based on the performance of the CNS systems supporting the operations or, alternatively, based on the performance of the CNS systems supporting the operations or, alternatively, the model can be applied to evaluate the CNS performance requirements to be enforced given a desired airspace capacity target. This two-way approach is conceptually illustrated in Figure 4. Aerospace 2020, 7, x FOR PEER REVIEW 9 of 25 Aerospace 2020, 7, x FOR PEER REVIEW 9 of 25 the model can be applied to evaluate the CNS performance requirements to be enforced given a Aerospace 2020, 7, 154 9 of 25 the model can be applied to evaluate the CNS performance requirements to be enforced given a desired airspace capacity target. This two-way approach is conceptually illustrated in Figure 4. desired airspace capacity target. This two-way approach is conceptually illustrated in Figure 4. Airspace capacity CNS Performance Airspace capacity CNS Performance CNS performance- CNS performance- based airspace model based airspace model Communication overhead Communication overhead Sensor and communication Sensor and communication and computational and computational error characteristics error characteristics complexity complexity Forward model application – Determination of airspace capacity Forward model application – Determination of airspace capacity Inverse model application – Determination of required CNS performance Inverse model application – Determination of required CNS performance Figure 4. Dual approach of performance-based airspace model. Figure 4. Dual approach of performance-based airspace model. Figure 4. Dual approach of performance-based airspace model. 3.1. CNS Risk Protection Volume Generation 3.1. CNS Risk Protection Volume Generation 3.1. CNS Risk Protection Volume Generation This sub-section covers the generation of the CNS risk protection volume. The methodology as This sub-section covers the generation of the CNS risk protection volume. The methodology as This sub-section covers the generation of the CNS risk protection volume. The methodology as illustrated in Figure 5 is adopted from the approach introduced in [36], which is used for solving the illustrated in Figure 5 is adopted from the approach introduced in [36], which is used for solving the illustrated in Figure 5 is adopted from the approach introduced in [36], which is used for solving the autonomous Detect-and-Avoid (DAA) problem. autonomous Detect-and-Avoid (DAA) problem. autonomous Detect-and-Avoid (DAA) problem. Surveillance system Surveillance system error error Augmentation Augmentation system(s) Measurement to system(s) Measurement to Transaction time position domain Transaction time position domain Sensor error Sensor error Combine residual Measurement to Communication modelling and error Combine residual Measurement to Surveillance inflation Communication modelling and error errors position domain Surveillance inflation inflation bound errors position domain inflation bound Figure 5. CNS risk protection volume determination methodology. Figure 5. CNS risk protection volume determination methodology. Figure 5. CNS risk protection volume determination methodology. The navigation component of the volume is generated by developing a sensor measurement The navigation component of the volume is generated by developing a sensor measurement The navigation component of the volume is generated by developing a sensor measurement model that accounts for both nominal and off-nominal performance of the sensors utilized by the model that accounts for both nominal and off-nominal performance of the sensors utilized by the model that accounts for both nominal and o-nominal performance of the sensors utilized by the UAS. UAS. These typically include a GNSS receiver, inertial sensors, visual sensors and altimeter. UAS. These typically include a GNSS receiver, inertial sensors, visual sensors and altimeter. These typically include a GNSS receiver, inertial sensors, visual sensors and altimeter. The distributions characterizing these errors from these sensors can be, in theory, completely The distributions characterizing these errors from these sensors can be, in theory, completely The distributions characterizing these errors from these sensors can be, in theory, completely arbitrary. However, the long-standing practise for certifying systems in aviation has been to replace arbitrary. However, the long-standing practise for certifying systems in aviation has been to replace arbitrary. However, the long-standing practise for certifying systems in aviation has been to replace the true arbitrary distributions with a simpler distribution. This simpler distribution, called the the true arbitrary distributions with a simpler distribution. This simpler distribution, called the the true arbitrary distributions with a simpler distribution. This simpler distribution, called the overbounding distribution, must be such that it overestimates the true error by a small margin (Since overbounding distribution, must be such that it overestimates the true error by a small margin (Since overbounding distribution, must be such that it overestimates the true error by a small margin this is preferable to underestimation of the true error, which would constitute hazardously this is preferable to underestimation of the true error, which would constitute hazardously (Since this is preferable to underestimation of the true error, which would constitute hazardously misleading information). misleading information). misleading information). Aerospace 2020, 7, 154 10 of 25 Aerospace 2020, 7, x FOR PEER REVIEW 10 of 25 Navigation systems in aviation have conventionally been modelled as Gaussian distributions. Navigation systems in aviation have conventionally been modelled as Gaussian distributions. The use of Gaussian distributions has been the basis for the design of all algorithms in safety-of-life The use of Gaussian distributions has been the basis for the design of all algorithms in safety-of-life GNSS-based systems, including SBAS and GBAS. In GNSS for instance, the characterization of GNSS-based systems, including SBAS and GBAS. In GNSS for instance, the characterization of positioning errors as Gaussians are well documented in current standards [37]. positioning errors as Gaussians are well documented in current standards [37]. Replacing the expected distribution by the overbounding Gaussian distribution is necessary for Replacing the expected distribution by the overbounding Gaussian distribution is necessary for at least two reasons. First, an arbitrary distribution would require an amount of data too large to at least two reasons. First, an arbitrary distribution would require an amount of data too large to be be sent through typically low bandwidth channels between dierent participants in the encounter sent through typically low bandwidth channels between different participants in the encounter whereas a Gaussian is described by only two parameters (a mean value and standard deviation). whereas a Gaussian is described by only two parameters (a mean value and standard deviation). Second, combining dierent error sources is accomplished through convolution of distributions. Second, combining different error sources is accomplished through convolution of distributions. Computing the convolution of many error sources characterized by arbitrary distributions is likely Computing the convolution of many error sources characterized by arbitrary distributions is likely to be prohibitive. However, Gaussian convolutions are simple to perform. The residual error from to be prohibitive. However, Gaussian convolutions are simple to perform. The residual error from each error source is combined and the resulting Gaussian distribution is a conservative error budget each error source is combined and the resulting Gaussian distribution is a conservative error budget of the sensor measurements. The conservative sensor errors are then translated from the sensor of the sensor measurements. The conservative sensor errors are then translated from the sensor measurement domain to the aircraft position domain and the Gaussian sphere bounding this position measurement domain to the aircraft position domain and the Gaussian sphere bounding this position error distribution is the navigation protection volume.A similar process is adopted for the surveillance error distribution is the navigation protection volume.A similar process is adopted for the sensor(s) involved in the scenario. Worst-case sensor measurement error overbounds are translated to surveillance sensor(s) involved in the scenario. Worst-case sensor measurement error overbounds are the position domain and constitute the surveillance layer of the risk protection volume. translated to the position domain and constitute the surveillance layer of the risk protection volume. The communication layer is essentially a buer which accounts for the transaction time taken to The communication layer is essentially a buffer which accounts for the transaction time taken to conduct a manoeuvre to maintain separation or avoid a collision. As detailed in Figure 5, this transaction conduct a manoeuvre to maintain separation or avoid a collision. As detailed in Figure 5, this time depends on a combination of human and technical factors. The overall risk protection volume is transaction time depends on a combination of human and technical factors. The overall risk illustrated in Figure 6. protection volume is illustrated in Figure 6. Navigation uncertainty volume Total uncertainty volume Communication buffer Surveillance buffer Figure 6. UAS fundamental protection volume. Figure 6. UAS fundamental protection volume. In addition to CNS performance, the volume can also be inflated to account for operational and In addition to CNS performance, the volume can also be inﬂated to account for operational and environmental factors such as relative aircraft dynamics and wake turbulence [36]. environmental factors such as relative aircraft dynamics and wake turbulence [36]. 3.1.1. Navigation Uncertainty 3.1.1. Navigation Uncertainty GNSS range measurement accuracy is well described by aviation standards [38], with Gaussian GNSS range measurement accuracy is well described by aviation standards [38], with Gaussian overbound distributions accounting for various sources of measurement error. These include: overbound distributions accounting for various sources of measurement error. These include: Satellite Clock and Ephemeris Errors (σ ) Satellite Clock and Ephemeris Errors ( ) clk&eph clk&eph Ionospheric Residual Error (σ ) Ionospheric Residual Error ( ) iono iono Tropospheric Residual Error (σ ) Tropospheric Residual Error (tropo ) tropo Thermal noise and interferences (σ ) Thermal noise and interferences (noise ) noise Multipath (σ ) multipath Multipath ( ) multipath The User Equivalent Range Error (UERE) σ is the root sum square of these independent UERE error sources: 2 2 2 2 2 2 σ = σ + σ + σ + σ + σ (2) UERE clk&eph iono tropo noise multipath Aerospace 2020, 7, 154 11 of 25 The User Equivalent Range Error (UERE) is the root sum square of these independent error UERE sources: 2 2 2 2 2 2 = + + + + (2) UERE clk&eph iono tropo noise multipath In addition to the errors in range measurement, accuracy is also dependent on the geometry of the satellites relative to the receiver. The eect of geometry on the accuracy of the solution is parameterized through a number of scalar factors termed the Dilution of Precision (DOP). The Position Dilution of Precision (PDOP) is the most commonly used factor and is presented in greater detail in [39]. The navigation position error is expressed in Equation (3), showing the relation of and PDOP. UERE The maximum errors are chosen to capture the worst performance, representing the conservative case. The covariance matrix of the navigation position uncertainty is presented in Equation (4). 2 2 = PDOP. (3) NAV,GNSS UERE 2 3 6 0 0 7 6 7 xNAV 6 7 6 7 2 6 7 6 0 0 7 = (4) 6 7 yNAV GNSS 6 7 6 7 4 5 0 0 zNAV In general, the components of are well described by long-standing models in the literature UERE that conservatively bound the error [39]. The only exception is the multipath error component , multipath which is a site-dependent source of error that is dicult to model. However, recent work in the domain is focussing on developing conservative multipath models for UAM given a terrain map and satellite ephemerides [40]. Threshold values on the positioning uncertainty are standardized in conventional aviation operations through so-called ‘Alert Limits’ [41]. These specify the maximum allowable positioning uncertainty. Since no standards are currently deﬁned for UAM navigation systems, this value is conservatively bounded based on urban canyon width dimensions (since this represents the most stringent performance requirements for UAM use-cases). Therefore, the navigation performance is declared insucient for W, where W is the smallest canyon width being traversed for a GNSS given planned operation (e.g., this corresponds to W = 25 m for most Australian city centres). 3.1.2. Surveillance Inﬂations As anticipated in Table 2, the uncertainty in the position of a detected intruder aircraft is dependent on two factors. First, the error in localizing the aircraft, and second, the latency between such detection and the position estimate in a downstream separation assurance system. A conventional surveillance scenario is illustrated in Figure 7. When Primary Surveillance Radar (PSR) is utilized, the estimated range, azimuth and elevation to the intruder aircraft is used to compute its position in a cartesian reference frame. Assuming the reference geometry depicted in Figure 7, target state vector information is measured relative to the radar site in a spherical coordinate system in range, elevation and azimuth (r , , respectively). The measurements in each of the elements are prone to an SNR RDR RDR RDR dependent random range measurement error, which can be calculated as: = (5) rRDR 2B 2(SNR) where B is waveform bandwidth, c is the speed of light and signal to noise ratio (SNR). Radar angular measurements are commonly made using monopulse receive antennas that provide a dierence pattern characterized by a deep null on boresight. The dierence pattern formed by these beams may be used to measure the target angular position with a single signal transmission. The measurement accuracy in each angular coordinate is characterized by the RMS of the SNR dependent random angular measurement error, angular bias, and random measurement error. As with the range error, we assume angular error to be normally distributed: Aerospace 2020, 7, 154 12 of 25 2 2 2 2 = + + (6) RDR AN AF AB e e e 2 2 2 2 = + + (7) Aerospace 2020, 7, x FOR PEER REVIEW RDR AN AF AB 12 of 25 ADS-B message • Position • Velocity • State uncertainty • Intents PSR station ADS-B Ground station Figure 7. Conventional surveillance scenario. Figure 7. Conventional surveillance scenario. As with the range errors, the SNR dependent error dominates the radar angle error: As with the range errors, the SNR dependent error dominates the radar angle error: σ = AN (8) = p (8) AN k 2(SNR) k 2(SNR) where: ν is the radar beamwidth in the angular coordinates and k is the monopulse pattern where: v is the radar beamwidth in the angular coordinates and k is the monopulse pattern dierence difference slope. The tracking covariance is then: slope. The tracking covariance is then: σ 0 0 TRK 2 3 SPH 2 0 0 6 7 Q = [ 0 σ 0 ] 6 r 7 (9) TRK TRK ϵ TRK 6 7 6 7 SPH 2 6 7 Q = 6 0 0 7 (9) 0 0 σ 6 7 TRK TRK η 6 TRK 7 4 5 0 0 TRK Errors in these measurements propagate to the Cartesian position uncertainty. A more detailed treatm Err en ors t of in th these ese er measur rors an ements d the re pr sul opagate ting ellto ips the oidCartesian is presented position in [42 uncertainty ]. On the o.th Aer mor han e d detailed , when A trDS eatment -B is empl of these oyed err as ors thand e prithe mary resulting means ellipsoid of surveill isapr nce esented , the acc in ur [42 acy ]. o On f lo the caliz other ing thand, he intrud when er a ADS-B ircraft i is s clo employed sely relaas tedthe to G primary NSS perf means orman of ce.surveillance, In addition, th the e tra accuracy cking un ofce localizing rtainty ellthe ipsointr id m uder ust a air lso craft acco is un closely t for th re elated error to due GNSS to latency. performance. The totaIn l tra addition, cking un the cert tracking ainty σ uncert is ainty then:ellipsoid must tracking also account for the error due to latency. The total tracking uncertainty is then: 2 2 2 tracking σ = σ + σ (10) tracking NAV com 2 2 2 = + (10) where σ is the uncertainty in GNSS-estimated position, and σ is the uncertainty due to com NAV tracking NAV com communication latency. where is the uncertainty in GNSS-estimated position, and is the uncertainty due to NAV com The overall latency of a SA/CA procedure depends not only on the surveillance latency but also communication latency. on the time required to assess a collision threat and to generate and execute a resolution. In The overall latency of a SA/CA procedure depends not only on the surveillance latency but also on conventional operations, these tasks would have been performed by the pilot and ATC operator. In the time required to assess a collision threat and to generate and execute a resolution. In conventional the UTM context, however, the bulk of these tasks are assumed to be performed autonomously, with operations, these tasks would have been performed by the pilot and ATC operator. In the UTM context, manual intervention required only in emergency conditions. however, the bulk of these tasks are assumed to be performed autonomously, with manual intervention The total deconfliction time (t ) is the result of the sequential tasks illustrated in Figure 8, SA/CA required only in emergency conditions. which includes the time taken to process the tracks of intruder aircraft, assess potential collision The total deconﬂiction time (t ) is the result of the sequential tasks illustrated in Figure 8, SA/CA threats, generate avoidance trajectories and execute them. which includes the time taken to process the tracks of intruder aircraft, assess potential collision threats, generate avoidance trajectories and execute them. Aerospace 2020, 7, 154 13 of 25 Aerospace 2020, 7, x FOR PEER REVIEW 13 of 25 SA & CA System Intruder Ownship data Navigation observations system Track processing Surveillance system Intruder tracks Threat assessment and alerting Alerts Trajectory generation Guidance and control System Figure 8. Overall SA/CA process for autonomous system. Figure 8. Overall SA/CA process for autonomous system. 3.1.3. Total CNS Protection Volume in a Cooperative ADS-B Encounter 3.1.3. Total CNS Protection Volume in a Cooperative ADS-B Encounter The volume generation process is illustrated here for a scenario in which ADS-B is used to monitor The volume generation process is illustrated here for a scenario in which ADS-B is used to an intruder aircraft. The volume in this instance bounds the GNSS error. This is then inﬂated to monitor an intruder aircraft. The volume in this instance bounds the GNSS error. This is then inflated account for the eect of latency. The radius of the protection volume r are then computed as: PRV to account for the effect of latency. The radius of the protection volume 𝑟 are then computed as: 𝑃𝑅𝑉 r = k + v 2 + (t ) + 2 (11) t v v,t PRV SA/CA SA/CA 2 GNSS (11) r = k ∙ √σ + [(v σ ) + (t σ ) + 2σ ] PRV GNSS t SA/CA v v,t SA/CA where is the variance in the GNSS solution; v and t are the estimated velocity and SA/CA GNSS wh deconﬂiction ere σ time; is th e va and ria nce in ar th ee the GNS standar S sold uti deviations on; v and in tthe estimated are the velocity estimated and v deconﬂiction elocity and GNSS v t SA/CA SA/CA d time econ respectively fliction tim . e; The σterm anin d the σ squar a ere par th entheses e stand[.] ard is d the evtranslation iations in of the the esti ADS-B mated latency veloci frty omathe nd v t SA/CA temporal to the spatial domain. This term is obtained by applying the law of propagation of uncertainty deconfliction time respectively. The term in the square parentheses [.] is the translation of the ADS-B to compute the distance covered by the vehicle over the duration of the SA/CA interval. This requires latency from the temporal to the spatial domain. This term is obtained by applying the law of the uncertainty in the estimated velocity and in the estimated deconﬂiction time. This term essentially propagation of uncertainty to compute the distance covered by the vehicle over the duration of the bounds the distance covered by the vehicle over the duration of the latency. k is an inﬂation factor SA/CA interval. This requires the uncertainty in the estimated velocity and in the estimated specifying the number of standard deviations to expand the volume. This is chosen as k = 3 i.e., a 3 deconfliction time. This term essentially bounds the distance covered by the vehicle over the duration or 99.7% bound. It is a commonly followed practice to employ a recursive state estimator such as the of the latency. k is an inflation factor specifying the number of standard deviations to expand the Extended Kalman Filter (EKF) as the GNSS navigation ﬁlter, to process GNSS pseudoranges and obtain volume. This is chosen as k = 3 i.e., a 3σ or 99.7% bound. It is a commonly followed practice to a position solution. Both and are inherent outputs of the GNSS EKF. can be obtained v t employ a recursive state est GNSS imator such as the Extended Kalman Filter (EKF) as the GNSS navigation SA/CA by ﬁtting an overbounding Gaussian distribution to empirical data of system computational times. filter, to process GNSS pseudoranges and obtain a position solution. Both σ and σ are inherent GNSS v outputs of the GNSS EKF. σ can be obtained by fitting an overbounding Gaussian distribution SA/CA 3.2. Performance-Based Elementary Cell Dimensions to empirical data of system computational times. The proposed performance-based urban airspace model is intended to enhance not only safety 3 and .2. P ee rforma ciency nce but -Base also d Eease lementa UTM ry C operator ell Dimens ’s io interpr ns etability by using CNS performance as the main driver in airspace structure design, determining spacing, sector dimensions and capacity altogether, The proposed performance-based urban airspace model is intended to enhance not only safety while also supporting an intuitive visualisation. Therefore, as previously indicated, the elementary cell and efficiency but also ease UTM operator’s interpretability by using CNS performance as the main dimensions shall be parameterized as functions of the CNS performance for a given airspace region. driver in airspace structure design, determining spacing, sector dimensions and capacity altogether, The reference UAS equipage is assumed as follows: while also supporting an intuitive visualisation. Therefore, as previously indicated, the elementary cell dimensions shall be parameterized as functions of the CNS performance for a given airspace Communication–datalink; region. The reference UAS equipage is assumed as follows: Navigation–GNSS/INS [43]; Communication–datalink; Navigation–GNSS/INS [43]; Aerospace 2020, 7, 154 14 of 25 Surveillance–Automatic Dependent Surveillance Broadcast (ADS-B) Within the scope of this study, the elementary cell dimensions are deﬁned as a function of GNSS performance. GNSS is in fact not only the most widely employed source of absolute positioning for UAS operations but, as recognised by ICAO, GNSS is also the key element of all Communications, Navigation, and Surveillance/ATM (CNS/ATM) systems [44]. The elementary cell dimension is set as a 3 (99.7%) bound on GNSS positioning errors, accounting for errors in pseudorange and unfavourable satellite geometry over the intended operational airspace region. The uncertainty in horizontal and vertical positioning is related to the pseudorange error and DOP as: = 3 HDOP (12) x,y UERE = 3 VDOP (13) z UERE where HDOP is the Horizontal Dilution of Precision; VDOP is the Vertical Dilution of Precision; UERE is the User Equivalent Range Error (m) as described by Equation (2) The elementary cell dimensions are then set as: r = r = (14) x y x,y r = (15) z z In order to calculate the DOP factor, the line-of-sight vector between the receiver and the satellites in the Earth-Centered Earth-Fixed (ECEF) coordinate system needs to be determined. This is readily obtained from the ephemera of the GNSS satellites. The elevation ( ) and azimuth (' ) of each mu mu satellite in the NED coordinate system is then: = arcsinu (16) mu as,D as,N ' = arctan (17) mu as,E where u , u , u are the Down, East, and North components of the line-of-sight unit as,D as,E as,N vector respectively. The geometric matrix is a concatenation of each of the visible satellite elevations and azimuths: 2 3 i i i i i cos cos' cos sin' sin' 1 6 7 mu mu mu mu mu 6 7 6 7 6 i+1 i+1 i+1 i+1 i+1 7 6 7 cos cos' cos sin' sin' 1 6 7 mu mu mu mu mu 6 7 6 7 H = 6 7 (18) . . . . G 6 7 6 . . . . 7 6 7 . . . . 6 7 6 7 4 5 n n n n n cos cos' cos sin' sin' 1 mu mu mu mu mu The Local-navigation-frame Cofactor Matrix is then: 2 3 D D D D 6 11 12 13 14 7 6 7 6 7 6 7 6 7 D D D D 6 21 22 23 24 7 n n T n 6 7 P = H H = (19) 6 7 6 7 G G 6 7 D D D D 6 31 32 33 34 7 6 7 4 5 D D D D 41 42 43 44 The HDOP and VDOP factors are calculated as: HDOP = D + D (20) 11 22 VDOP = D (21) 33 Aerospace 2020, 7, 154 15 of 25 Aerospace 2020, 7, x FOR PEER REVIEW 15 of 25 3.3. Occupancy Grid To facilitate ATM/UTM system processing and enhanced capability to balance airspace demand and To facilitate ATM/UTM system processing and enhanced capability to balance airspace demand capacity dynamically, the occupancy grid concept is proposed to approximate the volumetric and capacity dynamically, the occupancy grid concept is proposed to approximate the volumetric demand. A grid of elementary cells is constructed for the entire urban region adopting the models in demand. A grid of elementary cells is constructed for the entire urban region adopting the models in Section 3.2. The CNS protection volumes around each active aircraft are then superimposed to the Section 3.2. The CNS protection volumes around each active aircraft are then superimposed to the grid grid as illustrated in Figure 9. An elementary cell is designated as occupied when it contains at least as illustrated in Figure 9. An elementary cell is designated as occupied when it contains at least one one point of the CNS protection volume bounding each aircraft. Each occupied cell, therefore, is a point of the CNS protection volume bounding each aircraft. Each occupied cell, therefore, is a result result of the demand placed on the airspace by virtue of the CNS performance. For visualisation of the demand placed on the airspace by virtue of the CNS performance. For visualisation purposes, purposes, occupied and free elementary cells are designated different colours (red and grey occupied and free elementary cells are designated dierent colours (red and grey respectively). respectively). Figure Figure9. 9.UAS UASpr pr otection otectionvolume volume (left (left ) ) and and relative relative o occupancy ccupancygrid grid( right (righ). t). The occupancy grid approach allows for a rapid assessment of the occupied space against the The occupancy grid approach allows for a rapid assessment of the occupied space against the overall free space. A count performed over the cells supports the calculation of the theoretical remaining overall free space. A count performed over the cells supports the calculation of the theoretical capacity to service the airspace with the given CNS infrastructure. remaining capacity to service the airspace with the given CNS infrastructure. 3.4. Airspace Sectorisation 3.4. Airspace Sectorisation A sectorization scheme based on the previously deﬁned elementary cells is proposed here, which A sectorization scheme based on the previously defined elementary cells is proposed here, optimally balances computational complexity and communication overheads. which optimally balances computational complexity and communication overheads. A UAS operation spanning multiple sectors requires a handover to be performed at each sector A UAS operation spanning multiple sectors requires a handover to be performed at each sector boundary. The number of cells allocated to a given sector is determined so as to optimize two boundary. The number of cells allocated to a given sector is determined so as to optimize two key key factors: the communication overhead required to perform all the required handovers for all factors: the communication overhead required to perform all the required handovers for all involved involved trac (which increase proportionally with the density of sectors), and the computational traffic (which increase proportionally with the density of sectors), and the computational complexity complexity associated with computing an avoidance volume and avoidance trajectory (which increases associated with computing an avoidance volume and avoidance trajectory (which increases proportionally with sector dimensions). proportionally with sector dimensions). The communication overhead is parameterized as a dimensionless factor which is a function The communication overhead is parameterized as a dimensionless factor which is a function of of the time required to perform a ‘handshake’ or to establish a link between the UAS and the UTM the time required to perform a ‘handshake’ or to establish a link between the UAS and the UTM operator, and the average time spent in each sector: operator, and the average time spent in each sector: v ̅ x,y v ̅v x,y z OH = t + (22) OH = t HS ( + ) (22) HS d d x,y z d d x,y z where t is the handshake time (s); v is the average horizontal velocity of the vehicle over a sector where t is the handshake time (s); v ̅ is the average horizontal velocity of the vehicle over a HS x,y HS x,y (m/s); v is the average vertical velocity of the vehicle (m/s); d is the average horizontal distance sector ( zm/s); v ̅ is the average vertical velocity of the vehiclex,y (m/s); d is the average horizontal z x,y covered in a sector (m); d is the average horizontal distance covered in a sector (m). distance covered in a secto z r (m); d is the average horizontal distance covered in a sector (m). Equation (22) essentially represents the proportion of time spent in establishing a handover communication link between sectors. The sector horizontal and vertical dimensions are by definition Aerospace 2020, 7, 154 16 of 25 Equation (22) essentially represents the proportion of time spent in establishing a handover communication link between sectors. The sector horizontal and vertical dimensions are by deﬁnition set as multiples of the elementary cell dimensions. Thus, d and d are calculated from the x,y z multiplication of elementary cell numbers (k , k ) and elementary cell dimension (r , r ): x,y z x,y z d = k r (23) x,y x,y x,y d = k r (24) z z z where: HDOP k = k z x,y VDOP The computational complexity is approximated as a factor that captures the proportion of time spent in generating avoidance volumes and avoidance trajectories as well as the maximum potential number of aircraft simultaneously in the sector, which as a worst case can be assumed equal to the number of elementary cells in each sector: (n 1)[t + t ] CD AVG ATG CP = = (25) CoP CoP where CD is the computational demand (s); CoP is the available computing capacity (s); t is time to AVG generate avoidance volume; t is time to generate trajectory; n is the number of elementary cells in ATG a sector. n = k k (26) x,y 4. Veriﬁcation Case Study The proposed urban airspace model is veriﬁed in two simulation case studies; the ﬁrst case study is a UAM air-taxi scenario under open-sky conditions while the second case study is a sUAS delivery operation below the skyline, addressing both unmanned vehicle types: passenger carrying and small non-passenger carrying. 4.1. UAM Air-Taxi Scenario An on-demand air-taxi scenario in the Melbourne Central Business District (CBD) is simulated for this case study. The operational sequence is as described in Section 2.1. The test platform utilised in this study is the Electric Vertical Takeo and Landing (eVTOL) concept UAM vehicle developed in the Airbus Pop.Up One project [45]. An urban building height database is utilized in generating the scenario. The operational airspace in the vertical plane is assumed to extend from 330 m AGL upwards. This is to allow a minimum clearance of 30 m above the tallest building height (298 m as recorded in the height database). The urban airspace is assumed to originate at a speciﬁed point initialised at the corner of Flinders Street and Spencer Street (37.821058, 144.955217, 330) since it is the southmost point of Melbourne CBD. In this case, a GNSS constellation of 30 satellites with orbital radii of 26561.750 km is assumed. The inclination angle was set at 55 degrees, and longitude and timing osets are neglected. To obtain realistic and representative values of HDOP and VDOP, a ground-vehicle measurement campaign was conducted in an open sky area, a highway route proceeding towards a residential area. Although these values dier from the ones that will be observed in a UAM mission, the measurements obtained from this campaign are deemed to be more conservative since the satellite visibility is generally lower at ground level. The measured HDOP and VDOP values are therefore conservative. The time-series of these two parameters during the experiment are shown in Figure 10. The maximum recorded HDOP (0.703) and VDOP (1.103) are used to deﬁne the elementary cell dimensions according to Equations (12) and (13). The lengths of r and r for the grid above the skyline are computed to be x,y z 11.6 m and 18.2 m respectively. Aerospace 2020, 7, 154 17 of 25 Aerospace 2020, 7, x FOR PEER REVIEW 17 of 25 (a) (b) (c) (d) Figure 10. Cell dimensions from collected open-sky data: (a) HDOP (b) VDOP (c) horizontal cell Figure 10. Cell dimensions from collected open-sky data: (a) HDOP (b) VDOP (c) horizontal cell dimension (d) vertical cell dimension. dimension (d) vertical cell dimension. These inputs from Table 3 are fed into Equations (22)–(25). The values for the velocity are based These inputs from Table 3 are fed into Equations (22)–(25). The values for the velocity are based on on the specifications for the Airbus Pop.Up One [45]. Due to the dense population and small space in the speciﬁcations for the Airbus Pop.Up One [45]. Due to the dense population and small space in the the urban area, the strict speed restriction is applied at the maximum of 10 m/s. The airspace urban area, the strict speed restriction is applied at the maximum of 10 m/s. The airspace sectorization sectorization results are presented in Figure 11. The inefficiency surface for the communication results are presented in Figure 11. The ineciency surface for the communication overhead (OH) and overhead (OH) and Computational Complexity (CP) are plotted. The optimum points lie on the Computational Complexity (CP) are plotted. The optimum points lie on the intersection between these intersection between these two surfaces: OH and CP (Figure 11a). The vertical sector dimension is two surfaces: OH and CP (Figure 11a). The vertical sector dimension is derived to avoid overlap with derived to avoid overlap with the Terminal Area airspace (TMA) which begins from 457 m (1500 ft) the Terminal Area airspace (TMA) which begins from 457 m (1500 ft) and the base height of the over and the base height of the over skyline starts at 330 m. Hence, the preferred k is calculated from skyline starts at 330 m. Hence, the preferred k is calculated from dividing the totalz heights of low-level (457330) dividing the total heights of low-level airspace above skyline by the vertical dimension of the airspace above skyline by the vertical dimension of the elementary cell (r = 18.2 m), k = 7. z z 18.2 (457−330) elementary cell (r = 18.2 𝑚 ), k = ≈ 7. The optimal point is marked in red. Based on that, The optimal pointz is marked in r zed. Based on that, the value for k is 8. The result from Figure 11b is 18.2 the number of elementary cells in one sector as calculated from Equation (26): 448 cells. the value for k is 8. The result from Figure 11b is the number of elementary cells in one sector as calculated from Equation (26): 448 cells. Table 3. Airspace sectorisation input parameters. Hardware Software Velocity (m/s) 10 Time handshake (s) 0.3 Vertical rate (m/s) 4 Time avg (s) 112 Length (m) 4.4 Time moto (s) 12e5 Height (m) 0.85 Computational power 2.1e9 Width (m) 5 (a) (b) Figure 11. Airspace sectorisation: (a) overhead and computational complexity (b) optimal intersection point. Table 3. Airspace sectorisation input parameters. Hardware Software Velocity (m/s) 10 Time handshake (s) 0.3 Vertical rate (m/s) 4 Time avg (s) 112 Length (m) 4.4 Time moto (s) 12e5 Aerospace 2020, 7, x FOR PEER REVIEW 17 of 25 (a) (b) (c) (d) Figure 10. Cell dimensions from collected open-sky data: (a) HDOP (b) VDOP (c) horizontal cell dimension (d) vertical cell dimension. These inputs from Table 3 are fed into Equations (22)–(25). The values for the velocity are based on the specifications for the Airbus Pop.Up One [45]. Due to the dense population and small space in the urban area, the strict speed restriction is applied at the maximum of 10 m/s. The airspace sectorization results are presented in Figure 11. The inefficiency surface for the communication overhead (OH) and Computational Complexity (CP) are plotted. The optimum points lie on the intersection between these two surfaces: OH and CP (Figure 11a). The vertical sector dimension is derived to avoid overlap with the Terminal Area airspace (TMA) which begins from 457 m (1500 ft) and the base height of the over skyline starts at 330 m. Hence, the preferred k is calculated from dividing the total heights of low-level airspace above skyline by the vertical dimension of the ( ) 457−330 elementary cell (r = 18.2 𝑚 ), k = ≈ 7. The optimal point is marked in red. Based on that, z z 18.2 Aerospace 2020, 7, 154 18 of 25 the value for k is 8. The result from Figure 11b is the number of elementary cells in one sector as calculated from Equation (26): 448 cells. Aerospace 2020, 7, x FOR PEER REVIEW 18 of 25 (a) (b) Height (m) 0.85 Computational power 2.1e9 Figure 11. Airspace sectorisation: (a) overhead and computational complexity (b) optimal intersection Figure 11. Airspace sectorisation: (a) overhead and computational complexity (b) optimal Width (m) 5 point. intersection point. The resulting sectorisation of the skyline above the city centre is illustrated in Figure 12. The resulting sectorisation Table of 3the . Aiskyline rspace sec above torisatthe ion city inpucentr t parae m is etillustrated ers. in Figure 12. Hardware Software Velocity (m/s) 10 Time handshake (s) 0.3 Vertical rate (m/s) 4 Time avg (s) 112 Length (m) 4.4 Time moto (s) 12e5 (a) (b) Figure 12. Above skyline airspace sectorisation: (a) elementary cells (b) airspace sectors. Figure 12. Above skyline airspace sectorisation: (a) elementary cells (b) airspace sectors. In the case of an intruder aircraft, the required inputs for CNS protection volume generation are In the case of an intruder aircraft, the required inputs for CNS protection volume generation are shown in Table 4. The value of PDOP was set in a similar manner to the HDOP and VDOP from the shown in Table 4. The value of PDOP was set in a similar manner to the HDOP and VDOP from the measurement campaign. measurement campaign. A maximum allowable delay of 2 s for autonomous conﬂict detection and resolution is assumed. The maximum allowable tracking Table 4. I delay nputs accor for CN ding S prto otec the tion MOPS volum for e gen ADS-B eration and . TIS-B (RTCA-DO 260B) is 1.5 s. We assume a nominal tracking delay is assumed of 1 s. The total t is then 3 s with SA/CA Parameter Value Parameter Value the assumed value of at 0.5 s. The protection volume and occupancy grid are illustrated in SA/CA PDOP 1.3 σ 5.5 m UERE Figure 13. v 10 m/s 3 s SA/CA σ 0.058 m/s t 0.5 s v SA/CA σ 0.005 m r 26.2 m vd PRV A maximum allowable delay of 2 s for autonomous conflict detection and resolution is assumed. The maximum allowable tracking delay according to the MOPS for ADS-B and TIS-B (RTCA-DO 260B) is 1.5 s. We assume a nominal tracking delay is assumed of 1 s. The total t is then 3 s with SA/CA the assumed value of σ at 0.5 s. The protection volume and occupancy grid are illustrated in SA/CA Figure 13. Aerospace 2020, 7, 154 19 of 25 Table 4. Inputs for CNS protection volume generation. Parameter Value Parameter Value PDOP 1.3 5.5 m UERE v 10 m/s t 3 s SA/CA 0.058 m/s 0.5 s v t SA/CA 0.005 m r 26.2 m vd PRV Aerospace 2020, 7, x FOR PEER REVIEW 19 of 25 Aerospace 2020, 7, x FOR PEER REVIEW 19 of 25 Figure 13. The intruder’s CNS protection volume and associated occupancy grid above skyline. Figure 13. The intruder’s CNS protection volume and associated occupancy grid above skyline. Figure 13. The intruder ’s CNS protection volume and associated occupancy grid above skyline. 4.2. UAS Delivery Scenario 4.2. UAS Delivery Scenario 4.2. UAS Delivery Scenario This scenario investigates a sUAS delivery operation below the city skyline. To set the This scenario This scena invest rio in igates vestiga ates sUAS a sUA delivery S delioperati very ope on ra below tion bel the ow city the skyline. city skyl To inset e. To the set elementary the elementary cell dimensions for this low altitude region of airspace, the portion of the ground-vehicle elementary cell dimensions for this low altitude region of airspace, the portion of the ground-vehicle cell dimensions for this low altitude region of airspace, the portion of the ground-vehicle measurement measurement campaign within the city centre was extracted and the worst-case observed DOP measurement campaign within the city centre was extracted and the worst-case observed DOP campaign within the city centre was extracted and the worst-case observed DOP parameters were parameters were utilized. The DOP time series and correspondingly computed cell dimensions are parameters were utilized. The DOP time series and correspondingly computed cell dimensions are utilized. The DOP time series and correspondingly computed cell dimensions are shown in Figure 14. shown in Figure 14. As expected, the DOP is degraded (by nearly a factor of 2) from the open-sky shown in Figure 14. As expected, the DOP is degraded (by nearly a factor of 2) from the open-sky As expected, the DOP is degraded (by nearly a factor of 2) from the open-sky scenario. scenario. scenario. (a) (b) (a) (b) (c) (d) (c) (d) Figure 14. Cell dimensions for below skyline scenario: (a) HDOP (b) VDOP (c) horizontal cell Figure 14. Cell dimensions for below skyline scenario: (a) HDOP (b) VDOP (c) horizontal cell Figure 14. Cell dimensions for below skyline scenario: (a) HDOP (b) VDOP (c) horizontal cell dimension dimension (d) vertical cell dimension. dimension (d) vertical cell dimension. (d) vertical cell dimension. The worst HDOP and VDOP for this scenario are 1.337 and 2.067 respectively. The fundamental The worst HDOP and VDOP for this scenario are 1.337 and 2.067 respectively. The fundamental cell dimension, therefore, equates to r = 18.7 m and r = 34.1 m. Since airspace within urban x,y 𝑧 cell dimension, therefore, equates to r = 18.7 m and r = 34.1 m. Since airspace within urban x,y 𝑧 canyons is limited by buildings, the airspace generation procedure is different from the open-sky canyons is limited by buildings, the airspace generation procedure is different from the open-sky scenario. The algorithm pseudocode is presented in Table 5. scenario. The algorithm pseudocode is presented in Table 5. Aerospace 2020, 7, 154 20 of 25 The worst HDOP and VDOP for this scenario are 1.337 and 2.067 respectively. The fundamental cell dimension, therefore, equates to r = 18.7 m and r = 34.1 m. Since airspace within urban x,y z canyons is limited by buildings, the airspace generation procedure is dierent from the open-sky scenario. The algorithm pseudocode is presented in Table 5. Table 5. Pseudocode for within urban canyon airspace elementary cell generation. 1 Input: Width of Intersection W, Fundamental Cell Dimension r x, 2 Repeat 3 row = ﬂoor (W/r ) 4 If row = odd then 5 start = W/2 6 for I = 1:row do 7 i = i * r 8 Generateairspace (start, start + i, start i) 9 else 10 start1 = W/2 + r 11 start2 = W/2 r 12 for i = 1:row do 13 i = i * r 14 Generateairspace(start1,start1 + i,) 15 Generateairspace(start2,start2 + i,) 16 End 17 return Generateairspace 18 end procedure If the urban canyon width can only ﬁt a single row of cells, the width of the elementary cell is stretched to be equal to the canyon width. In this case study, the elementary cell width is 18.7 m, while the canon width is 25 m. Hence, the width of the cell is stretched to 25 m. For determining the optimal sector size, the same methodology presented in Section 3.4 is applied. The urban canyon airspace extends from 30m AGL to 330m AGL. It is desirable to have at least 2 altitude 300/34.1 levels. Accordingly, the vertical size of the elemental cells is set at k = 4. The optimal value for k is rounded down to 4 (Figure 15). The total number of elementary cells in one sector is 64 cells. x,y The resulting optimal sectorisation below the skyline is depicted in Figure 16. Similar to the ﬁrst case study, the required inputs for the generation of UAS protection volume are shown in Table 6. Currently, there is a lack of established standards for the sUAS speed limit in urban areas. The top speed of most sUAS ranges from 5 to 10 m/s [46]. However, a speed limit of 2 m/s is assumed, which is comparable to the top speed of several micro UAS [46], to allow sucient time for the system to respond in case of emergencies. The PDOP in this case is 2.6. The rest of the parameters remains the same as in Table 4. The protection volume and occupancy grid are illustrated in Figure 17. The resulting spherical CNS protection volume of an intruder and host UAS exceeds the below-skyline sector width (urban canyon width W = 25 m). Consistently with our proposed two-way approach (Figure 4), this is a clear indication that the observed CNS performance is insucient for safe navigation along the street and stricter CNS performance requirements are to be fulﬁlled to access such urban airspace portion. In future work, tailored directional inﬂations (in x, y, z) will be investigated, driven both by relative dynamics, weather, wake turbulence and other relevant factors. Aerospace 2020, 7, x FOR PEER REVIEW 20 of 25 Table 5. Pseudocode for within urban canyon airspace elementary cell generation. 1 Input: Width of Intersection W, Fundamental Cell Dimension rx, 2 Repeat 3 row = floor(W/rx) 4 If row = odd then 5 start = W/2 6 for I = 1:row do 7 i = i * rx 8 Generateairspace(start,start + i,start − i) 9 else 10 start1 = W/2 + rx 11 start2 = W/2 − rx 12 for i = 1:row do 13 i = i * rx 14 Generateairspace(start1,start1 + i,) 15 Generateairspace(start2,start2 + i,) 16 End 17 return Generateairspace 18 end procedure If the urban canyon width can only fit a single row of cells, the width of the elementary cell is stretched to be equal to the canyon width. In this case study, the elementary cell width is 18.7 m, while the canon width is 25 m. Hence, the width of the cell is stretched to 25 m. For determining the optimal sector size, the same methodology presented in Section 3.4 is applied. The urban canyon airspace extends from 30m AGL to 330m AGL. It is desirable to have at least 2 300/34.1 altitude levels. Accordingly, the vertical size of the elemental cells is set at k = ≈ 4. The Aerospace 2020, 7, 154 21 of 25 optimal value for k is rounded down to 4 (Figure 15). The total number of elementary cells in one x,y sector is 64 cells. Aerospace 2020, 7, x FOR PEER REVIEW 21 of 25 Table 6. Inputs for CNS protection volume generation. Parameter Value Parameter Value PDOP 2.6 5.5 m UERE v t 2 m/s 3 s UAS SA/CA (a) σ 0.058 m/s t 0.5 s (b) v SA/CA σ r 0.005 m 24.8 m vd PRV Figure 15. Airspace sectorisation study: (a) overall plot of overhead and computational complexity Figure 15. Airspace sectorisation study: (a) overall plot of overhead and computational complexity (b) optimal intersection point. (b) optimal intersection point. The resulting optimal sectorisation below the skyline is depicted in Figure 16. Similar to the first case study, the required inputs for the generation of UAS protection volume are shown in Table 6. Currently, there is a lack of established standards for the sUAS speed limit in urban areas. The top speed of most sUAS ranges from 5 to 10 m/s [46]. However, a speed limit of 2 m/s is assumed, which is comparable to the top speed of several micro UAS [46], to allow sufficient time for the system to respond in case of emergencies. The PDOP in this case is 2.6. The rest of the parameters remains the same as in Table 4. (a) (b) Figure 16. Airspace sectorisation below skyline: (a) elementary cells (b) airspace sectors. Figure 16. Airspace sectorisation below skyline: (a) elementary cells (b) airspace sectors. The protection volume and occupancy grid are illustrated in Figure 17. The resulting spherical CNS protection volume of an intruder and host UAS exceeds the below-skyline sector width (urban canyon width W = 25 m). Consistently with our proposed two-way approach (Figure 4), this is a clear indication that the observed CNS performance is insufficient for safe navigation along the street and stricter CNS performance requirements are to be fulfilled to access such urban airspace portion. In future work, tailored directional inflations (in x, y, z) will be investigated, driven both by relative dynamics, weather, wake turbulence and other relevant factors. Aerospace 2020, 7, 154 22 of 25 Table 6. Inputs for CNS protection volume generation. Parameter Value Parameter Value PDOP 2.6 5.5 m UERE v 2 m/s t 3 s UAS SA/CA 0.058 m/s 0.5 s v t SA/CA 0.005 m r 24.8 m vd PRV Aerospace 2020, 7, x FOR PEER REVIEW 22 of 25 (a) (b) Figure Figure 1 17. 7. UUAS AS PRV PRV vo volume lume wwith ith ococcupancy cupancy grid grid in in belbelow ow skyskyline line scen scenario ario (a) a (a x)on axonometric; ometric; (b) ( la bt )er lateral al view view . . 5. Conclusions and Future Research 5. Conclusions and Future Research The envisioned spread of Urban Air Mobility (UAM) and low altitude UAS services prompts the The envisioned spread of Urban Air Mobility (UAM) and low altitude UAS services prompts need to introduce new airspace structures and classiﬁcations, particularly in dense urban and suburban the need to introduce new airspace structures and classifications, particularly in dense urban and areas. This paper presented a new approach to the design of urban airspace based on the combined suburban areas. This paper presented a new approach to the design of urban airspace based on the performance of avionics systems and supporting trac management infrastructure. In particular, the combined performance of avionics systems and supporting traffic management infrastructure. In airspace was discretized into fundamental volume elements (elementary cells) as a function of the particular, the airspace was discretized into fundamental volume elements (elementary cells) as a Communication, Navigation and Surveillance (CNS) systems performance and a dedicated study was function of the Communication, Navigation and Surveillance (CNS) systems performance and a presented exploring the relevance of Global Navigation Satellite System (GNSS) performance to both dedicated study was presented exploring the relevance of Global Navigation Satellite System (GNSS) aircraft navigation and Automatic Dependent Surveillance-Broadcast (ADS-B) cooperative surveillance. performance to both aircraft navigation and Automatic Dependent Surveillance-Broadcast (ADS-B) Additionally, the paper introduced a methodology to determine the optimal sector dimension as a cooperative surveillance. Additionally, the paper introduced a methodology to determine the combination of multiple elementary cells, considering on one hand the trac complexity, and on the optimal sector dimension as a combination of multiple elementary cells, considering on one hand the other hand the overheads due to sector handovers. The proposed airspace structuring methodology traffic complexity, and on the other hand the overheads due to sector handovers. The proposed was then veriﬁed in two simulation case studies that utilized representative GNSS measurements airspace structuring methodology was then verified in two simulation case studies that utilized gathered from an urban measurement campaign. The case studies showed that the performance-based representative GNSS measurements gathered from an urban measurement campaign. The case discretization is numerically feasible and, based on the conservative assumptions and the principles studies showed that the performance-based discretization is numerically feasible and, based on the conservative assumptions and the principles adopted, it promotes operational safety while at the same time maximising the efficiency of airspace resource exploitation, particularly by efficiently accommodating unmanned aircraft with diverse avionics equipment. It is expected that the progressive transition from an RNP-based formulation to a full CNS performance-based approach, including other parameters from the communication system (e.g., signal to noise ratio, bit error rate, etc.), will offer significant benefits in all planning and operational UTM timeframes from strategic offline to tactical online scenarios. In this perspective, the proposed airspace sectorisation concept will contribute significantly to enhance decision support for demand-capacity balancing and dynamic airspace management in low-level ATM operations. Future work will also investigate the synergies between the proposed concept and adaptive/cognitive forms of Human–Machine Interfaces and Interactions (HMI ) to enhance the cooperation between human operators and the Aerospace 2020, 7, 154 23 of 25 adopted, it promotes operational safety while at the same time maximising the eciency of airspace resource exploitation, particularly by eciently accommodating unmanned aircraft with diverse avionics equipment. It is expected that the progressive transition from an RNP-based formulation to a full CNS performance-based approach, including other parameters from the communication system (e.g., signal to noise ratio, bit error rate, etc.), will oer signiﬁcant beneﬁts in all planning and operational UTM timeframes from strategic oine to tactical online scenarios. In this perspective, the proposed airspace sectorisation concept will contribute signiﬁcantly to enhance decision support for demand-capacity balancing and dynamic airspace management in low-level ATM operations. Future work will also investigate the synergies between the proposed concept and adaptive/cognitive forms of Human–Machine Interfaces and Interactions (HMI ) to enhance the cooperation between human operators and the increasingly automated (and trusted autonomous) system functions required in the UTM operational context. Application of directional volume inﬂations (as opposed to the conservative omnidirectional inﬂations adopted in this paper) to account for factors such as aircraft relative dynamics, weather eects, wake turbulence and downwash will also be investigated. Author Contributions: Conceptualization, N.P., S.B., A.G. and R.S.; Methodology, N.P, S.B, A.G., A.S. and Y.X.; Software, N.P., S.B. and A.G.; Formal Analysis, N.P. and S.B.; Investigation, R.S., N.P., S.B., A.S. and Y.X.; Resources, S.B. and A.G.; Writing-Original Draft Preparation, N.P, S.B, A.S. and Y.X.; Writing-Review & Editing, A.G., R.S. and T.K.; Supervision, R.S. and A.G.; Project Administration, R.S. and T.K. All authors have read and agreed to the published version of the manuscript. Funding: The authors wish to thank and acknowledge Thales Australia—Airspace Mobility Solutions for supporting this work under the collaborative research project RE-02544-0200315666. The authors also wish to thank and acknowledge the Australian Cooperative Research Centre on Intelligent Mobility and Vehicle Evolutions (iMOVE CRC) for supporting this work as part of the iMOVE industry-based student projects initiative. Conﬂicts of Interest: The authors declare no conﬂict of interest. References 1. Barr, L.C.; Newman, R.; Ancel, E.; Belcastro, C.M.; Foster, J.V.; Evans, J.; Klyde, D.H. Preliminary risk assessment for small unmanned aircraft systems. In Proceedings of the 17th AIAA Aviation Technology, Integration, and Operations Conference, Denver, CO, USA, 5–9 June 2017; p. 3272. 2. Prevot, T.; Homola, J.; Mercer, J. From rural to urban environments: Human/systems simulation research for low altitude UAS Trac Management (UTM). In Proceedings of the 16th AIAA Aviation Technology, Integration, and Operations Conference, Washington, DC, USA, 13–17 June 2016; p. 3291. 3. Cook, S.P.; Lacher, A.R.; Maroney, D.R.; Zeitlin, A.D. UAS sense and avoid development—The challenges of technology, standards, and certiﬁcation. In 50th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition; AIAA: Nashville, TN, USA, 2012. 4. Atkins, E.; Khalsa, A.; Groden, M. Commercial Low-Altitude UAS Operations in Population Centers. In Proceedings of the 9th AIAA Aviation Technology, Integration, and Operations Conference (ATIO) and Aircraft Noise and Emissions Reduction Symposium (ANERS), Hilton Head, SC, USA, 21–23 September 2009; p. 7070. 5. Heutger, M.; Kückelhaus, M. Unmanned Aerial Vehicles in Logistics—A DHL Perspective on Implications and Use Cases for the Logistics Industry; DHL Customer Solutions & Innovation: Troisdorf, Germany, 2014. 6. Kopardekar, P.H. Unmanned Aerial System (UAS) Trac Management (UTM): Enabling Low-Altitude Airspace and UAS Operations; TM—2014–218299; Ames Research Center: Moett Field, CA, USA, 2014. 7. FAA. 14 CFR: Aeronautics and Space. In Proceedings of Electronic Code of Federal Regulations; GPO: Washington, DC, USA, 2017. 8. Kopardekar, P.; Rios, J.; Prevot, T.; Johnson, M.; Jung, J.; Robinson, J.E. Unmanned Aircraft System Trac Management (UTM) Concept of Operations. In Proceedings of the 16th AIAA Aviation Technology, Integration, and Operations (ATIO), Washington, DC, USA, 13–17 June 2016. 9. Homola, J.; Martin, L.; Cencetti, M.; Aweiss, A. UAS Trac Management (UTM) Technical Capability Level 3 (TCL3) Flight Demonstration: Concept Tests and Results. In Proceedings of the 2019 IEEE/AIAA 38th Digital Avionics Systems Conference (DASC), San Diego, CA, USA, 8–12 September 2019; pp. 1–8. Aerospace 2020, 7, 154 24 of 25 10. Homola, J.; Dao, Q.; Martin, L.; Mercer, J.; Mohlenbrink, C.; Claudatos, L. Technical capability level 2 unmanned aircraft system trac management (UTM) ﬂight demonstration: Description and analysis. In Proceedings of the 2017 IEEE/AIAA 36th Digital Avionics Systems Conference (DASC), St. Petersburg, FL, USA, 17–21 September 2017; pp. 1–10. 11. Prevot, T.; Rios, J.; Kopardekar, P.; Robinson, J.E.; Johnson, M.; Jung, J. UAS trac management (UTM) concept of operations to safely enable low altitude ﬂight operations. In Proceedings of the 16th AIAA Aviation Technology, Integration, and Operations Conference, Washington, DC, USA, 13–17 June 2016; p. 3292. 12. Hackenberg, D. NASA Advanced Air Mobility (AAM) Mission [PowerPoint presentation]. In Proceedings of the AAM Ecosystem Working Groups (AEWG), National Aeronautical and Space Administration (NASA) AAM Mission Oce, Virtual Conference, 28 May 2020. 13. Jang, D.-S.; Ippolito, C.A.; Sankararaman, S.; Stepanyan, V. Concepts of airspace structures and system analysis for UAS trac ﬂows for urban areas. In Proceedings of the AIAA Information Systems—AIAA Infotech @ Aerospace, Grapevine, TX, USA, 9–13 January 2017; p. 0449. 14. Kistan, T.; Gardi, A.; Sabatini, R.; Ramasamy, S.; Batuwangala, E. An evolutionary outlook of air trac ﬂow management techniques. Progress Aerosp. Sci. 2017, 88, 15–42. [CrossRef] 15. Hoekstra, J.M.; Maas, J.; Tra, M.; Sunil, E. How do layered airspace design parameters aect airspace capacity and safety? In Proceedings of the Proceedings of the 7th International Conference on Research in Air Transportation, Philadelphia, PA, USA, 20–24 June 2016. 16. Majumdar, A.; Ochieng, W.; Polak, J. Estimation of European airspace capacity from a model of controller workload. J. Navig. 2002, 55, 381–403. [CrossRef] 17. Majumdar, A.; Ochieng, W.Y.; Bentham, J.; Richards, M. En-route sector capacity estimation methodologies: An international survey. J. Air Transp. Manag. 2005, 11, 375–387. [CrossRef] 18. Sunil, E.; Ellerbroek, J.; Hoekstra, J.; Vidosavljevic, A.; Arntzen, M.; Bussink, F.; Nieuwenhuisen, D. Analysis of airspace structure and capacity for decentralized separation using fast-time simulations. J. Guid. Control. Dyn. 2016, 40, 38–51. [CrossRef] 19. Sunil, E.; Hoekstra, J.; Ellerbroek, J.; Bussink, F.; Vidosavljevic, A.; Delahaye, D.; Aalmoes, R. The inﬂuence of trac structure on airspace capacity. In Proceedings of the 7th International Conference on Research in Air Transportation, Philadelphia, PA, USA, 20–24 June 2016. 20. Lim, Y.; Premlal, N.; Gardi, A.; Sabatini, R. Eulerian Optimal Control Formulation for Dynamic Morphing of Airspace Sectors. In Proceedings of the 31st Congress of the International Council of the Aeronautical Sciences, Belo Horizonte, Brazil, 9–14 September 2018. 21. Pongsakornsathien, N.; Bijjahalli, S.; Gardi, A.; Sabatini, R.; Kistan, T. A Novel Navigation Performance-based Airspace Model for Urban Air Mobility. In Proceedings of the 39th Digital Avionics System Conference (DASC), Virtual Conference, 11–16 October 2020. 22. Pongsakornsathien, N.; Gardi, A.; Sabatini, R.; Kistan, T.; Ezer, N. Human-Machine Interactions in Very-Low-Level UAS Operations and Trac Management. In Proceedings of the IEEE/AIAA 39th Digital Avionics Systems Conference (DASC 2020), Virtual Conference, 11–16 October 2020. 23. Sunil, E.; Hoekstra, J.; Ellerbroek, J.; Bussink, F.; Nieuwenhuisen, D.; Vidosavljevic, A.; Kern, S. Metropolis: Relating airspace structure and capacity for extreme trac densities. In Proceedings of the ATM Seminar 2015, 11th USA/EUROPE Air Trac Management R&D Seminar, Lisbon, Portugal, 23–26 June 2015. 24. Schneider, O.; Kern, S.; Knabe, F.; Gerdes, I.; Delahaye, D.; Vidosavljevic, A.; Leeuwen, P.; Nieuwenhuisen, D.; Sunil, E.; Hoekstra, J.; et al. METROPOLIS Concept Design Report; Delft University of Technology: Delft, The Netherlands, 2014. 25. Sacharny, D.; Henderson, T.C. A Lane-Based Approach for Large-Scale Strategic Conﬂict Management for UAS Service Suppliers. In Proceedings of the 2019 International Conference on Unmanned Aircraft Systems (ICUAS), Atlanta, GA, USA, 11–14 June 2019; pp. 937–945. 26. Sedov, L.; Polishchuk, V. Centralized and distributed UTM in layered airspace. In Proceedings of the 2018 8th International Conference on Research in Air Transportation (ICRAT), Catalonia, Spain, 26–29 June 2018. 27. Duchamp, V.; Sedov, L.; Polishchuk, V. Density-adapting layers towards PBN for UTM. In Proceedings of the 13th USA/Europe Air Trac Management Research and Development Seminar 2019, Vienna, Austria, 17–21 June 2019. Aerospace 2020, 7, 154 25 of 25 28. Amazon. Revising the Airspace Model for the Safe Integration of Small Unmanned Aircraft Systems. 2015. Available online: https://images-na.ssl-images-amazon.com/images/G/01/112715/download/Amazon_ Revising_the_Airspace_Model_for_the_Safe_Integration_of_sUAS.pdf (accessed on 15 September 2020). 29. Kopaderkar, P. Safely enabling UAS operations in low-altitude airspace. In Proceedings of the IEEE/AIAA 35th Digital Avionics Systems Conference (DASC), Sacramento, CA, USA, 25–29 September 2016; p. 33. 30. Mendonca, N.; Metcalfe, M.; Wiggins, S.; Grin, C.; DeCarme, D. AAM Ecosystem Working Groups (AEWG): Urban Air Mobility (UAM) Concept of Operations (ConOps) Airspace Breakout [PowerPoint presentation]. In Proceedings of the AAM Airspace Working Group Kicko, National Aeronautical and Space Administration (NASA) AAM Mission and Deloitte, Virtual Conference, 28 May 2020. 31. Tan, Q.; Wang, Z.; Ong, Y.-S.; Low, K.H. Evolutionary Optimization-based Mission Planning for UAS Trac Management (UTM). In Proceedings of the 2019 International Conference on Unmanned Aircraft Systems (ICUAS), Atlanta, GA, USA, 11–14 June 2019. 32. Nneji, V.C.; Stimpson, A.; Cummings, M.; Goodrich, K.H. Exploring concepts of operations for on-demand passenger air transportation. In Proceedings of the 17th AIAA Aviation Technology, Integration, and Operations Conference, Denver, CO, USA, 5–9 June 2017; p. 3085. 33. Samavati, S. Lilium Aviation Raises ¿10 million from Atomico for Vertical Take-O Electric Plane. Available online: https://tech.eu/brief/lilium-jet-secures-10-million/ (accessed on 27 November 2016). 34. Thipphavong, D.P.; Apaza, R.; Barmore, B.; Battiste, V.; Burian, B.; Dao, Q.; Feary, M.; Go, S.; Goodrich, K.H.; Homola, J. Urban air mobility airspace integration concepts and considerations. In Proceedings of the 2018 Aviation Technology, Integration, and Operations Conference, Atlanta, GA, USA, 25–29 June 2018; p. 3676. 35. Doc, I. 9613: Performance based Navigation (PBN) Manual; International Civil Aviation Organisation: Montréal, QC, Canada, 2008. 36. Ramasamy, S.; Sabatini, R.; Gardi, A. A Uniﬁed Analytical Framework for Aircraft Separation Assurance and UAS Sense-and-Avoid. J. Intell. Rob. Syst. Theor. Appl. 2018, 91, 735–754. [CrossRef] 37. RTCA. 229D—Minimum Operational Performance Standards for GPS WAAS Airborne Equipment; RTCA Inc.: Washington, DC, USA, 2006. 38. RTCA. DO-229D Minimum Operational Performance Standards for Global Positioning System/Wide Area Augmentation System Airborne Equipment; RTCA Inc.: Washington, DC, USA, 2006. 39. Sabatini, R.; Moore, T.; Ramasamy, S. Global navigation satellite systems performance analysis and augmentation strategies in aviation. Prog. Aerosp. Sci. 2017, 95, 45–98. [CrossRef] 40. Bijjahalli, S.; Sabatini, R.; Gardi, A. GNSS Performance Modelling and Augmentation for Urban Air Mobility. Sensors 2019, 19, 4209. [CrossRef] [PubMed] 41. RTCA. RTCA DO229E—Minimum Operational Performance Standards for Global Positioning System/Satellite-Based Augmentation System Airborne Equipment; RTCA, Inc.: Washington, DC, USA, 2016. 42. Hilton, S.; Cairola, F.; Gardi, A.; Sabatini, R.; Pongsakornsathien, N.; Ezer, N. Uncertainty Quantiﬁcation for Space Situational Awareness and Trac Management. Sensors 2019, 19, 4361. [CrossRef] [PubMed] 43. Bijjahalli, S.; Sabatini, R.; Gardi, A. Advances in intelligent and autonomous navigation systems for small UAS. Prog. Aerosp. Sci. 2020, 115, 100617. [CrossRef] 44. ICAO. Global Navigation Satellite System (GNSS) Manual. In Doc 9849 AN/457; International Civil Avaition Organization: Montreal, QC, Canada, 2005. 45. Airbus; Italdesign. Italdesign and Airbus unveil Pop.Up: A trailblazing modular ground and air passenger concept vehicle system. In Proceedings of the 87th Geneva International Motor Show, Geneva, Switzerland, 9–19 March 2017. 46. la Cour-Harbo, A. Mass threshold for ‘harmless’ drones. Int. J. Micro Air. Veh. 2017, 9, 77–92. [CrossRef] Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional aliations. © 2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Aerospace – Multidisciplinary Digital Publishing Institute

**Published: ** Oct 28, 2020

Loading...

You can share this free article with as many people as you like with the url below! We hope you enjoy this feature!

Read and print from thousands of top scholarly journals.

System error. Please try again!

Already have an account? Log in

Bookmark this article. You can see your Bookmarks on your DeepDyve Library.

To save an article, **log in** first, or **sign up** for a DeepDyve account if you don’t already have one.

Copy and paste the desired citation format or use the link below to download a file formatted for EndNote

Access the full text.

Sign up today, get DeepDyve free for 14 days.

All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.