Search NASA⌕ Search

SEARCH · Search NASA

Results for “Trajectory Clustering”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

Visualizing Corridors in Terminal Airspace using Trajectory Clustering

Context: Advances in battery and automation technology have made routine air taxi and cargo transport in urban areas a business model that can be attained by emerging aviation innovators. The community vision and work to enable these novel operations is discussed using the term ‘Urban Air Mobility’ or UAM. Small, piloted, airspace vehicles that fly with a few passengers do operate in urban areas today, and these vehicles can be studied as an early proxy for this future UAM traffic. Aim: We seek to identify corridors already in daily operation and their properties. Method: We applied DBSCAN and HDBSCAN to Dallas Forth-Worth TRACON flight data to identify corridors in use, their density, and devised a method to annotate landing sites used in these corridors with site metadata. Results: While DBSCAN was unable to group similar trajectories, we we were able to successfully identify corridors using HDBSCAN, measure their density and annotate them. Conclusion: The applied method can successfully identify corridors in daily operation with additional metadata to help domain expert understand the intent of UAM corridors.

UAM, Trajectory, TRACON, Clustering, DBSCAN, HDBSC↗

Visualizing Corridors in Terminal Airspace Using Trajectory Clustering

Context: Advances in battery and automation technology have made routine air taxi and cargo transport in urban areas a business model that can be attained by emerging aviation innovators. The community vision and work to enable these novel operations is discussed using the term ‘Urban Air Mobility’ or UAM. Small, piloted, airspace vehicles that fly with a few passengers do operate in urban areas today, and these vehicles can be studied as an early proxy for this future UAM traffic. Aim: We seek to identify corridors already in daily operation and their properties. Method: We applied DBSCAN and HDBSCAN to Dallas Forth-Worth TRACON flight data to identify corridors in use, their density, and devised a method to annotate landing sites used in these corridors with site metadata. Results: While DBSCAN was unable to group similar trajectories, we we were able to successfully identify corridors using HDBSCAN, measure their density and annotate them. Conclusion: The applied method can successfully identify corridors in daily operation with additional metadata to help domain expert understand the intent of UAM corridors.

UAM Trajectory, TRACON, Clustering, DBSCAN, HDBSCA↗

Monitoring Airspace Complexity and Determining Contributing Factors

The national airspace has evolved over many years to accommodate increased traffic demand [1] while simultaneously maintaining one of the safest forms of transportation [2], [3]. One of the reasons for this success is the ability of the system and the operators to adapt and accommodate to situations that routinely disrupt optimal operations. These situations may include: adverse weather, delays, early arrivals, equipment outages, and other factors that are outside the operators ability to control. These factors can lead to states where automation is unable to properly handle these issues and therefore air traffic controllers and pilots have to intervene, ultimately increasing communication between operators resulting in higher workload. As controller workload increases to handle sub-optimal operating conditions this can be viewed as an increase in complexity. The reasoning for this is because humans are now required to make tactical decisions in response to external factors, resulting in a departure from the strategic plan where operations would be more efficiently managed. Human operators control airspace complexity under rigid regulations that are constantly changing. The airspace is divided into sectors and the number of aircraft assigned to each controller is limited for safe handling. There has been past work that devised airspace complexity metrics in commercial aviation and related these metrics to controller workload (e.g., [4],[5]). The upper bounds on the system load are pre-determined. Such bounds on complexity make for a safe system, but the system cannot scale and adapt to autonomous, dense, and heterogeneous traffic, including the many types of Unmanned Aerial Vehicles (UAVs) envisioned to be added to the operations. We hypothesize that, as traffic density and heterogeneity grow, and other key metrics change, there will be phase transitions at which the way traffic should be managed changes significantly [6]. We offer a method for in-time detection of contributing factors that lead to phase transitions, characterized by increased complexity. To the best of our knowledge, there is no tool similar to our proposed effort that identifies such contributing factors or precursor patterns. To define the scope we are proposing to measure complexity from the viewpoint of the Terminal Radar Approach Control Facilities (TRACON) controller’s perspective. In particular we are analyzing arrivals into KSFO. With safety as the top concern for airspace operators, it is important to recognize that as density and heterogeneity grow, the focus of the system will change. Times of the day when the airspace has low density and heterogeneity, the flights will follow more efficient paths where the aircraft move on established routes that are more or less directly to the destination. However, when density and heterogeneity increases, the system will begin changing focus to avoiding conflicts and collisions and route the flights in a more flexible way. Higher flexibility requires more communication and coordination between controllers and pilots which the current automation is unable to handle. This paper proposes a novel approach that monitors airspace complexity at multiple scales, uses a Machine Learning-based tool that predicts when operations will transition to a regime of greater complexity, and identifies actions that can reduce the complexity while still maintaining efficient and safe operations. We demonstrate our proposed approach using data from multiple complementary sources. This includes, but is not limited to: historical aircraft surveillance data from NASA’s Sherlock Data Warehouse [7], METAR weather data, and airport configuration data from Aviation System Performance Metrics (ASPM). The surveillance data flight paths are sampled at a variable sample rate — increasing as the aircraft approaches the airport. This is due to how Sherlock manages flight track stitching between different radar facilities which have different sampling rates. The weather and performance data are logged at defined intervals throughout the day at a courser refresh rate. In addition to the logged data and metrics, we leverage pre-defined Standard Terminal Arrival Routes (STARs) procedures to characterize the path of each flight. Each flight files for one of these routes in the flight plan well before entering the terminal airspace, and approximately follows the route until it leaves the STAR, typically on the final fix of a runway transition. However, most flights do not always fly the full STAR procedure to completion [8], but the majority do adhere to the fixes within the common route of the procedure. Our approach leverages fixes in the common route of each of the STARs to build a reference path to the airport. This allows us to characterize the flight paths in what we are defining as the “maneuvering area” (the airspace between the STAR and before the flight is lined up on the runway’s final approach) to determine how off nominal the flights are to calculate its complexity score. Determining airspace complexity is a concept that does not have a concrete answer. In designing this metric, we consider what increases the workload for the air traffic controllers. Consequently more specialized vectoring maneuvers results in higher workload. Accordingly, we start with a theory: each flight has a direct path it takes from the STAR’s common route to the final approach’s outer marker fix for the flight’s landing runway. It is important to note that the direct path is only used as a reference. If the majority of the flights have a large consistent offset as compared to other routes it does not necessarily mean that those flights have higher complexity. We are merely building a distribution based on this direct path for that particular STAR and runway pair to determine the normal mode of operations for that route. Flights that are in the upper tail of these distributions will result in higher complexity scores and flights that fly in the median will represent the normal mode of operations and therefore will have lower complexity scores. Since flights following each STAR route take different paths to the airport, we have a different distribution for each STAR route and therefore can model these distributions to compute a complexity score from their respective normalized distributions. To evaluate the effectiveness of our proposed airspace complexity metric we will compare against an established approach based on trajectory clustering [9]. This unsupervised learning technique consists of the following steps: (1) identify the general maneuvering areas (waypoints) by performing $\kappa$-means or DBSCAN clustering on locations where aircraft frequently turn based on the surveillance radar track data, (2) map flight trajectories onto sequences of waypoints, and (3) cluster the sequences based on their common subsequences. From a high-level perspective, this baseline model learns nominal operations in the airspace through the sequence of waypoints that are representative of where aircraft change direction and defines deviations from the nominal operations as “complex.” Therefore, more deviations from the nominal operations correspond to higher complexity values. For our validation, we re-implemented this technique and tune model hyper-parameters to correctly detect waypoints for the arrival traffic into the San Francisco bay area. We will compute the complexity measure over a one-year period using our proposed technique as well as the baseline. Our validation will be based on each technique’s ability to detect a set of undesirable outcomes (e.g., go-arounds, holding patterns, average time in the airspace, etc.). Since our current complexity metric is derived from the offset from the direct reference path, it’s important to understand what causes these offsets. In many of the flights with high offset distance, flights performing holding patterns and S turns can be observed. These maneuvering tactics are utilized to add distance between the aircraft and the destination runway to prevent multiple flights from having conflicting arrival times. In order to predict a rise in complexity (or the precursor to complexity), it’s necessary to be able to identify these potential conflicts (which in turn, result in higher offsets). To do this, we define a “representative flight” for each STAR route and runway pair. This flight is approximately the path the flight would take if there was a clear path with no other flights in the airspace — including the time remaining to the airport. We first identify the flights for a given STAR runway pair using the offset to the reference path distributions that fall between the 44-55 percentiles. This yields the flights that conform to the most normal mode of operation. Each of these flights is partitioned based on the percent complete from the entry point into the maneuvering areas from 0\% – 100\% complete. Then for each percent “bin”, we take the median value of the flight’s latitude/longitude coordinates, airspeed, and (non causal) time remaining to the airport to construct a lookup table for each percent complete bin on a given route. As a flight enters the maneuvering area, we can find the estimated arrival time of a flight to the airport by finding the closest point to the representative path’s percent complete bin (relative to the flight’s current position at any snapshot in the airspace) and therefore retrieve the corresponding remaining time left on the “representative path”. We assume that the flight will follow the representative path to completion when deriving these estimates. We can then compare these estimated arrival times against other flights for the same snapshot in time to identify potential conflicts. If more flights are estimated to arrive within a tolerance window than there are runways available, then we have a potential conflict. We can use this derived measure along with other factors expected to add disruption to the operation such as weather and runway configuration changes as an input to machine learning tools to detect precursors that increases in our complexity measure. This novel method will assist in uncovering insights into the contributing factors that lead to increased complexity that may allow for in-time responses to avoid reaching a high complexity state in the airspace.

complexity↗

A Seasonal Air Transport Climatology for Kenya

A climatology of air transport to and from Kenya has been developed using kinematic trajectory modeling. Significant months for trajectory analysis have been determined from a classification of synoptic circulation fields. Five-point back and forward trajectory clusters to and from Kenya reveal that the transport corridors to Kenya are clearly bounded and well defined. Air reaching the country originates mainly from the Saharan region and northwestern Indian Ocean of the Arabian Sea in the northern hemisphere and from the Madagascan region of the Indian Ocean in the southern hemisphere. Transport from each of these source regions show distinctive annual cycles related to the northeasterly Asian monsoon and the southeasterly trade wind maximum over Kenya in May. The Saharan transport in the lower troposphere is at a maximum when the subtropical high over northern Africa is strongly developed in the boreal winter. Air reaching Kenya between 700 and 500 hPa is mainly from Sahara and northwest India Ocean flows in the months of January and March, which gives way to southwest Indian Ocean flow in May and November. In contrast, air reaching Kenya at 400 hPa is mainly from southwest Indian Ocean in January and March, which is replaced by Saharan transport in May and November. Transport of air from Kenya is invariant, both spatially and temporally, in the tropical easterlies to the Congo Basin and Atlantic Ocean in comparison to the transport to the country. Recirculation of air has also been observed, but on a limited and often local scale and not to the extent reported in southern Africa.

Gatebe, C. K.↗

Cluster-Based Flight Trajectory Outlier Detection

Given a set of flight trajectories, can we classify the trajectories that do not follow a normal path? By identifying abnormal trajectories, further analysis can be done to determine the reasoning for these actions. Addressing these scenarios can bring possible solutions for holding and rerouting problems when its time to incorporate UAM in the airspace.

Machine Learning↗

Aerosol Seasonal Variations over Urban-Industrial Regions in Ukraine According to AERONET and POLDER Measurements

The paper presents an investigation of aerosol seasonal variations in several urban-industrial regions in Ukraine. Our analysis of seasonal variations of optical and physical aerosol parameters is based on the sun-photometer 2008-2013 data from two urban ground-based AERONET (AErosol RObotic NETwork) sites in Ukraine (Kyiv, Lugansk) as well as on satellite POLDER instrument data for urban-industrial areas in Ukraine. We also analyzed the data from one AERONET site in Belarus (Minsk) in order to compare with the Ukrainian sites. Aerosol amount and optical depth (AOD) values in the atmosphere columns over the large urbanized areas like Kyiv and Minsk have maximum values in the spring (April-May) and late summer (August), whereas minimum values are observed in late autumn. The results show that fine-mode particles are most frequently detected during the spring and late summer seasons. The analysis of the seasonal AOD variations over the urban-industrial areas in the eastern and central parts of Ukraine according to both ground-based and POLDER data exhibits the similar traits. The seasonal variation similarity in the regions denotes the resemblance in basic aerosol sources that are closely related to properties of aerosol particles. The behavior of basic aerosol parameters in the western part of Ukraine is different from eastern and central regions and shows an earlier appearance of the spring and summer AOD maxima. Spectral single-scattering albedo, complex refractive index and size distribution of aerosol particles in the atmosphere column over Kyiv have different behavior for warm (April-October) and cold seasons. The seasonal features of fine and coarse aerosol particle behavior over the Kyiv site were analyzed. A prevailing influence of the fine-mode particles on the optical properties of the aerosol layer over the region has been established. The back-trajectory and cluster analysis techniques were applied to study the seasonal back trajectories and prevailing directions of the arrived air mass for the Kyiv and Minsk sites.

Ukraine↗

Clustering Approach To Identifying Low Lunar Frozen Orbits In A High-Fidelity Model

Low lunar frozen orbits are of continued interest in the astrodynamics community for trajectory design and space domain awareness. This paper presents a data-driven approach to analyzing a wide variety of numerically-generated lunar trajectories in a 100 × 100 lunar gravity model with the point mass gravity of the Earth and Sun. First, clustering is used to extract a summary of these trajectories: within each cluster, trajectories possess a geometrically similar evolution of perilune but varying drift and lifetimes. Within some clusters, trajectories with a bounded perilune evolution are also identified to produce candidates for lunar frozen orbits of distinct geometries.

lunar frozen orbit gravity model data clustering↗

Reactive Collision Avoidance Algorithm

The reactive collision avoidance (RCA) algorithm allows a spacecraft to find a fuel-optimal trajectory for avoiding an arbitrary number of colliding spacecraft in real time while accounting for acceleration limits. In addition to spacecraft, the technology can be used for vehicles that can accelerate in any direction, such as helicopters and submersibles. In contrast to existing, passive algorithms that simultaneously design trajectories for a cluster of vehicles working to achieve a common goal, RCA is implemented onboard spacecraft only when an imminent collision is detected, and then plans a collision avoidance maneuver for only that host vehicle, thus preventing a collision in an off-nominal situation for which passive algorithms cannot. An example scenario for such a situation might be when a spacecraft in the cluster is approaching another one, but enters safe mode and begins to drift. Functionally, the RCA detects colliding spacecraft, plans an evasion trajectory by solving the Evasion Trajectory Problem (ETP), and then recovers after the collision is avoided. A direct optimization approach was used to develop the algorithm so it can run in real time. In this innovation, a parameterized class of avoidance trajectories is specified, and then the optimal trajectory is found by searching over the parameters. The class of trajectories is selected as bang-off-bang as motivated by optimal control theory. That is, an avoiding spacecraft first applies full acceleration in a constant direction, then coasts, and finally applies full acceleration to stop. The parameter optimization problem can be solved offline and stored as a look-up table of values. Using a look-up table allows the algorithm to run in real time. Given a colliding spacecraft, the properties of the collision geometry serve as indices of the look-up table that gives the optimal trajectory. For multiple colliding spacecraft, the set of trajectories that avoid all spacecraft is rapidly searched on-line. The optimal avoidance trajectory is implemented as a receding-horizon model predictive control law. Therefore, at each time step, the optimal avoidance trajectory is found and the first time step of its acceleration is applied. At the next time step of the control computer, the problem is re-solved and the new first time step is again applied. This continual updating allows the RCA algorithm to adapt to a colliding spacecraft that is making erratic course changes.

Scharf, Daniel↗

Trajectory Design Considerations for Exploration Mission 1

Exploration Mission 1 (EM-1) will be the first mission to send an uncrewed Orion Multi-Purpose Crew Vehicle (MPCV) to cislunar space in the fall of 2019. EM-1 was originally conceived as a lunar free-return mission, but was later changed to a Distant Retrograde Orbit (DRO) mission as a precursor to the Asteroid Redirect Mission. To understand the required mission performance (i.e., propellant requirement), a series of trajectory optimization runs was conducted using JSC's Copernicus spacecraft trajectory optimization tool. In order for the runs to be done in a timely manner, it was necessary to employ a parallelization approach on a computing cluster using a new trajectory scan tool written in Python. Details of the scan tool are provided and how it is used to perform the scans and post-process the results. Initially, a scan of daily due east launched EM-1 DRO missions in 2018 was made. Valid mission opportunities are ones that do not exceed the useable propellant available to perform the required burns. The initial scan data showed the propellant and delta-V performance patterns for each launch period. As questions were raised from different subsystems (e.g., power, thermal, communications, flight operations, etc.), the mission parameters or data that were of interest to them were added to the scan output data file. The additional data includes: (1) local launch and landing times in relation to sunrise and sunset, (2) length of eclipse periods during the in-space portion of the mission, (3) Earth line of sight from cislunar space, (4) Deep Space Network field of view looking towards cislunar space, and (5) variation of the downrange distance from Earth entry interface to splashdown. Mission design trades can also be performed based on the information that the additional data shows. For example, if the landing is in darkness, but the recovery operations team desires a landing in daylight, then an analysis is performed to determine how to change the mission design to meet this request. Also, subsystems request feasibility of alternate or contingency mission designs, such as adding an Orion main engine checkout burn or Orion completing all of its burns using only its auxiliary thrusters. This paper examines and presents the evolving trade studies that incorporate subsystem feedback and demonstrate the feasibility of these constrained mission trajectory designs and contingencies.

Dawn, Timothy F.↗

A New Architecture for Parallelization of Complex Spacecraft Trajectory Optimization Scans

This paper describes CopScanner, a new component of the Copernicus ecosystem for spacecraft trajectory design and optimization. CopScanner is a Python library being developed at the NASA JSC which enables easy parallelization of Copernicus scans. CopScanner is currently being developed and implemented for production of Copernicus trajectory scans for upcoming Artemis Missions (Artemis II and beyond). On the backend, CopScanner utilizes Dask, an open-source Python library for parallel computing which enables parallelization over both multi-core local machines and large-scale distributed computing clusters. CopScanner abstracts the trajectory scanning process into a DAG which is constructed using a chain of individual subscans. Each node in the DAG executes a python module, called the callable, for which there are built-in defaults, or users may specify their own. Support for custom callables makes CopScanner a versatile trajectory optimization software. All output files and associated metadata from a CopScanner scan are compressed and stored in a two-file output, collectively called the FileStore, consisting of a SQLite database and a compressed JSON MessagePack file, for which CopScanner provides a Python class for interaction.

Quentin Moore↗

Origin, Transport, and Vertical Distribution of Atmospheric Polluntants over the Northern Sourth China Sea During the 7-SEAS-Dongsha Experiment

During the spring of 2010, comprehensive in situ measurements were made for the first time on a small atoll (Dongsha Island) in the northern South China Sea (SCS), a key region of the 7-SEAS (the Seven South East Asian Studies) program. This paper focuses on characterizing the source origins, transport processes, and vertical distributions of the Asian continental outflows over the region, using measurements including mass concentration, optical properties, hygroscopicity, and vertical distribution of the aerosol particles, as well as the trace gas composition. Cluster analysis of backward trajectories classified 52% of the air masses arriving at ground level of Dongsha Island as having a continental origin, mainly from northern China to the northern SCS, passing the coastal area and being confined in the marine boundary layer (0-0.5 km). Compared to aerosols of oceanic origin, the fine mode continental aerosols have a higher concentration, extinction coefficient, and single-scattering albedo at 550 nm (i.e., 19 vs. 14 microg per cubic meter in PM(sub 2.5); 77 vs. 59 M per meter in beta(sub e); and 0.94 vs. 0.90 in omega, respectively). These aerosols have a higher hygroscopicity (f at 85% RH = 2.1) than those in the upwind inland regions, suggesting that the aerosols transported to the northern SCS were modified by the marine environment. In addition to the near-surface aerosol transport, a significant upper-layer (3-4 km) transport of biomass-burning aerosols was observed. Our results suggest that emissions from both China and Southeast Asia could have a significant impact on the aerosol loading and other aerosol properties over the SCS. Furthermore, the complex vertical distribution of aerosols-coinciding-with-clouds has implications for remote-sensing observations and aerosol-cloud-radiation interactions.

pollutants↗

Radiative Effect of Springtime Biomass-Burning Aerosols over Northern Indochina During 7-SEAS Baseline 2013 Campaign

The direct aerosol radiative effects of biomass-burning (BB) aerosols over northern Indochina were estimated by using aerosol properties (physical, chemical, and optical) along with the vertical profile measurements from ground-based measurements with integration of an optical and a radiative transfer model during the Seven South East Asian Studies Biomass-Burning Aerosols Stratocumulus Environment: Lifecycles Interactions Experiment (7-SEASBASELInE) conducted in spring 2013. Cluster analysis of backward trajectories showed the air masses arriving at mountainous background site (Doi Ang Khang; 19.93degN, 99.05degE, 1536 m above mean sea level) in northern Indochina, mainly from near-source inland BB activities and being confined in the planetary boundary layer. The PM(sub10) and black carbon (BC)mass were 87 +/- 28 and 7 +/- 2 micrograms m(exp -3), respectively. The aerosol optical depth (AOD (sub 500) was found to be 0.26--1.13 (0.71 +/- 0.24). Finer (fine mode fraction is approximately or equal to 0.95, angstrom-exponent at 440-870 nm is approximately or equal to 1.77) and significantly absorbing aerosols(single scattering albedo is approximately or equal to 0.89, asymmetry-parameter is approximately or equal to 0.67, and absorption AOD 0.1 at 440 nm) dominated over this region. BB aerosols (water soluble and BC) were the main contributor to the aerosol radiative forcing (ARF), while others (water insoluble, sea salt and mineral dust) were negligible mainly due to their low extinction efficiency. BC contributed only 6 to the surface aerosol mass but its contribution to AOD was 12 (2 times higher). The overall mean ARF was 8.0 and -31.4 W m(exp -2) at top-of-atmosphere (TOA) and at the surface (SFC), respectively. Likely, ARF due to BC was +10.7 and -18.1 W m(exp -2) at TOA and SFC, respectively. BC imposed the heating rate of +1.4 K d(exp -1) within the atmosphere and highlighting its pivotal role in modifying the radiation budget. We propose that to upgrade our knowledge on BB aerosol radiative effects in BB source region, the long-term and extensive field measurements are needed.

aerosol optical properties↗

Machine Learning for Biological Trajectory Classification Applications

Machine-learning techniques, including clustering algorithms, support vector machines and hidden Markov models, are applied to the task of classifying trajectories of moving keratocyte cells. The different algorithms axe compared to each other as well as to expert and non-expert test persons, using concepts from signal-detection theory. The algorithms performed very well as compared to humans, suggesting a robust tool for trajectory classification in biological applications.

Sbalzarini, Ivo F.↗

Computational Approaches to Simulation and Optimization of Global Aircraft Trajectories

This study examines three possible approaches to improving the speed in generating wind-optimal routes for air traffic at the national or global level. They are: (a) using the resources of a supercomputer, (b) running the computations on multiple commercially available computers and (c) implementing those same algorithms into NASAs Future ATM Concepts Evaluation Tool (FACET) and compares those to a standard implementation run on a single CPU. Wind-optimal aircraft trajectories are computed using global air traffic schedules. The run time and wait time on the supercomputer for trajectory optimization using various numbers of CPUs ranging from 80 to 10,240 units are compared with the total computational time for running the same computation on a single desktop computer and on multiple commercially available computers for potential computational enhancement through parallel processing on the computer clusters. This study also re-implements the trajectory optimization algorithm for further reduction of computational time through algorithm modifications and integrates that with FACET to facilitate the use of the new features which calculate time-optimal routes between worldwide airport pairs in a wind field for use with existing FACET applications. The implementations of trajectory optimization algorithms use MATLAB, Python, and Java programming languages. The performance evaluations are done by comparing their computational efficiencies and based on the potential application of optimized trajectories. The paper shows that in the absence of special privileges on a supercomputer, a cluster of commercially available computers provides a feasible approach for national and global air traffic system studies.

global air traffic optimization↗

Towards Persistent Space Observations through Autonomous Multi-Agent Formations

Sensing platforms must advance in scale and sophistication in order to support increasingly ambitious missions across Earth and space science; intelligence, surveillance, and reconnaissance (ISR); and planetary exploration. Distributed, persistent observation platforms have the potential to play a pivotal role in next generation missions through improved area coverage, enhanced situational awareness, and faster identification of trends and changes. The Multi-Agent Clusters for Persistent Observations from Space (MACPOS) project at NASA Langley Research Center is developing key technologies for the autonomous, heterogeneous formations that will comprise such platforms. Research thrusts include dynamic formation negotiation for self-assembling clusters of agents, distributed motion planning, and coordinated trajectory execution. Adaptive leader-follower formation negotiation allows agents to cluster and break off as necessary to adapt to both nominal and new mission objectives. Coordinated motion planning and execution maintain the formation while ensuring safe separation distances among agents and obstacles in the environment. These capabilities align MACPOS with NASA’s initiative for space and surface in-situ assembly through fundamental technology development for autonomous multi-agent systems. This paper presents an overview and early progress for the MACPOS project. We describe the system architecture for both individual agents and the overall fleet. Design considerations are given for the planning, control, and metrology subsystems. Finally, we discuss planned project milestones and the expected course of development.

Matthew P Vaughan↗

Agricultural Burning and Air Quality over Northern India: A Synergistic Analysis Using NASA's A-Train Satellite Data and Ground Measurements

In the recent years, New Delhi, the capital city of India, has ranked among the most polluted cities in the world regarding its air quality related to the submicron Particulate Matter (PM(sub 2.5)). Using NASA's A-train satellite data (MODIS, OMI, and CALIOP), ground-level PM2.5 measured in New Delhi (2013-2016), and back-trajectory calculations, we show that the PM(sub 2.5). over New Delhi is strongly affected by the agricultural fires in the northwestern Indian states of Punjab and Haryana during the post-monsoon season (October and November). The mass concentration of PM2.5 escalates from ~50 μg m(exp -3) measured prior to the onset of residue burning in early October to as high as 300 μg m(exp -3) (24-hour averaged, 7- day running mean) during the peak burning period in early November. A linear regression analysis reveals that the variations in PM(sub 2.5) over New Delhi can be attributed to the concurrent changes in the satellite retrievals of fire counts and aerosols over the crop burning area. The back-trajectory analysis shows that most clusters (> 80%) of the northwesterly flow near the ground intercepted the crop burning region before arriving at the receptor location in New Delhi; this further corroborates the transport patterns inferred from the satellite data. A 15-year long satellite record (2002-2016) reveals an increasing trend in agricultural fires (~617 per year) and aerosol loading (0.031 and 0.04 per year in aerosol optical depth and UV aerosol index) in November. Increasing levels of crop residue burning and resulting particulate matter pollution at an alarming rate over northern India is a pressing concern demanding corrective measures to substantially reduce or completely diminish the crop burning through an effective residue management system.

Crop Residue Fires↗

Computational Approaches to Simulation and Optimization of Global Aircraft Trajectories

This study examines three possible approaches to improving the speed in generating wind-optimal routes for air traffic at the national or global level. They are: (a) using the resources of a supercomputer, (b) running the computations on multiple commercially available computers and (c) implementing those same algorithms into NASA’s Future ATM Concepts Evaluation Tool (FACET) and compares those to a standard implementation run on a single CPU. Wind-optimal aircraft trajectories are computed using global air traffic schedules. The run time and wait time on the supercomputer for trajectory optimization using various numbers of CPUs ranging from 80 to 10,240 units are compared with the total computational time for running the same computation on a single desktop computer and on multiple commercially available computers for potential computational enhancement through parallel processing on the computer clusters. This study also re-implements the trajectory optimization algorithm for further reduction of computational time through algorithm modifications and integrates that with FACET to facilitate the use of the new features which calculate time-optimal routes between worldwide airport pairs in a wind field for use with existing FACET applications. The implementations of trajectory optimization algorithms use MATLAB, Python, and Java programming languages. The performance evaluations are done by comparing their computational efficiencies and based on the potential application of optimized trajectories. The paper shows that in the absence of special privileges on a supercomputer, a cluster of commercially available computers provides a good option for computing wind-optimal trajectories for national and global air traffic system studies.

Ng, Hok K.↗