Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational inference”

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 1,207 records · Page 67

Development of a framework for sequential Bayesian design of experiments: Application to a pilot-scale solvent-based CO 2 capture process

In this paper, a methodology is developed for sequential design of experiments (SDoE) for process systems and applied to a solvent-based CO 2 capture system. In this approach, the prior knowledge of the system is used to prioritize process data collection at specific operating conditions. These data are then incorporated into a Bayesian inference methodology for updating a stochastic model by refining estimations of its underlying parameters, and the updated model is then used to generate the next set of test runs. Thus, the new knowledge obtained from the data is used to guide subsequent iterations of the experimental runs, ensuring that the overall data collection is maximally informative given that most experimental campaigns, especially at pilot or higher-scale plants, are costly, time-consuming, and resource-limited. The test run objective for this work was to minimize the maximum model prediction uncertainty for key output variables, but the methodology is generic and can be readily applied to other test run objectives. This methodology is applied to an aqueous monoethanolamine (MEA) pilot plant campaign at the National Carbon Capture Center (NCCC) in Wilsonville, Alabama, USA. The SDoE framework was utilized for two iterations, while collecting 18 sets of data representing different process conditions, and this resulted in an overall average reduction in uncertainty of approximately 50% in the prediction of CO 2 capture percentage. Moreover, 11 additional data sets were obtained with variation of absorber packing height for further model validation. This work shows the capability of the SDoE framework to maximize learning given limited resources, allowing for the reduction of model uncertainty, which is of great importance for many applications including reduction of technical risk associated with scale-up and economic analysis.

20 FOSSIL-FUELED POWER PLANTS↗

Polar Firn Properties in Greenland and Antarctica and Related Effects on Microwave Brightness Temperatures

In studying the mass balance of polar ice sheets, fluctuations of firn density near the surface is a major uncertainty. In this paper, we explore these variations at locations in the Greenland Ice Sheet and at the Dome C location in Antarctica. Borehole in situ measurements, Snow Radar echoes, microwave brightness temperatures, and modelling results from the Community firn model (CFM) are used. It is shown that firn density profiles can be represented using 3 processes: “long” and “short” length scale density variations and “refrozen layers”. Consistency with this description is observed in the dynamic range of airborne 0.5-2 GHz brightness temperatures and Snow Radar echo peaks in measurements performed in Greenland in 2017. Based on these insights, a new analytical partially coherent model is implemented to explain the microwave brightness temperatures using the three scale description of the firn. Short and long scale firn processes are modelled as a 3D continuous random medium with finite vertical and horizontal correlation lengths as opposed to past 1D random layered medium descriptions. Refrozen layers are described as deterministic sheets with planar interfaces, with the number of refrozen layer interfaces determined by radar observations. Firn density and correlation length parameters used in forward modelling to match measured 0.5-2 GHz brightness temperatures in Greenland show consistency with similar parameters in CFM predictions. Model predictions also are in good agreement with multi-angle 1.4 GHz vertically and horizontally polarized brightness temperature measured by the SMOS satellite at DOME C, Antarctica. This work shows that co-located active and passive microwave measurements can be used to infer polar firn properties that can be compared with predictions of the CFM. In particular, 0.5-2 GHz brightness temperature measurements are shown to be sensitive to long scale firn density fluctuations with density standard deviations in the range 0.01-0.06 g/cm^3 and vertical correlation lengths of 6-20 cm.

Haokui Xu↗

Unraveling the implications of finite specimen size on the interpretation of dynamic experiments for polycrystalline aluminum through direct numerical simulations

Normal and Pressure-shear plate impact (NPI and PSPI) experiments are popular experimental techniques for studying the mean-field macroscopic behavior of polycrystalline metals under high-rate dynamic loading. However, since both configurations rely upon geometry for subjecting the specimen to high strain rates, these experiments often involve a limited specimen size. Moreover, because of the inherent heterogeneities present within polycrystalline metals, it is difficult to ascertain if the size of the specimen and/or regions where measurements are made are sufficiently large for making representative inferences about the mean-field macroscopic properties from single-point velocity measurements. In the present study, we quantify the expected measurement variability on observable point measurements in NPI and PSPI experiments by carrying out direct numerical simulations (DNS) of statistically representative polycrystalline microstructures subjected to dynamic compression and compression-shear loading. In particular, we consider the role of specific material heterogeneities (e.g. the grain-to-grain difference in size, crystallographic orientation) on dispersion in the normal and transverse particle velocity records and on local fluctuations in key state variables (e.g. velocity, accumulated plastic strain) by incorporating these effects directly into a representative synthetic microstructure geometry and crystalline description of pure polycrystalline aluminum. The form of the present study is a large parametric investigation, consisting of ten ensembles of one hundred simulations. Each of the thousand simulations reflects a randomly realized synthetic microstructure in one of five cases of decreasing average grain size for the two loading configurations. Our analysis of the DNS results demonstrates that for both of these experimental configurations, the grain size directly correlates with the coefficient of variation (CV) in simulated point measurements, showing a convergent decrease in CV to zero (i.e. particle velocity record approaches the mean-field value) with decreasing grain size. Remarkably, the magnitude of variations in the particle velocity record is shown to be largest where the deviatoric stresses are most significant. In the case of NPI, this occurs at the elastic and plastic wavefront, whereas, in the case of PSPI, the magnitude of fluctuations are approximately constant throughout the experimental window time. The reasoning for the scatter in particle velocity due to the heterogeneous microstructure is demonstrated to be dependent on the mechanisms for accommodating deformation and on the interaction of reflection waves generated at sites of heterogeneities occurring at the scale of grains. Lastly, we develop a power-law description for the magnitude of scattering versus characteristic length, which provides a statistical framework for assessing the required number of grains per characteristic specimen dimension for minimizing scatter within these two experimental configurations (NPI, PSPI).

36 MATERIALS SCIENCE↗

The abundance and distribution of water vapor in Jupiter's atmosphere

The atmospheric transmission window between 1800 and 2250/cm in Jupiter's atmosphere was observed from the Kuiper Airborne Observatory and by the IR spectrometer (IRIS) on Voyager. The vertical distribution of H2O was derived for the 1-6 bar portion of Jupiter's troposphere. The spatial variation of H2O was measured using IRIS spectra of the Hot Spots in the North and South Equatorial Belts (NEB, SEB) and the Equatorial Zone and for an average of the North and South Tropical Zones. The H2O column abundance above the 4 bar level is the same in the zones as in the SEB Hot Spots, about 20 cm amagats. The NEB Hot Spots are desiccated by a factor of 3 with respect to the rest of Jupiter. For an average between -40 and +40 deg latitude, the H2O mole fraction, qH2O, is saturated for P less than 2 bars, qH2O = 4 millionths in the 2-4 bar range, and it increases to 3/100,000 at 6 bars. A similar vertical profile applies to the spatially resolved zone and belt spectra, except that H2O falls off more rapidly at P less than 4 bars in the NEB Hot Spots. A massive H2O cloud at 5 bars, T = 273 K is inconsistent with the observations. Instead, a thin H2O ice cloud would form at 2 bars, T = 200 K. The O/H ratio in Jupiter, inferred from H2O measurements in both belts and zones at 6 bars, is depleted by a factor of 50 with respect to the sun.

Bjoraker, Gordon L.↗

The footprints of visual attention in the Posner cueing paradigm revealed by classification images

In the Posner cueing paradigm, observers' performance in detecting a target is typically better in trials in which the target is present at the cued location than in trials in which the target appears at the uncued location. This effect can be explained in terms of a Bayesian observer where visual attention simply weights the information differently at the cued (attended) and uncued (unattended) locations without a change in the quality of processing at each location. Alternatively, it could also be explained in terms of visual attention changing the shape of the perceptual filter at the cued location. In this study, we use the classification image technique to compare the human perceptual filters at the cued and uncued locations in a contrast discrimination task. We did not find statistically significant differences between the shapes of the inferred perceptual filters across the two locations, nor did the observed differences account for the measured cueing effects in human observers. Instead, we found a difference in the magnitude of the classification images, supporting the idea that visual attention changes the weighting of information at the cued and uncued location, but does not change the quality of processing at each individual location.

Non-NASA Center↗

Inference of Ice Cloud Properties from High-spectral Resolution Infrared Observations

The theoretical basis is explored for inferring the microphysical properties of ice crystal from high-spectral resolution infrared observations. A radiative transfer model is employed to simulate spectral radiances to address relevant issues. The extinction and absorption efficiencies of individual ice crystals, assumed as hexagonal columns for large particles and droxtals for small particles, are computed from a combination of the finite- difference time-domain (FDTD) technique and a composite method. The corresponding phase functions are computed from a combination of FDTD and an improved geometric optics method (IGOM). Bulk scattering properties are derived by averaging the single- scattering properties of individual particles for 30 particle size distributions developed from in situ measurements and for additional four analytical Gamma size distributions for small particles. The non-sphericity of ice crystals is shown to have a significant impact on the radiative signatures in the infrared (IR) spectrum; the spherical particle approximation for inferring ice cloud properties may result in an overest&ation of the optical thickness and an inaccurate retrieval of effective particle size. Furthermore, we show that the error associated with the use of the Henyey-Greenstein phase function can be as larger as 1 K in terms of brightness temperature for larger particle effective size at some strong scattering wavenumbers. For small particles, the difference between the two phase functions is much less, with brightness temperatures generally differing by less than 0.4 K. The simulations undertaken in this study show that the slope of the IR brightness temperature spectrum between 790-960/cm is sensitive to the effective particle size. Furthermore, a strong sensitivity of IR brightness temperature to cloud optical thickness is noted within the l050-1250/cm region. Based on this spectral feature, a technique is presented for the simultaneous retrieval of the visible optical thickness and effective particle size from high spectral resolution infrared data under ice cloudy con&tion. The error analysis shows that the uncertainty of the retrieved optical thickness and effective particle size has a small range of variation. The error for retrieving particle size in conjunction with an uncertainty of 5 K in cloud'temperature, or a surface temperature uncertainty of 2.5 K, is less than 15%. The corresponding e m r in the uncertainty of optical thickness is within 5-2096, depending on the value of cloud optical thickness. The applicability of the technique is demonstrated using the aircraft-based High- resolution Interferometer Sounder (HIS) data from the Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) in 1996 and the First ISCCP Regional Experiment - Arctic Clouds Experiment (FIRE-ACE) in 1998.

Huang, Hung-Lung↗

Convective and stratiform rain: Multichannel microwave sensing over oceans

Measurements made by the Special Sensor Microwave/Imager (SSM/I) radiometer over the oceans, at 19, 37, and 85 GHz in dual polarization, are used to develop a model to classify rain into light-stratiform, moderately convective, and heavy convective types in the mesoscale convective systems (MCS). It is observed that the bulk of the 19- and 37-GHz data are linearly correlated with respect to one another, and generally increase together in brightness as the mean rain rate in the field of view (FOV) of the radiometer increases. However, a significant fraction of the data from these channels departs from this linear relationship, reflecting the nonuniform rain that is convective vs. the relatively light stratiform rain. It is inferred from the SSM/I data, in a MCS, when the slope dT sub 3/dT sub 19 is greater than unity there are optically thin clouds which produce light uniform rain. On the other hand, when dT sub 3/dT sub 19 is close to unity, the rain cells have an open structure and correspond to the convective type of rain. The openings between the cells are apparently a result of the downdrafts and/or entrainment. Relatively low values of 85-GHz brightness temperatures that are present when dT sub 37/dT sub 19 is close to unity support these views and, in addition, leads us to conclude that when the convection is heavy this brightness temperature decreases due to scattering by hydrometeors. On the basis of this explanation of the SSM/I data, an empirical rain retrieval algorithm is developed. Radar backscatter observations over the Atlantic Ocean next to Florida are used to demonstrate the applicability of this method. Three monthly mean maps of rainfall over the oceans from 50 degrees N to 50 degrees S, are presented to illustrate the ability of this method to sense seasonal and interannual variations of rain.

Prabhakara, C.↗

Oxygen and Magnesium Isotopic Compositions of Asteroidal Materials Returned from Itokawa by the Hayabusa Mission

The Hayabusa spacecraft made two touchdowns on the surface of Asteroid 25143 Itokawa on November 20th and 26th, 2005. The Asteroid 25143 Itokawa is classified as an S-type asteroid and inferred to consist of materials similar to ordinary chondrites or primitive achondrites [1]. Near-infrared spectroscopy by the Hayabusa spacecraft proposed that the surface of this body has an olivine-rich mineral assemblage potentially similar to that of LL5 or LL6 chondrites with different degrees of space weathering [2]. The spacecraft made the reentry into the Earth s atmosphere on June 12th, 2010 and the sample capsule was successfully recovered in Australia on June 13th, 2010. Although the sample collection processes on the Itokawa surface had not been made by the designed operations, more than 1,500 grains were identified as rocky particles in the sample curation facility of JAXA, and most of them were judged to be of extraterrestrial origin, and definitely from Asteroid Itokawa on November 17th, 2010 [3]. Although their sizes are mostly less than 10 microns, some larger grains of about 100 microns or larger were also included. The mineral assembly is olivine, pyroxene, plagioclase, iron sulfide and iron metal. The mean mineral compositions are consistent with the results of near-infrared spectroscopy from Hayabusa spacecraft [2], but the variations suggest that the petrologic type may be smaller than the spectroscopic results. Several tens of grains of relatively large sizes among the 1,500 grains will be selected by the Hayabusa sample curation team for preliminary examination [4]. Each grain will be subjected to one set of preliminary examinations, i.e., micro-tomography, XRD, XRF, TEM, SEM, EPMA and SIMS in this sequence. The preliminary examination will start from the last week of January 2011. Therefore, samples for isotope analyses in this study will start from the last week of February 2011. By the time of the LPSC meeting we will have measured the oxygen and magnesium isotopic composition of several grains. We will present the first results from the isotope analyses that will have been performed.

Yurimoto, H↗

Bladed Terrain on Pluto: Possible Origins and Evolutions

Pluto’s Bladed Terrain (centered roughly 20°N, 225°E) covers the flanks and crests of the informally named Tartarus Dorsa with numerous roughly aligned blade-like ridges oriented ~North-South; it may also stretch considerably farther east onto the non-close approach hemisphere but that inference is tentative. Individual ridges are typically several hundred meters high, and are spaced 5 to 10 km crest to crest, separated by Vshaped valleys. Many ridges merge at acute angles to form Y-shape junctions in plan view. The principle composition of the blades themselves we suspect is methane or a methane-rich mixture. (Methane is spectroscopically strongly observed on the optical surfaces of blades.) Nitrogen ice is very probably too soft to support their topography. Cemented mixtures of volatile and non-volatile ices may also provide a degradable but relief supporting “bedrock” for the blades, perhaps analogous to Callisto. Currently we are considering several hypotheses for the origins of the deposit from which Bladed Terrain has evolved, including aeolian disposition, atmospheric condensation, updoming and exhumation, volcanic intrusions or extrusions, crystal growth, among others. We are reviewing several processes as candidate creators or sculptors of the blades. Perhaps they are primary depositional patterns such as dunes, or differential condensation patterns (like on Callisto), or fissure extrusions. Or alternatively perhaps they are the consequence of differential erosion (such as sublimation erosion widening and deepening along cracks), variations in substrate properties, mass wasting into the subsurface, or sculpted by a combination of directional winds and solar isolation orientation. We will consider the roles of the long-term increasing solar flux and short periods of warm thick atmospheres. Hypotheses will be ordered based on observational constrains and modeling to be presented at the conference.

Nitrogen ice↗

Bladed Terrain on Pluto: Possible Origins and Evolutions

Pluto's Bladed Terrain (centered roughly 20 deg N, 225 deg E) covers the flanks and crests of the informally named Tartarus Dorsa with numerous roughly aligned blade-like ridges oriented approx. North-South; it may also stretch considerably farther east onto the non-close approach hemisphere but that inference is tentative. Individual ridges are typically several hundred meters high, and are spaced 5 to 10 km crest to crest, separated by V-shaped valleys. Many ridges merge at acute angles to form Y-shape junctions in plan view. The principle composition of the blades themselves we suspect is methane or a methane-rich mixture. (Methane is spectroscopically strongly observed on the optical surfaces of blades.) Nitrogen ice is very probably too soft to support their topography. Cemented mixtures of volatile and non-volatile ices may also provide a degradable but relief supporting "bedrock" for the blades, perhaps analogous to Callisto. Currently we are considering several hypotheses for the origins of the deposit from which Bladed Terrain has evolved, including aeolian disposition, atmospheric condensation, updoming and exhumation, volcanic intrusions or extrusions, crystal growth, among others. We are reviewing several processes as candidate creators or sculptors of the blades. Perhaps they are primary depositional patterns such as dunes, or differential condensation patterns (like on Callisto), or fissure extrusions. Or alternatively perhaps they are the consequence of differential erosion (such as sublimation erosion widening and deepening along cracks), variations in substrate properties, mass wasting into the subsurface, or sculpted by a combination of directional winds and solar isolation orientation. We will consider the roles of the long-term increasing solar flux and short periods of warm thick atmospheres. Hypotheses will be ordered based on observational constrains and modeling to be presented at the conference.

methane↗

The abundance and distribution of water vapor in Jupiter's atmosphere

The atmospheric transmission window between 1800 and 2250 cm(-1) in Jupiter's atmosphere was observed from the Kuiper Airborne Observatory (KAO) and by the infrared spectrometer (IRIS) on Voyager. The vertical distribution of H2O was derived for the 1 to 6 bar portion of Jupiter's troposphere. The spatial variation of H2O was measured using IRIS spectra of the Hot Spots in the North and South Equatorial Belts, the Equatorial Zone, and for an average of the North and South Tropical Zones. The H2O column abundance above the 4 bar level is the same in the zones as in the SEB Hot Spots, about 20 cm-amagat. The NEB Hot Spots are desiccated by a factor of 3 with respect to the rest of Jupiter. For an average between -40 to 40 deg latitude, the H2O mole fraction, qH2O, is saturated for P less than 2 bars, qH2O = 4x10 to the -6 in the 2 to 4 bar range and it increases to 3x10 to the -5 at 6 bars. A similar vertical profile applies to the spatially resolved zone and belt spectra, except that H2O falls off more rapidly at P less than 4 bars in the NEB Hot Spots. The massive H2O cloud at 5 bars, T = 273 K, proposed in solar composition models, is inconsistent with the observations. Instead, a thin H2O ice cloud would form at 2 bars, T = 200 K. The O/H ratio in Jupiter, inferred from H2O measurements in both belts and zones at 6 bars, is depleted by a factor of 50 with respect to the Sun. The implications for the origin of Jupiter of globally depleted O/H, but enhanced C/H and N/H, are discussed.

Bjoraker, Gordon L.↗

Internal Microstructure Dictates Interactions of Polymer-grafted Nanoparticles in Solution

Understanding the effects of polymer brush architecture on particle interactions in solution is requisite to enable the development of functional materials based on self-assembled polymer-grafted nanoparticles (GNPs). Static and dynamic light scattering of polystyrene-grafted silica particle solutions in toluene reveals that the pair interaction potential, inferred from the second virial coefficient, A 2 , is strongly affected by the grafting density, σ, and degree of polymerization, N, of tethered chains. In the limit of intermediate σ (~0.3 to 0.6 nm -2 ) and high N, A 2 is positive and increases with N. This confirms the good solvent conditions and can be qualitatively rationalized on the basis of a pair interaction potential derived for grafted (brush) particles. In contrast, for high σ > 0.6 nm -2 and low N, A 2 displays an unexpected reversal to negative values, thus indicating poor solvent conditions. These findings are rationalized by means of a simple analysis based on a coarse-grained brush potential, which balances the attractive core-core interactions and the excluded volume interactions imparted by the polymer grafts. The results suggest that the steric crowding of polymer ligands in dense GNP systems may fundamentally alter the interactions between brush particles in solution and highlight the crucial role of architecture (internal microstructure) on the behavior of hybrid materials. The effect of grafting density also illustrates the opportunity to tailor the physical properties of hybrid materials by altering geometry (or architecture) rather than a variation of the chemical composition.

36 MATERIALS SCIENCE↗

Comparing Two Approaches for Assessing Observation Impact

Langland and Baker introduced an approach to assess the impact of observations on the forecasts. In that approach, a state-space aspect of the forecast is defined and a procedure is derived ultimately relating changes in the aspect with changes in the observing system. Some features of the state-space approach are to be noted: the typical choice of forecast aspect is rather subjective and leads to incomplete assessment of the observing system, it requires availability of a verification state that is in practice correlated with the forecast, and it involves the adjoint operator of the entire data assimilation system and is thus constrained by the validity of this operator. This article revisits the topic of observation impacts from the perspective of estimation theory. An observation-space metric is used to allow inferring observation impact on the forecasts without the limitations just mentioned. Using differences of observation-minus-forecast residuals obtained from consecutive forecasts leads to the following advantages: (i) it suggests a rather natural choice of forecast aspect that directly links to the data assimilation procedure, (ii) it avoids introducing undesirable correlations in the forecast aspect since verification is done against the observations, and (iii) it does not involve linearization and use of adjoints. The observation-space approach has the additional advantage of being nearly cost free and very simple to implement. In its simplest form it reduces to evaluating the statistics of observationminus- background and observation-minus-analysis residuals with traditional methods. Illustrations comparing the approaches are given using the NASA Goddard Earth Observing System.

Variational Analysis↗

Spaceborne Lidar Retrievals of PM2.5 for Air Quality Studies and Applications

Fine particulate matter (PM2.5) substantially contributes to air pollution and negatively affects human health. While many studies have investigated the use of passive column-integrated aerosol optical depth to infer surface PM2.5, the use of lidar observations for air quality characterization is not nearly as extensive. Lidar measurements are critical, however, due to the vertical aerosol information they provide, including near the surface. In this presentation, we first provide an overview of various lidar-based approaches for estimating PM2.5 concentrations and then discuss how lidar measurements can assist other air quality applications. For example, estimates of PM2.5 have been obtained in a physics-based approach through CALIOP near-surface aerosol extinction retrievals, assumptions on the mass extinction efficiency, and incorporating other parameters (an aerosol hygroscopic growth factor and PM2.5/PM10 ratio). Application of this algorithm over the contiguous United States (CONUS) from 2006 to 2018 yielded larger PM2.5 values over the eastern and western CONUS (~10-15 μg/m³) and lower PM2.5 levels in the central CONUS (~5 μg/m³). These spatial patterns were similar to those from gridded PM2.5 concentrations obtained through in situ measurements at ground stations operated by the US Environmental Protection Agency. In another approach, the Cloud Aerosol Transport System (CATS) lidar was used with the Goddard Earth Observing System (GEOS) model in a 1D ensemble-based variational technique to obtain PM2.5 over the US and Europe, and the spatial patterns of the CATS/GEOS based PM2.5 concentrations generally captured those from surface stations (with corresponding hourly EPA PM2.5 vs CATS PM2.5 statistics of R=0.4 and bias=1.5 μg/m³). In our recent work, as part of the Models, In situ, and Remote sensing of Aerosols (MIRA) Working Group, we have applied both the CALIOP and CATS/GEOS based approaches over the highly polluted country of India during the post-monsoon season (September-October 2016). We derived elevated levels of two-month mean PM2.5 (~100 μg/m³) in northern India, especially near New Delhi. These high PM2.5 concentrations in the Indo-Gangetic plain are driven in large part from the seasonal burning of crop residue and meteorological conditions typical at this time of the year, such as low wind speeds and a shallow boundary layer. While the satellite-derived PM2.5 moderately replicates (R = ~0.7-0.9) the spatial variability in the two-month mean of surface in situ PM2.5 from monitoring sites operated by the Central and State Pollution Control Boards, we show results from specific scenes for which there are large deviations between the satellite-derived PM2.5 and in situ measurements. Other current work on this topic focuses on developing PM2.5 estimates using airborne high spectral resolution lidar measurements through machine learning regression algorithms and involves several parameters (e.g., aerosol extinction, color ratio, lidar ratio). Application of this method over major metropolitan areas in the US and Asia have resulted in high correlations (R = 0.93) with surface measurements. This airborne lidar approach can be adapted to spaceborne lidar measurements, and all three of these approaches can be applied to ESA’s EarthCARE Atmospheric Lidar instrument, setting the stage for the future Cloud Aerosol Lidar for Global Scale Observations of the Ocean-Land Atmosphere System (CALIGOLA) mission. Ultimately, beyond estimates of PM2.5, the aerosol vertical distribution from lidars can benefit studies involving passive sensor approaches for PM2.5 proxies (including from geostationary satellites), wildfire smoke plume injection heights, volcanic emissions (e.g., ash height retrievals), and aerosol/air quality model assimilation, evaluation, and forecasts.

Travis D Toth↗

Inferences About the Early Moon From Gravity and Topography

Recent spacecraft missions to the Moon have significantly improved our knowledge of the lunar gravity and topography fields, and have raised some new and old questions about the early lunar history. It has frequently been assumed that the shape of the Moon today reflects an earlier equilibrium state and that the Moon has retained some internal strength. Recent analysis indicating a superisostatic state of some lunar basins lends support to this hypothesis. On its simplest level, the present shape of the Moon is slightly flattened by 2.2 +/- 0.2 km while its gravity field, represented by an equipotential surface, is flattened only about 0.5 km. The hydrostatic component to the flattening arising from the Moon's present day rotation contributes only 7 m. This difference between the topographic shape of the MOon and the shape of its gravitational equipotential has frequently been explained as the "memory" of an earlier moon that was rotating faster and had a correspondingly larger hydrostatic flattening. To obtain this amount of hydrostatic flattening from rotation alone, and accounting for the contribution of the present-day gravity field, the Moon's rotation rate would need to be about 15x greater than at present, leading ot a period of < 2 days. Maintaining its synchronous rotation with Earth would require a radius for the Moon's orbit of approximately 9 Earth Radii. Unfortunately, our confidence in the observed lunar flattening is not as great as we would like. The uncertainty of .02 km may not properly reflect the limitations of the Clementine dataset, which did not sample poleward of latitudes 81 N and 79 S. Also, the large variation of topography +/- 8 km seen on the MOon dwarfs our estimate fo the flattening. Further the lunar south pole is on the edge of, or possibly inside the massive deep, South Pole-Aitken Basin. Thus, polar radii could be underestimated. This would yield a smaller flattening, which would imply a greater lunar rotation period and orbital radius. However, Basin compensation states and analyses of support and relaxation of topography at long wavelengths point to a lunar shape that has retained a flattening from an earlier faster rotation period.

Smith, D. E.↗

Surface property variations in Venusian fluidized ejecta blanket craters

A comprehensive study of Magellan Cycles 1 and 2 radar data from Venus reveals surface roughness and dielectric variations associated with fluidized ejecta blanket (FEB) craters that help illuminate styles of flow ejecta emplacement. This study develops new procedures of digital unit mapping and polygon-filling algorithms using Magellan synthetic aperture radar (SAR), altimetry, and radiometry data. These techniques allow the extraction of radiophysical information for FEB crater materials, nearby plains, and lava flows. Backscatter curve slopes of the FEBs studied here are consistent with surface textures that are transitional between a'a and pahoehoe-like. Average surface property values of ejecta units are relatively similar for a given crater, but are discernibly different from other craters. Individual crater ejecta reflectivity and emissivity values are relatively similar to those for the surrounding plains, which may suggest a link between plains material and ejecta dielectric properties. Increasing FEB roughness downflow are interpreted to be associated with more lava-like flows, while decreasing roughness are more similar to trends typical of gravity (pyroclastic-like or debris-like) flows. Most commonly, FEB crater flow materials exhibit transitions from proximal, lava/melt-like flow styles to distal, gravity flow-like styles. Some FEBs show more complicated behavior, however, or appear to be more dominated by dielectric differences downflow, as inferred from correlations between the data sets. Such transitions may result from changes in local topography or from overlapping of flow lobes during FEB emplacement.

Johnson, Jeffrey R.↗

Detection of Ice Polar Stratospheric Clouds from Assimilation of Atmospheric Infrared Sounder Data

A novel technique is presented for detection of ice polar stratospheric clouds (PSCs) that form at extremely low temperatures in the lower polar stratosphere during winter. Temperature is a major factor in determining abundance of PSCs, which in turn provide surfaces for heterogeneous chemical reactions leading to ozone loss and radiative cooling. The technique infers the presence of ice PSCs using radiances from the Atmospheric Infrared Sounder (AIRS) in the Goddard Earth Observing System version 5 (GEOS-5) data assimilation system. Brightness temperatures are computed from short-term GEOS-5 forecasts for several hundred AIRS channels, using a radiation transfer module. The differences between collocated AIRS observations and these computed values are the observed-minus-forecast (O-F) residuals in the assimilation system. Because the radiation model assumes clear-sky conditions, we hypothesize that these O-F residuals contain quantitative information about PSCs. This is confirmed using sparse data from the Polar Ozone and Aerosol Measurement (POAM) III occultation instrument. The analysis focuses on 0-F residuals for the 6.79pm AIRS moisture channel. At coincident locations, when POAM III detects ice clouds, the AIRS O-F residuals for this channel are lower than -2K. When no ice PSCs are evident in POAM III data, the AIRS 0-F residuals are larger. Given this relationship, the high spatial density of AIRS data is used to construct maps of regions where 0-F residuals are lower than -2K, as a proxy for ice PSCs. The spatial scales and spatio-temporal variations of these PSCs in the Antarctic and Arctic are discussed on the basis of these maps.

Stajner, Ivanka↗

Orbital Debris Shape Effect Investigations for Mitigating Risk

NASA’s Orbital Debris Program Office (ODPO) and Hypervelocity Impact Technology (HVIT) team have coordinated to better understand the risks to upper stages and spacecraft from non-spherical orbital debris. It is well understood that fragmentation (collision or explosion) events in orbit produce fragments of various materials, sizes, and shapes. To further characterize these parameters, the ODPO is developing the next-generation Orbital Debris Engineering Model (ORDEM) version 4.0 to include orbital debris shape distributions. Ground-based assets, such as radar and optical sensors, can provide size estimates and some insight into material based on radar return or optical filter photometry/spectroscopy, respectively. Characterizing an object’s shape requires more laboratory analyses to infer how shape affects these measurements. More importantly, in addition to size and material/density, the shape of fragments in orbit will alter the ballistic limit equations used in orbital debris risk assessments with NASA’s Bumper Code. The ODPO plans to release ORDEM 4.0 in the coming years. Performing ground-based laboratory impact tests on high-fidelity spacecraft mockups provides the means to directly measure size, mass, material/density, and shape of fragments, all key parameters needed to characterize real-world break up events. The DebriSat test, the results of which are provided, showcases the details of this type of experiment. The goal of this collaborative research between the ODPO and the HVIT team is to include a shape parameter in the environmental and breakup models used to assess risk for various space structures. This paper examines ground-based laboratory impact tests and the associated fragment shape categories. Provided these defined shapes, the approach is simplified by assuming a right circular cylinder (RCC) approximation with varying length-to-diameter ratios. Highlights of impact tests conducted by the HVIT team using non-spherical projectiles based on the RCC approximation are presented. Hydrocode simulations have also been performed to expand on the complexity of variations with non-spherical projectiles. Lastly, ray-tracing simulations of various RCCs of known material are provided to support the ongoing research on optical reflectance distributions with known shapes and to highlight how this may modify the current optical size estimation model. The status and plan forward are outlined for NASA's orbital debris shape effect investigation using a multidisciplinary approach by the ODPO and the HVIT team.

Heather Cowardin↗