Engineering PapersSearch

Engineering topics

Cuzzi, Jeffery N.

Publications and source records attributed to Cuzzi, Jeffery N..

Transport of Volcanic Clouds Injected in the Tropics: Case Studies of El Chichon and Mt. Pinatubo

The eruptions of Mt. Pinatubo and El Chichon at nearly the same tropical latitude produced two of the largest volcanic clouds this century. Climatic effects of any volcanic eruption depend on spatial distribution of the volcanic cloud. We report the first global simulations which reproduce the markedly different inter-hemispheric transport observed for these clouds in the months immediately following each eruption. Differences in stratospheric mean easterly winds and unequal radiative heating of the clouds explain the different latitudinal transport. It is found that radiative heating of tropical clouds is insensitive to particle size for typical size ranges. Radiative energy balance of optically thick tropical volcanic clouds tends to position the clouds near 25 km altitude.

Young, R. E.

Nonlinear Propagating Features in the Stratosphere of Jupiter Generated by the Impact of SL-9

Among the more intriguing features observed by HST after the impacts of the fragments of Shoemaker-Levy 9 were rings, encircling several of the impact sites, that expanded radially at constant velocities of order 450-500 m/s. One credible suggestion has been that the rings are stratospheric manifestations of linear gravity waves seated in the troposphere (Ingersoll and Kanamori, 1995). In this scenario the waves are made visible by a condensible material of unknown composition. Based on simulations of the comet impact by Zahnle and Mac Low (1995), we suggest an alternative hypothesis, namely that the rings are nonlinear stratospheric gravity waves which break and entrain impact debris material in the breaking wave front. The computations produce just such a feature, with the nonlinear breaking wave front in the stratosphere propagating radially away from the impact site with the correct speed. The computed ring of material is located in an altitude region where the static stability profile has a local maximum and therefore forms a waveguide. Because entrained impact debris forms the ring feature in the model, a condensate of some unknown composition is not required in order to view the wave front. However, there is one severe problem with the computed nonlinear wave: in contrast to a linear wave, it slows down. The computed feature stops within 1000 seconds, while the observed feature lasted at least 10,000 seconds. Since the nonlinear breaking wave front is a natural product of our impact simulations, we are now investigating the conditions that will enable the breaking wave front to propagate further, in accord with the observations. The results of these computations will be presented.

Young, R. E.

A Mars Rover Mission Simulation on Kilauea Volcano

A field experiment to simulate a rover mission on Mars was performed using the Russian Marsokhod rover deployed on Kilauea Volcano HI in February, 1995. A Russian Marsokhod rover chassis was equipped with American avionics equipment, stereo cameras on a pan and tilt platform, a digital high resolution body-mounted camera, and a manipulator arm on which was mounted a camera with a close-up lens. The six wheeled rover is 2 meters long and has a mass of 120 kg. The imaging system was designed to simulate that used on the planned "Mars Together" mission. The rover was deployed on Kilauea Volcano HI and operated from NASA Ames by a team of planetary geologists and exobiologists. Two modes of mission operations were simulated for three days each: (1) long time delay, low data bandwidth (simulating a Mars mission), and (2) live video, wide-bandwidth data (allowing active control simulating a Lunar rover mission or a Mars rover mission controlled from on or near the Martian surface). Simulated descent images (aerial photographs) were used to plan traverses to address a detailed set of science questions. The actual route taken was determined by the science team and the traverse path was frequently changed in response to the data acquired and to unforeseen operational issues. Traverses were thereby optimized to efficiently answer scientific questions. During the Mars simulation, the rover traversed a distance of 800 m. Based on the time delay between Earth and Mars, we estimate that the same operation would have taken 30 days to perform on Mars. This paper will describe the mission simulation and make recommendations about incorporating rovers into the Mars surveyor program.

Stoker, Carol

Heliospheric Termination Shock Motion Due to Fluctuations in the Solar Wind Upstream Conditions: Spherically Symmetric Model

Large-scale fluctuations in the solar wind plasma upstream of the heliospheric termination shock (TS) will cause inward and outward motions of the shock. Using numerical techniques, we extend an earlier strictly one-dimensional (planar) analytic gas dynamic model to spherical symmetry to investigate the features of global behavior of shock motion. Our starting point is to establish a steady numerical solution of the gasdynamic equations describing the interaction between the solar wind and the interstellar medium. We then introduce disturbances of the solar wind dynamic pressure at an inner boundary, and follow the subsequent evolution of the system, especially the motion of the termination shock. Our model solves spherically symmetric gasdynamic equations as an initial-boundary value problem. The equations in conservative form are solved using a fully implicit Total Variation Diminishing (TVD) upwind scheme with Roe-type Riemann solver. Boundary conditions are given by the solar wind parameters on an inner spherical boundary, where they are allowed to vary with time for unsteady calculations, and by a constant pressure (roughly simulating the effect of the local interstellar medium) on an outer boundary. We find that immediately after the interaction, the shock moves with speeds given by the earlier analogous analytic models. However, as the termination shock propagates it begins to slow down, seeking a new equilibrium position. In addition, the disturbance transmitted through the TS, either a shock or rarefaction wave, will encounter the heliopause boundary and be reflected back. The reflected signal will encounter the TS, causing it to oscillate. The phenomenon may be repeated for a number of reflections, resulting in a "ringing" of the outer heliosphere.

Ratkiewicz, R.

Why the SL9 Plumes Were All About the Same Height

Several of the SL9 ejecta plumes were observed by the HST to reach approximately the same height, about 3000 km above the jovian cloud tops. The duration of the infrared events, which were produced by the plume falling back on the atmosphere, measures time aloft and hence provides a second, more sensitive measure of plume height; the light curves indicate that the largest impacts produced modestly higher plumes. Evidently these plumes were launched with about the same vertical velocity, roughly 10-13 kilometers per second. As the impactors themselves were not all the same, nor the impacts equally luminous, nor the plumes equally opaque, the similar plume heights has been seen as a puzzle needing explanation. A second, closely related matter that needs to addressed quantitatively is the popular contention that a big plume requires a big impact. This view is misleading at best, yet plume heights can be used to constrain impact parameters. Dimensional analysis indicates that plume height goes as z alpha v (sup 2) (sub ej) alpha E/pH (sup 2), where v (sub ej) is the ejection velocity, E the explosion energy, and p and H the ambient pressure and scale height at termination. Using a semi-analytic model for the deceleration, disintegration, and destruction of intruding bodies by an ever-vigilant atmosphere, we find that the ratio E/pH(sup 2) is roughly constant for fragments with diameters of order 100 m to 1000 m. Constancy of v(sub ej) is in part due to the greater role of radiative ablation on the flight of smaller objects. We conclude that similar plume heights is a direct consequence of smaller impactors exploding at higher altitudes, in such a way that the different explosions were geometrically similar.

Zahnle, K.

Profiles of Gamma-Ray Bursts and Their Component Pulses

One physically informative regularity of their otherwise heterogeneous ensemble, is that many Gamma-Ray Bursts consist of well defined pulses. To objectively quantify the temporal structure of BATSE bursts, we have developed an automatic modeling procedure that separates overlapping pulses and determines the energy-dependence of the pulse-shape parameters. No binning of photon arrival times is needed, so when applied to time-tagged events (TTE) the procedure captures variability information down to the shortest time scales present in the raw data. Maximizing the Bayesian likelihood function Pr(data/model) yields estimates of the model parameters, including the number of pulses present, and allows intercomparison of models of different forms. As with any nonlinear optimization, good initial guesses are crucial to avoid convergence to undesirable local minima. We find excellent initial pulse decompositions by wavelet-denoising a cumulative distribution of the raw photon arrival data; differentiation then gives a time profile mostly free of the systematic effects of degraded resolution (as in ordinary Fourier smoothing) and binning. We present statistical information on pulse rise-time, decay-time, peakedness, and amplitudes, plus their energy dependences - both within a single burst and for a large ensemble of bursts.

Scargle, Jeff D.

Tropospheric Lapse Rate and Methane on Titan

We have reanalyzed the Voyager radio occultation data for Titan with two alternative approaches to methane condensation. In one approach, methane condensation is enhanced by the presence of nitrogen. In the other approach, methane condensation does not occur. As pointed out by Thompson, nitrogen lowers the condensation level for a methane/nitrogen Mixture and we find that the upper limit on surface relative humidity of methane obtained from the Voyager occultation data is lowered from 0.7 to 0.6. However, 140% supersaturation of methane in the troposphere, suggested by Courtin et al., allows all surface humidities to be consistent with the Voyager occultation data and the upper limit is set by other considerations. We conclude that if supersaturation is not included then the surface relative humidity of methane is between 0.08 and 0.6, with values close to 0.6 indicated. If supersaturation is included then the surface relative humidity of methane is between 0.08 and 0.85, again, with values close to 0.6 indicated. The tropospheric lapse rate on Titan appears to be determined by radiative equilibrium. It is everywhere stable against dry convection but is unstable to moist convection. This is consistent with a supersaturated atmosphere in which condensation - and hence moist convection - is inhibited. The absence of dry convection in the troposphere of Titan can be explained by a simple grey model which shows that the radiative profile of any gas for which the ratio of the gas constant to the specific heat at constant pressure is greater than 0.25 never becomes unstable to dry convection.

McKay, C. P.

Interstellar Dust Connections to Interplanetary Dust

The organic component of the interstellar medium, which has revealed itself through the ubiquitous 3.4 micrometer hydrocarbon absorption feature, is widespread throughout the disk of our galaxy and has been attributed to dust grains residing in the diffuse interstellar medium. The absorption band positions near 3.4 micrometers are characteristic of C-H stretching vibrations in the -CH3 and -CH2- groups of saturated aliphatic hydrocarbons associated with perturbing chemical groups. The production of complex molecules is thought to occur within dense molecular clouds when ice- mantled grains are processed by various energetic mechanisms. Studies of the processing of interstellar ices and the subsequent production of organic residues have relevance to studies of interplanetary dust, because primitive solar system bodies are likely reservoirs of organic material, either preserved from the interstellar medium or produced in situ. A comparison of the interstellar organics and the Murchison meteorite illustrates the importance of probing the interstellar connection to the solar system, because although the carbonaceous meteorites are undoubtedly highly processed, they do retain specific interstellar signatures (such as diamonds, SiC grains, graphite and enriched D/H). The organic component, while not proven interstellar, has a remarkable similarity to the interstellar organics observed in over a dozen sightlines through our galaxy. This paper compares spectra from laboratory organics produced through the processing of interstellar ice analog materials with the high resolution infrared observations of the interstellar medium in order to investigate the mechanisms (such as ion bombardment, plasma processing, and UV photolysis) that may be producing the organics in the ISM. Comparisons to interplanetary dust particles, comets, and meteorites are discussed.

Pendleton, Y. J.

Calibration of Tests for Time-Dilation in GRB Pulse Structures

If gamma-ray bursts are at cosmological distances, then not only their constituent pulses but also the intervals between pulses should be time-dilated. Unlike time-dilation measures of pulse emission, intervals would appear to require less "K-correction" for redshift of narrower temporal structure from higher energy into the band of observation. However, stretching of pulse intervals is intrinsically difficult to measure without incurring a timescale-dependent bias since, as time profiles are stretched, more structure can appear at the limit of resolution. This problem is compounded in dimmer bursts since identification of significant structures is problematic. We have attempted to minimize brightness bias problems by equalizing the signal-to-noise level of all bursts. Then, analyzing wavelet-denoised profiles binned to several resolutions, we identify significant (3 - 4 sigma) fluctuations between pulse structures and interjacent valleys. When bursts are ranked by peak flux, an interval time-dilation signature is evident, but its measure is dependent upon choice of temporal resolution. This result appears to tentatively and qualitatively confirm the work of Davis, who found such a signature using a pulse fitting approach.

Norris, J. P.

The 5-8 Micron Infrared Spectrum of the Galactic Center and the Composition of Interstellar Dust

Interstellar dust is an important component of the interstellar medium which dominates the opacity and, hence, regulates radiative transfer, molecule formation, and thermal balance of the ISM. Much of our knowledge on the composition of interstellar dust results from infrared spect,oscopy. The extinction along the line of sight towards the galactic center is believed to be dominated by dust in the diffuse ISM. Because of the high extinction and high IR flux, IR spectra of galactic center sources have been a prime sampling ground for the characteristics of interstellar dust. We have obtained 5-8 micrometer spectra of the galactic center using the KAO. These spectra show absorption features at 5.5, 5.8, 6.1, and 6.8 micrometers. Together with features in the 3 micrometer region previously observed by us using the IRTF, these features are compared to laboratory spectra of candidate materials. We conclude that the 3.0 and 6.1 micrometer feature are carried by H2O, likely in the form of water of hydration in interstellar silicates. The 3.4, 5.5, 5.8, and 6.8 micrometer features are due to CH2, CH3 and C=O stretching and deformation modes in a hydrocarbon grain component. Comparing derived dust abundances, we conclude that silicates dominate the interstellar dust volume. Hydrocarbon and (small) graphite grains contribute each about 0.1. The remainder of the interstellar dust volume does not show strong IR absorption features and is likely in the form of large graphite, amorphous, carbon, or diamond grains.

Tielens, A. G. G. M.

Formation Of the Giant Planets By Concurrent Accretion Of Solids And Gas

New numerical simulations of the formation of the giant planets are presented, in which for the first time both the gas and planetesimal accretion rates are calculated in a self-consistent, interactive fashion. The simulations combine three elements: 1) three-body accretion cross-sections of solids onto an isolated planetary embryo, 2) a stellar evolution code for the planet's gaseous envelope, and 3) a planetesimal dissolution code within the envelope, used to evaluate the planet's effective capture radius and the energy deposition profile of accreted material. Major assumptions include: The planet is embedded in a disk of gas and small planetesimals with locally uniform initial surface mass density, and planetesimals are not allowed to migrate into or out of the planet's feeding zone. All simulations are characterized by three major phases. During the first phase, the planet's mass consists primarily of solid material. The planetesimal accretion rate, which dominates that of gas, rapidly increases owing to runaway accretion, then decreases as the planet's feeding zone is depleted. During the second phase, both solid and gas accretion rates are small and nearly independent of time. The third phase, marked by runaway gas accretion, starts when the solid and gas masses are about equal. It is engendered by a strong positive feedback on the gas accretion rates, driven by the rapid contraction of the gaseous envelope and the rapid expansion of the outer boundary, which depends on the planet's total mass. The overall evolutionary time scale is generally determined by the length of the second phase. The actual rates at which the giant planets accreted small planetesimals is probably intermediate between the constant rates assumed in most previous studies and the highly variable rates that we have used. Within the context, of the adopted model of planetesimal accretion, the joint constraints of the time scale for dissipation of the solar nebula and the current high-Z masses of the giant planets lead to estimates of the initial surface density (sigma(sub init)) of planetesimals in the outer region of the solar nebula. The results show sigma(sub init) approx. = 10 g/sq cm near Jupiter's orbit and that sigma(sub init) proportional to alpha(sup -2), where alpha is the distance from the Sun. These values are a factor of 3 - 4 times as high as that of the "minimum mass" solar nebula at Jupiter's distance and a factor of 2 - 3 times as high it Saturn's distance. Our estimates for the formation time of Jupiter and Saturn are 1 - 10 million years while those for Uranus fall in the range of 2 - 16 million years. These estimates follow from the properties of our Solar System and do not necessarily apply to giant planets in other planetary systems.

Pollack, James B.

Obliquity Experiments with a Mars General Circulation Model

We have simulated the seasonal variation of the general circulation on Mars for obliquities of 0deg and 60deg. These obliquities represent the minimum and maximum values the planet has experienced during the past 10(exp 7) years (e.g., Laskar and Robutel, 1993, Nature, 361, 608-614). The model we use is the NASA/Ames Mars General Circulation Model (Pollack et al., 1993, J. Geophys. Res. 98, 3149-3181). We vary only the obliquity; all other model parameters are as in Pollack et al. At high obliquity, the model shows dramatic seasonal variations in the polar caps and in the structure and intensity of the circulation. At the solstices the winter cap extends to the equator. Thus, surface temperatures throughout the entire winter hemisphere are fixed at the CO2 frost point. During summer surface temperatures at the poles reach 269K in the north and 295K in the south. The most notable changes to the circulation at solstice compared to our standard runs are a general weakening of the winter westerlies, a Hadley cell of greater latitudinal extent, and the development of very strong, possibly unstable, low-level jets in midlatitudes of the summer hemisphere. Surface stresses associated with these jets are sufficient to raise dust continuously. Thus, dust storms should be frequent features of the high obliquity climate. This result is independent of any desorbed regolith CO2 which would raise mean surface pressures. At zero obliquity the structure of the circulation resembles that of present day equinox conditions modulated by the varying insolation associated with orbital eccentricity. Notable features include equatorial superrotation, asymmetric Hadley cells, and stronger poleward heat fluxes in the northern hemisphere. Since the poles do not receive solar energy at any time of year, permanent caps form which extend to about 70deg in each hemisphere. However, the north permanent cap is growing at a rate 40% faster than the south cap. This is due to the differences in topography, albedo, emissivity, and heat transport in each hemisphere. Thus the zero obliquity climate system appears to be evolving toward a state in which the north cap is the dominant cap.

Harberle, R. M.