Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “large deviations”

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 757 records · Page 42

Measuring the Hubble constant and our Virgo-infall velocity independently

A sample of spiral galaxies with B(sub T) less than 14.5 located in two local volumes, one in the direction of, but behind, the Virgo Cluster (behind-Virgo volume (BV)) and the other in the opposite direction (anti-Virgo volume (AV)), were used via a Tully-Fisher (TF) relation to derive the following two parameters: H(sub AB), the mean Hubble ratio between AV and BV, and delta v(sub parallel), the peculiar velocity of the Local Group in the direction of the Virgo Cluster (VC) with respect to a uniformly expanding reference system defined by our AV and BV sub-samples. The two sampled volumes, separated by a velocity interval of 5600 km/s, form an antipodal pair. This particular geometry not only allows us to derive the two parameters independently but also reduces the dynamical effect of the Local Supercluster on H(sub AB) without increasing the Malmquist bias. By limiting our sample to spiral galaxies having large velocity widths W(sub R), we effectively reduce the TF scatter and Malmquist bias in our sample. The TF zero point and dispersion were then determined by further correcting for the small residual Malmquist bias. An additional sample of fainter galaxies was used to test for a non-Gaussian tail to the TF disperison. We found no evidence for such a tail and formally give an upper limit of about 18% for the fractional contribution of an unseen tail. The average intrinsic TF dispersion for the dominant Gaussian component is sigma(sub TF)(sup 0) approximately 0.33 mag for W(sub R) approximately equal to or greater than 180 km/s. Our numerical results are delta v(sub parallel) approximately equals 414 +/- 82 km/s and H(sub AB) approximately equals (84.0 +/- 2.4)(1 + epsilon) km/s Mpc, where (1 + epsilon) accounts for any systematic error between the calibrators and the sample galaxies. Various dynamical models were tested to explore the effect on H(sub AB) of the uncertainties in the local velocity field. Constrained by our observed delta v(sub parallel) as well as other observational quantities, we found that the rms deviation from unity of H(sub AB)/H(sub 0) (where H(sub 0) is the Hubble constant for each model) is 5%, making H(sub AB) a good indicator for H(sub 0). Taking this variation as an additional error, our formal estimate for the Hubble constant is H(sub 0) approximately equals (84 +/- 5)(1 + epsilon) km/s Mpc.

Lu, Nanyao Y.↗

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

Accurate age determinations of hyper-velocity impact and cratering events remains difficult and often controversial, while less than half of all known impact craters are regarded as accurately and precisely dated. Besides 40Ar/39Ar and U-Pb methods, zircon (U-Th)/He dating of impactites is a burgeoning technique to date large- to medium-sized impact structures. Zircon (U-Th)/He ages can be fully reset in minutes at 1000°C, T commonly reached in and directly adjacent to impact melt domains, whereas complete resetting of zircon (U-Th)/He at <300°C, which might be encountered near the crater margins or persist in post-impact hydrothermal systems, may take >103-4 years. However, there is a critical need to test the reliability of (U-Th)/He impact dating in shock deformed zircon, and to quantify helium diffusion kinetics in well-characterized grains with a broad spectrum of shock-induced defect substructures. For this purpose, we investigated samples from two impact structures, the 66 Ma Chicxulub multi-ring basin and the 15 Ma Ries complex crater, to compare zircon diffusion kinetics from impact structures with varying parameters, including size, age, and hydrothermal system longevity. Shock microstructures were characterized by backscattered-electron (BSE) imaging prior to determination of diffusion step-heating fractional release experiments using light-bulb furnace with prograde and retrograde incrementally 10°C steps from 300°C to 600°C. We find that zircon with low-level shock microstructures exhibit no significant deviation from helium diffusion kinetics of undamaged zircon. In contrast, zircon grains with planar deformation features and granular textures classified by SEM are characterized by a dramatic decrease in helium retentivity, similar to radiation damage, due to the reduction in the effective domain size and the introduction of fast diffusion pathways. This likely renders shocked grains more susceptible to impact-induced hydrothermal resetting. Hence, characterization of impact microstructure is critical for determining accurate impact ages, but also offers the opportunity to determine the magnitude and duration of post-impact hydrothermal circulation.

Catherine Ross↗

Precipitation and Latent Heating Distributions from Satellite Passive Microwave Radiometry: Method and Uncertainties - Part 1

A revised Bayesian algorithm for estimating surface rain rate, convective rain proportion, and latent heating/drying profiles from satellite-borne passive microwave radiometer observations over ocean backgrounds is described. The algorithm searches a large database of cloud-radiative model simulations to find cloud profiles that are radiatively consistent with a given set of microwave radiance measurements. The properties of these radiatively consistent profiles are then composited to obtain best estimates of the observed properties. The revised algorithm is supported by an expanded and more physically consistent database of cloud-radiative model simulations. The algorithm also features a better quantification of the convective and non-convective contributions to total rainfall, a new geographic database, and an improved representation of background radiances in rain-free regions. Bias and random error estimates are derived from applications of the algorithm to synthetic radiance data, based upon a subset of cloud resolving model simulations, and from the Bayesian formulation itself. Synthetic rain rate and latent heating estimates exhibit a trend of high (low) bias for low (high) retrieved values. The Bayesian estimates of random error are propagated to represent errors at coarser time and space resolutions, based upon applications of the algorithm to TRMM Microwave Imager (TMI) data. Errors in instantaneous rain rate estimates at 0.5 deg resolution range from approximately 50% at 1 mm/h to 20% at 14 mm/h. These errors represent about 70-90% of the mean random deviation between collocated passive microwave and spaceborne radar rain rate estimates. The cumulative algorithm error in TMI estimates at monthly, 2.5 deg resolution is relatively small (less than 6% at 5 mm/day) compared to the random error due to infrequent satellite temporal sampling (8-35% at the same rain rate).

Olson, William S.↗

Multi-wavelength Radio Continuum Emission Studies of Dust-free Red Giants

Multi-wavelength centimeter continuum observations of non-dusty, non-pulsating K spectral-type red giants directly sample their chromospheres and wind acceleration zones. Such stars are feeble emitters at these wavelengths, however, and previous observations have provided only a small number of modest signal-to-noise measurements slowly accumulated over three decades. We present multi-wavelength Karl G. Jansky Very Large Array thermal continuum observations of the wind acceleration zones of two dust-free red giants, Arcturus (alpha Boo: K2 III) and Aldebaran (alpha Tau: K5 III). Importantly, most of our observations of each star were carried out over just a few days, so that we obtained a snapshot of the different stellar atmospheric layers sampled at different wavelengths, independent of any long-term variability. We report the first detections at several wavelengths for each star including a detection at 10 cm (3.0 GHz: S band) for both stars and a 20 cm (1.5 GHz: L band) detection for alpha Boo. This is the first time single (non-binary) luminosity class III red giants have been detected at these continuum wavelengths. Our long-wavelength data sample the outer layers of alpha Boo's atmosphere where its wind velocity is approaching (or possibly has reached) its terminal value and the ionization balance is becoming frozen-in. For alpha Tau, however, our long-wavelength data are still sampling its inner atmosphere, where the wind is still accelerating probably due to its lower mass-loss rate. We compare our data with published semi-empirical models based on ultraviolet data, and the marked deviations highlight the need for new atmospheric models to be developed. Spectral indices are used to discuss the possible properties of the stellar atmospheres, and we find evidence for a rapidly cooling wind in the case of alpha Boo. Finally, we develop a simple analytical wind model for alpha Boo based on our new long-wavelength flux measurements.

chromospheres↗

Contribution of vertebral deformities to chronic back pain and disability. The Study of Osteoporotic Fractures Research Group

Among 2992 white women aged 65-70 years recruited from population-based listings, we measured radiographic vertebral dimensions of T5-L4 and calculated ratios of heights: anterior/posterior, mid/posterior, and posterior/posterior of either adjacent vertebra. The degree of deformity for each vertebra was analyzed in terms of the number of standard deviations (SD) that ratio differed from the mean ratio calculated for the same vertebral level in this population. We correlated the severity of each woman's worst vertebral deformity with back pain, back disability in six activities of daily living, and height loss since age 25. Only 39.4% of the cohort had no vertebral deformity; 10.2% had a deformity greater than or equal to 4 SD. Vertebral deformities less than 4 SD below the mean were not associated with increased back pain, disability, or loss of height. In contrast, women whose deformity was greater than or equal to 4 SD had a 1.9 (95% CI, 1.5-2.4) times higher risk of moderate to severe back pain and a 2.6 (95% CI, 1.7-3.9) times higher risk of disability involving the back; they were also 2.5 (95% CI, 2.0-3.2) times more likely to have lost greater than or equal to 4 cm in height. All three types of vertebral deformity (wedge, end plate, and crush) were equally associated with these outcomes. Multiple deformities less than 4 SD did not increase the likelihood of these three outcomes, but multiple deformities greater than or equal to 4 SD tended to be associated with increased back pain, disability, and height loss. This large cross-sectional study suggests that vertebral deformities cause substantial pain, disability, or loss of height only if vertebral height ratios fall 4 SD below the normal mean. Much back pain could not be attributed to vertebral deformities, suggesting other causes.

Non-NASA Center↗

An Evaluation of the Liquid Cloud Droplet Effective Radius Derived from MODIS, Airborne Remote Sensing, and in Situ Measurements from CAMP2Ex

The cloud drop effective radius (Re) of the drop size distribution derived from passive satellite sensors is a key variable used in climate research. Validation of these satellite products has often taken place under stratiform cloud conditions that favor the assumption of cloud horizontal homogeneity used by the retrieval techniques. However, many studies have noted concerns with respect to significant biases in retrieved Re arising from cloud heterogeneity, for example, in cumulus cloud fields. Here, we examine data collected during the 2019 “Cloud, Aerosol and Monsoon Processes Philippines Experiment” (CAMP2Ex), which, in part, targeted the objective of providing the first detailed evaluation of Re retrieved across multiple platforms and techniques in a cumulus and congestus cloud region. Our evaluation consists of cross-comparisons of Re between the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the Terra satellite, the Research Scanning Polarimeter (RSP) onboard the NASA P-3 aircraft, and in situ measurements from both the NASA P-3 and Learjet aircraft that are all taken in close spatiotemporal proximity to the same cloud fields. A particular advantage of our approach lies in the capability of the RSP to retrieve Re using a bi-spectral MODIS approach and a polarimetric approach, which allows for the evaluation of bi-spectral and polarimetric Re retrievals from an airborne perspective using the same samples. Averaged over all P-3 flight segments examined here for warm clouds, the RSP polarimetric method, the in situ method, and the bias-adjusted MODIS method of Fu et al. (2019) show a comparable median (mean ± standard deviation) for the Re samples of 9.6 (10.2 ± 4.0) µm, 11.0 (13.6 ± 11.3) µm, and 10.4 (10.8 ± 3.8) µm, respectively. These values are far lower than the values of 15.1 (16.2 ± 5.5) µm and 17.2 (17.7 ± 5.7) µm from the bi-spectral retrievals of RSP and MODIS, respectively. Similar results are observed when Re is segregated by cloud-top height and in detailed case studies. The clouds sampled during CAMP2Ex consist of mostly small (mean transect length ∼ 1.4 km) and low clouds (mean cloud-top height ∼ 1 km), which had more numerous small clouds than the trade wind cumuli sampled in past field campaigns such as Rain in Shallow Cumulus over the Ocean (RICO) and the Indian Ocean Experiment (INDOEX). The overestimates of Re from the RSP bi-spectral technique compared with the polarimetric technique increased as cloud size and cloud optical depth decreased. Drizzle, cloud-top bumpiness, and solar zenith angle, however, are not closely correlated with the overestimate of bi-spectral Re. For shallow clouds that dominated the liquid cloud cover for the CAMP2Ex region and period, we show that 3-D radiative transfer and cloud heterogeneity, particularly for the optically thin and small clouds, appear to be the leading cause of the large positive biases in bi-spectral retrievals. Because this bias varies with the underlying structure of the cloud field, caution continues to be warranted in studies that use bi-spectral Re retrievals in cumulus cloud fields.

Cloud microphysics↗

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↗

Arctic sea ice leads from advanced very high resolution radiometer images

A large number of advanced very high resolution radiometer (AVHRR) images from throughout 1989 are analyzed to determine lead characteristics. The units of analysis are square 200-km cells, and there are 270 such cells in the data set. Clouds are masked manually. Leads determine from images of the potential open water delta, a scaled version of the surface temperature or albedo that weights thin ice by its thermal or brightness impact. The lead fraction is determined as the mean delta, the monthly mean lead fraction ranges from 0.02 in winter to 0.06 in summer in the central Arctic and is near 0.08 in the winter in the peripheral seas. A method of accounting for lead width sampling errors due to the finite sample areas is introduced. In the central Arctic the observed mean lead width for a threshold of delta = 0.1 ranges from 2 or 3 km (near the resolution of the instrument) in the winter to 6 km in the summer. In the peripheral seas it is about 5 km in the winter. Width distributions are often more heavily weighted in the tail than exponential distributions and are well approximated by a power law. The along-track, number density power law N = aw(exp -6) has a mean exponent of b = 1.60 (standard deviation 0.18) and shows some seasonal variability. Mean floe widths in the central Arctic are 40 to 50 km in the winter, dropping to about 10 km in the summer. For floes the power law has a mean exponent of 0.93 and exhibits a clearer annual cycle. Lead orientation is determined with a method based on the direction of maximum extent.

Lindsay, R. W.↗

Characterization of the Ryugu Surface By Means Of the Variability of the Near-Infrared Spectral Slope in NIRS3 Data

The Near-Earth Asteroid 162173 Ryugu (1999 JU3) was investigated by the JAXA Hayabusa2 mission from June 2018 to November 2019. The data acquired by NIRS3 spectrometer revealed a dark surface with a positive near-infrared spectral slope. In this work we investigated the spectral slope variations across the Ryugu surface, providing information about physical/chemical properties of the surface. We analysed the calibrated, thermally and photometrically corrected NIRS3 data and we estimated the mean value of spectral slope between 1.9 μm and 2.5 μm, corresponding to 0.163. Starting from the mean value of slope and moving in step of 1 standard deviation (0.022), we defined 9 “families of slope”, the Low-Red-Slope families (LR1, LR2 and LR3) and the High-Red-Sloped families (HR1, HR2, HR3, HR4, HR5, HR6). The mean values of some spectral parameters were estimated for each family, such as the reflectance factor at 1.9 μm, the spectral slope, the depth of bands at 2.7 μm and at 2.8 μm. A progressive spectral reddening, darkening and weakening/narrowing of OH bands is observed moving from the LR families to the HR families. We concluded that the spectral variability observed among families is the resulting contribution of the thermal metamorphism experienced by Ryugu after the catastrophic disruption of its parent body and space weathering processes that occurred on airless bodies as Ryugu, such as impact cratering and solar wind irradiation. As a consequence, the HR1, LR1, LR2 and LR3 families, corresponding to equatorial ridge and crater rims, are the less altered regions on Ryugu surface, which experienced the minor alteration and OH devolatilization; the HR2, HR3, HR4, HR5 families, coincident with floors and walls of impact craters, are the most altered areas, result of the three processes occurred on Ryugu. The strong reddening of the HR6 family (coincident with Ejima Saxum) is likely due to the fine-sized material covering the large boulder.

Ryugu↗

Diurnal Variability of Vertical Structure from a TRMM Passive Microwave "Virtual Radar" Retrieval

Robust description of the diurnal cycle from TRMM observations is complicated by the limitations of Low Earth Orbit (LEO) sampling; from a 'climatological' perspective, sufficient sampling must exist to control for both spatial and seasonal variability, before tackling an additional diurnal component (e.g., with 8 additional 3-hourly or 24 1-hourly bins). For documentation of vertical structure, the narrow sample swath of the TRMM Precipitation Radar limits the resolution of any of these components. A neural-network based 'virtual radar" retrieval has been trained and internally validated, using multifrequency / multipolarization passive microwave(TM1) brightness temperatures and textures parameters and lightning (LIS) observations, as inputs, and PR volumetric reflectivity as targets (outputs). By training the algorithms (essentially highly multivariate, nonlinear regressions) on a very large sample of high-quality co-located data from the center of the TRMM swath, 3D radar reflectivity and derived parameters (VIL, IWC, Echo Tops, etc.) can be retrieved across the entire TMI swath, good to 8-9% over the dynamic range of parameters. As a step in the retrieval (and as an output of the process), each TMI multifrequency pixel (at 85 GHz resolution) is classified into one of the 25 archetypal radar profile vertical structure "types", previously identified using cluster analysis. The dynamic range of retrieved vertical structure appears to have higher fidelity than the current (Version 6) experimental GPROF hydrometeor vertical structure retrievals. This is attributable to correct representation of the prior probabilities of vertical structure variability in the neural network training data, unlike the GPROF cloud-resolving model training dataset used in the V6 algorithms. The LIS lightning inputs are supplementary inputs, and a separate offline neural network has been trained to impute (predict) LIS lightning from passive-microwave-only data. The virtual radar retrieval is thus, in principle, extensible to Aqua/AMSR-E and NPOESS/CMIS passive microwave instruments. The virtual radar approach yields a threefold increase in effective sampling from the mission, albeit of lower-quality "retrieved" data, reducing the variance of local estimates by one third (or the standard deviation by-0.57). In this talk, the variance reduction is leveraged to more finely resolve global diurnal variability in both space and time (local hour).

Boccippio, Dennis J.↗

Gravitational collapse of small-scale structure as the origin of the Lyman-alpha forest

If gravitational clustering is a hierarchical process, the present large-scale structure of the galaxy distribution implies that structures on smaller scales must have formed at high redshift. We simulate the formation of small-scale structure (average cell mass: (delta) bar m(sub b) = 10(sup 4.2) solar mass) and the evolution of photoionized gas, in the specific case of the CDM + Lambda model. The photoionized gas has a natural minimal scale of collapse, the Jeans scale (m(sub b,J) is approximately equal to 10(exp 9) solar mass). We find that low column density (N(sub HI) is less than or equal to 10(exp 14)/sq cm) lines originate in regions resembling Zel'dovich pancakes, where gas with overdensities in the range 3-30 is enclosed by two shocks but is typically reexpanding at approximately the Hubble velocity. However, higher column density lines stem from more overdense regions where the shocked gas is cooling. We show that this model can probably account for the observed number of lines, their distribution in column density and b-parameters, as well as the cloud physical sizes as observed in gravitationally lensed quasars. We find a redshift evolution that is too steep; however, this may be due to insufficient dynamical range in the simulation or because the specific model is incorrect. The model predicts that high signal-to-noise observations should find systematic deviations from Voigt profiles, mainly in the form of broad wings in the line profiles, and that a fluctuating Gunn-Peterson effect will be detected, which can be modeled as a superposition of weak lines with a wide range of b-parameters.

Cen, Renyue↗

Monitoring Airspace Complexity and Determining Contributing Factors

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

complexity↗

A Peculiar Class of Debris Disks from Herschel/DUNES: A Steep Fall Off in the Far Infrared

Context. The existence of debris disks around old main sequence stars is usually explained by continuous replenishment of small dust grains through collisions from a reservoir of larger objects. Aims. We present photometric data of debris disks around HIP 103389 (HD199260), HIP 100350 (HN Peg, HD206860), and HIP 114948 (HD 219482), obtained in the context of our Herschel Open TIme Key Program DUNES (DUst around NEarby Stars). Methods. We used Herschel/PACS to detect the thermal emission of the three debris disks with a 30 sigma sensitivity of a few mJy at l00 micron and 160 micron. In addition, we obtained Herschel/PACS photometric data at 70 micron for HIP 103389. These observations are complemented by a large variety of optical to far-infrared photometric data. Two different approaches are applied to reduce the Herschel data to investigate the impact of data reduction on the photometry. We fit analytical models to the available spectral energy distribution (SED) data using the fitting method of simulated therma1 annealing as well as a classical grid search method. Results. The SEDs of the three disks potentially exhibit an unusually steep decrease at wavelengths >= 70 micron. We investigate the significance of the peculiar shape of these SEDs and the impact on models of the disks provided it is real. Using grain compositions that have been applied successfully for modeling of many other debris disks, our modeling reveals that such a steep decrease of the SEDs in the long wavelength regime is inconsistent with a power-law exponent of the grain size distribution -3.5 expected from a standard equilibrium collisional cascade. In contrast, a steep grain size distribution or, alternatively an upper grain size in the range of few tens of micrometers are implied. This suggests that a very distinct range of grain sizes would dominate the thermal. emission of such disks. However, we demonstrate that the understanding of the data of faint sources obtained with Herschel is still incomplete and that the significance of our results depends on the version of the data reduction pipeline used. Conclusions. A new mechanism to produce the dust in the presented debris disks, deviations from the conditions required for a standard equilibrium collisional cascade (grain size exponent of -3.5), and/or significantly different dust properties would be necessary to explain the potentially steep SED shape of the three debris disks presented.

Ertel, S.↗

An Upper Bound on Orbital Debris Collision Probability When Only One Object has Position Uncertainty Information

Upper bounds on high speed satellite collision probability, P (sub c), have been investigated. Previous methods assume an individual position error covariance matrix is available for each object. The two matrices being combined into a single, relative position error covariance matrix. Components of the combined error covariance are then varied to obtain a maximum P (sub c). If error covariance information for only one of the two objects was available, either some default shape has been used or nothing could be done. An alternative is presented that uses the known covariance information along with a critical value of the missing covariance to obtain an approximate but useful P (sub c) upper bound. There are various avenues along which an upper bound on the high speed satellite collision probability has been pursued. Typically, for the collision plane representation of the high speed collision probability problem, the predicted miss position in the collision plane is assumed fixed. Then the shape (aspect ratio of ellipse), the size (scaling of standard deviations) or the orientation (rotation of ellipse principal axes) of the combined position error ellipse is varied to obtain a maximum P (sub c). Regardless as to the exact details of the approach, previously presented methods all assume that an individual position error covariance matrix is available for each object and the two are combined into a single, relative position error covariance matrix. This combined position error covariance matrix is then modified according to the chosen scheme to arrive at a maximum P (sub c). But what if error covariance information for one of the two objects is not available? When error covariance information for one of the objects is not available the analyst has commonly defaulted to the situation in which only the relative miss position and velocity are known without any corresponding state error covariance information. The various usual methods of finding a maximum P (sub c) do no good because the analyst defaults to no knowledge of the combined, relative position error covariance matrix. It is reasonable to think, given an assumption of no covariance information, an analyst might still attempt to determine the error covariance matrix that results in an upper bound on the P (sub c). Without some guidance on limits to the shape, size and orientation of the unknown covariance matrix, the limiting case is a degenerate ellipse lying along the relative miss vector in the collision plane. Unless the miss position is exceptionally large or the at-risk object is exceptionally small, this method results in a maximum P (sub c) too large to be of practical use. For example, assuming that the miss distance is equal to the current ISS alert volume along-track (+ or -) distance of 25 kilometers and that the at-risk area has a 70 meter radius. The maximum (degenerate ellipse) P (sub c) is about 0.00136. At 40 kilometers, the maximum P (sub c) would be 0.00085 which is still almost an order of magnitude larger than the ISS maneuver threshold of 0.0001. In fact, a miss distance of almost 340 kilometers is necessary to reduce the maximum P (sub c) associated with this degenerate ellipse to the ISS maneuver threshold value. Such a result is frequently of no practical value to the analyst. Some improvement may be made with respect to this problem by realizing that while the position error covariance matrix of one of the objects (usually the debris object) may not be known the position error covariance matrix of the other object (usually the asset) is almost always available. Making use of the position error covariance information for the one object provides an improvement in finding a maximum P (sub c) which, in some cases, may offer real utility. The equations to be used are presented and their use discussed.

Frisbee, Joseph H., Jr.↗

Evaluating Mass Analyzers as Candidates for Small, Portable, Rugged Single Point Mass Spectrometers for Analysis of Permanent Gases

For Space Shuttle launch safety, there is a need to monitor the concentration of H2, He, O2 and Ar around the launch vehicle. Currently a large mass spectrometry system performs this task, using long transport lines to draw in samples. There is great interest in replacing this stationary system with several miniature, portable, rugged mass spectrometers which act as point sensors which can be placed at the sampling point. Five commercial and two non-commercial analyzers are evaluated. The five commercial systems include the Leybold Inficon XPR-2 linear quadrupole, the Stanford Research (SRS-100) linear quadrupole, the Ferran linear quadrupole array, the ThermoQuest Polaris-Q quadrupole ion trap, and the IonWerks Time-of-Flight (TOF). The non-commercial systems include a compact double focusing sector (CDFMS) developed at the University of Minnesota, and a quadrupole ion trap (UF-IT) developed at the University of Florida. The System Volume is determined by measuring the entire system volume including the mass analyzer, its associated electronics, the associated vacuum system, the high vacuum pump and rough pump. Also measured are any ion gauge controllers or other required equipment. Computers are not included. Scan Time is the time required for one scan to be acquired and the data to be transferred. It is determined by measuring the time required acquiring a known number of scans and dividing by said number of scans. Limit of Detection is determined first by performing a zero-span calibration (using a 10-point data set). Then the limit of detection (LOD) is defined as 3 times the standard deviation of the zero data set. (An LOD of 10 ppm or less is considered acceptable.)

Arkin, C. Richard↗

Development of New Lunar Highland Regolith Simulant, NUW-LHT-5M

Introduction: NASA has a need for large quantities of lunar simulants that closely match future manned lunar missions at landing sites near the Lunar South Pole. The current simulant was designed to approximate NU-LHT-2M, and -4M, except using natural minerals and a fully synthetic, non-basaltic, high calcium glass. Washington Mills in Niagara Falls, NY was contracted to produce this new simulant due to their electric arc furnace (EAF) fusion technology and capabilities for crushing and sizing ceramics. The primary objectives of this work were: (1) to create a commercially produced simulant that contains a high fidelity, high calcium (An90+), low Mg/Fe glass without basaltic constituents; (2) to create a simulant that uses existing, terrestrially mined, and readily available minerals blended with a synthetic glass, and (3) accomplish objectives (1) and (2) using conventional processing methods capable of production levels that can meet substantial current and future demands of NASA projects. Another consideration was that highland simulants closest to meeting these requirements, NU-LHT-2M and -4M types, are out of production and largely out of stock. Design: The design chosen for development and production for this work is NU-LHT-2M, de-signed by Doug Stoeser and Douglas Rickman, which uses Stillwater anorthosite and norite rocks and synthetic glass. The composition targeted the average composition, glass content and particle size distribution of Apollo 16 samples [1]. The design utilizes two Stillwater minerals, anorthosite (37.7 wt.%), norite (17.6 wt.%), a commercially sourced olivine (4.7 wt.%), and a synthetic glass (40 wt.%). Processing: Stillwater rocks were hand collected and initially crushed by the USGS. Further crushing and milling was performed at Washing-ton Mills. For the glass, oxides were blended in the ap-propriate ratios, blended in a v-blender, and placed in a graphite lined, water cooled pot. The mixture was fused in Washington Mills’ pilot scale, 500 kW EAF. The molten glass was then poured into water and quenched to inhibit crystallization. Several pour/quench cycles were required to produce the total quantity of glass needed. Anorthosite, norite and olivine were milled together to the target particle size in a ball mill. The glass was ball milled separately. The complete simulant was produced by blending the milled materials. Particle size distribution, phase and chemical analyses of the materials was were per-formed at Washington Mills. Results: The glass composition was found to match the Apollo 16 oxide content average [1] quite well, with the exception of SiO2 and Fe2O3. The deviation of SiO2 and Fe2O3 from batched content was believed to be due to limitations inherent in the use of the arc furnace technique for these types of materials. A particle size distribution of the full simulant formulation closely matching the average for the Apollo 16 samples was also achieved. Fig. 1 shows the morphology of the final simulant particles. The size distribution is compared to the target in Fig. 2. Conclusions: A new lunar highlands type simulant containing a high-quality synthetic glass, closely matching an average composition of the Apollo 16 regolith was produced. Processing utilized easily scalable mineral and glass making methods.

simulant↗

Transparent magnetic state in single crystal Nd(1.85)Ce(0.15)CuO(4-y) superconductors

Several experimental studies have been reported as evidence of Josephson coupling between the superconducting layers in the highly anisotropic oxide such as the Bi2Sr2CaCu2O8 and Tl2Ba2CuO6 systems. These include the large penetration depth of 100 mu m measured, ac and dc Josephson effects. Recently two critical temperatures corresponding to Josephson coupling in between the layers and the Berezinskii-Kosterlitz-Thouless transition in the ab-plane have been directly observed in the transport measurements. If the field is applied parallel to the superconducting layers, the magnetic excitation is not the conventional Abrikosov vortices, but the Josephson vortices which extend lambda(sub ab) in the c-axis direction and lambda(sub J) = gamma s in the plane (s is the interlayer distance, gamma is the anisotropy constant). Because of the weak screening effect associated with the Josephson vortices, there have been predictions of magnetic transparent states at magnetic field above a characteristic field H(sub J), a behavior distinctively different from that of the type-II superconductors. In this paper, we report an experimental result which illustrates a transition from the Meissner state to the magnetic transparent state in single crystal of Nd(1.85)Ce(0.15)CuO(4-y). Magnetization has been measured as a function of temperature and field in the magnetic field parallel or close to ab-plane geometry. For a fixed magnetic field, the magnetization shows a two-step transition in M(T); for a fixed temperature, the magnetization shows an abrupt change to almost zero value above a characteristic field H(sub J), an indication of magnetic transparent state. The data of magnetization as a function of field clearly deviates from the behavior predicted by the Abrikosov theory for type-II superconductors. Instead, the data fit well into the picture of Josephson decoupling between the CuO2 layers.

Zuo, F.↗

Advancements in the Aerosol Robotic Network (AERONET) Version 3 database – automated near-real-time quality control algorithm with improved cloud screening for Sun photometer aerosol optical depth (AOD) measurements

The Aerosol Robotic Network (AERONET) has provided highly accurate, ground-truth measurements of the aerosol optical depth (AOD) using Cimel Electronique Sun–sky radiometers for more than 25 years. In Version 2 (V2) of the AERONET database, the near-real-time AOD was semiautomatically quality controlled utilizing mainly cloud-screening methodology, while additional AOD data contaminated by clouds or affected by instrument anomalies were removed manually before attaining quality-assured status (Level 2.0). The large growth in the number of AERONET sites over the past 25 years resulted in significant burden to the manual quality control of millions of measurements in a consistent manner. The AERONET Version 3 (V3) algorithm provides fully automatic cloud screening and instrument anomaly quality controls. All of these new algorithm updates apply to near-real-time data as well as post-field-deployment processed data, and AERONET reprocessed the database in 2018. A full algorithm redevelopment provided the opportunity to improve data inputs and corrections such as unique filter-specific temperature characterizations for all visible and near-infrared wavelengths, updated gaseous and water vapor absorption coefficients, and ancillary data sets. The Level 2.0 AOD quality-assured data set is now available within a month after post-field calibration, reducing the lag time from up to several months. Near-real-time estimated uncertainty is determined using data qualified as V3 Level 2.0 AOD and considering the difference between the AOD computed with the pre-field calibration and AOD computed with pre-field and post-field calibration. This assessment provides a near-real-time uncertainty estimate for which average differences of AOD suggest a +0.02 bias and one sigma uncertainty of 0.02, spectrally, but the bias and uncertainty can be significantly larger for specific instrument deployments. Long-term monthly averages analyzed for the entire V3 and V2 databases produced average differences (V3–V2) of +0.002 with a ±0.02 SD (standard deviation), yet monthly averages calculated using time-matched observations in both databases were analyzed to compute an average difference of −0.002 with a ±0.004 SD. The high statistical agreement in multiyear monthly averaged AOD validates the advanced automatic data quality control algorithms and suggests that migrating research to the V3 database will corroborate most V2 research conclusions and likely lead to more accurate results in some cases.

David M. Giles↗