Search NASA⌕ Search

SEARCH · Search NASA

Results for “uncertainty in complex systems”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Ar-39-Ar-40 Ages of Euerites and the Thermal History of Asteroid 4-Vesta

Eucrite meteorites are igneous rocks that derive from a large asteroid, probably 4 Vesta. Prior studies have shown that after eucrites formed, most were subsequently metamorphosed to temperatures up to equal to or greater than 800 C, and much later many were brecciated and heated by large impacts into the parent body surface. The uncommon basaltic, unbrecciated eucrites also formed near the surface but presumably escaped later brecciation, whereas the cumulate eucrites formed at depth where metamorphism may have persisted for a considerable period. To further understand the complex HED parent body thermal history, we determined new Ar-39-Ar-40 ages for nine eucrites classified as basaltic but unbrecciated, six eucrites classified as cumulate, and several basaltic-brecciated eucrites. Relatively precise Ar-Ar ages of two cumulate eucrites (Moama and EET87520) and four unbrecciated eucrites give a tight cluster at 4.48 +/1 0.01 Gyr. Ar-Ar ages of six additional unbrecciated eucrites are consistent with this age, within their larger age uncertainties. In contrast, available literature data on Pb-Pb isochron ages of four cumulate eucrites and one unbrecciated eucrite vary over 4.4-4.515 Gyr, and Sm-147 - Nd-143 isochron ages of four cumulate and three unbrecciated eucrites vary over 4.41-4.55 Gyr. Similar Ar-Ar ages for cumulate and unbrecciated eucrites imply that cumulate eucrites do not have a younger formation age than basaltic eucrites, as previously proposed. Rather, we suggest that these cumulate and unbrecciated eucrites resided at depth where parent body temperatures were sufficiently high to cause the K-Ar and some other chronometers to remain open diffusion systems. From the strong clustering of Ar-Ar ages at approximately 4.48 Gyr, we propose that these meteorites were excavated from depth in a single large impact event approximately 4.48 Gyr ago, which quickly cooled the samples and started the K-Ar chronometer. A large (approximately 460 km) crater postulated to exist on Vesta may be the source of these eucrites and of many smaller asteroids thought to be spectrally or physically associated with Vesta. Some Pb-Pb and Sm-Nd ages of cumulate and unbrecciated eucrites are consistent with the 4.48 Gyr Ar-Ar age, and the few older Pb-Pb and Sm-Nd ages may reflect isotopic closure prior to the large cratering event. One cumulate eucrite gives an Ar-Ar age of 4.25 Gyr; three additional cumulate eucrites give Ar-Ar ages of 3.4-3.7 Gyr; and two unbrecciated eucrites give Ar-Ar ages of approximately 3.55 Gyr. We attribute these younger ages to later impact heating. In addition, we find Ar-Ar impact-reset ages of several brecciated eucrites and eucritic clasts in howardites to fall in the range of 3.5-4.1 Gyr. Among these, Piplia Kalan, the first eucrite to show evidence for extinct 26 Al, was strongly impact heated approximately3.5 Gyr ago. When these data are combined with eucrite Ar-Ar ages in the literature, they confirm the previous suggestion that several large impact heating events occurred on Vesta over the time period approximately 4.1-3.4 Gyr ago. The onset of major impact heating may have occurred at similar times for both Vesta and the Moon, but impact heating appears to have persisted to a somewhat later time on Vesta compared to the Moon.

Bogard, Donald D.↗
NASA NTRSDate not supplied

Zircon (U-Th)/He Impact Crater Thermochronometry and the Effects of Shock Microstructures on Helium Diffusion Kinetics

Absolute age determination of impact cratering events remains difficult and often controversial; a challenge that has resulted in < 50% of known impact craters regarded as accurately and precisely dated. Besides conventional 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He (ZHe) dating of impactites has been applied to large- to medium-sized impact structures. ZHe dates can be fully reset in minutes at 1000°C, which is commonly reached in central sections of the melt sheet, whereas resetting of ZHe at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-6 years. There is a critical need to quantify the effects of shock-induced microstructures and impact metamorphism on helium diffusion kinetics in well-characterized, variably shocked zircon to further establish the reliability of (U-Th)/He for dating impacts. For this purpose, we investigated suevite and impact melt samples from two impact structures, the Chicxulub multi-ring basin and the Ries complex crater, which enables us to compare zircon helium diffusion kinetics from impact structures with differing sizes, ages, and hydrothermal system longevities. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incremental 10°C steps from 250°C to 600°C. Afterward, we characterize the diffusion domain sizes and their interconnectivity within the shocked zircon grains using electron backscatter diffraction (EBSD). We find that zircon with few shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features (PDFs) and granular textures classified by BSE and EBSD are characterized by a dramatic decrease in helium retentivity, similar to radiation damaged grains, due to a reduction in the effective domain size and the introduction of interconnected fast diffusion pathways created by shock microstructures. A subset of grains were dated by ZHe after the external morphology of the grains was determined by BSE imaging. The euhedral grains yielded a weighted mean age within the uncertainty of the accepted impact ages, whereas the grains with PDFs or granular textures exhibited younger ages. Thus, these diffusion experiments and ZHe dates suggest that the dramatic decrease in domain size likely renders shocked grains more susceptible to impact-induced hydrothermal resetting and subsequent overprinting. Hence, characterization of shock microstructures is critical for determining accurate impact ages using ZHe methods especially when applied to previously unconstrained craters. The thermochronometer also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

zircon↗
NASA NTRS2020

Application of Sparse Identification of Nonlinear Dynamics for Physics-Informed Learning

Advances in machine learning and deep neural networks has enabled complex engineering tasks like image recognition, anomaly detection, regression, and multi-objective optimization, to name but a few. The complexity of the algorithm architecture, e.g., the number of hidden layers in a deep neural network, typically grows with the complexity of the problems they are required to solve, leaving little room for interpreting (or explaining) the path that results in a specific solution. This drawback is particularly relevant for autonomous aerospace and aviation systems, where certifications require a complete understanding of the algorithm behavior in all possible scenarios. Including physics knowledge in such data-driven tools may improve the interpretability of the algorithms, thus enhancing model validation against events with low probability but relevant for system certification. Such events include, for example, spacecraft or aircraft sub-system failures, for which data may not be available in the training phase. This paper investigates a recent physics-informed learning algorithm for identification of system dynamics, and shows how the governing equations of a system can be extracted from data using sparse regression. The learned relationships can be utilized as a surrogate model which, unlike typical data-driven surrogate models, relies on the learned underlying dynamics of the system rather than large number of fitting parameters. The work shows that the algorithm can reconstruct the differential equations underlying the observed dynamics using a single trajectory when no uncertainty is involved. However, the training set size must increase when dealing with stochastic systems, e.g., nonlinear dynamics with random initial conditions.

Corbetta, Matteo↗
NASA NTRS2014

A Multi-Parameter Approach for Calculating Crack Instability

An accurate fracture control analysis of spacecraft pressure systems, boosters, rocket hardware and other critical low-cycle fatigue cases where the fracture toughness highly impacts cycles to failure requires accurate knowledge of the material fracture toughness. However, applicability of the measured fracture toughness values using standard specimens and transferability of the values to crack instability analysis of the realistically complex structures is refutable. The commonly used single parameter Linear Elastic Fracture Mechanics (LEFM) approach which relies on the key assumption that the fracture toughness is a material property would result in inaccurate crack instability predictions. In the past years extensive studies have been conducted to improve the single parameter (K-controlled) LEFM by introducing parameters accounting for the geometry or in-plane constraint effects]. Despite the importance of the thickness (out-of-plane constraint) effects in fracture control problems, the literature is mainly limited to some empirical equations for scaling the fracture toughness data] and only few theoretically based developments can be found. In aerospace hardware where the structure might have only one life cycle and weight reduction is crucial, reducing the design margin of safety by decreasing the uncertainty involved in fracture toughness evaluations would result in lighter hardware. In such conditions LEFM would not suffice and an elastic-plastic analysis would be vital. Multi-parameter elastic plastic crack tip field quantifying developments combined with statistical methods] have been shown to have the potential to be used as a powerful tool for tackling such problems. However, these approaches have not been comprehensively scrutinized using experimental tests. Therefore, in this paper a multi-parameter elastic-plastic approach has been used to study the crack instability problem and the transferability issue by considering the effects of geometrical constraints as well as the thickness. The feasibility of the approach has been examined using a wide range of specimen geometries and thicknesses manufactured from 7075-T7351 aluminum alloy.

Zanganeh, M.↗
NASA NTRS1990

Situationally Driven Local Navigation for Mobile Robots

For mobile robots to autonomously accommodate dynamically changing navigation tasks in a goal-directed fashion, they must employ navigation plans. Any such plan must provide for the robot's immediate and continuous need for guidance while remaining highly flexible in order to avoid costly computation each time the robot's perception of the world changes. Due to the world's uncertainties, creation and maintenance of navigation plans cannot involve arbitrarily complex processes, as the robot's perception of the world will be in constant flux, requiring modifications to be made quickly if they are to be of any use. This work introduces navigation templates (NaT's) which are building blocks for the construction and maintenance of rough navigation plans which capture the relationship that objects in the world have to the current navigation task. By encoding only the critical relationship between the objects in the world and the navigation task, a NaT-based navigation plan is highly flexible; allowing new constraints to be quickly incorporated into the plan and existing constraints to be updated or deleted from the plan. To satisfy the robot's need for immediate local guidance, the NaT's forming the current navigation plan are passed to a transformation function. The transformation function analyzes the plan with respect to the robot's current location to quickly determine (a few times a second) the locally preferred direction of travel. This dissertation presents NaT's and the transformation function as well as the needed support systems to demonstrate the usefulness of the technique for controlling the actions of a mobile robot operating in an uncertain world.

Slack, Marc Glenn↗
NASA NTRS2010

Simplified Vicarious Radiometric Calibration

A measurement-based radiance estimation approach for vicarious radiometric calibration of spaceborne multispectral remote sensing systems has been developed. This simplified process eliminates the use of radiative transfer codes and reduces the number of atmospheric assumptions required to perform sensor calibrations. Like prior approaches, the simplified method involves the collection of ground truth data coincident with the overpass of the remote sensing system being calibrated, but this approach differs from the prior techniques in both the nature of the data collected and the manner in which the data are processed. In traditional vicarious radiometric calibration, ground truth data are gathered using ground-viewing spectroradiometers and one or more sun photometer( s), among other instruments, located at a ground target area. The measured data from the ground-based instruments are used in radiative transfer models to estimate the top-of-atmosphere (TOA) target radiances at the time of satellite overpass. These TOA radiances are compared with the satellite sensor readings to radiometrically calibrate the sensor. Traditional vicarious radiometric calibration methods require that an atmospheric model be defined such that the ground-based observations of solar transmission and diffuse-to-global ratios are in close agreement with the radiative transfer code estimation of these parameters. This process is labor-intensive and complex, and can be prone to errors. The errors can be compounded because of approximations in the model and inaccurate assumptions about the radiative coupling between the atmosphere and the terrain. The errors can increase the uncertainty of the TOA radiance estimates used to perform the radiometric calibration. In comparison, the simplified approach does not use atmospheric radiative transfer models and involves fewer assumptions concerning the radiative transfer properties of the atmosphere. This new technique uses two neighboring uniform ground target areas having different reflectance values. The target areas can be natural or artificial and must be large enough to minimize adjacent-pixel contamination effects. The radiative coupling between the atmosphere and the terrain needs to be approximately the same for the two targets. This condition can be met for relatively uniform backgrounds when the distance between the targets is within a few hundred meters. For each target area, the radiance leaving the ground in the direction of the satellite is measured with a radiometrically calibrated spectroradiometer. Using the radiance measurements from the two targets, atmospheric adjacency and atmospheric scattering effects can be subtracted, thereby eliminating many assumptions about the atmosphere and the radiative interaction between the atmosphere and the terrain. In addition, the radiometrically calibrated spectroradiometer can be used with a known reflectance target to estimate atmospheric transmission and diffuse- to-global ratios without the need for ancillary sun photometers. Several comparisons between the simplified method and traditional techniques were found to agree within a few percent. Hence, the simplified method reduces the overall complexity of performing vicarious calibrations and can serve as a method for validating traditional radiative transfer models

Stanley, Thomas↗
NASA NTRS2024

Evaluating Differences Among Crop Models in Simulating Soybean in-Season Growth

Crop models are useful tools for simulating agricultural systems that require continued model development and testing to increase their robustness and improve how they describe our current understanding of processes. Coordinated and “blind” evaluation of multiple models using same protocols and experimental datasets provides unique opportunities to further improve models and enhance their reliability. For soybean [Glycine max (L.) Merr.], there has been limited coordinated multi-model evaluations for the simulation of in-season plant growth dynamics. We evaluated ten dynamic soybean crop models for their simulation of in-season plant growth using data from five experiments conducted in Argentina, Brazil, France, and USA. We evaluated models after a Blind (using only phenology data) and a Full calibration (with in-season and end-of-season variables). Calibration reduced model uncertainty by reducing standard bias for the simulation of in-season variables (biomass, leaf, pod, and stem weights, and leaf area index, LAI). However, we found that most models had difficulty in reproducing leaf growth dynamics, with normalized root mean squared error (nRMSE) of 56% for leaf weight and 43% for LAI (across locations and models after Full calibration). Models with different levels of complexity and experience were capable of simulating final seed yield at maturity with reasonable accuracy (nRMSE of 8–31% after Full calibration). However, the nRMSE for pod weight (of 17–64% after Full calibration) was two-fold larger than that of seed yield. Moreover, the models differed in how they simulated timing from sowing to beginning seed growth (47–93 days) and effective seed filling period (18–54 days), owing to model structural differences in defining the reproductive developmental stages. Overall, we identified the following processes that can benefit from further model improvement: leaf expansion and senescence, reproductive phenology, and partitioning to reproductive growth. Simulation of pod wall tissue and individual seed cohorts is another aspect that many models currently lack. Model improvement can benefit from high-temporal resolution experimental datasets that concurrently account for phenology, plant growth, and partitioning. Further, we recommend collecting reproductive phenology in the field consistent with actual dry matter allocation to organs in the models and collecting multiple observations of seed and pod weight to aid model improvement for simulation of seed growth and yield formation.

Agricultural Model Intercomparison and Improvement↗
NASA NTRS2020

150 Shades of Green: Using the Full Spectrum of Remote Sensing Reflectance to Elucidate Color Shifts in the Ocean

This article proposes a simple and intuitive classification system by which to define full spectral remote sensing reflectance (Rrs(λ)) data with a quantitative output that enables a more manageable handling of spectral information for aquatic science applications. The weighted harmonic mean of the Rrs(λ) wavelengths outputs an Apparent Visible Wavelength (in units of nanometers), representing a one-dimensional geophysical metric of color that is inherently correlated to spectral shape. This dimensionality reduction of spectral information combined with the output along a continuum of wavelength values offers a robust and user-friendly means to describe and analyze spectral Rrs(λ) in terms of spatial and temporal trends and variability. The uncertainty in the algorithm's estimation of spectral shape is demonstrated on a global scale, in addition to the utility of the algorithm to discern spectral-spatial-temporal trends in the ocean, on a per-pixel basis for the entire 22 year continuous ocean color (SeaWiFS and MODIS-Aqua) time-series. This technique can be applied to datasets of varying multi- and hyper-spectral resolutions, providing continuity between heritage and future satellite sensors, and further enabling an effective means of elucidating similarities or differences in complex spectral signatures within the constraints of two dimensions. This straightforward means of conceptualizing multi-dimensional variability can help maximize the potential of the spectral information embedded in remote sensing data.

ocean color↗
NASA NTRS2008

AERONET-OC: Strengths and Weaknesses of a Network for the Validation of Satellite Coastal Radiometric Products

The Ocean Color component of the Aerosol Robotic Network (AERONET-OC) has been implemented to support long-term satellite ocean color investigations through cross-site consistent and accurate measurements collected by autonomous radiometer systems deployed on offshore fixed platforms. The ultimate purpose of AERONET-OC is the production of standardized measurements performed at different sites with identical measuring systems and protocols, calibrated using a single reference source and method, and processed with the same code. The AERONET-OC primary data product is the normalized water leaving radiance determined at center-wavelengths of interest for satellite ocean color applications, with an uncertainty lower than 5% in the blue-green spectral regions and higher than 8% in the red. Measurements collected at 6 sites counting the northern Adriatic Sea, the Baltic Proper, the Gulf of Finland, the Persian Gulf, and, the northern and southern margins of the Middle Atlantic Bay, have shown the capability of producing quality assured data over a wide range of bio-optical conditions including Case-2 yellow substance- and sedimentdominated waters. This work briefly introduces network elements like: deployment sites, measurement method, instrument calibration, processing scheme, quality-assurance, uncertainties, data archive and products accessibility. Emphases is given to those elements which underline the network strengths (i.e., mostly standardization of any network element) and its weaknesses (i.e., the use of consolidated, but old-fashioned technology). The work also addresses the application of AERONET-OC data to the validation of primary satellite radiometric products over a variety of complex coastal waters and finally provides elements for the identification of new deployment sites most suitable to support satellite ocean color missions.

Zibordi, Giuseppe↗
NASA NTRS2012

Anomaly Detection in Test Equipment via Sliding Mode Observers

Nonlinear observers were originally developed based on the ideas of variable structure control, and for the purpose of detecting disturbances in complex systems. In this anomaly detection application, these observers were designed for estimating the distributed state of fluid flow in a pipe described by a class of advection equations. The observer algorithm uses collected data in a piping system to estimate the distributed system state (pressure and velocity along a pipe containing liquid gas propellant flow) using only boundary measurements. These estimates are then used to further estimate and localize possible anomalies such as leaks or foreign objects, and instrumentation metering problems such as incorrect flow meter orifice plate size. The observer algorithm has the following parts: a mathematical model of the fluid flow, observer control algorithm, and an anomaly identification algorithm. The main functional operation of the algorithm is in creating the sliding mode in the observer system implemented as software. Once the sliding mode starts in the system, the equivalent value of the discontinuous function in sliding mode can be obtained by filtering out the high-frequency chattering component. In control theory, "observers" are dynamic algorithms for the online estimation of the current state of a dynamic system by measurements of an output of the system. Classical linear observers can provide optimal estimates of a system state in case of uncertainty modeled by white noise. For nonlinear cases, the theory of nonlinear observers has been developed and its success is mainly due to the sliding mode approach. Using the mathematical theory of variable structure systems with sliding modes, the observer algorithm is designed in such a way that it steers the output of the model to the output of the system obtained via a variety of sensors, in spite of possible mismatches between the assumed model and actual system. The unique properties of sliding mode control allow not only control of the model internal states to the states of the real-life system, but also identification of the disturbance or anomaly that may occur.

Solano, Wanda M.↗
NASA NTRS1995

Multiparameter linear least-squares fitting to Poisson data one count at a time

A standard problem in gamma-ray astronomy data analysis is the decomposition of a set of observed counts, described by Poisson statistics, according to a given multicomponent linear model, with underlying physical count rates or fluxes which are to be estimated from the data. Despite its conceptual simplicity, the linear least-squares (LLSQ) method for solving this problem has generally been limited to situations in which the number n(sub i) of counts in each bin i is not too small, conventionally more than 5-30. It seems to be widely believed that the failure of the LLSQ method for small counts is due to the failure of the Poisson distribution to be even approximately normal for small numbers. The cause is more accurately the strong anticorrelation between the data and the wieghts w(sub i) in the weighted LLSQ method when square root of n(sub i) instead of square root of bar-n(sub i) is used to approximate the uncertainties, sigma(sub i), in the data, where bar-n(sub i) = E(n(sub i)), the expected value of N(sub i). We show in an appendix that, avoiding this approximation, the correct equations for the Poisson LLSQ (PLLSQ) problems are actually identical to those for the maximum likelihood estimate using the exact Poisson distribution. We apply the method to solve a problem in high-resolution gamma-ray spectroscopy for the JPL High-Resolution Gamma-Ray Spectrometer flown on HEAO 3. Systematic error in subtracting the strong, highly variable background encountered in the low-energy gamma-ray region can be significantly reduced by closely pairing source and background data in short segments. Significant results can be built up by weighted averaging of the net fluxes obtained from the subtraction of many individual source/background pairs. Extension of the approach to complex situations, with multiple cosmic sources and realistic background parameterizations, requires a means of efficiently fitting to data from single scans in the narrow (approximately = 1.2 keV, HEAO 3) energy channels of a Ge spectrometer, where the expected number of counts obtained per scan may be very low. Such an analysis system is discussed and compared to the method previously used.

Wheaton, Wm. A.↗
NASA NTRS2015

Comparison of Spares Logistics Analysis Techniques for Long Duration Human Spaceflight

As the durations and distances involved in human exploration missions increase, the logistics associated with the repair and maintenance becomes more challenging. Whereas the operation of the International Space Station (ISS) depends upon regular resupply from the Earth, this paradigm may not be feasible for future missions. Longer mission durations result in higher probabilities of component failures as well as higher uncertainty regarding which components may fail, and longer distances from Earth increase the cost of resupply as well as the speed at which the crew can abort to Earth in the event of an emergency. As such, mission development efforts must take into account the logistics requirements associated with maintenance and spares. Accurate prediction of the spare parts demand for a given mission plan and how that demand changes as a result of changes to the system architecture enables full consideration of the lifecycle cost associated with different options. In this paper, we utilize a range of analysis techniques - Monte Carlo, semi-Markov, binomial, and heuristic - to examine the relationship between the mass of spares and probability of loss of function related to the Carbon Dioxide Removal System (CRS) for a notional, simplified mission profile. The Exploration Maintainability Analysis Tool (EMAT), developed at NASA Langley Research Center, is utilized for the Monte Carlo analysis. We discuss the implications of these results and the features and drawbacks of each method. In particular, we identify the limitations of heuristic methods for logistics analysis, and the additional insights provided by more in-depth techniques. We discuss the potential impact of system complexity on each technique, as well as their respective abilities to examine dynamic events. This work is the first step in an effort that will quantitatively examine how well these techniques handle increasingly more complex systems by gradually expanding the system boundary.

Owens, Andrew↗
NASA NTRSDate not supplied

A Machine Learning Approach to Improve Air Traffic Management Initiatives

Collaborating closely with commercial air carriers and related organizations, the Federal Aviation Administration(FAA) regulates air traffic and ensures the safety and efficiency of air operations. Air traffic controllers make strategic decisions, such as delaying, rerouting, or canceling flights, partly based on guidance provided by the FAA’s Air TrafficControl System Command Center (ATCSCC). The guidance includes, among other things, control measures known asTraffic Management Initiatives (TMIs) designed to enhance safety and improve operational efficiency. TMIs play a crucial role in managing the demand and capacity within the U.S. National Airspace System (NAS). Two major TMIs that are routinely used (primarily to mitigate the adverse effects of bad weather) are Ground Delay Programs (GDPs) andGround Stops (GSs). In a GDP, flights destined for airports facing thunderstorm activity experience delays at their origin airports. This proactive approach minimizes the risk of routing aircraft through hazardous weather conditions and also replaces (fuel burning) airborne delays with ground delays. In a GS, a temporary restriction is imposed on the departure or arrival of aircraft at a specific airport or within a designated airspace. Although other TMIs (e.g., miles-in-trail) are also implemented as part of (air) traffic flow management in the NAS, the focus of this work is on GDPs and GSs. Since TMIs, by design, lead to flight delays or cancellations, it is crucial to put in place the right set of parameters(e.g., scope and duration of the GDP). For example, when the end time of a GDP extends beyond what is necessary, it imposes unnecessary delays on departing flights. This situation could occur as a result of inaccurate prediction of the(required) duration of the GDP based on the weather forecast. On the other hand, if a GDP ends prematurely before the underlying capacity constraints are resolved at the destination airport, it may result in airborne holding. The delicate balance lies in matching the termination of the GDP precisely with the resolution of capacity constraints, avoiding both the imposition of unnecessary ground delays and the need for airborne holding due to premature program termination.Failing to specify the right parameters for TMIs also leads to flight delays, creating a significant obstacle in managing the increasing traffic volumes causing increased work load for the controllers. To address this issue, we propose the integration of Machine Learning (ML) models in the traffic flow management(TFM) pipeline. In current operations, decisions are made by human experts based on extensive training, historical patterns, available traffic and weather data. Since we have an abundance of data from past events that tell us the likely impact of various TMIs, by ingesting historical data, properly trained ML models can offer valuable insights and aid human decision-making. With the FAA increasingly exploring advanced analytics, ML emerges as a focal point for enhancing TFM within the National Airspace System (NAS). As a first step, this study aims to provide traffic controllers with decision-making support for the issuance and adjustment of TMIs. Data analytics and machine learning have been previously employed to address some of the challenges associated with TMIs. Numerous studies have concentrated on various facets of TMI issuance, exploring factors influencing TMI parameters, including arrival rate, airport capacity, and delay prediction. For example, using weather forecasts, several statistical methods were used to produce probabilistic capacity profiles which in conjunction with deterministic models provided insights into the GDP planning process [1–4]. The downside of using deterministic models is that they rely on fixed inputs and predetermined rules, which lack the ability to account for the inherent uncertainty and variability present in real-world scenarios. In a separate series of studies, researchers aimed to predict the occurrences of GDPs and GSs. The majority of these studies utilized various supervised learning methods, including Decision Trees, Naive Bayes, Support VectorMachines, and Random Forests to analyze the influence of weather conditions and arrival demand on TMI incidents[5–8]. However, these studies primarily focused on predicting the incidence of TMIs without explicitly addressing the scope of TMIs, including their duration and their geographical coverage. Furthermore, the emphasis of these studies was largely on GDPs, given their higher frequency and longer duration when compared to GSs. A limited number of studies focused on predicting the parameters of TMIs, specifically addressing their duration and extent. In one such study focusing on optimizing the TMI parameters at San Francisco International Airport (SFO),the authors utilized a probabilistic forecast of fog [9]. They simulated various capacity scenarios based on the (fog)burn-off forecasts, selecting GDP parameters that minimized airborne and overall ground delays. However, this approach exclusively emphasizes stratus (fog) burn-off as the primary determinant of GDP and GS, neglecting other influential factors like severe weather events, runway closures, lower capacity than traffic demand, and other important variables. Given the complexity of predicting the TMI and determining its scope, we seek a more holistic approach. We aim to consider all significant factors that could impact TMIs and their parameters. What sets this research apart is the fusion of all data sources relevant to the issuance and adjustment of TMIs and it represents the first comprehensive attempt to optimize TMIs in this manner. Since this comprehensive solution involves various aspects, we break down the problem into smaller components and input all parameters into a unified model called the “TMI Adjuster”. Figure 1 shows the overall framework and the list of datasets used in each model. The objective of the TMI Adjuster module is to deliver reliable, consistent and expedited recommendations for the progression, adjustment, and termination of TMIs. The ML solution entails developing a pipeline capable of predicting the necessity of a TMI (e.g., GS or GDP) along with its various parameters. For example, in the case of a GS, this includes the scope of the GS either in terms of distance from the destination airport or based on pre-defined airspace sectors. Here, scope refers to those regions and departing airports that are subject to the GS. In this paper, we concentrate on the issuance of GSs in the three major airports in the New York area — LaGuardia(LGA), John F. Kennedy International (JFK), and Newark Liberty International (EWR). We fuse traffic, weather and other relevant aviation data from years 2017 to 2019 to train and validate the ML models. In particular, we use the following datasets: •Terminal Aerodrome Forecast (TAF): meteorological forecasts specific to each airport, issued four times a day, covering predefined time periods. •TMI data: includes all GSs and GDPs along with their respective parameters. •Aviation System Performance Metrics (ASPM): includes traffic related data such as aircraft delays, arrival, and departure rates. •Notices to Airmen (NOTAMs): utilized to extract runway closure data and manage interdependencies between terminals in close proximity. •Flight cancellation data •Airspace Flow Programs (AFP): includes information on flight airborne holdings caused by TMIs. The data preprocessing entails transforming ASPM, TMI, AFP, NOTAMs, and weather data into an hourly format and consolidating all datasets by merging them based on date and time as the primary key. The TMI Adjuster framework comprises two parallel models: one dedicated to GS and a second model focused on GDP. As previously mentioned, our specific focus is on the GS model as a multi-classification problem. In this framework, each data point of the GS model input summarizes ten hours of data. Specifically, the data loader for the GS model generates the input and output of the model as follows: at a given time step, the input includes the actual traffic, weather, and TMI data from the two-hour window before the time step, alongside the weather forecast and scheduled traffic for the next 8 hours starting from the time step. Based on this information, the output of the GS model for each time interval consists of three dimensions. The first dimension represents a binary decision on whether there should be a GS in place for the next hour or not. The second dimension is related to the scope of the GS in the United States, and the third dimension is related to the scope of the GS in Canada (i.e., to determine if the GS impacts airports in Canada).One of the challenges with TMI modeling is the sparsity of TMI events, particularly regarding its scope. To address this challenge in the scope of the GS model output, we implement grouping. The GS scope for the US region is defined based on a list of centers that should be included when the GS is in place. With 20 centers in the US, we utilized historical data to group them into 4 categories. In particular, we summarized our historical data in a graph format where nodes represent centers, and link weights are defined based on the co-occurrence of centers in the scope parameter ofTMIs. By identified strongly connected components in this graph, we were able to partition the centers into four groups. We consider two model structures for the GS Model. Firstly, a hierarchical classification model [10], where the human decision-making for a GS is of hierarchical nature. The decision-maker first decides whether there is a need fora GS, and if the answer is yes, determines the scope. A hierarchical classification model organizes the problem into a class hierarchy, typically a tree or a Directed Acyclic Graph (DAG) structure, and considers the dependency of the decision in the previous step to the next component [10]. Here, we employ the local classifier per level approach, which involves training one multi-class classifier for each level of the class hierarchy. The second structure is the independent structure. In this setting, as the name suggests, we do not consider the dependency of the decisions in the different dimensions of the output of the model. Instead, for each dimension, we train a multi-class classifier independently. Table 1 summarizes GS model statistics for training, validation and testing. The table documents the effect of limiting data to the time steps when there was actually a TMI in place or when a TMI had just terminated. This resulted in a more balanced distribution of the GS class(GS positive class)versus “No GS”(GS negative class), which might help the training process. While JFK and LGA follow very similar distributions, with 40% and 42% GS positive class respectively, EWR has proportionally fewer GS incidents at 28%. Our subsequent phase involves evaluating the performance of both hierarchical structure and independent structure using different state-of-the-art multi-class classifier models such as Random Forest, Decision Trees, K-nearest Neighbors, and Logistic Regression and forecast the duration and scope of the GSs.

Farzan Masrour Shalmani↗
NASA NTRS2016

Predicting Real-Time Safety of the National Airspace System

To gain the situational awareness necessary for informed decision making regarding avoidance of airspace hazards, each operator must consolidate operations-relevant information from disparate sources and apply extensive domain knowledge to correctly interpret not just the current state of the NAS but forecast its (combined) evolution over the duration of the operation. This time- and workload-intensive process is periodically repeated throughout the operation so that changes can be managed in a timely manner.The imprecision, inaccuracies, inconsistency, and incompleteness of the incoming data further challenges the process. To facilitate informed decision making, this paper presents a model-based framework for the textitautomated real-time monitoring and prediction of possible effects of airspace hazards on the safety of the National Airspace System (NAS). First, hazards to flight are identified and transformed into sms, that is, quantities of interest that could be evaluated based on available data and are predictive of an unsafe event. The sms and associated thresholds that specify when an event transitions from emphsafe to emphunsafe are combined with models of airspace operations and aircraft dynamics. The framework can include any hazard to flight that can be modeled quantitatively. Models can be detailed and complex, or they can be considerably simplifed, as appropriate to the application. Real-time NAS safety monitoring and prediction begins with an estimate of the state of the NAS using the dynamic models. Given the state estimate and a probability distribution of future inputs to the NAS, we can then predict the evolution of the NAS - the future state - and the occurrence of hazards and unsafe events. The entire probability distribution of airspace sms is computed, not just point estimates, without significant assumptions regarding the distribution type andor parameters. We demonstrate our overall approach through a simulated scenario in which we predict the occurrence of some unsafe events and show how these predictions evolve in time as flight operations progress. Predictions accounting for common sources of uncertainty are included and it is shown how the predictions improve in time, become more confident, and change dynamically as new information is made available to the prediction algorithm.

safety↗
NASA NTRS2018

Characterizing Differences in the Aerosol Plume and Cloud Structure over Ascension Island During the 2016 and 2017 Biomass Burning Seasons

Marine boundary layer clouds, including the transition from stratocumulus to cumulus, are poorly represented in numerical weather prediction and general circulation models. In many cases, the complex physical relationships between marine boundary cloud morphology and the environmental conditions in which the clouds exist are not well understood. Such uncertainties arise in the presence of biomass burning carbonaceous aerosol, as is the case over the southeast Atlantic Ocean. It is likely that the absorbing and heating properties of these aerosols influence the microphysical composition and macrophysical arrangement of marine stratocumulus and trade cumulus in this region; however, this has yet to be quantified. The deployment of the Atmospheric Radiation Measurement Mobile Facility #1 (AMF1) in support of LASIC (Layered Atlantic Smoke Interactions with Clouds) provided a unique opportunity to collect observations of cloud and aerosol properties during two consecutive biomass burning seasons during July through October of 2016 and 2017 over Ascension Island (7.96 S, 14.35 W). Through the use of AMF1 observations, the Modern Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), and back trajectories from the Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT), it will be demonstrated that differences in the atmospheric circulation during the two years result in varying aerosol conditions over Ascension Island. When the aerosol plume is overhead, the aerosol loading is higher during the 2016 season as a result of a weaker subtropical high-pressure system. Furthermore, the aerosol plume originates from central Africa in 2016, but further south in 2017. Contrasts in the season-to-season and day-to-day aerosol loading are used to categorize boundary layer cloud and sub-cloud turbulence measurements above Ascension Island using the AMF1 Doppler lidar and cloud radar.

Collow, Allison B. Marquardt↗
NASA NTRSDate not supplied

How Well Satellite Remote Sensing Can Inform Surface Level Ozone Production Sensitivity to Precursor Trace Gas Emissions?

Surface-level ozone (O3) pollution mitigation strategies rely on a clear understanding of the atmospheric chemical processes involving O3, NOx (nitric oxide (NO) + nitrogen dioxide (NO2)) and volatile organic compounds (VOCs). A robust spatiotemporal classification of O3 photochemical regimes with relative chemical sensitivity of local O3 formation to emission reductions of NOx (NOx-limited regime) versus VOCs (radical-limited regime) is required to design anthropogenic emission reduction policies in order to improve surface air quality. Numerous studies found that the ratio of satellite-retrieved Vertical Column Densities (VCDs) of formaldehyde (HCHO) to NO2 can diagnose O3 production sensitivity. However, uncertainty still remains about the accuracy of such estimates using satellite-based retrievals. This study conducts an investigation to identify how well satellite data can inform surface-level O3 sensitivity by validating satellite-based estimates against aircraft-based retrievals taken during the Long Island Sound Tropospheric Ozone Study (LISTOS-2018) and Ozone Water-Land Environmental Transition Study (OWLETS-2) field campaigns. For this study, we use HCHO and NO2 retrievals from the Ozone Monitoring Instrument (OMI) onboard the Aura satellite, and high spatiotemporal resolution TROPOspheric Monitoring Instrument (TROPOMI) onboard the Sentinel-5 Precursor (S5P) satellite. The Community Multiscale Air Quality (CMAQ) modelling system is used to provide trace gas vertical profiles to recalculate air mass factor for deriving consistent satellite-based VCDs. To demonstrate the ability of satellites to diagnose O3 production regimes, we validate satellite-derived O3 sensitivity regimes, along with the OMI and TROPOMI products applied, using HCHO and NO2 retrievals from the aircraft-based instruments: Geostationary Trace gas and Aerosol Sensor Optimization (GeoTASO) spectrometer and Geostationary Coastal and Air Pollution Events Airborne Simulator (GCAS) instruments deployed during the LISTOS-2018 and OWLETS-2 campaigns. The results of this work will advance our fundamental knowledge of the spatiotemporal evolution of the complex O3-NOx-VOC chemistry, and will advance our understanding about the quality of remote-sensing, suborbital, and satellite data for accurately observing surface-level O3 production sensitivity.

Satellite↗
NASA NTRS2022

Evaluation of IMERG Satellite Precipitation Over the Land-Coast-Ocean Continuum – Part Ii: Quantification

To understand and manage water systems under a changing climate and meet an increasing demand for water, a quantitative understanding of precipitation is most important in coastal regions. The capabilities of the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG) V06B product for precipitation quantification are examined over three coastal regions of the United States: the West Coast, the Gulf of Mexico, and the East Coast, all of which are characterized by different topographies and precipitation climatologies. A novel uncertainty analysis of IMERG is proposed that considers environmental and physical parameters such as elevation and distance to the coastline. The IMERG performance is traced back to its components, i.e., passive microwave (PMW), infrared (IR), and morphing-based estimates. The analysis is performed using high-resolution, high-quality Ground Validation Multi-Radar/Multi-Sensor (GV-MRMS) rainfall estimates as ground reference at the native resolution of IMERG of 30 min and 0.1°. IMERG Final (IM-F) quantification performance heavily depends on the respective contribution of PMW, IR, and morph components. IM-F and its components overestimate the contribution of light rainfall (<1 mm h−1) and underestimate the contribution of high rainfall rates (>10 mm h−1) to the total rainfall volume. Strong regional dependencies are highlighted, especially over the West Coast, where the proximity of complex terrain to the coastline challenges precipitation estimates. Other major drivers are the distance from the coastline, elevation, and precipitation types, especially over the land and coast surface types, that highlight the impact of precipitation regimes.

Atmosphere↗
NASA NTRS2024

Pliocene Model Intercomparison Project Phase 3 (PlioMIP3) – Science Plan and Experimental Design

The Pliocene Model Intercomparison Project (PlioMIP) was initiated in 2008. Over two phases PlioMIP has helped co-ordinate the experimental design and publication strategy of the community, which has included an increasing number of climate models and modelling groups from around the world. It has engaged with palaeoenvironmental scientists to foster new data synthesis supporting the construction of new model boundary conditions, as well as to facilitate new data-model comparisons. The work has advanced our understanding of Pliocene climates and environments, enhanced our knowledge regarding the ability of complex climate and Earth System models to accurately simulate climate change, and helped to refine our estimates of how sensitive the climate system is to forcing conditions. In this community protocol paper, we outline the scientific plan for PlioMIP Phase 3 (PlioMIP3). This plan provides the required guidance to participating modelling groups from around the world to successfully set up and perform PlioMIP3 climate model experiments. The project is open to new participants from the scientific community (both from the climate modelling and geosciences communities). In PlioMIP3, we retain the PlioMIP2 Core experiments (Eoi 400 , E 280 ) and extend the Core requirements to include either an experiment focussed on the Early Pliocene or an alternative Late Pliocene simulation (or both). These additions (a) allow a comparison of Early and Late Pliocene warm intervals and help build research connections and synergy with the MioMIP (Miocene Model Intercomparison Project - also known as DeepMIP-Miocene) and PlioMioVAR projects (Pliocene-Miocene Variability Working Group), and (b) create an alternative time slice simulation for 3.205 Ma (MIS KM5c) through removal of some of the largest palaeogeographic differences introduced between PlioMIP1 and 2 resulting in minimal land-sea mask variations from the modern. In addition, we present ten optional experiments designed to enhance our assessment of climate sensitivity and to explore the uncertainty in greenhouse gas-related forcing. For the first time, we introduce orbital sensitivity experiments into the science plan, as well as simulations incorporating dynamic vegetation-climate feedbacks and an experiment designed to examine the potential significance of East Antarctic Ice Sheet boundary condition uncertainty. These changes enhance palaeo-to-future scientific connections and enable an exploration of the significance of palaeogeographic uncertainties on climate simulations.

Climate modelling↗