SEARCH · Engineering Papers

Results for “distributed 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,225 records · Page 68

All sources NASA NTRS DOE OSTI

NASA NTRS2021

Contribution of the World's Main Dust Source Regions to the Global Cycle of Desert Dust

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, the relative contributions of the world's major source regions to the global dust cycle remain poorly constrained. This problem hinders accounting for the potentially large impact of regional differences in dust properties on clouds, the Earth's energy balance, and terrestrial and marine biogeochemical cycles. Here, we constrain the contribution of each of the world's main dust source regions to the global dust cycle. We use an analytical framework that integrates an ensemble of global aerosol model simulations with observationally informed constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth (DAOD). We obtain a dataset that constrains the relative contribution of nine major source regions to size-resolved dust emission, atmospheric loading, DAOD, concentration, and deposition flux. We find that the 22–29 Tg (1 standard error range) global loading of dust with a geometric diameter up to 20 µm is partitioned as follows: North African source regions contribute ∼ 50 % (11–15 Tg), Asian source regions contribute ∼ 40 % (8–13 Tg), and North American and Southern Hemisphere regions contribute ∼ 10 % (1.8–3.2 Tg). These results suggest that current models on average overestimate the contribution of North African sources to atmospheric dust loading at ∼ 65 %, while underestimating the contribution of Asian dust at ∼ 30 %. Our results further show that each source region's dust loading peaks in local spring and summer, which is partially driven by increased dust lifetime in those seasons. We also quantify the dust deposition flux to the Amazon rainforest to be ∼ 10 Tg yr−1, which is a factor of 2–3 less than inferred from satellite data by previous work that likely overestimated dust deposition by underestimating the dust mass extinction efficiency. The data obtained in this paper can be used to obtain improved constraints on dust impacts on clouds, climate, biogeochemical cycles, and other parts of the Earth system.

Dust↗

NASA NTRS2012

The Role of Lightning in Controlling Interannual Variability of Tropical Tropospheric Ozone and OH and its Implications for Climate

Nitrogen oxides (NO(x) = NO + NO2) produced by lightning make a major contribution to the production of the dominant tropospheric oxidants (OH and ozone). These oxidants control the lifetime of many trace gases including long-lived greenhouse gases, and control the source-receptor relationship of inter-hemispheric pollutant transport. Lightning is affected by meteorological variability, and therefore represents a potentially important tropospheric chemistry-climate feedback. Understanding how interannual variability (IAV) in lightning affects IAV in ozone and OH in the recent past is important if we are to predict how oxidant levels may change in a future warmer climate. However, lightning parameterizations for chemical transport models (CTMs) show low skill in reproducing even climatological distributions of flash rates from the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) satellite instruments. We present an optimized regional scaling algorithm for CTMs that enables sufficient sampling of spatiotemporally sparse satellite lightning data from LIS to constrain the spatial, seasonal, and interannual variability of tropical lightning. We construct a monthly time series of lightning flash rates for 1998-2010 and 35degS-35degN, and find a correlation of IAV in total tropical lightning with El Nino. We use the IAV-constraint to drive a 9-year hindcast (1998-2006) of the GEOS-Chem 3D chemical transport model, and find the increased IAV in LNO(x) drives increased IAV in ozone and OH, improving the model fs ability to simulate both. Although lightning contributes more than any other emission source to IAV in ozone, we find ozone more sensitive to meteorology, particularly convective transport. However, we find IAV in OH to be highly sensitive to lightning NO(x), and the constraint improves the ability of the model to capture the temporal behavior of OH anomalies inferred from observations of methyl chloroform and other gases. The sensitivity of OH is explained using photochemical reaction rates which show a "magnification" effect of the initial lightning NO perturbation on OH primary production, HO(x) recycling, and OH loss frequencies. This influence on OH may represent a negative feedback, if lightning increases in a warming world..

Murray, Lee T.↗

NASA NTRS2004

Physical Characterization of the Near-Earth Object Population

Many pieces of the puzzle must be brought together in order to have a clear picture of the near-Earth object (NEO) population. Four of the pieces that can be described include: i) the taxonomic distribution of the population as measured by observational sampling, ii) the determination of albedos that can be associated with the taxonomic distribution, iii) discovery statistics for the NE0 population, and iv) the debiasing of the discovery statistics using the taxonomic and albedo information. Support from this grant enables us to address three of these four pieces. Binzel et al. (2004, submitted) presents the first piece, detailing the observations and observed characteristics of the NE0 and Mars-crossing (MC) population. For the second piece, a complementary program of albedo measurements is pursued at the Keck Observatory (Binzel, P. I.) with first results published in Delbo et al. (2003). For the third piece, the most extensive NE0 discovery statistics are provided by the LINEAR survey. Binzel has supervised the MIT Ph. D. thesis work of Stuart (2003) to bring the fourth piece, submitted for publication by Stuart and Binzel (2004). Our results provide new constraints for the NE0 population and progress for the Spaceguard Survey, illuminate asteroid and comet source regions for the NEOs, and provide new evidence for space weathering processes linking asteroids and meteorites. Further, we are identifying top priority near-Earth spacecraft mission candidates based on their spectral properties and inferred compositions.

Binzel, Richard P.↗

NASA NTRS1997

Pioneer 10 and 11 Data Analysis

This report finishes the work of NASA Grant NAS2-153, which supported data analysis for the UCSD instruments on Pioneers 10 and 11. The data analyzed under this grant span 22 years of interplanetary measurements in the inner and outer heliosphere. The UCSD instruments made their mark in cosmic ray research based upon their high energy thresholds, directional responses, and reliable data streams. one of their primary scientific objectives concerns the size, configuration, and time behavior of the heliosphere. The size scale is inferred from the radial intensity gradient, which is measured between the two spacecraft and extrapolated to interstellar intensity levels at the cosmic ray modulation boundary. This boundary still eludes us, and its position, motion, and the best method of extrapolation are ongoing problems. Current projections place the boundary beyond 100 AU, which may be beyond the termination shock, and raises the question of possible modulation in the heliosheath. Probably our only hope of seeing this region in the immediate future rides on the possibility that the boundary will move inward. Our instruments have recorded many Forbush, or transient, decreases in the outer heliosphere. These observations led us to a model that attributes many of the decreases to solar wind stream-stream interactions, and relates the cosmic ray variations to the locally observed magnetic field magnitude. As the cosmic ray variations in this model result only from topological changes in the modulation integral, the model is a tool for studying the possibility that the 11 year cosmic ray modulation cycle can be accounted for by a superposition of Forbush decreases. The cosmic ray angular distribution function is measurable, given a good telemetry rate, by the UCSD Cerenkov detector which counts particles of energy greater than 500 MeV/n. We obtained statistically significant samples from 1 to 9 AU, at 13 AU, and at 34 AU. The anisotropy tends to be a few tenths of a per cent at all radial distances. A quasiperiodic variation in the east-west anisotropy with period of about 50 days remains unexplained.

Fillius, Walker↗

NASA NTRS2006

Formaldehyde Distribution over North America: Implications for Satellite Retrievals of Formaldehyde Columns and Isoprene Emission

Formaldehyde (HCHO) columns measured from space provide constraints on emissions of volatile organic compounds (VOCs). Quantitative interpretation requires characterization of errors in HCHO column retrievals and relating these columns to VOC emissions. Retrieval error is mainly in the air mass factor (AMF) which relates fitted backscattered radiances to vertical columns and requires external information on HCHO, aerosols, and clouds. Here we use aircraft data collected over North America and the Atlantic to determine the local relationships between HCHO columns and VOC emissions, calculate AMFs for HCHO retrievals, assess the errors in deriving AMFs with a chemical transport model (GEOS-Chem), and draw conclusions regarding space-based mapping of VOC emissions. We show that isoprene drives observed HCHO column variability over North America; HCHO column data from space can thus be used effectively as a proxy for isoprene emission. From observed HCHO and isoprene profiles we find an HCHO molar yield from isoprene oxidation of 1.6 +/- 0.5, consistent with current chemical mechanisms. Clouds are the primary error source in the AMF calculation; errors in the HCHO vertical profile and aerosols have comparatively little effect. The mean bias and 1Q uncertainty in the GEOS-Chem AMF calculation increase from <1% and 15% for clear skies to 17% and 24% for half-cloudy scenes. With fitting errors, this gives an overall 1 Q error in HCHO satellite measurements of 25-31%. Retrieval errors, combined with uncertainties in the HCHO yield from isoprene oxidation, result in a 40% (1sigma) error in inferring isoprene emissions from HCHO satellite measurements.

Millet, Dylan B.↗

NASA NTRS1999

A TRMM Microwave Radiometer Rain Rate Estimation Method with Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima (T85min) in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of surface rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85(sub min), produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer empirically from TMI data the presence of three different kinds of thunderstorms or Cbs. These Cbs are classified as young, mature, and decaying types, and are assumed to have a scale of about 20 km on the average. Two parameters are used to classify these three kinds of Cbs based on the T85 data: a) the magnitude of scattering depression deduced from local T85(sub min) and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85min to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (>= 20 mm/hr) rain rates are retrieved satisfactorily. Accuracy in the estimates of the light, moderate, and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15 %. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗

NASA NTRS2005

Short Gamma-Ray Bursts with Extended Emission

The recent association of several short gamma-ray bursts (GRBs) with early type galaxies with low star formation rate demonstrates that short bursts arise from a different progenitor mechanism than long bursts. However, since the duration distributions of the two classes overlap, membership is not always easily established. The picture is complicated by the occasional presence of softer, extended emission lasting tens of seconds after the initial spike- like emission comprising an otherwise short burst. Using the large BATSE sample with time-tagged event (TTE) data, we show that the fundamental defining characteristic of the short burst class is that the initial spike exhibits negligible spectral evolution at energies above approx. 25 keV. This is behavior is nearly ubiquitous for the 260 bursts with T(sub 90) less than 2s where the BATSE TTE data type completely included the initial spike: Their spectral lags measured between the 25-50 keV and 100-300 energy ranges are consistent with zero in 90-95% of the cases, with most outliers probably representing the tail of the long burst class. We also analyze a small sample of "short" BATSE bursts - those with the most fluent, intense extended emission. The same lack of evolution on the pulse timescale obtains for the extended emission in the brighter bursts where significant measurements can be made. One possible inference is that both emission components may arise in the same region. We also show that the dynamic range in the ratio of peak intensities, spike : extended, is at least approx. l0(exp 3), and that for some bursts, the extended emission is only a factor of 2-5 lower. However, for our whole sample the total counts fluence of the extended component equals or exceeds that in the spike by a factor of several.

Norris, J. P.↗

NASA NTRS1994

Colors and color gradients in bulges of galaxies

We have obtained surface photometry in U, B, R, and I for a complete optically selected sample of 45 early-type spiral galaxies, to investigate the colors and color gradients of spiral bulges. Color profiles in U-R, B-R, U-B, and R-I have been determined in wedges opening on the semiminor axes. Based on several criteria, like the smoothness of the color profiles, the absence of dust lanes, and the central colors, we have defined a subsample of 18 objects whose colors are largely unaffected by dust. We believe such colors are suitable for inferring properties of the stellar populations of bulges. We find that the colors of bulges are predominantly bluer than those of ellipticals. This result holds even when bulges are compared to ellipticals of the same luminosity, and indicates that bulges are younger and/or more metal poor than old elliptical galaxies. Most bulges do not reach solar metallicities. Bulges show predominantly negative color gradients (bluer outward). For bright bulges (M(sup Bulge sub R) is less than -20.0), the magnitude of the gradient increases with bulge luminosity. For fainter bulges, gradients scatter around large negative values. The behavior of color gradients as a function of bulge luminosity suggests different formation mechanisms for faint and bright spheroids. For bright bulges, the scaling of gradients with luminosity suggests a formation process involving dissipation. The similarity with ellipticals suggests that the formation of the disk did not affect the stellar populations of the bulge in a major way. For small bulges (M(sub R) is greater than -20), the existence of pronounced color gradients suggests a different formation mechanism. For these objects, the presence of the disk may have severely affected the radial population distribution in the bulge.

Balcells, Marc↗

NASA NTRS1994

Landau damping of auroral hiss

Auroral hiss is observed to propagate over distances comparable to an Earth radius from its source in the auroral oval. The role of Landau damping is investigated for upward propagating auroral hiss. By using a ray tracing code and a simplified model of the distribution function, the effect of Landau damping is calculated for auroral hiss propagation through the environment around the auroral oval. Landau damping is found to be the likely mechanism for explaining some of the one-sided auroral hiss funnels observed by Dynamics Explorer 1. It is also found that Landau damping puts a lower limit on the wavelength of auroral hiss. Poleward of the auroral oval, Landau damping is found in a typical case to limit omega/k(sub parallel) to values of 3.4 x 10(exp 4) km/s or greater, corresponding to resonance energies of 3.2 keV or greater and wavelengths of 2 km or greater. For equatorward propagation, omega/k(sub parallel) is limited to values greater than 6.8 x 10(exp 4) km/s, corresponding to resonance energies greater than 13 keV and wavelengths greater than 3 km. Independent estimates based on measured ratios of the magnetic to electric field intensity also show that omega/k(sub parallel) corresponds to resonance energies greater than 1 keV and wavelengths greater than 1 km. These results lead to the difficulty that upgoing electron beams sufficiently energetic to directly generate auroral hiss of the inferred wavelength are not usually observed. A partial transmission mechanism utilizing density discontinuities oblique to the magnetic field is proposed for converting auroral hiss to wavelengths long enough to avoid damping of the wave over long distances. Numerous reflections of the wave in an upwardly flared density cavity could convert waves to significantly increased wavelengths and resonance velocities.

Morgan, D. D.↗

NASA NTRS2011

Spectral and Vertical Distribution Properties of Titan's Particulates from Thermal-IR CIRS Data: Physical Implications

Analyses of far-IR spectra between 20 and 560/cm (500 and 18 micron) recorded by the Cassini Composite Infrared Spectrometer (CIRS) yield the spectral dependence and the vertical distribution of Titan's photochemical aerosol and stratospheric ice clouds. Below the stratopause (approx. 300 km) the aerosol appears to be incompletely mixed for the following reasons: 1) the altitude dependence of the aerosol mass absorption coefficient is larger at higher altitudes than at lower altitudes, 2} the aerosol scale height varies with altitude, which implies some kind of layering effect, and 3) the aerosol abundance varies with latitude. The spectral shape of the aerosol opacity appears to be independent in altitude and latitude below the stratopause, even though inhomogeneities in the abundance appear to be prevalent throughout this altitude region. This implies that aerosol chemistry is restricted to altitude regions above the stratopause, where pressures are less than approx 0.1 mbar. The aerosol exhibits an extremely broad emisSion feature with a spectral peak at 140/cm (71 micron), which is not evident in laboratory simulated Titan aerosols (tholin) that are created at pressures greater than 0.1 mbar. A strong broad emission feature centered roughly around 160 cm-1 corresponds very closely to those found in nitrile ice spectra. This feature is pervasive throughout the region from high northern to high southern latitudes. The inference of nitrile ices is consistent with the highly restricted altitude ranges over which these features are observed, and appear to be dominated by HCN and HC3N. At low and moderate latitudes these clouds are observed to be located between 60 and 100 km, whereas at high northern latitudes during northern winter these clouds are observed at altitudes between 150 and 165 km. The ubiquitous nature of these nitrile ice clouds is inconsistent with a simple meridional circulation concept, suggesting that the true dynamical situation is more complex.

Anderson, Carrie M.↗

NASA NTRS1993

Angular momentum regulation in low-mass young stars surrounded by accretion disks

From study of a sample of 34 T Tauri stars with photometrically derived rotation periods and spectral types later than KS, we find that the observed periods appear to be related to the presence or absence of an accretion disk. Those stars which we infer to be surrounded by accretion disks have rotation periods P(rot) over 4 days with a most probable P(rot) of about 8.5 days, while those stars which lack accretion disk signatures cover a wide range of P(rot) from 1.5 to 16 days, including a significant number of objects with P(rot) less than 4 days. This suggests the possibility that the 'initial' angular momentum of a star is not established until it dissipates its circumstellar accretion disk. During the disk accretion phase, the stellar angular velocity appears to be regulated at a low value, countering the tendency of the star to spin up both from contraction toward the main sequence and from the accretion of inner disk material of high specific angular momentum. When the accretion disk is dissipated, this regulation mechanism will cease to function. At this point, the star is no longer maintained at a low angular velocity, but is 'free' to conserve its angular momentum, and thus to increase its angular velocity in response to contraction and changes in moment of inertia. This hypothesis, combined with a spread in disk dispersal time scales, provides a context for explaining the observed distribution of stellar rotational velocities for stars on the ZAMS in young clusters.

Edwards, Suzan↗

NASA NTRS2018

Evolution of H2O Production in Comet C/2012 S1 (ISON) as Inferred from Forbidden Oxygen and OH Emission

We present H2O production rates for comet C/2012 S1 (ISON) derived from observations of [O I] and OH emission during its inbound leg, covering a heliocentric distance range of 1.8–0.44 AU. Our production rates are in agreement with previous measurements using a variety of instruments and techniques and with data from the various observatories greatly differing in their projected fields of view. The consistent results across all data suggest the absence of an extended source of H2O production, for example sublimation of icy grains in the coma, or a source with spatial extent confined to the dimensions of the smallest projected field of view (in this case < 1000 km). We find that ISON had an active area of around 10 sq km for heliocentric distances Rh > 1.2 AU, which then decreased to about half this value from Rh = 1.2–0.9 AU. This was followed by a rapid increase in active area at about Rh = 0.6 AU, corresponding to the first of three major outbursts ISON experienced inside of 1 AU. The combination of a detected outburst in the light curve and rapid increase in active area likely indicates a major nucleus fragmentation event. The 5–10 sq km active area observed outside of Rh = 0.6 AU is consistent with a 50–100% active fraction for the nucleus, larger than typically observed for cometary nuclei. Although the absolute value of the active area is somewhat dependent on the thermal model employed, the changes in active area observed are consistent among models. The conclusion of a 50–100+% active fraction is robust for realistic thermal models of the nucleus. However the possibility of a contribution of a spatially unresolved distribution of icy grains cannot be discounted. As our [OI]-derived H2O production rates are consistent with values derived using other methods, we conclude that the contribution of O2 photodissociation to the observed [O I] emission is at most 5–10% that of the contribution of H2O for ISON. This is consistent with the expected contribution of O2 photodissociation if O2/H2O ∼ 4%, meaning [O I] emission can still be utilized as a reliable proxy for H2O production in comets as long as O2/H2O ≲ 4%, similar to the abundance measured by the ROSINA instrument on Rosetta at comet 67P/Churyumov–Gerasimenko.

Coma↗

NASA NTRS1986

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.↗

NASA NTRS2011

A Fractal Model for the Capacitance of Lunar Dust and Lunar Dust Aggregates

Lunar dust grains and dust aggregates exhibit clumping, with an uneven mass distribution, as well as features that span many spatial scales. It has been observed that these aggregates display an almost fractal repetition of geometry with scale. Furthermore, lunar dust grains typically have sharp protrusions and jagged features that result from the lack of aeolian weathering (as opposed to space weathering) on the Moon. A perfectly spherical geometry, frequently used as a model for lunar dust grains, has none of these characteristics (although a sphere may be a reasonable proxy for the very smallest grains and some glasses). We present a fractal model for a lunar dust grain or aggregate of grains that reproduces (1) the irregular clumpy nature of lunar dust, (2) the presence of sharp points, and (3) dust features that span multiple scale lengths. We calculate the capacitance of the fractal lunar dust analytically assuming fixed dust mass (i.e. volume) for an arbitrary number of fractal levels and compare the capacitance to that of a non-fractal object with the same volume, surface area, and characteristic width. The fractal capacitance is larger than that of the equivalent non-fractal object suggesting that for a given potential, electrostatic forces on lunar dust grains and aggregates are greater than one might infer from assuming dust grains are sphericaL Consequently, electrostatic transport of lunar dust grains, for example lofting, appears more plausible than might be inferred by calculations based on less realistic assumptions about dust shape and associated capacitance.

Collier, Michael R.↗

NASA NTRS2022

UV Spectropolarimetry with Polstar: Massive Star Binary Colliding Winds

The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mapping from that Doppler shift to the radial distance from the source is ambiguous. Binary systems can reduce this ambiguity by providing a second light source at a known radius in the wind, seen from orbitally modulated directions. From the nature of the collision between the winds, a massive companion also provides unique additional information about wind momentum fluxes. Since massive stars are strong ultraviolet (UV) sources, and UV resonance line opacity in the wind is strong, UV instruments with a high resolution spectroscopic capability are essential for extracting this dynamical information. Polarimetric capability also helps to further resolve ambiguities in aspects of the wind geometry that are not axisymmetric about the line of sight, because of its unique access to scattering direction information. We review how the proposed MIDEX-scale mission Polstar can use UV spectropolarimetric observations to critically constrain the physics of colliding winds, and hence radiatively-driven winds in general. We propose a sample of 20 binary targets, capitalizing on this unique combination of illumination by companion starlight, and collision with a companion wind, to probe wind attributes over a range in wind strengths. Of particular interest is the hypothesis that the radial distribution of the wind acceleration is altered significantly, when the radiative transfer within the winds becomes optically thick to resonance scattering in multiple overlapping UV lines.

Ultraviolet Astronomy↗

NASA NTRS1986

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.↗

NASA NTRSDate not supplied

Adding GPU Support to the Markov Chain Monte Carlo Code Catmip

In geophysics, we are confronted with many under-determined inverse problems. For example, all of our observations of earthquakes are made at the Earth’s surface. So, when we try to infer how slip during an earthquake evolves in space and time, we find that there are many potential slip histories that are consistent with our limited observations and our understanding of earthquake physics. One way to approach these problems is with Bayesian analysis which allows us to infer the ensemble of all potential slip models that satisfy the observations and our prior knowledge of earthquake physics. In Bayesian analysis, our prior knowledge is known as the prior probability density function or prior PDF, the fit to the data is known as the data likelihood, and the target PDF that satisfies both the prior PDF and data likelihood is known as the posterior PDF. However, simulating the posterior PDF typically requires using Markov Chain Monte Carlo (MCMC) to draw tens of billions of random realizations of earthquake slip models, which may not be computationally feasible. To make this and similar geophysical inversions computationally tractable, we developed the Cascading Adaptive Transitional Metropolis In Parallel (CATMIP) algorithm. CATMIP is an efficient parallel Markov Chain Monte Carlo (MCMC) sampler that is used for model fitting and uncertainty quantification in geophysics. Example use cases are earthquake rupture modeling, determining mineral composition on Mars, reconstructing the history of ocean salinity, and historical earthquake relocation. CATMIP employs many parallel instances of the Metropolis algorithm for sampling in a transitioning framework. Transitioning is a process in which a set of random samples at equilibrium with a known probability density function (PDF) are used as seeds for the Markov chains to sample successive target PDFs that incrementally move the distribution from the starting seeds to the final desired PDF that describes the relative plausibility of potential values for the model parameters. The algorithm is implemented as a Master-Worker model employing MPI for communication. The worker processes are loosely coupled with global parameters periodically optimized by the master process. This provides a very high amount of parallelism with little communication between updates. During the presentation we will discuss the history of the algorithm and elaborate the earthquake rupture modeling use case for the CATMIP package. Our first step toward GPU optimization was to optimize the code for the CPU. CPU profiling revealed that most of the compute time is spent in calls to level 2 BLAS routines and calls to GSL random number generators. We revised the algorithm to employ level 3 BLAS routines instead. In our presentation we will describe how this was accomplished. Adding GPU support to CATMIP consisted mostly of replacing the calls to GSL with calls to GPU vendor-provided library routines. A small number of loops were directly implemented in CUDA. In the presentation will provide implementation details. Finally, we will discuss methods for profiling and opportunities for further optimizing GPU execution. By creating a code with the flexibility to run on either a CPU or GPU architecture, CATMIP can be used on systems ranging from large CPU-based HPC environments to single servers with GPU acceleration and everything in between.

HECC↗

NASA NTRS1990

Mechanism of transport and distribution of organic solvents in blood

Little is known about the mechanism of transport and distribution of volatile organic compounds in blood. Studies were conducted on five typical organic solvents to investigate how these compounds are transported and distributed in blood. Groups of four to five rats were exposed for 2 hr to 500 ppm of n-hexane, toluene, chloroform, methyl isobutyl ketone (MIBK), or diethyl ether vapor; 94, 66, 90, 51, or 49%, respectively, of these solvents in the blood were found in the red blood cells (RBCs). Very similar results were obtained in vitro when aqueous solutions of these solvents were added to rat blood. In vitro studies were also conducted on human blood with these solvents; 66, 43, 65, 49, or 46%, respectively, of the added solvent was taken up by the RBCs. These results indicate that RBCs from humans and rats exhibited substantial differences in affinity for the three more hydrophobic solvents studied. When solutions of these solvents were added to human plasma and RBC samples, large fractions (51-96%) of the solvents were recovered from ammonium sulfate-precipitated plasma proteins and hemoglobin. Smaller fractions were recovered from plasma water and red cell water. Less than 10% of each of the added solvents in RBC samples was found in the red cell membrane ghosts. These results indicate that RBCs play an important role in the uptake and transport of these solvents. Proteins, chiefly hemoglobin, are the major carriers of these compounds in blood. It can be inferred from the results of the present study that volatile lipophilic organic solvents are probably taken up by the hydrophobic sites of blood proteins.

Solvents/administration & dosage/metabolism/pharma↗