Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “deviated”

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 775 records · Page 43

Validation of Special Sensor Microwave Imager Monthly-Mean Wind Speed From July 1987 to December 1989

Since July 1987 the Remote Sensing Systems, Santa Rosa, California, routinely computes from first principles the wind speed from Special Sensor Microwave Imager (SSM/I) measurements of the intensity of microwave radiation emitted at the ocean surface. The accuracy of monthly-mean SSM/I wind speeds is determined by comparisons with moored-buoy wind measurements, which were recorded by four different institutions at 44 sites in the Atlantic and Pacific Oceans during July 1987 to December 1989. All results for 1988 were virtually identical with 1989. The range of monthly mean moored-buoy wind speeds was 2-10 m s^(-1). During 1987, the equatorial matchups were not equivalent with 1998 and 1989, and the cause remains unknown. The root-mean-square (rms) difference of 697 monthly-mean matchups of the composite 1988 and 1989 data set was 1.2 m s^(-1). The rms differences were smaller in the equatorial zone and higher in middle latitudes, in accord with the monthly standard deviation. At middle latitudes the time series of rms differences displayed an annual cycle. In the equatorial zone the agreement between SSM/I and in situ data was better in regions with a lesser amount of clouds, and vice versa. For SSM/I monthly standard deviations of 1-2, 2-3, and 3-4 m s^(-1), the average absolute values of the monthly-mean difference between SSM/I and moored-buoy wind speeds were 0.6, 0.9, and 1.4 m s^(-1), respectively.

Halpern, D.↗

Radiography in High Mass X-Ray Binaries: Micro-Structure of the Stellar Wind Through Variability of the Column Density

Context. In high mass X-ray binaries, an accreting compact object orbits a high mass star, which loses mass through a dense and in homogeneous wind. Aims. Using the compact object as an X-ray backlight, the time variability of the absorbing column density in the wind can be exploited in order to shed light on the micro-structure of the wind and obtain unbiased stellar mass-loss rates for high mass stars. Methods. We developed a simplified representation of the stellar wind where all the matter is gathered in spherical “clumps” that are radially advected away from the star. This model enables us to explore the connections between the stochastic properties of the wind and the variability of the column density for a comprehensive set of parameters related to the orbit and to the wind micro-structure, such as the size of the clumps and their individual mass. In particular, we focus on the evolution with the orbital phase of the standard deviation of the column density and of the characteristic duration of enhanced absorption episodes. Using the porosity length, we derive analytical predictions and compare them to the standard deviations and coherence time scales that were obtained. Results. We identified the favorable systems and orbital phases to determine the wind micro-structure. The coherence time scale of the column density is shown to be the self-crossing time of a single clump in front of the compact object. We thus provide a procedure to get accurate measurements of the size and of the mass of the clumps, purely based on the observable time variability of the column density. Conclusions. The coherence time scale grants direct access to the size of the clumps, while their mass can be deduced separately from the amplitude of the variability. We further show how monitoring the variability at superior conjunctions can probe the onset of the clump-forming region above the stellar photosphere. If the high column density variations in some high mass X-ray binaries are due to unaccreted clumps which are passing by the line-of-sight, this would require high mass clumps to reproduce the observed peak-to-peak amplitude and coherence time scales. These clump properties are marginally compatible with the ones derived from radiative-hydrodynamics simulations. Alternatively, the following components could contribute to the variability of the column density: larger orbital scale structures produced by a mechanism that has yet to be identified or a dense environment in the immediate vicinity of the accretor, such as an accretion disk, an outflow, or a spherical shell surrounding the magnetosphere of the accreting neutron star.

stars: mass loss↗

Development of the WFIRST CGI Integral Field Spectrograph

The principal science goal of directly imaging exoplanets is spectroscopy. Future mission concepts such as LU- VOIR and HabEx have baselined their science instruments to use an integral field spectrograph (IFS) to maximize observing efficiency. Coronagraphic imaging with an IFS has strong heritage at ground observatories, having produced planet and brown dwarf near-infrared spectra at contrast floors as low as high 1 × 10-6 levels. As a technology demonstration for future missions the WFIRST coronagraph instrument (CGI) will demonstrate the ability to measure spectra in visible light with at contrast levels sufficient for detecting Earth-like planets. The spectrograph was originally meant to match that of future missions; an independent IFS channel for wavefront control and imaging spectroscopy. Due to a series of de-scopes as the project progressed, the baseline spectroscopic capability of CGI is no longer an IFS but a zero deviation optical prism and slit used in conjunction with the direct imaging camera. Regardless, CGI advanced a lenslet IFS through the preliminary design phase of a flight project. Here we present that development as the challenges and requirements of an IFS for coronagraphy are quite unique and relevant to future missions. With that development in place, we also present the initial design of the zero optical deviation prism, highlighting the relative challenges and performance of each concept.

Tyler D Groff↗

Ephemeris refinement of 21 hot Jupiter exoplanets with high timing uncertainties

Transit events of extrasolar planets offer a wealth of information for planetary characterization. However, for many known targets, the uncertainty of their predicted transit windows prohibits an accurate scheduling of follow-up observations. In this work, we refine the ephemerides of 21 hot Jupiter exoplanets with the largest timing uncertainties. We collected 120 professional and amateur transit light curves of the targets of interest, observed with a range of telescopes of 0.3 m–2.2 m, and analyzed them along with the timing information of the planets discovery papers. In the case of WASP-117b, we measured a timing deviation compared to the known ephemeris of about 3.5 h, and for HAT-P-29b and HAT-P-31bthe deviation amounted to about 2 h and more. For all targets, the new ephemeris predicts transit timings with uncertainties of less than 6 min in the year 2018 and less than 13 min until 2025. Thus, our results allow for an accurate scheduling of follow-up observations in the next decade.

photometric↗

Precipitable water vapor over oceans from the Maritime Aerosol Network: Evaluation of global models and satellite products under clear sky conditions

We present results from an evaluation of precipitable water vapor (W) over remote oceanic areas as derived from global reanalysis models and from satellites against observations from the Maritime Aerosol Network (MAN) for cloudless skies during the period of 2004–2017. They cover polar, mid latitude and tropical oceanic regions and represent a first effort to use MAN observations for such evaluation. The global reanalysis model products evaluated in this study are from the Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA-2), the European Centre for Medium-Range Weather Forecasts (ECMWF) Interim Reanalysis (ERA I), and the Climate Forecast System Reanalysis (CFSR) model. The satellite products evaluated are from the Moderate Resolution Imaging Spectroradiometer (MODIS), the Polarization and Directionality of the Earth's Reflectances (POLDER), the Global Ozone Monitoring Experiment (GOME-2), the Scanning Imaging Absorption Spectrometer for Atmospheric Chartography (SCIAMACHY), and the Atmospheric Infra-red Sounder (AIRS). Satellite retrievals of W are based on the attenuation of solar reflected light by water vapor absorption bands, except those from AIRS that rely on brightness temperature measurements. A very good agreement is observed between the model estimates and MAN, with mean differences of ~5% and standard deviations of ~15%. These results are within the uncertainties associated with the models and the measurements, indicating the skill of the reanalysis models to estimate W over oceans under clear sky conditions. Mean differences of W between the satellite and MAN products are ~11, 6.7, 12, −7, and 3% for MODIS, POLDER, GOME-2, SCIAMACHY and AIRS respectively, while their standard deviations are 31, 29, 28, 20 and 17%. These differences reveal the need to address inconsistencies among different satellite sensors and ground-based measurements to reduce the uncertainties associated with the retrievals.

Daniel Pérez-Ramírez↗

Tropospheric Ozone Assessment Report: Tropospheric Ozone From 1877 to 2016, Observed Levels, Trends and Uncertainties

From the earliest observations of ozone in the lower atmosphere in the 19th century, both measurement methods and the portion of the globe observed have evolved and changed. These methods have different uncertainties and biases, and the data records differ with respect to coverage (space and time), information content, and representativeness. In this study, various ozone measurement methods and ozone datasets are reviewed and selected for inclusion in the historical record of background ozone levels, based on relationship of the measurement technique to the modern UV absorption standard, absence of interfering pollutants, representativeness of the well-mixed boundary layer and expert judgement of their credibility. There are significant uncertainties with the 19th and early 20th-century measurements related to interference of other gases. Spectroscopic methods applied before 1960 have likely underestimated ozone by as much as 11% at the surface and by about 24% in the free troposphere, due to the use of differing ozone absorption coefficients. There is no unambiguous evidence in the measurement record back to 1896 that typical mid-latitude background surface ozone values were below about 20 nmol mol–1, but there is robust evidence for increases in the temperate and polar regions of the northern hemisphere of 30–70%, with large uncertainty, between the period of historic observations, 1896–1975, and the modern period (1990–2014). Independent historical observations from balloons and aircraft indicate similar changes in the free troposphere. Changes in the southern hemisphere are much less. Regional representativeness of the available observations remains a potential source of large errors, which are difficult to quantify. The great majority of validation and intercomparison studies of free tropospheric ozone measurement methods use ECC ozonesondes as reference. Compared to UV-absorption measurements they show a modest (~1–5% ±5%) high bias in the troposphere, but no evidence of a change with time. Umkehr, lidar, and FTIR methods all show modest low biases relative to ECCs, and so, using ECC sondes as a transfer standard, all appear to agree to within one standard deviation with the modern UV-absorption standard. Other sonde types show an increase of 5–20% in sensitivity to tropospheric ozone from 1970–1995. Biases and standard deviations of satellite retrieval comparisons are often 2–3 times larger than those of other free tropospheric measurements. The lack of information on temporal changes of bias for satellite measurements of tropospheric ozone is an area of concern for long-term trend studies.

UV absorption standard↗

Uncertainty Estimates for Sonic-Boom Pressure Signatures and Loudness Carpets

A non-intrusive uncertainty quantification method is applied to computational analysis of supersonic, low-boom aircraft. The mean and standard deviation statistics of the pressure waveforms and loudness metrics are evaluated through use of numerical quadrature. The probability density function (p.d.f.) of these outputs is evaluated via kernel density estimation. The simulations use an inviscid, embedded-boundary Cartesian-mesh flow solver in the nearfield combined with an augmented Burgers’ equation solver for propagation in the farfield. The results show that the p.d.f. of the waveform is bimodal at shocks, which makes the mean and standard deviation statistics inappropriate. Despite this limitation, we show that the moment statistics can provide effective assessment of discrepancies when comparing with experimental data. This is demonstrated by presenting uncertainty analysis of a wind-tunnel test and showing that we significantly improve the predictions when we include the test uncertainties in the simulation. Normal distributions are obtained for the ground signature and loudness metrics, which is primarily due to the careful shaping of the low-boom waveform. Separation of variables and error control are used to reduce computational cost. We demonstrate that this is an efficient approach in the sense of balancing numerical errors in the statistics quadrature with discretization errors in the solvers.

ARMD↗

Quantifying the Observational Requirements of a Space-borne LiDAR Snow Mission

This study quantifies the level of observational accuracy required from a spaceborne light detection and ranging (LiDAR) snow depth retrieval mission for enabling beneficial impacts for snow estimation. The study is conducted over a region in Western Colorado using a suite of observing system simulation experiments(OSSEs).The Joint UK Land Environment Simulator, version 5.0 (JULES v5.0) is employed to simulate a suite of idealized LiDAR observations, considering a range of LiDAR snow depth retrieval errors, different hypothetical sensor swath widths, and the impact of cloud cover on observability. These simulated observations are then assimilated into the Noah land surface modelwith multi-parameterization options, version 3.6 (Noah-MP v3.6) model. This data assimilation setup is used to systematically evaluate the potential utility of LiDAR observations for improving modeled snow water equivalent (SWE)estimates and water budget variables such as runoff. Results from the OSSE runsshow that, in general, assimilation of synthetic LiDAR observations provide beneficial impacts when theLiDAR snow depth retrieval error standard deviation (σerror) is below 60 cm.Based on comparisons between the realistic (i.e., swath-limited and cloud-attenuated) case and the idealized (i.e., infinite swath width in the absence of cloud cover) case,this study concludes that observations with a conservative error standard deviation threshold of 40 cm (i.e., upper limit of the snow depth retrieval error that adds value to the SWE estimates viaassimilation) are needed for improving modeled snow estimates. More than a 33% reduction in SWEroot mean square errors and more than a 15% increase in correlation coefficientsare achieved when σ error ≤ 40 cmusing a 170-km sensor swath width in the presence of cloud attenuation effects. Further, the integrated hydrologic response, as represented by total (surface and subsurface) runoff estimates during the snow ablation season, are also enhanced when assimilating synthetic LiDAR snow depth retrievals with errors below this level.

Yonghwan Kwon↗

A Study of Precipitation Mean Particle Size with the GPM Combined Radar-Radiometer Algorithm

The Global Precipitation Measurement (GPM) mission provides quantitative estimates of the mass-weighted mean diameter, Dm, a characteristic parameter of the precipitation particle size distribution (PSD). This study investigates the three-dimensional distribution of Dm retrieved by the GPM Combined Radar-Radiometer Algorithm (CORRA) using observations from GPM’s dual-frequency precipitation radar and microwave imager. We examine Dm in individual precipitation systems in different climate regimes and investigate a six-year global climatology. The vertical structure of Dm is demonstrated with cases of deep convection, frontal rain and snow, and stratocumulus light rain. The Dm values, detectable by GPM, range from appx 0.7 mm in stratocumulus precipitation to >3.5 mm in the ice layers of intense convection. Within the constraint of the 12-dBZ detectability threshold, the smallest climatological mean Dm (~ 0.8 mm) are found in the eastern oceans, and the largest values (~ 2 mm) occur above the melting levels in convection over land in summer. The standard deviation is generally < 0.45 mm below 6 km. Climate regimes are characterized with Dm annual/seasonal variations, its convective/stratiform components, and vertical variabilities (2-10 km). The US Central Plains and Argentina are associated with the largest Dm in a deep layer. Tropical Africa has larger Dm and standard deviation than Amazon. Large convective Dm occurs at high latitudes of Eurasia and North America in summer; the southern hemisphere high latitudes have shallower systems with smaller Dm. Oceanic storm tracks in both hemispheres have relatively large Dm, particularly for convective Dm in winter.

Mei Han↗

Future Climate Change Significantly Alters Interannual Wheat Yield Variability Over Half of Harvested Areas

Climate change affects the spatial and temporal distribution of crop yields, which can critically impair food security across scales. A number of previous studies have assessed the impact of climate change on mean crop yield and future food availability, but much less is known about potential future changes in interannual yield variability. Here, we evaluate future changes in relative interannual global wheat yield variability (the coefficient of variation (CV)) at 0.25° spatial resolution for two representative concentration pathways (RCP4.5 and RCP8.5). A multi-model ensemble of crop model emulators based on global process-based models is used to evaluate responses to changes in temperature, precipitation, and CO2. The results indicate that over 60% of harvested areas could experience significant changes in interannual yield variability under a high-emission scenario by the end of the 21st century (2066–2095). About 31% and 44% of harvested areas are projected to undergo significant reductions of relative yield variability under RCP4.5 and RCP8.5, respectively. In turn, wheat yield is projected to become more unstable across 23% (RCP4.5) and 18% (RCP8.5) of global harvested areas—mostly in hot or low fertilizer input regions, including some of the major breadbasket countries. The major driver of increasing yield CV change is the increase in yield standard deviation, whereas declining yield CV is mostly caused by stronger increases in mean yield than in the standard deviation. Changes in temperature are the dominant cause of change in wheat yield CVs, having a greater influence than changes in precipitation in 53% and 72% of global harvested areas by the end of the century under RCP4.5 and RCP8.5, respectively. This research highlights the potential challenges posed by increased yield variability and the need for tailored regional adaptation strategies.

Climate change↗

KMT-2019-BLG-1715: Planetary Microlensing Event with Three Lens Masses and Two Source Stars

We investigate the gravitational microlensing event KMT-2019-BLG-1715, the light curve of which shows two short-term anomalies from a caustic-crossing binary-lensing light curve: one with a large deviation and the other with a small deviation. We identify five pairs of solutions, in which the anomalies are explained by adding an extra lens or source component in addition to the base binary-lens model. We resolve the degeneracies by applying a method in which the measured flux ratio between the first and second source stars is compared with the flux ratio deduced from the ratio of the source radii. Applying this method leaves a single pair of viable solutions, in both of which the major anomaly is generated by a planetary-mass third body of the lens, and the minor anomaly is generated by a faint second source. A Bayesian analysis indicates that the lens comprises three masses: a planet-mass object with ∼2.6 MJ and binary stars of K and M dwarfs lying in the galactic disk. We point out the possibility that the lens is the blend, and this can be verified by conducting high-resolution follow-up imaging for the resolution of the lens from the source.

Cheongho Han↗

Understanding the Global Three-Dimensional Distribution of Precipitation Mean Particle Size with the Global Precipitation Measurement Mission

This study addresses the global distribution of precipitation mean particle size using data from the Global Precipitation Measurement (GPM) mission. The mass-weighted mean diameter Dm is a characteristic parameter of the precipitation particle size distribution (PSD), estimated from the GPM Combined Radar–Radiometer Algorithm (CORRA) using data from GPM’s dual-frequency precipitation radar and microwave imager. We examine Dm in individual precipitation systems in different climate regimes and investigate a 6-yr (2014–20) global climatology within 70°N–70°S. The vertical structure of Dm is demonstrated with cases of deep convection, frontal rain and snow, and stratocumulus light rain. The Dm values, detectable by GPM, range from ~0.7 mm in stratocumulus precipitation to >3.5 mm in the ice layers of intense convection. Within the constraint of the 12-dBZ detectability threshold, the smallest annual mean Dm (~0.8 mm) are found in the eastern oceans, and the largest values (~2 mm) occur above the melting levels in convection over land in summer. The standard deviation of the annual mean is generally <0.45 mm below 6 km. Climate regimes are characterized with Dm annual/seasonal variations, its convective/stratiform components, and vertical variabilities (2–10 km). The U.S. Central Plains and Argentina are associated with the largest Dm in a deep layer. Tropical Africa has larger Dm and standard deviation than Amazon. Large convective Dm occurs at high latitudes of Eurasia and North America in summer; the Southern Hemisphere high latitudes have shallower systems with smaller Dm. Oceanic storm tracks in both hemispheres have relatively large Dm, particularly for convective Dm in winter. Relatively small Dm occurs over tropical oceans, including ITCZ, requiring further investigation.

Mei Han↗

VESIcal: An Open-source Thermodynamic Model Engine for Mixed Volatile (H2O-CO2) Solubility in Silicate Melts

Modeling the solubility of volatiles in silicate melts is fundamental to the interpretation of volcanic systems and has implications for magma dynamics, eruption style, and material transport between the mantle, crust, and atmosphere. Recent advancements in computational capabilities and access to computing tools has outpaced the functionality and extensibility of previously available modeling platforms. Here we present VESIcal (Volatile Equilibria and Saturation Index calculator), the first comprehensive modeling tool for H2O, CO2, and mixed (H2O-CO2) solubility in silicate melts that: a) allows users access to seven popular models, with easy inter-comparison between models; b) provides universal functionality for all models (e.g., functions for calculating saturation pressures, degassing paths, etc.); c) can process large datasets (1,000’s of samples) automatically; d) can output computed data into an Excel spreadsheet or CSV file for post-modeling analysis; e) integrates plotting capabilities directly within the tool; and f) provides all of this within the framework of a python library, making the tool extensible by the user and allowing any of the model functions to be incorporated into any other code capable of calling python.Here we will provide a demonstration of VESIcal and its capabilities with applications to various volcanic processes affected by volatiles. VESIcal represents the first tool capable of directly comparing multiple solubility models and equations of state. We find that commonly used models predict surprisingly different volatile solubilities, particularly for pure CO2 or mixed CO2-H2O fluids. Even for melt compositions that are well represented in the calibration datasets of multiple models (e.g., MORBs), calculated solubilities for pure CO2 and pure H2O can deviate from one another by factors of >2 leading to 2x deviations in calculated saturation pressures (e.g., 5 to 10 kbar). The solubility of CO2 predicted by different rhyolitic models also differs substantially, overwhelming other sources of uncertainty such as analytical errors on measurements of volatile contents or uncertainties in crustal density profiles. This highlights the importance of model choice when drawing geological conclusions based on volatiles in magmas.

Kayla Iacovino↗

Predicting Two-Dimensional Airfoil Performance Using Graph Neural Networks

Computer simulations require the use of meshes to simulate geometries. These meshes capture important geometric features of the design and can be used in machine learning modeling. This report explores the use of graph neural networks (GNNs) to learn features from two-dimensional (2D) airfoil designs represented as a set of nodes connected using edges. This type of network is common in aerospace applications: most geometries are represented as a mesh in order to perform analysis. The objective of this work is to use GNNs to predict the performance of 2D airfoils generated using the program XFOIL. The predicted performance parameters include bulk quantities such as coefficients of lift (C L ), drag (C d , C dp ), moment (C m ), and node-specific quantities such as coefficient of pressure (C p ). In this report, a spline convolutional graph-based neural network is compared with deep learning neural networks to predict both bulk and node-specific quantities. The findings indicate the GNNs are able to predict bulk quantities quite well; however, when the number of outputs is increased, the deep neural network (DNN) proves to be better in its prediction capability. Two different normalization strategies were compared in the training of both GNNs and DNNs: minmax and standard deviation. In both types of networks, standard deviation scaling proved to be the best.

machine learning↗

Impact of Recent Climate Extremes on Mosquito-borne Disease Transmission in Kenya

Climate change and variability influence temperature and rainfall, which impact vector abundance and the dynamics of vector-borne disease transmission. Climate change is projected to increase the frequency and intensity of extreme climate events. Mosquito-borne diseases, such as dengue fever, are primarily transmitted by Aedes aegypti mosquitoes. Freshwater availability and temperature affect dengue vector populations via a variety of biological processes and thus influence the ability of mosquitoes to effectively transmit disease. However, the effect of droughts, floods, heat waves, and cold waves is not well understood. Using vector, climate, and dengue disease data collected between 2013 and 2019 in Kenya, this retrospective cohort study aims to elucidate the impact of extreme rainfall and temperature on mosquito abundance and the risk of arboviral infections. To define extreme periods of rainfall and land surface temperature (LST), we calculated monthly anomalies as deviations from long-term means (1983–2019 for rainfall, 2000–2019 for LST) across four study locations in Kenya. We classified extreme climate events as the upper and lower 10% of these calculated LST or rainfall deviations. Monthly Ae. aegypti abundance was recorded in Kenya using four trapping methods. Blood samples were also collected from children with febrile illness presenting to four field sites and tested for dengue virus using an IgG enzyme-linked immunosorbent assay (ELISA) and polymerase chain reaction (PCR). We found that mosquito eggs and adults were significantly more abundant one month following an abnormally wet month. The relationship between mosquito abundance and dengue risk follows a non-linear association. Our findings suggest that early warnings and targeted interventions during periods of abnormal rainfall and temperature, especially flooding, can potentially contribute to reductions in risk of viral transmission

Aedes aegypti mosquitoes↗

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↗

M2-SCREAM: A Stratospheric Composition Reanalysis of Aura MLS Data With MERRA-2 Transport

MERRA-2 Stratospheric Composition Reanalysis of Aura Microwave Limb Sounder (M2-SCREAM) is a new reanalysis of stratospheric ozone, water vapor, hydrogen chloride (HCl), nitric acid (HNO 3 ) and nitrous oxide (N 2 O) between 2004 and the present (with a latency of several months). The assimilated fields are provided at a 50-km horizontal resolution and at a three-hourly frequency. M2-SCREAM assimilates version 4.2 Microwave Limb Sounder (MLS) profiles of the five constituents alongside total ozone column from the Ozone Monitoring Instrument. Dynamics and tropospheric water vapor are constrained by the MERRA-2 reanalysis. The assimilated species are in excellent agreement with the MLS observations, except for HNO 3 in polar night, where data are not assimilated. Comparisons against independent observations show that the reanalysis realistically captures the spatial and temporal variability of all the assimilated constituents. In particular, the standard deviations of the differences between M2-SCREAM and constituent mixing ratio data from The Atmospheric Chemistry Experiment Fourier Transform Spectrometer are much smaller than the standard deviations of the measured constituents. Evaluation of the reanalysis against aircraft data and balloon-borne frost point hygrometers indicates a faithful representation of small-scale structures in the assimilated water vapor, HNO 3 and ozone fields near the tropopause. Comparisons with independent observations and a process-based analysis of the consistency of the assimilated constituent fields with the MERRA-2 dynamics and with large-scale stratospheric processes demonstrate the utility of M2-SCREAM for scientific studies of chemical and transport variability on time scales ranging from hours to decades. Analysis uncertainties and guidelines for data usage are provided.

MERRA-2↗

Ames-2021 CO 2 Dipole Moment Surface and IR Line Lists: Toward 0.1% Uncertainty for CO 2 IR Intensities

A highly accurate CO 2 ab initio dipole moment surface (DMS), Ames-2021, is reported along with 12 C 16 O 2 infrared (IR) intensity comparisons approaching a 1–4‰ level of agreement and uncertainty. The Ames-2021 DMS was accurately fitted from CCSD(T) finite-field dipoles computed with the aug-cc-pVXZ (X = T, Q, 5) basis for C atom and the d-aug-cc-pVXZ (X = T, Q, 5) basis for O atoms, and extrapolated to the one particle basis set limit. Fitting σ(rms) is 3.8 × 10^(–7) au for 4443 geometries below 15 000 per cm . The corresponding IR intensity, S Ames-2021 , are computed using the Ames-2 potential energy surface (PES), which is the best PES available for CO 2 . Compared to high accuracy IR studies for 2001i–00001 and 3001i–00001 bands, S Ames-2021 matches NIST experiment-based intensities [S NIST-HIT16 or S HIT20 ] to −1.0 ± 1.3‰, or matches DLR experiment-based intensities [S DLR-HIT16/UCL/Ames ] to 1.9 ± 3.7‰. This indicates the systematic deviations and uncertainties have been significantly reduced in S Ames-2021 . The S UCL2015 (or S HITRAN2016 ) have larger deviations (vs S DLR ) and uncertainties (vs S DLR , S NIST ) which are attributed to the less accurate Ames-1 PES adopted in UCL-296 line list calculation. The S Ames-2021 intensity of 12 C 16 O 2 and 13 C 16 O 2 is utilized to derive new absolute 13 C/ 12 C ratios for Vienna PeeDee Belemnite (VPDB) with uncertainty reduced by 1/3 or 2/3. Further evaluation of S Ames-2021 are carried out on those CO 2 bands discussed in the HITRAN2020 update paper. Consistent improvements and better accuracies are found in band-by-band analysis, except for those bands strongly affected by Coriolis couplings, or very weak bands measured with relatively larger experimental uncertainties. The Ames-2021 296 K IR line lists are generated for 13 CO 2 isotopologues, with 18,000 per and S 296K > 1E–31) cm/molecule cutoff and then combined with CDSD line positions (except 14 C 16 O 2 ). The Ames-2021 DMS and 296 K IR line lists represent a major improvement over previous CO 2 theoretical IR intensity studies, including Ames-2016, UCL-296, and recent UCL DMS 2021 update. A real 1 permille level of agreement and uncertainty will definitely require both more accurate PES and more accurate DMS.

Chemical calculations↗