Search NASA⌕ Search

SEARCH · Search NASA

Results for “Random Uncertainty”

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.

430 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Computational Model Prediction and Biological Validation Using Simplified Mixed Field Exposures for the Development of a GCR Reference Field

The yield of chromosomal aberrations has been shown to increase in the lymphocytes of astronauts after long-duration missions of several months in space. Chromosome exchanges, especially translocations, are positively correlated with many cancers and are therefore a potential biomarker of cancer risk associated with radiation exposure. Although extensive studies have been carried out on the induction of chromosomal aberrations by low- and high-LET radiation in human lymphocytes, fibroblasts, and epithelial cells exposed in vitro, there is a lack of data on chromosome aberrations induced by low dose-rate chronic exposure and mixed field beams such as those expected in space. Chromosome aberration studies at NSRL will provide the biological validation needed to extend the computational models over a broader range of experimental conditions (more complicated mixed fields leading up to the galactic cosmic rays (GCR) simulator), helping to reduce uncertainties in radiation quality effects and dose-rate dependence in cancer risk models. These models can then be used to answer some of the open questions regarding requirements for a full GCR reference field, including particle type and number, energy, dose rate, and delivery order. In this study, we designed a simplified mixed field beam with a combination of proton, helium, oxygen, and iron ions with shielding or proton, helium, oxygen, and titanium without shielding. Human fibroblasts cells were irradiated with these mixed field beam as well as each single beam with acute and chronic dose rate, and chromosome aberrations (CA) were measured with 3-color fluorescent in situ hybridization (FISH) chromosome painting methods. Frequency and type of CA induced with acute dose rate and chronic dose rates with single and mixed field beam will be discussed. A computational chromosome and radiation-induced DNA damage model, BDSTRACKS (Biological Damage by Stochastic Tracks), was updated to simulate various types of CA induced by acute exposures of the mixed field beams used for the experiments. The chromosomes were simulated by a polymer random walk algorithm with restrictions to their respective domains in the nucleus [1]. The stochastic dose to the nucleus was calculated with the code RITRACKS [2]. Irradiation of a target volume by a mixed field of ions was implemented within RITRACKs, and the fields of ions can be delivered over specific periods of time, allowing the simulation of dose-rate effects. Similarly, particles of various types and energies extracted from a pre-calculated spectra of galactic cosmic rays (GCR) can be used in RITRACKS. The number and spatial location of DSBs (DNA double-strand breaks) were calculated in BDSTRACKS using the simulated chromosomes and local (voxel) dose. Assuming that DSBs led to chromosome breaks, and simulating the rejoining of damaged chromosomes occurring during repair, BDSTRACKS produces the yield of various types of chromosome aberrations as a function of time (only final yields are presented). A comparison between experimental and simulation results will be shown.

Hada, M.↗
NASA NTRS2024

Summary and Annotated Bibliography of Measurement Error Corrections with Potential Application in Future Quesst Mission Community Noise Studies

This document is motivated by likely needs of the Quesst mission community response tests, which will culminate in data collection and estimation of dose-response regression relationships for consideration by domestic and international aviation regulators. Furthermore, basic research questions evaluating interactions between rates of community annoyance, dose levels, and indicators of the presence of rattle, vibration, and startle hinge on hypothesis testing in the context of regression models. For a variety of reasons, noise doses may be known only imprecisely and may not reflect the actual level experienced by responding subjects. These differences between true dose and estimated dose, be they systematic or random, constitute covariate measurement error. Available statistics literature speaks to the impacts of measurement error on regression models, both in terms of bias in estimated coefficients and predicted values, and in terms of the loss of statistical power for hypothesis testing. Given the particulars of a categorical annoyance response variable and a continuous noise dose predictor variable subject to measurement error during testing, the emphasis of this report is on findings and methods pertinent to generalized linear (and mixed) models likely to be employed during the Quesst mission community tests. We reach the following conclusions: 1. Of four reviewed methods, structural Bayesian measurement error models and simulation extrapolation (SIMEX) may be the most readily applicable to Quesst mission community noise study objectives. 2. If warranted, a linear measurement model can help model systematic sources of measurement error that the classical measurement error does not. 3. For its ready implementation and small additional input requirements, simulation extrapolation may be ideally suited for addressing secondary research questions involving interactions between annoyance, noise dose, and other factors through hypothesis testing. 4. For their flexibility and ability to propagate uncertainty, structural Bayesian hierarchical models have great appeal for mission purposes; some care may be needed in developing appropriate probability models describing actual noise exposure during testing. An annotated bibliography logs additional papers and resources that may be of value to analysts in other projects and disciplines.

Dose-Response Model↗
NASA NTRS2007

Os-186 and Os-187 Enrichments and High-He-3/He-4 sources in the Earth's Mantle: Evidence from Icelandic Picrites

Picrites from the neovolcanic zones in Iceland display a range in Os-187/Os-188O from 0.1297 to 0.1381 ((gamma)Os = 0.0 to 6.5) and uniform Os-186/Os-188 of 0.1198375+/-32 (2 (sigma)). The value for Os-186/Os-188 is within uncertainty of the present-day value for the primitive upper mantle of 0.1198398+/-16. These Os isotope systematics are best explained by ancient recycled crust or melt enrichment in the mantle source region. If so, then the coupled enrichments displayed in Os-186/Os-188 and Os-187/Os-188 from lavas of other plume systems must result from an independent process, the most viable candidate at present remains core-mantle interaction. While some plumes with high He-3/He-4, such as Hawaii, appear to have been subjected to detectable addition of Os (and possibly He) from the outer core, others such as Iceland do not. A positive correlation between Os-187/Os-188 and He-3/He-4 from 9.6 to 19 RA in Iceland picrites is best modeled as mixtures of 500 Ma or older ancient recycled crust mixed with primitive mantle, creating a hybrid source region that subsequently mixes with the convecting MORB mantle during ascent and melting. This multistage mechanism to explain these isotope systematics is consistent with ancient recycled crust juxtaposed with more primitive, relatively He-rich mantle, in convective isolation from the upper mantle, most likely in the lowermost mantle. This is inconsistent with models that propose random mixing between heterogeneities in the convecting upper mantle as a mechanism to explain the observed isotopic variation in oceanic lavas or models that produce a high He-3/He-4 signature in melt depleted and strongly outgassed, He-poor mantle. Instead these systematics require a deep mantle source to explain the 3He/4He signature in Iceland lavas. The He-3/He-4 of lavas derived from the Iceland plume changed over time, from a maximum of 50 RA at 60 Ma, to approximately 25-27 RA at present. The changes are coupled with distinct compositional gaps between the different aged lavas when H-3/He-4 is plotted versus various geochemical parameters such as Nd-143/Nd-144 and La/Sm. These relationships can be interpreted as an increase in the proportion of ancient recycled crust in the upwelling plume over this time period.

Brandon, Alan D.↗
NASA NTRS2020

Effects of electromagnetic wave interference on observations of the Earth radiation budget

This paper investigates conditions necessary to match the irradiance derived by integrating radiances measured by a narrow field of view scanning radiometer with the irradiance measured by a hemispherical radiometer, both placed at a satellite altitude for Earth radiation budget estimates. When all sources are similar and they are spatially distributed randomly, then integrating radiance for the irradiance does not introduce a bias. Although the exact magnitude of the bias in other conditions is unknown, a finite area of the aperture that is much larger than the coherence area of radiation contributing to the Earth radiation budget, and a finite time to take a single measurement that is longer than the coherence time are likely to make the difference of the irradiance integrated from radiances and the irradiance measured by a hemispherical instrument insignificant. This conclusion does not contradict the existence of spatial coherence of light from incoherent sources. Therefore, electromagnetic energy absorbed by Earth is derivable from radiances measured by a scanning radiometer integrated over the Earth-viewing hemisphere and then averaging across all locations on the satellite orbital sphere when combined with solar irradiance measurements. Comparisons made in earlier studies show that the difference is less than 1%. In addition, when surface irradiances computed by a radiative transfer model constrained by top-of-atmosphere irradiances derived from radiance measurements are compared with downward shortwave irradiances taken by combinations of a pyreheliometer and a shaded pyranometer, or pyranometers, and with longwave ir- radiances taken by pyrgeometers, the biases in monthly mean irradiances are less than the uncertainties in the surface observations.

Seiji Kato↗
NASA NTRS2000

A Probabilistic Framework for the Validation and Certification of Computer Simulations

The paper presents a methodology for quantifying, propagating, and managing the uncertainty in the data required to initialize computer simulations of complex phenomena. The purpose of the methodology is to permit the quantitative assessment of a certification level to be associated with the predictions from the simulations, as well as the design of a data acquisition strategy to achieve a target level of certification. The value of a methodology that can address the above issues is obvious, specially in light of the trend in the availability of computational resources, as well as the trend in sensor technology. These two trends make it possible to probe physical phenomena both with physical sensors, as well as with complex models, at previously inconceivable levels. With these new abilities arises the need to develop the knowledge to integrate the information from sensors and computer simulations. This is achieved in the present work by tracing both activities back to a level of abstraction that highlights their commonalities, thus allowing them to be manipulated in a mathematically consistent fashion. In particular, the mathematical theory underlying computer simulations has long been associated with partial differential equations and functional analysis concepts such as Hilbert spares and orthogonal projections. By relying on a probabilistic framework for the modeling of data, a Hilbert space framework emerges that permits the modeling of coefficients in the governing equations as random variables, or equivalently, as elements in a Hilbert space. This permits the development of an approximation theory for probabilistic problems that parallels that of deterministic approximation theory. According to this formalism, the solution of the problem is identified by its projection on a basis in the Hilbert space of random variables, as opposed to more traditional techniques where the solution is approximated by its first or second-order statistics. The present representation, in addition to capturing significantly more information than the traditional approach, facilitates the linkage between different interacting stochastic systems as is typically observed in real-life situations.

Ghanem, Roger↗
NASA NTRS2024

Synergistic Retrievals of Ice in High Clouds From Elastic Backscatter Lidar, Ku-band Radar and Submillimeter Wave Radiometer Observations

In this study, we investigate the synergy of elastic backscatter lidar, Ku-band radar, and sub-millimeter-wave radiometer measurements in the retrieval of ice from satellite observations. The synergy is analyzed through the generation of a large dataset of IceWater Content (IWC) profiles and simulated lidar, radar and radiometer observations. The characteristics of the instruments e.g. frequencies, sensitivities, etc. are set based on the expected characteristics of instruments of the Atmosphere Observing System (AOS) mission. A hold-out validation methodology is used to assess the accuracy of the IWC profiles retrieved from various combinations of observations from the three instruments. Specifically, the IWC and associated observations are randomly divided into two datasets, one for training and the other for evaluation. The training dataset is used to train the retrieval algorithm, while the evaluation dataset is used to assess the retrieval performance. The dataset of IWC profiles is derived from CloudSat reflectivity and CALIOP lidar observations. The retrieval of the ice water content IWC profiles from the computed observations is achieved in two steps. In the first step, a class, out of 18 potential classes characterized by different vertical distribution of IWC, is estimated from the observations. The 18 classes are predetermined based on the k-Means clustering algorithm. In the second step, the IWC profile is estimated using an Ensemble Kalman Smoother (EKS) algorithm that uses the estimated class as a priori information. The results of the study show that the synergy of lidar, radar, and radiometer observations is significant in the retrieval of the IWC profiles. Nevertheless, it should be mentioned that this synergy was found under idealized conditions, and additional work might be required to materialize it in practice. The inclusion of the lidar backscatter observations in the retrieval process has a larger impact on the retrieval performance than the inclusion of the radar observations. As ice clouds have a significant impact on atmospheric radiative processes, this work is relevant to ongoing efforts to reduce uncertainties in climate analyses and projections.

Mircea Grecu↗
NASA NTRS2004

Large Alluvial Fans on Mars

Several dozen distinct alluvial fans, 10 to greater than 40 km long downslope are observed exclusively in highlands craters. Within a search region between 0 deg. and 30 deg. S, alluvial fan-containing craters were only found between 18 and 29 S, and they all occur at around plus or minus 1 km of the MOLA-defined Martian datum. Within the study area they are not randomly distributed but instead form three distinct clusters. Fans typically descend greater than 1 km from where they disgorge from their alcoves. Longitudinal profiles show that their surfaces are very slightly concave with a mean slope of 2 degrees. Many fans exhibit very long, narrow low-relief ridges radially oriented down-slope, often branching at their distal ends, suggestive of distributaries. Morphometric data for 31 fans was derived from MOLA data and compared with terrestrial fans with high-relief source areas, terrestrial low gradient alluvial ramps in inactive tectonic settings, and older Martian alluvial ramps along crater floors. The Martian alluvial fans generally fall on the same trends as the terrestrial alluvial fans, whereas the gentler Martian crater floor ramps are similar in gradient to the low relief terrestrial alluvial surfaces. For a given fan gradient, Martian alluvial fans generally have greater source basin relief than terrestrial fans in active tectonic settings. This suggests that the terrestrial source basins either yield coarser debris or have higher sediment concentrations than their Martian counterpoints. Martian fans and Basin and Range fans have steeper gradients than the older Martian alluvial ramps and terrestrial low relief alluvial surfaces, which is consistent with a supply of coarse sediment. Martian fans are relatively large and of low gradient, similar to terrestrial fluvial fans rather than debris flow fans. However, gravity scaling uncertainties make the flow regime forming Martian fans uncertain. Martian fans, at least those in Holden crater, apparently formed around the time of the Noachian-Hesperian boundary. We infer that these fans formed during an episode of enhanced precipitation (probably snow) and runoff, which exhibited both sudden onset and termination.

Moore, Jeffrey M.↗
NASA NTRS2017

Mars2020 Entry, Descent, and Landing Instrumentation 2 (MEDLI2) Do No Harm Test Series

A total of seventeen instrumented thermal sensor plugs, eight pressure transducers, two heat flux sensors, and one radiometer are planned to be utilized on the Mars 2020 missions thermal protection system (TPS) as part of the Mars Entry, Descent, and Landing Instrumentation II (MEDLI2) project. Of the MEDLI2 instrumentation, eleven instrumented thermal plugs and seven pressure transducers will be installed on the heatshield of the Mars 2020 vehicle while the rest will be installed on the backshell. The goal of the MEDLI2 instrumentation is to directly inform the large performance uncertainties that contribute to the design and validation of a Mars entry system. A better understanding of the entry environment and TPS performance could lead to reduced design margins enabling a greater payload mass-fraction and smaller landing ellipses. To prove that the MEDLI2 system will not degrade the performance of the Mars 2020 TPS, an Aerothermal Do No Harm (DNH) test series was designed and conducted. Like Mars 2020s predecessor, Mars Science Laboratory (MSL), the heatshield material will be Phenolic Impregnated Carbon Ablator (PICA); the Mars 2020 entry conditions are enveloped by the MSL design environments, therefore the development and qualification testing performed during MEDLI is sufficient to show that the similar MEDLI2 heatshield instrumentation will not degrade PICA performance. However, given that MEDLI did not include any backshell instrumentation, the MEDLI2 team was required to design and execute a DNH test series utilizing the backshell TPS material (SLA-561V) with the intended flight sensor suite. To meet the requirements handed down from Mars 2020, the MEDLI2 DNH test series emphasized the interaction between the MEDLI2 sensors and sensing locations with the surrounding backshell TPS and substrucutre. These interactions were characterized by performing environmental testing of four 12 by 12 test panels, which mimicked the construction of the backshell TPS and the integration of the MEDLI2 sensors as seen in Figure 1. The testing included thermal vacuumcycling, random vibration, shock, and arc jet testing. The test panels were fabricated by Lockheed Martin, establishing techniques that will be utilized during the Mars 2020 vehicle installation. Each test panel included one thermal sensor plug (two embedded thermocouples), one heat flux sensor, and multiple pressure port holes for evaluation.This presentation will discuss the planning and execution of the MEDLI2 DNH test series. Selected highlights and results of each environmental test will be presented, and lessons learned will be addressed that will feed forward into the planning for the MEDLI2 flight system certification testing.

backshell↗
NASA NTRS2017

MELDI2 Do No Harm Test Series

Mars 2020 will fly the Mars Entry, Descent, and Landing Instrumentation II (MEDLI2) sensor suite consisting of a total of seventeen instrumented thermal sensor plugs, eight pressure transducers, two heat flux sensors, and one radiometer embedded in the thermal protection system (TPS). Of the MEDLI2 instrumentation, eleven instrumented thermal plugs and seven pressure transducers will be installed on the heatshield of the Mars 2020 vehicle while the rest will be installed on the backshell. The goal of the MEDLI2 instrumentation is to directly inform the large performance uncertainties that contribute to the design and validation of a Mars entry system. A better understanding of the entry environment and TPS performance could lead to reduced design margins enabling a greater payload mass-fraction and smaller landing ellipses. To prove that the MEDLI2 system will not degrade the performance of the Mars 2020 TPS, an Aerothermal Do No Harm (DNH) test series was designed and conducted. Like Mars 2020's predecessor, Mars Science Laboratory (MSL), the heatshield material will be Phenolic Impregnated Carbon Ablator (PICA); the Mars 2020 entry conditions are enveloped by the MSL design environments, therefore the development and qualification testing performed during MEDLI is sufficient to show that the similar MEDLI2 heatshield instrumentation will not degrade PICA performance. However, given that MEDLI did not include any backshell instrumentation, the MEDLI2 team was required to design and execute a DNH test series utilizing the backshell TPS material (SLA-561V) with the intended flight sensor suite. To meet the requirements handed down from Mars 2020, the MEDLI2 DNH test series emphasized the interaction between the MEDLI2 sensors and sensing locations with the surrounding backshell TPS and substrucutre. These interactions were characterized by performing environmental testing of four 12" by 12" test panels, which mimicked the construction of the backshell TPS and the integration of the MEDLI2 sensors as seen in Figure 1. The testing included thermal vacuum/ cycling, random vibration, shock, and arc jet testing. The test panels were fabricated by Lockheed Martin, establishing techniques that will be utilized during the Mars 2020 vehicle installation. Each test panel included one thermal sensor plug (two embedded thermocouples), one heat flux sensor, and multiple pressure port holes for evaluation. This presentation will discuss the planning and execution of the MEDLI2 DNH test series. Selected highlights and results of each environmental test will be presented, and lessons learned will be addressed that will feed forward into the planning for the MEDLI2 flight system certification testing.

Heatshield↗
NASA NTRS2020

Recommendations for Using Noise Monitors to Estimate Noise Exposure During X-59 Community Tests

A low fidelity simulation approach is used to explore how to place and use noise monitors during X-59 QueSST community tests, where people’s annoyance to the noise produced by the X-59 aircraft will be gathered. Several recommendations are provided including: 1) the desired number of sparsely spaced noise monitor sites within the survey area, 2) whether to group and average measurements across multiple noise monitors located at a site, 3) what spacing should be used if grouped noise monitors are used, 4) an approach to mitigate ambient noise contamination at the measurement sites, 5) a method to combine empirical and predicted dose estimates to provide a single dose estimate for respondents, and 6) assessing how changes in turbulence intensity and array configuration affect dose uncertainty. To make these recommendations, the error that is expected when fitting contrived, smoothly varying sonic boom “reference exposure surfaces” is studied when a spatially sparse and scattered set of samples is used as responses for the fit. The reference exposure surfaces mimic the sonic boom exposure at ground level that might be expected in the X-59 survey area in the absence of atmospheric turbulence, ambient noise, and other localized effects. The spatial extent of these surfaces varies and is representative of the different survey area sizes that might be expected during future X-59 community overflight tests. These contrived reference surfaces are sampled, and those reference samples are then perturbed to mimic atmospheric turbulence, ambient noise and other localized effects that might affect noise monitor measurements within overflown communities. Two different surface fitting methods are investigated when fitting these perturbed samples to approximate the reference surface. The first method uses interpolation between the perturbed data at the scattered sites to compute the fit. The second method fits a polynomial surface model to the perturbed data using ordinary least squares regression analysis. For both fitting methods, the root mean square fit error is computed from the pointwise difference between the fit surface and the reference surface as the count and configuration of the sites is varied while also averaging the error across many different realizations of both the smooth variation of the reference exposure surface and the random, localized perturbations at the sample sites. Different site configurations are compared using this error statistic to make the recommendations noted above. Additionally, the two fitting approaches (interpolation vs linear regression) are compared based on the fit error observed in these simulations. These analyses, comparisons, and recommendations should inform future decisions on the noise monitor placement and the methods used to analyze the noise monitor data that is collected during X-59 community overflights.

sonic boom↗
NASA NTRS2011

Thoughts Without Content are Empty, Intuitions Without Concepts are Blind - Determinism and Contingency Revisited

Was the emergence of life a predictable outcome of chemical evolution on earth? Could evolution produce life very different from ours? These are one of the oldest questions in the field of the origin of life that not only have broad philosophical implications but also impact how we approach the problem from the methodological standpoint. Framing the issue in terms of the dichotomy between contingency and determinism is not a fortunate because these two terms in their conventional meaning are neither mutually exclusive nor jointly exhaustive. Determinism, represented in natural sciences by Newtonian physics, relies on the assumption that every event is causally determined by a chain of previous events. In the context of the origin of life it means that once the initial conditions on the early earth have been specified further evolution follows inevitably. Considering uncertainties about conditions on the prebiotic earth, many plausible sets of initial conditions can be defined, each followed by a separate deterministic trajectory. This conventional understanding of determinism does not admit contingency. Further, it has no implications for evaluating how many sets of initial conditions lead to the emergence of life. It appears that a better framing of the problem is as follows: given plausible sets of initial conditions on the early earth how probable and broadly spread are evolutionary trajectories that lead to life? Instead of undertaking an impossible task of specifying microscopic initial conditions for all components of the system one uses a reduced representation of this system and specify only a small set of essential macroscopic parameters, values (or ranges of values) of which can be identified, inferred or estimated from experiment, theory or historical record. The following evolutionary trajectories are still governed by laws of physics and chemistry but become probabilistic and "contingency" is admitted as variations in other variables in the system. A similar reasoning is common in other fields of science, for example in statistical mechanics. Some trajectories lead to life, perhaps in different forms, whereas others do not. Of our true interest is the ratio of these two outcomes. The issue of determinism does not directly enter the picture. The debate about the likelihood of the emergence of life is quite old. One view holds that the origin of life is an event governed by chance, and the result of so many random events (contingencies) is unpredictable. This view was eloquently expressed by Monod. In his book "Chance or Necessity" he argued that life was a product of "nature's roulette." In an alternative view, expressed in particular by deDuve and Morowitz, the origin of life is considered a highly probable or even inevitable event (although its details need not be determined in every respect). Only in this sense the origin of life can be considered a "deterministic event".

Pohorille, Andrew↗
NASA NTRS2018

Snowex 2017 Community Snow Depth Measurements: A Quality-Controlled, Georeferenced Product

Snow depth was one of the core ground measurements required to validate remotely-sensed data collected during SnowEx Year 1, which occurred in Colorado. The use of a single, common protocol was fundamental to produce a community reference dataset of high quality. Most of the nearly 100 Grand Mesa and Senator Beck Basin SnowEx ground crew participants contributed to this crucial dataset during 6-25 February 2017. Snow depths were measured along ~300 m transects, whose locations were determined according to a random-stratified approach using snowfall and tree-density gradients. Two-person teams used snowmobiles, skis, or snowshoes to travel to staked transect locations and to conduct measurements. Depths were measured with a 1-cm incremented probe every 3 meters along transects. In shallow areas of Grand Mesa, depth measurements were also collected with GPS snow-depth probes (a.k.a. MagnaProbes) at ~1-m intervals. During summer 2017, all reference stake positions were surveyed with <10 cm accuracy to improve overall snow depth location accuracy. During the campaign, 193 transects were measured over three weeks at Grand Mesa and 40 were collected over two weeks in Senator Beck Basin, representing more than 27,000 depth values. Each day of the campaign depth measurements were written in waterproof field books and photographed by National Snow and Ice Data Center (NSIDC) participants. The data were later transcribed and prepared for extensive quality assessment and control. Common issues such as protocol errors (e.g., survey in reverse direction), notebook image issues (e.g., halo in the center of digitized picture), and data-entry errors (sloppy writing and transcription errors) were identified and fixed on a point-by-point basis. In addition, we strove to produce a georeferenced product of fine quality, so we calculated and interpolated coordinates for every depth measurement based on surveyed stakes and the number of measurements made per transect. The product has been submitted to NSIDC in csv format. To educate data users, we present the study design and processing steps that have improved the quality and usability of this product. Also, we will address measurement and design uncertainties, which are different in open vs. forest areas.

Brucker, L.↗
NASA NTRS1997

The Age of the Surface of Venus

Impact craters on Venus appear to be uniformly and randomly scattered over a once, but no longer, geologically active planet. To first approximation, the planet shows a single surface of a single age. Here we use Monte Carlo cratering simulations to estimate the age of the surface of Venus. The simulations are based on the present populations of Earth-approaching asteroids, Jupiter-family, Halley-family, and long period comets; they use standard Schmidt-Housen crater scalings in the gravity regime; and they describe interaction with the atmosphere using a semi-analytic 'pancake' model that is calibrated to detailed numerical simulations of impactors striking Venus. The lunar and terrestrial cratering records are also simulated. Both of these records suffer from poor statistics. The Moon has few young large craters and fewer still whose ages are known, and the record is biased because small craters tend to look old and large craters tend to look young. The craters of the Earth provide the only reliable ages, but these craters are few, eroded, of uncertain diameter, and statistically incomplete. Together the three cratering records can be inverted to constrain the flux of impacting bodies, crater diameters given impact parameters, and the calibration of atmospheric interactions. The surface age of Venus that results is relatively young. Alternatively, we can use our best estimates for these three input parameters to derive a best estimate for the age of the surface of Venus. Our tentative conclusions are that comets are unimportant, that the lunar and terrestrial crater records are both subject to strong biases, that there is no strong evidence for an increasing cratering flux in recent years, and that that the nominal age of the surface of Venus is about 600 Ma, although the uncertainty is about a factor of two. The chief difference between our estimate and earlier, somewhat younger estimates is that we find that the venusian atmosphere is less permeable to impacting bodies than supposed by earlier studies. An older surface increases the likelihood that Venus is dead.

Zahnle, K. J.↗
NASA NTRS2017

Nominal 30-M Cropland Extent Map of Continental Africa by Integrating Pixel-Based and Object-Based Algorithms Using Sentinel-2 and Landsat-8 Data on Google Earth Engine

A satellite-derived cropland extent map at high spatial resolution (30-m or better) is a must for food and water security analysis. Precise and accurate global cropland extent maps, indicating cropland and non-cropland areas, is a starting point to develop high-level products such as crop watering methods (irrigated or rainfed), cropping intensities (e.g., single, double, or continuous cropping), crop types, cropland fallows, as well as assessment of cropland productivity (productivity per unit of land), and crop water productivity (productivity per unit of water). Uncertainties associated with the cropland extent map have cascading effects on all higher-level cropland products. However, precise and accurate cropland extent maps at high spatial resolution over large areas (e.g., continents or the globe) are challenging to produce due to the small-holder dominant agricultural systems like those found in most of Africa and Asia. Cloud-based Geospatial computing platforms and multi-date, multi-sensor satellite image inventories on Google Earth Engine offer opportunities for mapping croplands with precision and accuracy over large areas that satisfy the requirements of broad range of applications. Such maps are expected to provide highly significant improvements compared to existing products, which tend to be coarser in resolution, and often fail to capture fragmented small-holder farms especially in regions with high dynamic change within and across years. To overcome these limitations, in this research we present an approach for cropland extent mapping at high spatial resolution (30-m or better) using the 10-day, 10 to 20-m, Sentinel-2 data in combination with 16-day, 30-m, Landsat-8 data on Google Earth Engine (GEE). First, nominal 30-m resolution satellite imagery composites were created from 36,924 scenes of Sentinel-2 and Landsat-8 images for the entire African continent in 2015-2016. These composites were generated using a median-mosaic of five bands (blue, green, red, near-infrared, NDVI) during each of the two periods (period 1: January-June 2016 and period 2: July-December 2015) plus a 30-m slope layer derived from the Shuttle Radar Topographic Mission (SRTM) elevation dataset. Second, we selected Cropland/Non-cropland training samples (sample size 9791) from various sources in GEE to create pixel-based classifications. As supervised classification algorithm, Random Forest (RF) was used as the primary classifier because of its efficiency, and when over-fitting issues of RF happened due to the noise of input training data, Support Vector Machine (SVM) was applied to compensate for such defects in specific areas. Third, the Recursive Hierarchical Segmentation (RHSeg) algorithm was employed to generate an object-oriented segmentation layer based on spectral and spatial properties from the same input data. This layer was merged with the pixel-based classification to improve segmentation accuracy. Accuracies of the merged 30-m crop extent product were computed using an error matrix approach in which 1754 independent validation samples were used. In addition, a comparison was performed with other available cropland maps as well as with LULC maps to show spatial similarity. Finally, the cropland area results derived from the map were compared with UN FAO statistics. The independent accuracy assessment showed a weighted overall accuracy of 94, with a producers accuracy of 85.9 (or omission error of 14.1), and users accuracy of 68.5 (commission error of 31.5) for the cropland class. The total net cropland area (TNCA) of Africa was estimated as 313 Mha for the nominal year 2015.

Cropland mapping; cropland areas; 30-m; Landsat-8;↗
NASA NTRSDate not supplied

Influence of Thermophysical Property Variability on Thermal Predictions of Ti-6Al-4V

The additive manufacturing (AM) industry has experienced rapid growth in recent decades as industrial interest in the process has grown. Because of this, there are many groups interested in simulations of the AM process. However, the influence of thermophysical property variability on the melt pool geometry during processing of various AM alloys is unclear. The goal of my work at NASA Langley Research Center (LaRC) was to characterize this influence on Ti-6Al-4V during AM processing alongside developing a tool to enable equivalent studies on other AM materials. To facilitate this process, a database tool was developed to contain and organize thermophysical data previously reported by primary sources. The specific thermophysical properties of interest for this study were density, specific heat capacity, and conductivity. Thermal diffusivity was calculated using the other three properties. This data was then fit with a Gaussian distribution and then randomly sampled from using Monte Carlo random value sampling. Using the Rosenthal equation, this sample data was used to simulate the temperature field during additive manufacturing and extract the melt pool geometry. Specifically, the melt pool’s length, width, and depth. This process was repeated an arbitrary number of times, with the default being 1000. Once this information was obtained, histograms were made showing the distributions of the sizes of the melt pool’s length, width, and depth. Representative statistical metrics of the distributions were calculated (mean, standard deviation, and coefficient of variance). This work found that at the simulated processing values, the melt pool’s width and depth had a 5.7% variation while the length had a 3.2% variation. Additionally, a tool was constructed to graph a thermal color map representation of the melt pool for specific, arbitrary values of density, specific heat, and conductivity. This tool allowed for more efficient plotting of individual queries of the thermophysical properties. The results that were observed in the research were that values reported in literature vary and this variance can have a significant impact on the simulation. Predictions of the melt pool’s dimensions show that length has a standard deviation of approximately 6 μm, width has a standard deviation of approximately 5 μm, and depth can vary by approximately 3 μm. Considering the mean sizes of length, width, and depth are 185 μm, 88 μm, and 44 μm respectively, such a deviation is significant. This shows that depth, for example, could be more than 10% larger or smaller than expected. This demonstrates that variability in the reported thermophysical properties are not negligible and should be expected to have an influence on the results of the laser powder bed fusion additive manufacturing process. When simulating this process in the future, measures should be taken to account for this discrepancy and the uncertainty involved

Justin Martin↗
NASA NTRS2015

Rao-Blackwellization for Adaptive Gaussian Sum Nonlinear Model Propagation

When dealing with imperfect data and general models of dynamic systems, the best estimate is always sought in the presence of uncertainty or unknown parameters. In many cases, as the first attempt, the Extended Kalman filter (EKF) provides sufficient solutions to handling issues arising from nonlinear and non-Gaussian estimation problems. But these issues may lead unacceptable performance and even divergence. In order to accurately capture the nonlinearities of most real-world dynamic systems, advanced filtering methods have been created to reduce filter divergence while enhancing performance. Approaches, such as Gaussian sum filtering, grid based Bayesian methods and particle filters are well-known examples of advanced methods used to represent and recursively reproduce an approximation to the state probability density function (pdf). Some of these filtering methods were conceptually developed years before their widespread uses were realized. Advanced nonlinear filtering methods currently benefit from the computing advancements in computational speeds, memory, and parallel processing. Grid based methods, multiple-model approaches and Gaussian sum filtering are numerical solutions that take advantage of different state coordinates or multiple-model methods that reduced the amount of approximations used. Choosing an efficient grid is very difficult for multi-dimensional state spaces, and oftentimes expensive computations must be done at each point. For the original Gaussian sum filter, a weighted sum of Gaussian density functions approximates the pdf but suffers at the update step for the individual component weight selections. In order to improve upon the original Gaussian sum filter, Ref. [2] introduces a weight update approach at the filter propagation stage instead of the measurement update stage. This weight update is performed by minimizing the integral square difference between the true forecast pdf and its Gaussian sum approximation. By adaptively updating each component weight during the nonlinear propagation stage an approximation of the true pdf can be successfully reconstructed. Particle filtering (PF) methods have gained popularity recently for solving nonlinear estimation problems due to their straightforward approach and the processing capabilities mentioned above. The basic concept behind PF is to represent any pdf as a set of random samples. As the number of samples increases, they will theoretically converge to the exact, equivalent representation of the desired pdf. When the estimated qth moment is needed, the samples are used for its construction allowing further analysis of the pdf characteristics. However, filter performance deteriorates as the dimension of the state vector increases. To overcome this problem Ref. [5] applies a marginalization technique for PF methods, decreasing complexity of the system to one linear and another nonlinear state estimation problem. The marginalization theory was originally developed by Rao and Blackwell independently. According to Ref. [6] it improves any given estimator under every convex loss function. The improvement comes from calculating a conditional expected value, often involving integrating out a supportive statistic. In other words, Rao-Blackwellization allows for smaller but separate computations to be carried out while reaching the main objective of the estimator. In the case of improving an estimator's variance, any supporting statistic can be removed and its variance determined. Next, any other information that dependents on the supporting statistic is found along with its respective variance. A new approach is developed here by utilizing the strengths of the adaptive Gaussian sum propagation in Ref. [2] and a marginalization approach used for PF methods found in Ref. [7]. In the following sections a modified filtering approach is presented based on a special state-space model within nonlinear systems to reduce the dimensionality of the optimization problem in Ref. [2]. First, the adaptive Gaussian sum propagation is explained and then the new marginalized adaptive Gaussian sum propagation is derived. Finally, an example simulation is presented.

state estimation↗