ALBERT

Description: This paper describes a computer program developed for structural dynamic analysis of horizontal axis wind turbines (HAWT's). It is based on the finite element method through its reliance on NASTRAN for the development of mass, stiffness, and damping matrices of the tower end rotor, which are treated in NASTRAN as separate structures. The tower is modeled in a stationary frame and the rotor in one rotating at a constant angular velocity. The two structures are subsequently joined together (external to NASTRAN) using a time-dependent transformation consistent with the hub configuration. Aerodynamic loads are computed with an established flow model based on strip theory. Aeroelastic effects are included by incorporating the local velocity and twisting deformation of the blade in the load computation. The turbulent nature of the wind, both in space and time, is modeled by adding in stochastic wind increments. The resulting equations of motion are solved in the time domain using the implicit Newmark-Beta integrator. Preliminary comparisons with data from the Boeing/NASA MOD2 HAWT indicate that the code is capable of accurately and efficiently predicting the response of HAWT's driven by turbulent winds.

Description: Underintegrated methods are investigated with respect to their stability and convergence properties. The focus was on identifying regions where they work and regions where techniques such as hourglass viscosity and hourglass control can be used. Results obtained show that underintegrated methods typically lead to finite element stiffness with spurious modes in the solution. However, problems exist (scalar elliptic boundary value problems) where underintegrated with hourglass control yield convergent solutions. Also, stress averaging in underintegrated stiffness calculations does not necessarily lead to stable or convergent stress states.

Description: The Magellan radar-mapping mission collected a large amount of science and engineering data. Now available to the general scientific community, this data set can be overwhelming to someone who is unfamiliar with the mission. This user guide outlines the mission operations and data set so that someone working with the data can understand the mapping and data-processing techniques used in the mission. Radar-mapping parameters as well as data acquisition issues are discussed. In addition, this user guide provides information on how the data set is organized and where specific elements of the set can be located.

Description: A statistically generated weighting function for a second-order polynomial curve fit of residual functions has been developed. The residual flexibility test method, from which a residual function is generated, is a procedure for modal testing large structures in an external constraint-free environment to measure the effects of higher order modes and interface stiffness. This test method is applicable to structures with distinct degree-of-freedom interfaces to other system components. A theoretical residual function in the displacement/force domain has the characteristics of a relatively flat line in the lower frequencies and a slight upward curvature in the higher frequency range. In the test residual function, the above-mentioned characteristics can be seen in the data, but due to the present limitations in the modal parameter evaluation (natural frequencies and mode shapes) of test data, the residual function has regions of ragged data. A second order polynomial curve fit is required to obtain the residual flexibility term. A weighting function of the data is generated by examining the variances between neighboring data points. From a weighted second-order polynomial curve fit, an accurate residual flexibility value can be obtained. The residual flexibility value and free-free modes from testing are used to improve a mathematical model of the structure. The residual flexibility modal test method is applied to a straight beam with a trunnion appendage and a space shuttle payload pallet simulator.

Description: In elasticity, the method of forces, wherein stress parameters are considered as the primary unknowns, is known as the Beltrami-Michell formulation (BMF). The existing BMF can only solve stress boundary value problems; it cannot handle the more prevalent displacement of mixed boundary value problems of elasticity. Therefore, this formulation, which has restricted application, could not become a true alternative to the Navier's displacement method, which can solve all three types of boundary value problems. The restrictions in the BMF have been alleviated by augmenting the classical formulation with a novel set of conditions identified as the boundary compatibility conditions. This new method, which completes the classical force formulation, has been termed the completed Beltrami-Michell formulation (CBMF). The CBMF can solve general elasticity problems with stress, displacement, and mixed boundary conditions in terms of stresses as the primary unknowns. The CBMF is derived from the stationary condition of the variational functional of the integrated force method. In the CBMF, stresses for kinematically stable structures can be obtained without any reference to the displacements either in the field or on the boundary. This paper presents the CBMF and its derivation from the variational functional of the integrated force method. Several examples are presented to demonstrate the applicability of the completed formulation for analyzing mixed boundary value problems under thermomechanical loads. Selected example problems include a cylindrical shell wherein membrane and bending responses are coupled, and a composite circular plate.

Description: The Martian regolith is the most substantial volatile reservoir on the planet; estimates of its adsorbed inventory have been based on simple measurements of the adsorption of either water or CO2 in isolation. Under some conditions, H2O can poison adsorbate surfaces, such that CO2 uptake is greatly reduced. We have made the first measurements of the simultaneous adsorption of CO2 and H2O under conditions appropriate to the Martian regolith and have found that at H2O monolayer coverage above about 0.5, CO2 begins to be displaced into the gas phase. We have developed an empirical expression that describes our co-adsorption data and have applied it to standard models of the Martian regolith. We find that currently, H2O does not substantially displace CO2, implying that the adsorbate inventories previously derived may be accurate, not more than 3-4 kPa (30-40 mbar). No substantial increase in atmospheric pressure is predicted at higher obliquities because high-latitude ground ice buffers the partial pressure of H2O in the pores, preventing high monolayer coverages of H2O from displacing CO2. The peak atmospheric pressure at high obliquity does increase as the total inventory of exchangeable CO2 increases.

Description: We present plasma, magnetic field, and electric field data of the Pioneer Venus Orbiter (PVO) showing that the shocked solar wind in the Venus inner ionosheath exhibits flow conditions substantially different from those in the outer ionosheath. In particular, the plasma density is seen to drop significantly to low values within a layer adjacent to, and downstream from, the planet's ionopause. This change is not seen to develop gradually as the PVO moves into that region of space but occurs abruptly across a well-defined transition which extends downstream along the flanks of the Venus ionosheath. We explore the implications that these observations have in regard to the character of the interaction process between the shocked solar wind and the ionospheric plasma. It is argued that the existence of a sharply bounded region in the inner ionosheath within which the plasma density is severely depressed is consistent with the existence of friction at and near the ionopause. Plasma perturbations generated at this latter boundary, and distributed downstream through the ionosheath flow, may be responsible for the change of properties exhibited by the solar wind plasma in the inner ionosheath.

Description: Statistical studies of the VLF bursts detected in the nightside ionosphere of Venus show that the bursts fall into two classes. The first consists of signals detected when vertical propagation within the whistler-mode resonance cone is allowed. The second consists of signals whose burst rate decreases rapidly as a function of increasing altitude, with a scale height of about 20 km. These non-whistler-mode signals also display a strong dependence on local time, with the burst rates being largest in the post-dusk local time sector. Since these signals are not propagating we assume that they correspond to a 'near-field' or prompt response to a lightning stroke. As such we can use these signals to estimate the planetary lightning rate, and we find that the rate at Venus is comparable to or greater than the terrestrial planetary rate of 100 flashes/sec.

Description: In this study, we sought to characterize variations in the abundance and distribution of subcloud H2SO4(g) in the Venus atmosphere by using a number of 13cm radio occultation measurements conducted with the Pioneer Venus Orbiter near the inferior conjunction of 1991. A total of ten data sets were examined and analyzed, producing vertical profiles of temperature and pressure in the neutral atmosphere, and sulfuric acid vapor abundance below the main cloud layer. Two of the vertical profiles of the abundance of H2SO4(g) were correlated with NIR images of the night side of Venus made during the same period of time by Boris Ragent (under a separate PVO Guest Investigator Grant). Initially, we had hoped that the combination of these two different types of data would make it possible to constrain or identify the composition of the large particles causing the features observed in the NIR images. However, the sparseness of the radio occultation data set, combined with the sparseness of the NIR data set (one image per day over an 8 day period) made it impossible to draw strong conclusions. Considered on their own, however, the parameters retrieved from the radio occultation experiments are valuable science products.

Description: A workshop was convened at Ames Research Center on September 28 and 29, 1993, to address the unexplained electrical anomalies experienced in December 1978 by the four Pioneer Venus probes below a Venus altitude of 12.5 km. These anomalies caused the loss of valuable data in the deep atmosphere, and, if their cause were to remain unexplained, could reoccur on future Venus missions. The workshop participants reviewed the evidence and studied all identified mechanisms that could consistently account for all observed anomalies. Both hardware problems and atmospheric interactions were considered. Based on a workshop recommendation, subsequent testing identified the cause as being an insulation failure of the external harness. All anomalous events are now explained.

Description: Structural analyses are developed to determine the linear elastic and the geometrically nonlinear elastic response of an internally pressurized, orthogonally stiffened, composite material cylindrical shell. The configuration is a long circular cylindrical shell stiffened on the inside by a regular arrangement of identical stringers and identical rings. Periodicity permits the analysis of a unit cell model consisting of a portion of the shell wall centered over one stringer-ring joint. The stringer-ring-shell joint is modeled in an idealized manner; the stiffeners are mathematically permitted to pass through one another without contact, but do interact indirectly through their mutual contact with the shell at the joint. Discrete beams models of the stiffeners include a stringer with a symmetrical cross section and a ring with either a symmetrical or an asymmetrical open section. Mathematical formulations presented for the linear response include the effect of transverse shear deformations and the effect of warping of the ring's cross section due to torsion. These effects are important when the ring has an asymmetrical cross section because the loss of symmetry in the problem results in torsion and out-of-plane bending of the ring, and a concomitant rotation of the joint at the stiffener intersection about the circumferential axis. Data from a composite material crown panel typical of a large transport fuselage structure are used for two numerical examples. Although the inclusion of geometric nonlinearity reduces the 'pillowing' of the shell, it is found that bending is localized to a narrow region near the stiffener. Including warping deformation of the ring into the analysis changes the sense of the joint rotation. Transverse shear deformation models result in increased joint flexibility.

Description: The general objective for this work was to develop a theoretically and experimentally consistent explanation for the diffuse component of radar backscatter from Mars. The strength, variability, and wavelength independence of Mars' diffuse backscatter are unique among our Moon and the terrestrial planets. This diffuse backscatter is generally attributed to wavelength-scale surface roughness and to rock clasts within the Martian regolith. Through the combination of theory and experiment, the authors attempted to bound the range of surface characteristics that could produce the observed diffuse backscatter. Through these bounds they gained a limited capability for data inversion. Within this umbrella, specific objectives were: (1) To better define the statistical roughness parameters of Mars' surface so that they are consistent with observed radar backscatter data, and with the physical and chemical characteristics of Mars' surface as inferred from Mariner 9, the Viking probes, and Earth-based spectroscopy; (2) To better understand the partitioning between surface and volume scattering in the Mars regolith; (3) To develop computational models of Mars' radio emission that incorporate frequency dependent, surface and volume scattering.

Description: A methodology to compute cumulative probability distribution functions (CDF) of fatigue life for different ratios, r of applied stress to the laminate strength based on first ply failure criteria has been developed and demonstrated. Degradation effects due to long term environmental exposure and mechanical cyclic loads are considered in the simulation process. A unified time-stress dependent multi-factor interaction equation model developed at NASA Lewis Research Center has been used to account for the degradation/aging of material properties due to cyclic loads. Fast probability integration method is used to perform probabilistic simulation of uncertainties. Sensitivity of fatigue life reliability to uncertainties in the primitive random variables are computed and their significance in the reliability based design for maximum life is discussed. The results show that the graphite/epoxy (0/+45/90) deg laminate with ply thickness 0.125 in. has 500,000 cycles life for applied stress to laminate strength ratio of 0.6 and a reliability of 0.999. Also, the fatigue life reliability has been found to be most sensitive to the ply thickness and matrix tensile strength. Tighter quality controls must therefore be enforced on ply thickness and matrix strength in order to achieve high reliability of the structure.

Description: Specific forms for both the Gibb's and complementary dissipation potentials are chosen such that a complete (i.e., fully associative) potential base multiaxial, nonisothermal unified viscoplastic model is obtained. This model possesses one tensorial internal state variable (that is, associated with dislocation substructure) and an evolutionary law that has nonlinear kinematic hardening and both thermal and strain induced recovery mechanisms. A unique aspect of the present model is the inclusion of nonlinear hardening through the use of a compliance operator, derived from the Gibb's potential, in the evolution law for the back stress. This nonlinear tensorial operator is significant in that it allows both the flow and evolutionary laws to be fully associative (and therefore easily integrated), greatly influences the multiaxial response under non-proportional loading paths, and in the case of nonisothermal histories, introduces an instantaneous thermal softening mechanism proportional to the rate of change in temperature. In addition to this nonlinear compliance operator, a new consistent, potential preserving, internal strain unloading criterion has been introduced to prevent abnormalities in the predicted stress-strain curves, which are present with nonlinear hardening formulations, during unloading and reversed loading of the external variables. The specific model proposed is characterized for a representative titanium alloy commonly used as the matrix material in SiC fiber reinforced composites, i.e., TIMETAL 21S. Verification of the proposed model is shown using 'specialized' non-standard isothermal and thermomechanical deformation tests.

Description: Solving for the displacements of free-free coupled systems acted upon by static loads is commonly performed throughout the aerospace industry. Many times, these problems are solved using static analysis with inertia relief. This solution technique allows for a free-free static analysis by balancing the applied loads with inertia loads generated by the applied loads. For some engineering applications, the displacements of the free-free coupled system induce additional static loads. Hence, the applied loads are equal to the original loads plus displacement-dependent loads. Solving for the final displacements of such systems is commonly performed using iterative solution techniques. Unfortunately, these techniques can be time-consuming and labor-intensive. Since the coupled system equations for free-free systems with displacement-dependent loads can be written in closed-form, it is advantageous to solve for the displacements in this manner. Implementing closed-form equations in static analysis with inertia relief is analogous to implementing transfer functions in dynamic analysis. Using a MSC/NASTRAN DMAP Alter, displacement-dependent loads have been included in static analysis with inertia relief. Such an Alter has been used successfully to solve efficiently a common aerospace problem typically solved using an iterative technique.

Description: The objectives of the research are to measure low temperature laboratory rate coefficients for key reactions relevant to the atmospheres of Titan and Saturn. These reactions are, for example, C2H + H2, CH4, C2H2, and other hydrocarbons which need to be measured at low temperatures, down to approximately 150 K. The results of this work are provided to NASA specialists who study modeling of the hydrocarbon chemistry of the outer planets. The apparatus for this work consists of a pulsed laser photolysis system and a tunable F-center probe laser to monitor the disappearance of C2H. A low temperature cell with a cryogenic circulating fluid in the outer jacket provides the gas handling system for this work. These elements have been described in detail in previous reports. Several new results are completed and the publications are just being prepared. The reaction of C2H with C2H2 has been measured with an improved apparatus down to 154 K. An Arrhenius plot indicates a clear increase in the rate coefficient at the lowest temperatures, most likely because of the long-lived (C4H3) intermediate. The capability to achieve the lowest temperatures in this work was made possible by construction of a new cell and addition of a multipass arrangement for the probe laser, as well as improvements to the laser system.

Description: We report on further studies of radio wave bursts detected by the Orbiting Electric Field Detector (OEFD) on the Pioneer Venus Orbiter (PVO) in the nightside ionosphere of Venus. We have tested a total of 25 cases of wave burst activity for evidence of whistler-mode propagation to the spacecraft from impulsive subionospheric sources. As in a previous study of 11 of these cases (Sonwalkar et al., 1991) we find at least two distinct classes of events, one, mostly involving bursts at 100 Hz only, that passes certain tests for whistler-mode propagation, and another, mostly involving bursts in two or more of the four PVO narrowband channels (at 100 Hz, 730 Hz, 5.4 kHz, and 30 kHz), that fails to pass the tests. The subionospheric lightning hypothesis continues to be tenable as a candidate explanation for many of the 100 Hz-only events, but its number of 100 Hz-only cases that do no pass all the applicable whistler-mode tests, as well as the existence at a wide range of altitudes of multichannel cases that are clearly not propagating whistler-mode waves. The wideband bursts are often observed at altitudes above 1000 km and frequently occur in regions of locally reduced electron density. Those observed at high altitude (and possibly low altitude as well) are believed to be generated near the spacecraft, possibly by an as yet unknown mechanism responsible for similar burst observations made near Earth and other planets.

Description: Plasma wave data from the Pioneer Venus Orbiter provide the largest body of data cited as evidence for lightning on Venus. These data are also the most controversial, mainly because of the ambiguity in mode identification due to limited spectral information. We review some of the more recent studies of the plasma wave data at Venus, and we demonstrate that the characteristics of the 100 Hz waves are consistent with whistler-mode waves propagating vertically from below the ionosphere. We further show that in situ instabilities are too weak to generate whistler-mode waves, mainly because the thermal pressure is comparable with the magnetic field pressure in the ionosphere of Venus. The lower hybrid drift instability has also been suggested as an alternative source for the 100 Hz waves. However, the wave properties are more consistent with whistler-mode propagation; the lower hybrid dirft instability requires very short gradient scale lengths to overcome damping due to collisions. We also note that an apparent association between Langmuir probe anomalies and 100 Hz waves is much lower than previously reported, once we apply a consistent intensity threshold for identifying wave bursts. The lightning hyposthesis remains the most probable explanation of the plasma waves detected at low altitudes in the nightside ionosphere of Venus.

Description: Very low frequency (VLF) signals detected in the nightside ionosphere of Venus have generally been attributed to atmospheric lightning. However it has recently been suggested that these bursts could be generated by either whistler-mode or lower hybrid drift instabilities. It has previously been shown that the growth rate for whistler-mode instabilities in the nightside ionosphere is too small for appreciable growth at altitudes less than 200 km, where the VLF burst rate is highest. We show that the bursts are usually observed in regions of low electron beta, where whistler-mode attenuation is small. We further show that many of the bursts are detected in regions of high collision frequency, which stabilizes the lower hybrid drift instability. Lastly, the waves are also detected in regions where the wavelength required for Doppler-shift of lower hybrid waves to 100 Hz is shorter than the electron Larmor radius, which also argues against a lower hybrid drift instability. Planetary lightning is consequently a more likely source for the VLF bursts.

Description: The basic objectives of this NASA Grant are to develop theoretical understandings (tested with spacecraft data) of the generation and characteristics of electron plasma waves, commonly known as Langmuir-like waves, and associated radiation near f(sub p) and 2f(sub p) in planetary foreshocks. (Here f(sub p) is plasma frequency.) Related waves and radiation in the source regions of interplanetary type III solar radio bursts provide a simpler observational and theoretical context for developing and testing such understandings. Accordingly, applications to type III bursts constitute a significant fraction of the research effort. The testing of the new Stochastic Growth Theory (SGT) for type III bursts, and its extension and testing for foreshock waves and radiation, constitutes a major longterm strategic goal of the research effort.

Description: The objective of the proposal was to model the rototranslational and rotovibrational collision induced absorption spectral bands of importance for the radiative transfer analysis of the atmosphere of Venus. Our main task has involved CO2 pairs. The approach is not straightforward: whereas computational techniques to compute CIA spectra of small linear molecules exist, and were successfully applied to molecules like H2 or N2, they fail when applied to large molecules like CO2. For small molecules one can safely assume that the interaction potential is isotropic. The same approximation does not work for CO2, and when employed, it gives an incorrect band shape and only 50 percent of the CIA intensity.

Description: Progress in the development of the model for the circumplanetary distribution of atomic hydrogen in the Saturn system produced by a Titan source is discussed. Because of the action of the solar radiation acceleration and the obliquity of Saturn, the hydrogen distribution is shown to undergo seasonal changes as the planet moves about the Sun. Preliminary model calculations show that for a continuous Titan source, the H distribution is highly asymmetric about the planet and has a density maximum near the dusk side of Saturn, qualitatively similar to the pattern recently deduced by Shemansky and Hall from observations acquired by the UVS instruments aboard the Voyager spacecrafts. The investigation of these Voyager data will be undertaken in the next project year.

Description: The fast fracture strength distribution of uniaxially ground, alpha silicon carbide was investigated as a function of grinding angle relative to the principal stress direction in flexure. Both as-ground and ground/annealed surfaces were investigated. The resulting flexural strength distributions were used to verify reliability models and predict the strength distribution of larger plate specimens tested in biaxial flexure. Complete fractography was done on the specimens. Failures occurred from agglomerates, machining cracks, or hybrid flaws that consisted of a machining crack located at a processing agglomerate. Annealing eliminated failures due to machining damage. Reliability analyses were performed using two and three parameter Weibull and Batdorf methodologies. The Weibull size effect was demonstrated for machining flaws. Mixed mode reliability models reasonably predicted the strength distributions of uniaxial flexure and biaxial plate specimens.

Description: Probabilistic predictions based on the IPACS code are presented for the material and structural response of unnotched and notched, IM6/3501-6 Gr/Ep laminates. Comparisons of predicted and measured modulus and strength distributions are given for unnotched unidirectional, cross-ply and quasi-isotropic laminates. The predicted modulus distributions were found to correlate well with the test results for all three unnotched laminates. Correlations of strength distributions for the unnotched laminates are judged good for the unidirectional laminate and fair for the cross-ply laminate, whereas the strength correlation for the quasi-isotropic laminate is judged poor because IPACS did not have a progressive failure capability at the time this work was performed. The report also presents probabilistic and structural reliability analysis predictions for the strain concentration factor (SCF) for an open-hole, quasi-isotropic laminate subjected to longitudinal tension. A special procedure was developed to adapt IPACS for the structural reliability analysis. The reliability results show the importance of identifying the most significant random variables upon which the SCF depends, and of having accurate scatter values for these variables.

Description: The cyclic deformation behavior of a wrought cobalt-base superalloy, Haynes 188, has been investigated over a range of temperatures between 25 and 1000 C under isothermal and in-phase thermomechanical fatigue (TMF) conditions. Constant mechanical strain rates (epsilon-dot) of 10(exp -3)/s and 10(exp -4)/s were examined with a fully reversed strain range of 0.8%. Particular attention was given to the effects of dynamic strain aging (DSA) on the stress-strain response and low cycle fatigue life. A correlation between cyclic deformation behavior and microstructural substructure was made through detailed transmission electron microscopy. Although DSA was found to occur over a wide temperature range between approximately 300 and 750 C the microstructural characteristics and the deformation mechanisms responsible for DSA varied considerably and were dependent upon temperature. In general, the operation of DSA processes led to a maximum of the cyclic stress amplitude at 650 C and was accompanied by pronounced planar slip, relatively high dislocation density, and the generation of stacking faults. DSA was evidenced through a combination of phenomena, including serrated yielding, an inverse dependence of the maximum cyclic hardening with epsilon-dot, and an instantaneous inverse epsilon-dot sensitivity verified by specialized epsilon-dot -change tests. The TMF cyclic hardening behavior of the alloy appeared to be dictated by the substructural changes occuring at the maximum temperature in the TMF cycle.

Description: Support is requested for continuation of a program of dynamic impact (harpoon) coring of planetary, comet, or asteroid surface materials. We have previously demonstrated that good quality cores are obtainable for planetary materials with compressive strengths less than 200 MPa. Since the dynamics of penetration are observable on a Discovery class spacecraft, which images the sampling operation, these data can be used with a model developed under this project, to measure in-situ strength and frictional strength of the crust of the object. During the last year we have developed a detailed analytic model of penetrator mechanics. Progress is reported for the solid penetrators experiments, the CIT penetrator model, and the impact spall sampling apparatus.

Description: In spite of great improvements in accuracy through the use of computers, design methods, which can be equally critical in establishing the commercial success of a material, have been treated as afterthoughts. Early investment in design and development technologies can easily reduce manufacturing costs later in the product cycle. To avoid lengthy product development times for ceramic composites, funding agencies for materials research must commit resources to support design and development technologies early in the material life cycle. These technologies need not focus on designing the material, rather, the technology must focus on designing with the material, i. e., developing methods to design components fabricated from the new material. Thus a basic tenet that motivated this research effort is that a persistent need exists for improvements in the analysis of components fabricated from CMC material systems. From an aerospace design engineer's perspective the new generation of ceramic composites offers a significant potential for raising the thrust/weight ratio and reducing NOx emissions of gas turbine engines. Continuous ceramic fiber composites exhibit an increase in work of fracture, which allows for 'graceful' rather than catastrophic failure. When loaded in the fiber direction, these composites retain substantial strength capacity beyond the initiation of transverse matrix cracking despite the fact that neither of its constituents would exhibit such behavior if tested alone. As additional load is applied beyond first matrix cracking, the matrix tends to break in a series of cracks bridged by the ceramic fibers. Thus any additional load is born increasingly by the fibers until the ultimate strength of the composite is reached. Establishing design protocols that enable the engineer to analyze and predict this type of behavior in ceramic composites was the general goal of this project.

Description: We have carried out instrumental neutron activation analysis of 11 enstatite chondrites and electron microprobe analyses of 17 enstatite chondrites, most of which were previously little described. We report here the third known EH5 chondrite (LEW 88180) and an unusual EL6 chondrite (LEW 87119), new data on four EL3 chondrites (ALH 85119, EET 90299, PCA 91020, and MAC 88136, which is paired with MAC 88180 and MAC 88184), the second EL5 chondrite (TIL 91714), and an unusual metal-rich and sulfide-poor EL3 chondrite (LEW 87223). The often discussed differences in mineral composition displayed by the EH and EL chondrites are not as marked after the inclusion of the new samples in the database, and the two classes apparently experienced a similar range of equilibrium temperatures. However, texturally the EL chondrites appear to have experienced much higher levels of metamorphic alteration than EH chondrites of similar equilibration temperatures. Most of the petrologic type criteria are not applicable to enstatite chondrites and, unlike the ordinary chondrites, texture and mineralogy reflect different aspects of the meteorite history. We therefore propose that the existing petrologic type scheme not be used for enstatite chondrites. We suggest that while 'textural type' reflects peak metamorphic temperatures, the 'mineralogical type' reflects equilibration during postmetamorphic (probably regolith) processes. Unlike the ordinary chondrites and EH chondrites, EL chondrites experienced an extensive low-temperature metamorphic episode. There are now a large number of enstatite meteorite breccias and impact melts, and apparently surface processes were important in determining the present nature of the enstatite chondrites.

Description: The thermal emision of two palagonitic soils, common visible and near infrared spectral analogs for bright soils on Mars, was measured over the wavelength range of 5 to 25 micrometers (2000 to 400/cm) for several partical size separates. All spectra exhibit emissivity features due to vibrations associated with H2O and SiO. The maximum variability of emissivity is approximately 20% in the short wavelength region (5 to 6.5 mirometers, 2000 to 1500/cm), and is more subdued, less than 4%, at longer wavelengths. The strengths of features present in the infrared spectra of Mars cannot be solely provided by emissivity variations of palagonite; some other material or mechanism must provide additional absorptions(s).

Description: We have reanalyzed three sets of Viking Lander 1 and 2 (VL1 and VL2) images of the Martian atmosphere to better evaluate the radiative properties of the atmospheric dust particles. The properties of interest are the first two moments of the size distribution, the single-scattering albedo, the dust single-scattering phase function, and the imaginary index of refraction. These properties provide a good definition of the influence that the atmospheric dust has on heating of the atmosphere. Our analysis represents a significant improvement over past analyses (Pollack et al. 1977, 1979) by deriving more accurate brightness closer to the sun, by carrying out more precise analyses of the data to acquire the quantities of interest, and by using a better representation of scattering by nonspherical particles. The improvements allow us to better define the diffraction peak and hence the size distribution of the particles. For a lognormal particle size distribution, the first two moments of the size distribution, weighted by the geometric cross section, are found. The geometric cross-section weighted mean radius r(sub eff) is found to be 1.85 +/- 0.3 micrometers at VL2 during northern summer when dust loading was low and 1.52 +/- 0.3 micrometers at VL1 during the first dust storm. In both cases the best cross-section weighted mean variance nu(sub eff) of the size distribution is equal to 0.5 +/- 0.2 micrometers. The changes in size distribution, and thus radiative properties, do not represent a substantial change in solar energy deposition in the atmosphere over the Pollak et al. (1977, 1979) estimates.

Description: A time-domain study of the random response of a laminated plate subjected to combined acoustic and thermal loads is carried out. The features of this problem also include given uniform static inplane forces. The formulation takes into consideration a possible initial imperfection in the flatness of the plate. High decibel sound pressure levels along with high thermal gradients across thickness drive the plate response into nonlinear regimes. This calls for the analysis to use von Karman large deflection strain-displacement relationships. A finite element model that combines the von Karman strains with the first-order shear deformation plate theory is developed. The development of the analytical model can accommodate an anisotropic composite laminate built up of uniformly thick layers of orthotropic, linearly elastic laminae. The global system of finite element equations is then reduced to a modal system of equations. Numerical simulation using a single-step algorithm in the time-domain is then carried out to solve for the modal coordinates. Nonlinear algebraic equations within each time-step are solved by the Newton-Raphson method. The random gaussian filtered white noise load is generated using Monte Carlo simulation. The acoustic pressure distribution over the plate is capable of accounting for a grazing incidence wavefront. Numerical results are presented to study a variety of cases.

Description: Contrary to recent assertions in the literature, lunar emission spectra obtained at an altitude of 32 km with a balloon-borne telescope are virtually undistorted by atmospheric absorption, except in the ozone region. These spectra have been found to correlate closely with laboratory emission spectra of returned lunar samples measured in a simulated lunar environment. Thus, lunar spectra obtained with the balloon-borne telescope system are the standards against which groundbased spectral measurements of the Moon should be compared. Comparison of balloon-borne measurements with recent groundbased spectral measurements suggests that the latter suffer from some source of systematic error. When high-quality lunar spectra are available, they will be best interpreted by comparison with returned lunar samples, rather than terrestrial minerals or rocks, because of subtle spectral differences between lunar and terrestrial minerals perhaps associated with incipient alteration of the latter.

Description: The quantum yield and reaction threshold for the photochemical dissociation of cyanoacetylene into a hydrogen atom and the cyanoethynyl radical have been determined. The quantum yield at 185 nm is approximately 0.09. The threshold is approximately 240 nm. Combination of this data with literature values shows that production of excited-state cyanoacetylene is the major primary process resulting from irradiation between 185 and 254 nm. Also determined are the relative rate constants for the abstraction of a hydrogen atom from hydrogen, methane, and ethane by the cyanoethynyl radical (k(H2):k(CH4):k(C2H6) = 1:9.3:63). Implications of these results for the proposal that hydrogen abstraction plays an important role in the conversion of methane to ethane and in the protection of unsaturated compounds from photoconsumption in the atmosphere of Titan are discussed.

Description: The atmospheric flow on Io is numerically computed in a flat 2-D axisymmetric geometry for a sublimation atmosphere on the trailing hemisphere subjected to plasma bombardment, UV heating, and IR cooling. Calculations are performed for subsolar vapor pressures of approximately 6.5 x 10(exp -3) Pa (approximately 3 x 10(exp 18) SO2/sq cm) and 6.8 x 10(exp -4) Pa (approximately 4 x 10(exp 17) SO2/sq cm); the latter approximates the vapor pressure of F. P. Fanale et al. (1982). The amount of plasma energy deposited in the atmosphere is 20% of the plasma flow energy due to corotation (J. A. Linker et al., 1988). It is found that plasma heating significantly inflates the upper atmosphere, increasing both the exobase altitude and the amount of surface covered by more than an exospheric column of gas. This in turn controls the supply of the Io plasma torus (M. A. McGrath and R. E. Johnson, 1987). The horizontal flow of mass and energy is also important in determining the exobase altitude; and it is shown that IR cooling can be important, although our use of the equilibrium, cool-to-space approximation for a pure SO2 gas (E. Lellouch et al., 1992) may overestimate this effect. The calculated exobase altitudes are somewhat lower than those suggested by McGrath and Johnson (1987) for supplying the torus, indicating the details of the plasma energy deposition and sputter ejection rate near the exobase, as well as the IR emission from this region need to be examined. In addition, the molecules sublimed (or sputtered) from the surface are transported to the exobase in times short compared to the molecular photodissociation time. Therefore, the exobase is dominated by molecular species and the exobase is supplied by a small region of the surface.

Description: Iapetus shows a greater hemispheric albedo asymmetry than any other body in the solar system. Hapke scattering theory and optical constants measured in the laboratory are used to identify possible compositions for the dark material on the leading hemisphere of Iapetus. The materials considered are poly-HCN, kerogen, Murchison organic residue, Titan tholin, ice tholin, and water ice. Three-component mixtures of these materials are modeled in intraparticle mixture of 25% poly-HCN, 10% Murchison residue, and 65% water ice is found to best fit the spectrum, albedo, and phase behavior of the dark material. The Murchison residue and/or water ice can be replaced by kerogen and ice tholin, respectively, and still produce very good fits. Areal and particle mixtures of poly-HCN, Titan tholin, and either ice tholin or Murchison residue are also possible models. Poly-HCN is a necessary component in almost all good models. The presence of poly-HCN can be further tested by high-resolution observations near 4.5 micrometers.

Description: We prospose a new center and ring assignment for the original Chryse impact basin based upon photogeologic mapping of Noachian outcrops and re-examination of the published geology using orthographic projections. While others have centered the Chryse impact on the topographic low associated with Hesperian volcanic and fluvial deposits, we suggest that the center of the Noachian-age excavation cavity was located approximately 800 km to the north, and that the basin topography was modified significantly from the Noachian into the Hesperian. Evolution of the topographic low included structural modification by a later impact centered in Acidalia, restricted volcanic deposition and loading, localized subsidence, and restricted deposition from the circum-Chryse outlfow channels.

Description: This report describes a general plan and the pertinent technological requirements for TOPS (Toward Other Planetary Systems), a staged program to ascertain the prevalence and character of other planetary systems and to construct a definitive picture of the formation of stars and their planets. The first stages focus on discovering and studying a significant number of fully formed planetary systems, as well as expanding current studies of protoplanetary systems. As the TOPS Program evolves, emphasis will shift toward intensive study of the discovered systems and of individual planets. Early stages of the TOPS Program can be undertaken with ground-based observations and space missions comparable in scale to those now being performed. In the long term, however, TOPS will become an ambitious program that challenges our capabilities and provides impetus for major space initiatives and new technologies.

Description: The idea of searching for evidence of life on Mars may strike some as far-fetched, even fanciful. But there is a compelling logic to such a quest, as well as an equally compelling excitement. Early environments were apparently sufficiently similar on Mars and Earth, and life arose so rapidly on Earth once conditions became clement, that emergence of life on both planets at that time is scarcely less plausible than emergence on only one. Furthermore, although a fossil on Mars might seem at first like a proverbial needle in a haystack, experience on Earth tell us that if we know where to look, finding evidence of ancient life is not particularly difficult, especially when one considers that such evidence can be relatively widely disseminated in the form of chemical or isotopic signatures. The key is to recognize that the search for ancient life on Mars will involve a logically designed sequence of missions, each of which will focus on defining ever more closely where and how biosignatures may be found. Although one can never rule out a chance discovery, this quest should not be approached as one that will yield to a single, expeditious mission. (In fact, the proposed strategy lends itself particularly well to the use of a series of relatively small, inexpensive spacecraft, rather than a single flagship-class mission). The search for life on Mars will take time and commitment, but the reward could be a discovery of inestimable importance, not just to science, but to humanity as a whole.

Description: The spectral characterization of small asteroids is important for understanding the evolution of their compositional and mineralogical properties. We report the results of a CCD spectroscopic survey of small main-belt asteroids which we call the Small Main-belt Asteroid Spectroscopic Survey (SMASS). Spectra of 316 asteroids were obtained, with wavelength coverage ranging from 4000 to 10000 A (0.4 to 1 micrometers). More than half of the objects in our survey have diameters less than 20 km. Survey results include the identification of the first object resembling ordinary chondrite meteorites among the main-belt asteroids (Binzel, R. P., et al, 1993) and observations of more than 20 asteroids showing basaltic achondrite spectral absorption features that strongly link Vesta as the parent body for the basaltic achondrite meteorites (Binzel, R. P., and S. Xu 1993). A potential Mars-crossing asteroid analog to ordinary chondrite meteorites (H chondrites), 2078 Nanking, is reported here. Through a principal component analysis, we have assigned classifications to the members of our sample. The majority of the small main-belt asteroids belong to S and C classes, similar to large asteroids. Our analysis shows that two new classes are justified which we label as J and O. Small asteroids display more diversity in spectral absorption features than the larger ones, which may indicate a greater variation of compositions in the small asteroid population. We found a few candidates for olivine-rich asteroids within the S class. Although the total number of olivine-rich candidates is relatively small, we present evidence suggesting that such objects are more prevalent at smaller sizes.

Description: We report the results of a detailed experiment study of the kinetics and mechanism of pyrite (FeS2) chemical weathering under Venus surface conditions. Pyrite is thermodynamically unstable on the surface of Venus and will spontaneously decompose to pyrrhotite (Fe7S8) because the observed S2 partial pressure in the lower atmosphere of Venus is lower than the S2 vapor pressure over coexisting pyrite and pyrrhotite. Pyrite decomposition kinetics were studied in pure CO2 and CO2 gas mixtures along five isotherms in the temperature range 390-531 C. In all gas mixtures studied, pyrite thermally decomposes to pyrrhotite (Fe7S8), which on continued heating loses sulfur to form more Fe-rich pyrrhotites. During this process the pyrrhotites are also being oxidized to form magnetite (Fe3O4), which converts to maghemite (gamma-Fe2O3), and then to hematite (alpha-Fe2O3). The reaction rates for pyrite thermal decomposition to pyrrhotite were determined by measuring the weight loss. The thickness of the unreacted pyrite in the samples provided a second independent reaction rate measurement. Finally, Mossbauer spectra done on 42 of the 115 experimental samples provided a third set of independent reaction rate data. Pyrite decomposition follows zero-order kinetics and is independent of the amount of pyrite present. The rate of pyrite decomposition is apparently independent of the gas compositions used and of the CO2 number density over a range of a factor of 40. The derived activation energy of approximately 150 kJ/mole is the same in pure CO2, two different CO-CO2 mixtures, and a ternary CO-SO2-CO2 mixture. Based on data for a CO-CO2-SO2 gas mixture with a CO number density approximately 10 times higher than at the surface of Venus and a SO2 number density approximately equal to that at the surface of Venus, the rate of pyrite destruction on the surface of Venus varies from about 1225 +/- 238 days/cm at the top of Maxwell Montes (approximately 660 K) to about 233 +/- 133 days/cm in the plains of Venus (approximately 740 K). These lifetimes are very short on a geological time scale and show that pyrite cannot exist on the surface of Venus for any appreciable length of time.

Description: Observations of CO at 2.6 mm (115.27 GHz) were made with the Owens Valley Radio Observatory (OVRO) millimeter interferometer in 1986 and 1988, yielding high-quality disk-resolved spectra which were inverted to determine the CO mixing ratio profile from distinct regions on the disk, allowing us to map the distribution of CO in the upper mesosphere of Venus both horizontally and vertically. The 1986 observations were of the morning terminator and were particularly useful in searching for a suspected CO maximum ('bulge') on the nightside. The resulting CO mixing ratio profiles were mapped for various altitudes as functions of latitude and local time, and we report that we have resolved the previously inferred CO bulge. The bulge increases in magnitude from a small day-night variation at 90 km to an extensive nightside peak at 100 km, the upper limit of our observations. The peak bulge-to-dayside ratio approached 20-30 at 100 km in 1986 and may have been as large as 50-100, assuming late-afternoon CO abundances found in 1988 were similar to those in 1986. Three-dimensional mapping shows that in the upper mesosphere the bulge was displaced from local midnight toward the morning equator, centered at 3:30 AM local time. Using the qualitative model of mesospheric circulation on Venus proposed by Clancy and Muhleman (1985, 1991), we explain this shift in terms of strong retrograde zonal winds throughout the mesosphere, matching the directly detected mesospheric circulation (Shah, K., D. O. Muhleman, and G. L. Berge 1991) observed with the same dataset in 1988.

Description: Glassy objects reportedly found in Pima County, Arizona, have been identified as tektites. A field survey of the area where they were said to occur, however, did not produce any other tektites, nor did it reveal any other geologic features that might indicate a nearby impact crater. The major-, minor-, and trace-element composition of one specimen is similar to those measured in indochinites, which suggests the objects reportedly from Pima County were instead transported to southern Arizona from Indochina by people.

Description: The abundances of alkali elements in the Earth's core are predicted by assuming that accretion of the Earth started from material similar in composition to enstatite chondrites and that enstatite achondrites (aubrites) provide a natural laboratory to study core-mantle differentiation under extremely reducing conditions. If core formation on the aubrite parent body is comparable with core formation on the early Earth, it is found that 2600 (+/- 1000) ppm Na, 550 (+/- 260) ppm K, 3.4 (+/- 2.1) ppm Rb, and 0.31 (+/- 0.24) ppm Cs can reside in the Earth's core. The alkali-element abundances are consistent with those predicted by independent estimates based on nebula condensation calculations and heat flow data.

Description: A structurally ordered cubic metal grain containing approximately 39 wt% Co, 61 wt% Fe, and less than 0.6 wt% Ni (approximately Fe3Co2) was found associated with troilite and pentlandite in the matrix of the Ningqiang carbonaceous chondrite. This mineral is similar to terrestrial wairauite. Experimental data in the Fe-Co system indicate that this CsCl-type Co-rich metal is stable below 700 C. Phase relations in the Fe-Co-Ni system show that Co cannot fractionate from Ni above 500 C. The dominant opaque minerals of awaruite, magnetite, and pentlandite in Ningqiang suggest relatively oxidizing conditions.

Description: A spherical 220-micrometer diameter, spinel-hibonite-perovskite inclusion from the Allende C3V meteorite contains a central hibonite cluster with an angular coundary. This central hibonite is enclosed within spinel that is zoned from Mg-rich at the hibonite boundary to more Fe-rich at the inclusion boundary. This spinel zone includes lath-shaped hibonites usually subradial to the central hibonites. Two textural types of perovskites are present as exsolution from the central hibonite and as equidimensional grains within both the central hibonite and spinel. These second perovskites have exsolution lamellae of Al2O3. Within the central hibonite and adjacent to some equidimensional perovskites, a fine porous phase interpreted as alteration has a composition of nearly pure Al2O3 with minor amounts of Na and Si. This is possibly either an intergrowth of corundum and nepheline or a modified Al2O3, Beta-alumina.

Description: Some fraction of Zn, Cu, Se, Ga and Ge in chondritic interplanetary dust particles (IDPs) collected in the lower stratosphere between 1981 May and 1984 June has a volcanic origin. I present a method to evaluate the extent of this unavoidable type of stratospheric contamination for individual particles. The mass-normalized abundances for Cu and Ge as a function of mass-normalized stratospheric residence time show their time-integrated stratospheric aerosol abundances. The Zn, Se and Ga abundances show a subdivision into two groups that span approximately two-year periods following the eruptions of the Mount St. Helens (1980 May) and El Chichon (1982 April) volcanoes. Elemental abundances in particles collected at the end of each two-year period indicate low, but not necessarily ambient, volcanic stratospheric abundances. Using this time-integrated baseline, I calculate the straospheric contaminant fractions in nine IDPs and show that Zn, SE and Ga abundances in chondritic IDPs derive in part from stratospheric aerosol contaminants. Post-entry elemental abundances (i.e., the amount that survived atmospheric entry heating of the IDP) show enrichments relative to the CI abundances but in a smaller number of particles than previously suggested.

Description: Emissivity spectra of particulate mineral samples are highly dependent on particle size when that size is comparable to the wavelength of light emitted (5-50 micrometers for the midinfrared). Proper geologic interpretation of data from planetary infrared spectrometers will require that these particle size effects be well understood. To address this issue, samples of quartz powders were produced with narrow, well-characterized particle size distributions. Mean particle diameters in these samples ranged from 15 to 227 micrometers. Emission spectra of these powders allow the first detailed comparison of the complex spectral variations with particle size observed in laboratory data with the predictions of radiative transfer models. Four such models are considered here. Hapke's relectance theory (converted to emissivity via Kirchoff's law) is the first model tested. Hapke's more recently published emission theory is also employed. The third model, the 'Mie/Conel' model, uses Mie single scattering with a two-stream approximation for multiple scattering. This model, like the first, is a converted reflec- tance model. Mie scattering assumes particles are both spherical and well separated, which is not true for the quartz powders, but includes diffraction effects. The fourth model uses the Mie solution for single scattering by spheres and inputs those results into the multiple scattering formalism of Hapke's emission theory. The results of the four models are considered in relation to the values of the optical constants n and k. We have grouped these as class 1 (k large), class 2 (k moderate, n is approximately 2), class 3 (k small, n is approximately 2), and class 4 (k small, n is approximately 1). In general, the Mie/Hapke hybrid model does best at predicting variations with grain size. In particular, it predicts changes of the correct pattern, although incorrect magnitude, for class 1 bands, where large increases in emissivity with decreasing grain size are observed. This model also does an excellent job on moderate (class 2) and very weak and intraband (class 3) regions, and correctly predicts the emission maximum and its invariance with grain size near the Christiansen frequency (class 4). The Mie/Hapke hybrid model also has the fewest free parameters of the four models examined, while maintaining the most physical treatment of the radiative transfer. The Mie/Conel model performs as well as the Mie/Hapke hybrid model in strong bands (class 1) but does not accurately model the behavior of moderate (class 2) and very weak (class 3) bands.

Description: The identification of the mechanism(s) responsible for heating asteroids is among the major problems in planetary science. Because of difficulties with models of electromagnetic induction and the decay of short-lived radionuclides, it is worthwhile to evaluate the evidence for collisional heating. New evidence for localized impact heating comes from the high proportion of relict type-6 material among impact-melt-bearing ordinary chondrites (OC). This relict material was probably metamorphosed by residual heat within large craters. Olivine aggregates composed of faceted crystals with 120 deg triple junctions occur within the melted regions of the Chico and Rose City OC melt rocks; the olivine aggregates formed from shocked, mosaicized olivine grains that underwent contact metamorphism. Large-scale collisional heating is supoorted by the correlation in OC between petrologic type and shock stage; no other heating mechanism can readily account for this correlation. The occurrence of impact-melt-rock clasts in OC that have been metamorphosed along with their whole rocks indicates that some impact events preceded or accompanied thermal metamorphism. Such impacts events, occurring during or shortly after accretion, are probably responsible for substantially melting approximately 0.5% of OC. These events must have heated a larger percentage of OC to subsolidus temperatures sufficient to have caused significant metamorphism. If collisional heating is viable, then OC parent asteroids must have been large; large OC asteroids in the main belt may include those of the S(IV) spectral subtype. Collisional heating is inconsistent with layered ('onion-shell') structures in OC asteroids (wherein the degree of metamorphism increases with depth), but the evidence for such structures is weak. It seems likely that collisional heating played an important role in metamorphosing chondritic asteroids.

Description: We report data for the trace elements Au, Co, Sb, Ga, Rb, Ag, Se, Cs, Te, Zn, Cd, Bi, T1, and In (ordered by putative volatility during nebular condensation and accretion) determined by radiochemical neutron activation analysis of 14 additional H5 and H6 chondrite falls. Data for the 10 most volatile elements (Rb to In) treated by the multivariate techniques of linear discriminant analysis and logistic regression in these and 44 other falls are compared with those of 59 H4-6 chondrites from Antarctica. Various populations are tested by the multivariate techniques, using the previously developed method of randomization-simulation to assess significance levels. An earlier conclusion, based on fewer examples, that H4-6 chondrite falls are compositionally distinguishable from the Antarctic suite is verified by the additional data. This distinctiveness is highly significant because of the presence of samples from Victoria Land in the Antarctic population, which differ compositionally from falls beyond any reasonable doubt. However, it cannot be proven unequivocally that falls and Antarctic samples from Queen Maud Land are compositionally distinguishable. Trivial causes (e.g., analyst bias, weathering) cannot explain the Victoria Land (Antarctic)/non-Antarctic compositional difference for paradigmatic H4-6 chondrites. This seems to reflect a time-dependent variation of near-Earth meteoroid source regions differing in average thermal history.

Description: We report Cl-36 (301-kyr half-life) data obtained by accelerator mass spectrometry allowing nominal terrestrial ages to be determined for 39 Antarctic H4-6 chondrites for which contents of volatile trace elements are known. The compositional difference between these Antarctic meteorites and 58 non-Antarctic falls increases with terrestrial age and, using multivariate statistical techniques, becomes highly significant for Antarctic samples with ages greater than 50 kyr. The compositional difference is inconsistent with trivial causes such as weathering and seems to reflect differences in thermal histories of parent sources. Temporal source variations for the H chondrite flux on Earth thus exist not only on a short-term, 40 years, basis (Dodd et al., 1993) but also on a long-term, greater than 50 kyr, basis.

Description: We report data for the trace elements Au, Co, Sb, Ga, Rb, Ag, Se, Cs, Te, Zn, Cd, Bi, Ti, and In (ordered by putative volatility during nebular condensation and accretion) determined by neutron activation analysis in 13 H5 chondrites from Victoria Land and 20 H4-6 chondrites from Queen Maud Land, Antarctica. These and earlier results provide Antarctic sample suites of 34 chondrites from Victoria Land and 25 from Queen Maud Land. Treatment of data for the most volatile 10 elements (Rb to In) in these studies by multivariate statistical techniques more robust, as well as more conservative, than conventional linear discriminant analysis and logistic regression demonstrates that compositions differ at marginally significant levels. This difference cannot be explained by trivial (terrestrial) causes and becomes more significant, despite the smaller size of the database, when comparisons are limited to data from a single analyst and when all upper limits are eliminated from consideration. The Victoria Land and Queen Maud Land suites have different mean terrestrial ages (approximately 300 kyr and approximately 100 kyr, respectively) and age distributions, suggesting that a time-dependent variation of chondritic sources with different thermal histories is responsible. As a result, these two Antarctic suites are, on average, chemically distinguishable from each other. Since H chondrites serve as a paradigm for other meteorite classes, these results indicate that the near-Earth populations of planetary materials varied with time on the 10(exp 5)-year timescale.

Description: Models of planet formation and of the orbital stability of planetary systems are used to predict the variety of planetary and satellite systems that may be present within our galaxy. A new approximate global criterion for orbital stability of planetary systems based on an extension of the local resonance overlap criterion is proposed. This criterion implies that at least some of Uranus' small inner moons are significantly less massive than predicted by estimates based on Voyager volumes and densities assumed to equal that of Miranda. Simple calculations (neglecting planetary gravity) suggest that giant planets which acrete substantial amounts of gas while their envelopes are extremely distended ultimately rotate rapidly in the prgrade direction.

Description: The approximately concentric arrangement of layered deposits and dune fields at the two Martian poles may reflect a nearly steady state dispersal of material from the polar deposits. Data on effective surface winds from high resolution Viking Images combined with theory of local winds suggest that the northern dunes are in part confined to a latitude band by winds generated by their own low albedo. Dispersal of the dark sand from the southern polar region is not subject to this kind of feedback because the irregular topography prevents areal accumulations sufficiently extensive to produce winds.

Description: The circulation of the Martian atmosphere during late southern summer is derived from atmospheric temperature and dust distributions retrieved from a subset of the Mariner 9 infrared interferometer spectrometer (IRIS) thermal emission spectra (LS = 343-348 deg) (Santee and Crisp, 1933). Zonal-mean zonal winds are calculated by assuming gradient wind balance and zero surface zonal wind. Both hemispheres have intense midlatitude westerly jets with velocities of 80-90 m/s near 50 km; in the southern tropics the winds are easterly with velocities of 40 m/s near 50 km. The net effect of the zonal-mean meridional circulation and large-scale waves can be approximated by the diabatic ciculation, which is defined from the atmospheric thermal structure and net radiative heating rates. The radiative transfer model described by Crisp (1990) and Santee (1993) is used to compute solar heating and thermal cooling rates from diurnal averages of the retrieved IRIS temperature and dust distributions. At pressures below 4 mbar, there are large net radiative heating rates (up to 5 K/d) in the equatorial region and large net radiative cooling rates (up to 12 K/d) in the polar regions. These net radiative heating rates are used in a diagnostic stream function model which solves for the meridonal and vertical components of the diabatic circulation simultaneously. We find a two cell circulation, with rising motion over the equator, poleward flow in both hemispheres, sinking motion over both polar regions, and return flow in the lowest atmospheric levels. The maximum poleward velocity is 3 m/s in the tropics at approximately 55 km altitude, and the maximum vertical velocity is 2.5 cm/s downward over the north pole at approximately 60 km altitude. If these large transport rates are sustained for an entire season, the Martian atmosphere above the 1-mbar level is overturned in about 38 days. This diabatic circulation is qualitatively similar to the terrestial diabatic circulation at the comparable season, but is more vigorous.

Description: Infrared reflectance spectra of carefully selected Mars soil analog materials have been measured under low atmospheric pressures and temperatures. Chemically altered montmorillonites containing ferrihydrite and hydrated ferric sulfate complexes are examined, as well as synthetic ferrihydrate and a palagonitic soil from Haleakala, Maui. Reflectance spectra of these analog materials exhibit subtle visible to near-infrared features, which are indicative of nanophase ferric oxides or oxyhydroxides and are similar to features observed in the spectra of the bright regions of Mars. Infrared reflectance spectra of these analogs include hydration features due to structural OH, bound H2O and adsorbed H2O. The spectal character of these hydration features is highly dependent on the sample environment and on the nature of the H2O/OH in the analogs. The behavior of the hydration features near 1.9 micrometers, 2.2 micrometers, 2.7 micrometers, 3 micrometers, and 6 micrometers are reported here in spetra measured under Marslike atmospheric environment. In spectra of these analogs measured under dry Earth atmospheric conditions the 1.9-micrometer band depth is 8-17%; this band is much stonger under moist conditions. Under Marslike atmospheric conditions the 1.9-micrometer feature is broad and barely discernible (1-3% band depth) in spectra of the ferrihydrite and palagonitic soil samples. In comparable spectra of the ferric sulfate-bearing montmorillonite the 1.9-micrometer feature is also broad, but stronger (6% band depth). In the low atmospheric pressure and temperature spectra of the ferrihydrite-bearing montmorillonite this feature is sharper than the other analogs and relatively stronger (6% band depth). Although the intensity of the 3- micrometer band is weaker in spectra of each of the analogs when measured under Marslike conditions, the 3-micromter band remains a dominant feature and is especially broad in spectra of the ferrihydrite and palagonitic soil. The structural OH features observed in these materials at 2.2-2.3 micrometers and 2.27 micrometers remain largely unaffected by the environmental conditions. A shift in the Christiansen feature towards shorter wavelengths has also been observed with decreasing atmospheric pressure and temperature in the midinfrared spectra of these samples.

Description: The first thermal images of Phobos' shadow on the surface of Mars, in addition to simultaneous visible images, were obtained by the Phobus '88 Termoskan instrument. The best observed shadow occurence was on the flanks of Arsia Mons. For this occurence, we combined the observed decrease in visible illumination of the surface with the observed decrease in brightness temperature to calculate thermal inertias of the Martian surface. The most realistic of our three models of eclipse cooling improves upon our preliminary model by including nonisothermal initial conditions and downward atmospheric flux. Most of our derived inertias fall within the range 38 to 59 J/Sq m/S(exp 0.5)K (0.9 to 1.4 10(exp -3)Cal/Sq m/S(exp 0.5)/K), corresponding to dust-sized particles (for a homogeneous surface), consistent with previous theories of Tharsis as a currrent area of dust deposition. Viking infrared thermal mapper (IRTM) inertias are diurnally derived and are sensitive to centimeter depths, whereas the shadow-derived inertias sample the upper tenths of a millimeter of the surface. The shadow-derived inertias are lower than those derived from Viking IRTM measurements (84 to 147), however, uncertainties in both sets of derived inertias make conclusions about layering tenuous. Thus, near-surface millimeter versus centimeter layering may exist in this region, but if it does, it is likely not very significant. Both eclipse and diurnal inertias appear to increase near the eastern end of the shadow occurence. We also analyzed a shadow occurence near the crater Herschel that showed no observed cooling. This analysis was limited by cool morning temperatures and instrument sensitivity, but yielded a lower bound of 80 on eclipse inertias in that region. Based upon our results, we strongly recommend future spacecraft thermal observations of Phobus' shadow, and suggest that they will be most useful if they improve upon Termoskan's geographic and temporal coverage and its accuracy.

Description: The transection and superposition relationships among channels, chaos, surface materials units, and other features in the circum-Chryse region of Mars were used to evaluate relative age relationships and evolution of flood events. Channels and chaos in contact (with one another) were treated as single discrete flood-carved systems. Some outflow channel systems form networks and are inferred to have been created by multiple flood events. Within some outflow channel networks, several separate individual channel systems can be traced to a specific chaos which acted as flood-source area to that specific flood channel. Individual flood-carved systems were related to widespread materials units or other surface features that served as stratigraphic horizons. Chryse outflow channels are inferred to have formed over most of the perceivable history of Mars. Outflow channels are inferred to become younger with increasing proximity to the Chryse basin. The relationship of subsequent outflow channel sources to the sources of earlier floods is inferred to disfavor episodic flooding due to the progresssive tapping of a juvenile near-surface water supply. Instead, we propose the circum-Chryse region as a candidate site of past hydrological recycling. The discharge rates necessary to carve the circum-Chryse outflow channels would have inevitably formed temporary standing bodies of H2O on the Martian surface where the flood-waters stagnated and pooled (the Chryse basin is topographically enclosed). These observations and inferences have led us to formulate and evaluate two hypotheses. Our numerical evaluations indicate that of these two hypotheses formulated, the groundwater seep cycle seems by far the more viable. Further observations from forthcoming missions may permit the determination of which mechanisms may have operated to recycle the Chryse flood-waters.

Description: We propose key modifications to the Toon et al. (1977) model of the particle size distribution and composition of Mars atmospheric dust, based on a variety of spacecraft and wavelength observations of the dust. A much broader (r(sub eff) variance approximately 0.8 micrometers), smaller particle size (r(sub mode) approximately 0.02 micrometers) distribution coupled with a 'palagonite-like' composition is argued to fit the complete ultraviolet-to-30-micrometer absorption properties of the dust better than the montmorillonite-basalt, r(sub eff) variance = 0.4 micrometers, r(sub mode) = 0.40 dust model of Toon et al. Mariner 9 (infrared interferometer spectrometer) IRIS spectra of high atmospheric dust opacities during the 1971-1972 Mars global dust storm are analyzed in terms of the Toon et al. dust model, and a Hawaiian palagonite sample (Rousch et al., 1991) with two different size distribution models incorporating smaller dust particle sizes. Viking Infrared Thermal Mapper (IRTM) emmission-phase-function (EPF) observations at 9 micrometers are analyzed to retrieve 9-micrometer dust opacities coincident with solar band dust opacities obtained from the same EPF sequences (Clancy and Lee, 1991). These EPF dust opacities provide an independent measurement of the visible/9-micrometer extinction opacity ratio (greater than or = 2) for Mars atmospheric dust, which is consistent with a previous measurement by Martin (1986). Model values for the visible/9-micrometer opacity ratio and the ultraviolet and visible single-scattering albedos are calculated for the palagonite model with the smaller particle size distributions compared to the same properties for the Toon et al. model of dust. The montmorillonite model of the dust is found to fit the detailed shape of the dust 9-micrometer absorption well. However, it predicts structured, deep aborptions at 20 micrometers which are not observed and requires a separate ultraviolet-visible absorbing component to match the observed behavior of the dust in this wavelength region. The modeled palagonite does not match the 8-to 9-micrometer absorption presented by the dust in the IRIS spectra, probably due to its low SiO2 content (31%). However, it does provide consistent levels of ultraviolet/visible absorption, 9-to 12-micrometer absorption, and a lack of structured absorption at 20 micrometers. The ratios of dust extinction opacities at visible, 9 micrometers, and 30 micrometers are strongly affected by the dust particle size distribution. The Toon et al. dust size distribution (r(sub mode) = 0.40,r(sub eff) variance = 0.4 micrometers, r(sub cwmu) = 2.7 micrometers) predicts the correct ratio of the 9- to 30-micrometer opacity, but underpredicts the visible/9-micrometer opacity ratio considerably (1 versus greater than or = 2). A similar particle distribution width with smaller particle sizes (r(sub mode) = 0.17, r(sub eff) variance = 0.4 micrometers, r(sub cwmu) = 1.2 micrometers) will fit the observed visible/9-micrometer opacity ratio, but overpredicts the observed 9-micrometer/30-micrometer opacity ratio. A smaller and much broader particle size distribution (r(sub mode) = 0.002, r(sub eff) variance = 0.8 micrometers, r(sub cwmu) = 1.8 micrometers) can fit both dust opacity ratios. Overall, the nanocrystalline structure of palagonite coupled with a smaller, broader distribution of dust particle sizes provides a more consistent fit than the Toon et al. model of the dust to the IRIS spectra, the observed visible/9-micrometer dust opacity ratio, the Phobos occulation measurements of the dust particle sizes (Chassefiere et al., 1992), and the weakness of surface near IR absorptions expected for clay minerals (Clark, 1992; Bell and Crisp, 1993).

Description: Recent calculations of the Martian obliquity suggests that it varies chaotically on timescales longer than about 10(exp 7) years and varies between about 0 and 60 deg. We examine the seasonal water behavior at obliquities between 40 and 60 deg. Up to several tens of centimeters of water may sublime from the polar caps each year, and possibly move to the equator, where it is more stable. The CO2 frost and CO2-H2O clathrate hydrate are stable in thepolar deposits below a few tens of meters depth, so that the polar cap could contain a significant CO2 reservoir. If CO2 is present, it could be left over from the early history of Mars; also, it could be released into the atmosphere during periods of high obliquity, causing occasional periods of more-clement climate.

Description: The high vapor pressure of volatile metal halides and chalcogenides (e.g., of Cu, Zn, Sn, Pb, As, Sb, Bi) at typical Venus surface temperatures, coupled with the altitude-dependent temperature gradient of approximately 8.5 K/km, is calculated to transport volatile metal vapors to the highlands of Venus, where condensation and accumulation will occur. The predicted geochemistry of volatile metals on Venus is supported by observations of CuCl in volcanic gases at Kilauea and Nyiragongo, and large enrichments of these and other volatile elements in terrestrial volcanic aerosols. A one-dimensional finite difference vapor transport model shows the diffusive migration of a thickness of 0.01 to greater than 10 microns/yr of moderately to highly volatile phases (e.g., metal halides and chalcogenides) from the hot lowlands (740 K) to the cold highlands (660 K) on Venus. The diffusive transport of volatile phases on Venus may explain the observed low emissivity of the Venusian highlands, hazes at 6-km altitude observed by two Pioneer Venus entry probes, and the Pioneer Venus entry probe anomalies at 12.5 km.

Description: A rift over 6000 km in total length runs along the border of Lada Terra, a highland in the southern hemisphere of Venus, and Lavinia Planitia, a basin that has been interpreted as a site of early-stage mantle downwelling. Along the length of the rift are a number of volcanic centers of widely varying morphology and volcanic output. These include coronae, radially fractured domes, and large flow fields similar in scale to terrestrial flood basalts. We develop a model for the origin of extension related to passive rifting in response to stresses created by the adjacent downwelling. Volcanism and extension at other rifts on Venus, such as Devana Chasma, have been attributed to deep-seated mantle plume activity. In contrast, we interpret the origin of extension and volcanism along the Lada rift to be linked to upwelling and decompression melting of mantle material due to rifting and, possibly, to counterflow associated with downwelling. Extension occurred generally prior to the formation of volcanic centers and the eruption of large-scale flow fields, although most of the volcanic centers have been fractured by continued extension along the rift. Current debate over the formation of terrestrial flood basalts centers on the necessity of preexisting extension and stretched and thinned lithosphere to produce enhanced decompression melting within a large plume head or mantle thermal anomaly. Our studies of large-scale flow fields associated with the Lada rift and coronae on Venus indicate that extension is a prerequisite for the formation of the majority of large-scale flow units on Venus.

Description: Although the Phobos-2 spacecraft recently obtained important results relevant to some of the major remaining questions in Mars aeronomy, much remains to be done. In particular, not since the Viking Landers have we made in-situ measurements of aeronomical quantities such as atmospheric and ionospheric densities and temperatures below 400 km altitude. We have never made magnetic field measurements at these altitudes. Without such measurements we cannot unambiguously resolve arguments concerning issues such as the significance of the planetary magnetic field in the solar wind interaction, or understand the atmospheric cycle that leads to escape to space. With the trio of future orbiters including Mars Observer, Mars-94, and Planet-B we should see a veritable explosion of new knowledge, but some gaps in aeronomical science coverage will still remain. This paper briefly reviews some of the major unsolved problems in Mars aeronomy, and points out which are expected to remain outstanding after this flotilla of missions.

Description: A general theoretical overview for the sources, sinks, gas-surface interactions, and transport dynamics of sodium and potassium in the exospheric atmsophere of Mercury is given. Information for these four factors, which control the spatial distribution of these two alkali-group gases about the planet, is incorporated in numerical models. The spatial nature and relative importance of the initial source atom atmosphere and the ambient (ballistic hopping) atom atmosphere are then examined and are shown to be controlled and coupled to a great extent by the extremely large and variable solar radiation acceleration experienced by sodium and potassium as they resonantly scatter solar photons. The lateral (antisunward) transport rate of thermally accommodated sodium and potassium ambient atoms is shown to be driven by the solar radiation acceleration and, over a significant portion of Mercury's orbit about the Sun, is sufficiently rapid to be competitive with the short photoionization lifetimes for these atoms when they are located on the summit surface near or within about 30 deg of the terminator. The lateral transport rate is characterized by a migration time determined by model calculations for an ensemble of atoms initially starting at a point source on the surface (i.e., a numerical spacetime dependent Green's function). Four animations for the spacetime evolution of the sodium (or potassium) atmosphere produced by a point source on the surface are presented on a videotape format. For extended surface sources for sodium and potassium, the local column density is determined by competition between the photoionization lifetimes and the lateral transport times of atoms originating from different surface source locations. Sodium surface source fluxes (referenced to Mercury at perihelion) that are required on the sunlit hemisphere to reproduce the typically observed several megarayleighs of D2 emission-line brightness and the inferred column densities of 1-2 x 10(exp 11) atoms per sq cm range from approximately 2-5 x 10(exp 7) atoms/sq cm/sec. The sodium model is applied to study observational data that document an anticorrelation in the average sodium column density and solar radiation acceleration. Lateral transport driven by the solar radiation acceleration is shown to produce this behavior for combinations of different sources and surface accomodation coefficients. The best fit model fits to the observational data require a significant degree of thermal accommodation of the ambient sodium atoms to the surface and a source rate that decreases as an inverse power of 1.5 to 2 in heliocentric distance.

Description: The spatial variation of the geoid/topography ratio over the large Venusian volcanic highland Beta Regio is suggestive of thermal compensation, i.e., support of the highland's topography by lithospheric thinning. Both the thickness of the lithosphere and the density contrast at its base can be inferred from a quadratic regression of suitably filtered (600 km less than wavelength less than 4000 km) geoid vs. topography data. The regression yields a mean lithospheric thickness of 270 km and a density contrast of magnitude 2.5% to 3.0%. Simple isostatic balance of the long-wavelength topography at Beta Regio requires thinning of the lithosphere by 50-60% beneath the rise.

Description: The earliest evolution of the Moon likely included the formation of a magma ocean and the subsequent development of anorthositic flotation cumulates. This primary anorthositic crust was then intruded by mafic magmas which crystallized to form the lunar highlands magnesian suite. The present study is a compilation of petrologic, mineral-chemical, and geochemical information on all pristine magnesian-suite plutonic rocks and the interpretation of this data in light of 18 'new' samples. Of these 18 clasts taken from Apollo 14 breccias, 12 are probably pristine and include four dunites, two norites, four troctolites, and two anorthosites. Radiogenic isotopic whole rock data also are reported for one of the 'probably pristine' anorthositic troctolites, sample 14303,347. The relatively low Rb content and high Sm and Nd abundances of 14303,347 suggest that this cumulate rock was derived from a parental magma which had these chemical characteristics. Trace element, isotopic, and mineral-chemical data are used to interpret the total highlands magnesian suite as crustal precipitates of a primitive KREEP (possessing a K-, rare earth element (REE)-, and P-enriched chemical signature) basalt magma. This KREEP basalt was created by the mixing of ascending ultramafic melts from the lunar interior with urKREEP (the late, K-, REE-, and P-enriched residuum of the lunar magma ocean). A few samples of the magnesian suite with extremely elevated large-ion lithophile elements (5-10x other magnesian-suite rocks) cannot be explained by this model or any other model of autometasomatism, equilibrium crystallization, or 'local melt-pocket equilibrium' without recourse to an extremely large-ion lithophile element-enriched parent liquid. It is difficult to generate parental liquids which are 2-4 x higher in the REE than average lunar KREEP, unless the liquids are the basic complement of a liquid-liquid pair, i.e., the so-called 'REEP-fraction,' from the silicate liquid immiscibility of urKREEP. Scarce age information on lunar rocks suggests that magnesian-suite magmatism was initiated at progressively more recent time from the northeast to the southwest on the lunar nearside from 4.45 to 4.25 Ga.

Description: A technique for estimating the ground and near-ground atmospheric temperatures within a Martian local dust storm is presented. It is applied to soundings taken by the Viking orbiter infrared thermal mapper (IRTM) instrument at four times-of-day for one storm. Essentially, a comparison is made between infrared radiances emerging from the storm interior and those from the region surrounding the storm. Particle extinction properties are assumed to be independent of position in the storm region, and scattering properties must be selected arbitrarily. For the storm studied here, the ground temperature in the interior is at least 6 K cooler, whereas the near-ground atmospheric temperature may be less than or comparable to, those of the surroundings. The thermal structure of the storm interior did not change measurably between 11.5 and 16.6 hours local time. These observations favor the theories of dust storm development in which regional winds rather than local, dust driven convection initiate the mobilization of dust from the surface. It also concluded that the optical properties of dust particles in this local storm differ from those observed by Mariner 9 during the 1971-1972 global dust storm.

Description: We present the first results of a new technique for imaging the lunar Na atmosphere. The technique employs high resolution, a narrow bandpass, and specific observing geometry to suppress scattered light and image lunar atmospheric Na I emission down to approximately 50 km altitude. Analysis of four latitudinally dispersed images shows that the lunar Na atmosphere exhibits intersting latitudinal and radial dependencies. Application of a simple Maxwellian collisionless exosphere model indicates that: (1) at least two thermal populations are required to adequately fit the soldium's radial intensity behavior, and (2) the fractional abundances and temperatures of the two components vary systematically with latitude. We conclude that both cold (barometric) and hot (suprathermal) Na may coexist in the lunar atmosphere, either as distinct components or as elements of a continuum of populations ranging in temperature from the local surface temperature up to or exceeding escape energies.

Description: We have measured significant concentrations of Cl-36, Ca-41, Ar-36 from decay of Cl-36, and Sm-150 produced from the capture of thermalized neutrons in the large Chico L6 chondrite. Activities of Cl-36 and Ca-41, corrected for a high-energy spallogenic component and a terrestrial age of approximately 50 ka, give average neutron-capture production rates of 208 atoms/min/g-Cl and 1525 atoms/min/kg-Ca, which correspond to thermal neutron (n) fluxes of 6.2 n/sq cm/s and 4.3 n/sq cm/s, respectively. If sustained for the approximately 65 Ma single-stage, cosmic ray exposure age of Chico, these values correspond to thermal neutron fluences of approximately 1.3 x 10(exp 16) and 0.8 x 10(exp 16) n/sq cm for Cl-36 and Ca-41, respectively. Stepwise temperature extraction of Ar in Chico impact melt shows Ar-36/Ar-38 ratios as large as approximately 9. The correlation of high Ar-36/Ar-38 with high Cl/Ca phases in neutron-irradiated Chico indicates that the excess Ar-36 above that expected from spallation is due to decay of neutron-produced Cl-36. Excess Ar-36 in Chico requires a thermal neutron fluence of 0.9-1.7 x 10(exp 16) n/sq cm. Decreases in Sm-149/Sm-152 due to neutron-capture by Sm-149 correlate with increases in Sm-150/Sm-152 for three samples of Chico, and one of the Torino H-chondrite. The 0.08% decrease in Sm-149 shown by Chico corresponds to a neutron fluence of 1.23 x 10(exp 16) n/sq cm. This fluence derived from Sm considers capture of epithermal neutrons and effects of chemical composition on the neutron energy distribution. Excess Ar-36 identified in the Arapahoe, Bruderheim, and Torino chondrites and the Shallowater aubrite suggest exposure to neutron fluences of approximately 0.2-0.2 x 10(exp 16) n/sq cm. Depletion of Sm-149 in Torino and the LEW86010 angrite suggest neutron fluences of 0.8 x 10(exp 16) n/sq cm and 0.25 x 10(exp 16) n/sq cm, respectively. Neutron fluences of approximately 10(exp 16) n/sq cm in Chico are almost as large as those previously observed for some lunar soils. Consideration of exposure ages suggests that the neutron flux in Chico may have been greater than that in many lunar soils.

Description: All the available nightside electron data obtained during circular orbits at Mars from the Phobos-2 Hyperbolic Retarded Potential Analyzer (HARP) instrument have been examined in detail and are summarized in this paper. An electron flux component with energies exceeding that of the unperturbed solar wind was observed inside the magnetosheath, indicating the presence of acceleration mechanism(s). The character of the electron fluxes measured in the magnetotail cannot be classified in any simple manner, however, there is a correlation between the electron fluxes measured well inside this region and the unperturbed solar wind ram pressure.

Description: Dust grains in planetary rings acquire stochastically fluctuating electric charges as they orbit through any corotating magnetospheric plasma. Here we investigate the nature of this stochastic charging and calculate its effect on the Lorentz resonance (LR). First we model grain charging as a Markov process, where the transition probabilities are identified as the ensemble-averaged charging fluxes due to plasma pickup and photoemission. We determine the distribution function P(t;N), giving the probability that a grain has N excess charges at time t. The autocorrelation function tau(sub q) for the strochastic charge process can be approximated by a Fokker-Planck treatment of the evolution equations for P(t; N). We calculate the mean square response to the stochastic fluctuations in the Lorentz force. We find that transport in phase space is very small compared to the resonant increase in amplitudes due to the mean charge, over the timescale that the oscillator is resonantly pumped up. Therefore the stochastic charge variations cannot break the resonant interaction; locally, the Lorentz resonance is a robust mechanism for the shaping of etheral dust ring systems. Slightly stronger bounds on plasma parameters are required when we consider the longer transit times between Lorentz resonances.

Description: We calculated the vapor flux of the icy components in the surface layer of a porous, short-period, Jupiter-class comet, in order to investigate the relationship of the observed relative molecular abundances in the coma with those in the nucleus. The model assumes a body containing one major ice component (H20) and up to three minor components of higher volatility (e.g., CO, CO2, CH3OH). The body's porous structure is modeled as a bundle of tubes with a given tortuosity and initially a constant pore diameter. The mass and energy equations for the different volatiles are solved simultaneously under appropiate boundary conditions. Heat is conducted by the matrix and carried by the vapors. The one-dimensional model includes radially inward and outward flowing vapor within the body, complete depletion of less volatile ices in outer layers, and recondensation of vapor in deeper, coller layers. As a result, we obtain the temperature and abundance distribution in the nucleus and the gas flux into the interior and into the coma for each of the volatiles at various positions in the orbit. The ratio of the gas flux of minor volatiles to that of H2) in the coma varies by several orders of magnitude throughout the orbit. Thus, the relative abundances of species observed in the coma are in most cases not the same as those in the nucleus. Results also indicate that it will be impossible to determine the relative abundances of ices more volatile than water from samples taken a few meters below the surface during a comet rendezvous mission. We made calculations for a wide range of different parameters, such as porosity, pore radius, and thermal conductivity of the matrix. To introduce the model we present typical results for a dust-free comet.

Description: Longitudinal dunes are extremely rare on Mars, but constitute a substantial fraction of terrestrial desert dunes. We report finding isolated examples of longitudinal dunes on Mars and relate their occurence to expected sand transport regimes. Terrestrial longitudinal dunes form in bimodal and multimodal transport regimes. General circulation models and streak data indicate that bimodal and multimodal transport of sand should be very rare on Mars. Thus the dearth of longitudinal dunes on Mars is consistant with their apparent formation conditions on Earth.

Description: We examine the effects of a dusty CO2 atmosphere on the thermal inertia and thermally derived albedo of Mars and we present a new map of thermal inertias. This new map was produced using a coupled surface atmosphere (CSA) model, dust opacities from Viking infrared thermal mapper (IRTM) data, and CO2 columns based on topography. The CSA model thermal inertias are smaller than the 2% model thermal inertias, with the difference largest at large thermal inertia. Although the difference between the thermal inertias obtained with the two models is moderate for much of the region studied, it is largest in regions of either high dust opacity or of topographic lows, including the Viking Lander 1 site and some geologically interesting regions. The CSA model thermally derived albedos do not acurately predict the IRTM measured albedos and are very similar to the thermally derived albedos obtained with models making the 2% assumption.

Description: The Pioneer 10 ultraviolet measurements obtained during the Jupiter encounter in 1973 have been further examined by using improved data handling and analysis techniques. The Pioneer 10 observations of Jupiter and its satellites during the encounter have been carefully reviewed in order to improve our understanding of the morphology of the Io plasma torus and Jupiter's upper atmosphere and to investigate the possible existence of other emission source such as Europa. In addition, the morphology of Io's bimodal torus observed during the Pioneer 10 encounter has been compared with the Voyager observations obtained approximately 6 years after the Pioneer 10 flyby and significant differences in the torus characteristics are found. The Io torus in 1973 was more similar to the 1992 Ulysses observations of a longitudinally asymmetric ring than to the complete ring observed by Voyager. Pioneer 10 observed a significantly dimmer Io torus and Jupiter upper atmosphere in the EUV compared to the Voyager observations. Apart from the torus and Jupiter, Pioneer 10 observed additional emissions which have been attributed to Io itself. Three distinct possibilities have been discussed to explain these additional emissions. The most likely is that Pioneer 10 observed volcanism on Io. There is also evidence of Pioneer 10 observing emissions from Europa. The present analysis clearly shows that the Jovian system in 1973 was significantly different from that observed in 1979.

Description: A new meteorite find from the Nullarbor Plain in Australia was studied using optical, scanning electron microscopy (SEM), and electron microprobe techniques. The meteorite, Nullarbor 018, is an orthodox L6 chondrite that experienced minor-to-moderate alteration of metal during terrestrial weathering (grade A-B to B). During weathering, trolite was preferentially altered, and roughly 20% of the original complement of S in the meteorite was removed. Shock metamorphic effects corresponding to shock stage S4 (or shock facies d) are found, including the presence of some diaplectic feldspar (maskelynite). The meteorite is not obviously paired with other finds from the Nullarbor region, but the possibility that it is paired cannot be excluded.

Description: A large (greater than 4.5 x 7 x 4 mm), igneous-textured clast in the Bovedy (L3) chondrite is notable for its high bulk SiO2 content (is approximately equal to 57.5 wt%). The clast consists of normally zoned orthopyroxene (83.8 vol%), tridymite (6.2%), an intergrowth of feldspar (5.8%) and sodic glass (3.1%), pigeonite (1.0%), and small amounts of chromite (0.2%), augite, and Fe,Ni-metal; it is best described as a silica-rich orthopyroxenite. The oxygen-isotopic composition of the clast is similar, but not identical, to Bovedy and other ordinary chondrites. The clast has a superchondritic Si/Mg ratio, but has Mg/(Mg + Fe) and Fe/Mn ratios that are similar to ordinary chondrite silicate. The closest chemical analogues to the clast are radial-pyroxene chondrules, diogenites, pyroxene-silica objects in ordinary chondrites, and silicates in the IIE iron meteorite Weekeroo Station. The clast crystallized from a siliceous melt that cooled fast enough to prevent complete attainment of equilibrium but slow enough to allow nearly complete crystallization. The texture, form, size and composition of the clast suggestion that it is an igneous differentiate from an asteroid or planetesimal that formed in the vicinity of ordinary chondrites. The melt probably cooled in the near-surface region of the parent object. It appears that in the source region of the clast, metallic and silicate partial melt were largely-to-completely lost during a relatively low degree of melting, and that during a higher degree of melting, olivine and low-Ca pyroxene separated from the remaining liquid, which ultimately solidified to form the clast. While these fractionation steps could not have all occurred at the same temperature, they could have been accomplished in a single melting episode, possibly as a result of heating by radionuclides or by electromagnetic induction. Fractionated magmas can also account for other Si-rich objects in chondrites.

Description: We describe a previously unreported meteorite found in Axtell, Texas, in 1943. Based on the mineralogical composition and texture of its matrix and the sizes and abundance of chondrules, we classify it as a CV3 carbonaceous chondrite. The dominant opaque phase in the chondrules is magnetite, and that in refractory inclusions is Ni-rich metal (awaruite). Axtell, therefore, belongs to the oxidized subgroup of CV3 chondrites, although unlike Allende it escaped strong sulfidation. The meteorite bears a strong textural resemblance to Allende, and its chondrule population and matrix appear to be quite similar to those of Allende, but its refractory inclusions, thermoluminescence properties, and cosmogenic Co-60 abundances are not. Our data are consistent with a terrestrial age for Axtell of approximately 100 years and a metamorphic grade slightly lower than that of Allende.

Description: The 3.0-micrometers water of hydration absorption feature observed in the IR photometry of many low-albedo and some medium-albedo asteroids strongly correlates with the 0.7-micrometers Fe(+2) to Fe(+3) oxidized iron absorption feature observed in narrowband spectrophotometry of these asteroids. Using this relationship, an empirical algorithm for predicting the presence of water of hydration in the surface material of a Solar System body using photometry obtained through the Eight-Color Asteroid Survey nu (0.550 micrometers), w (0.701 micrometers), and x (0.853 micrometers) filters was developed and applied to the ECAS photometry of asteroids and outer planet satellites. The percentage of objects in low-albedo, outer main-belt asteroid classes that test positively for water of hydration increases from P to B to C to G class and correlates linearly with the increasing mean albedos of those objects testing positively. The medium-albedo M-class asteroids do not test positively in large number using this algorithm. Aqueously altered asteroids dominate the Solar System population between heliocentric distances of 2.6 to 3.5 AU, bracketing the Solar System region where the aqueous alteration mechanism operated most strongly. One jovian satellite, J VI Himalia, and one saturnian satellite. Phoebe, tested positively for water of hydration, supporting the hypothesis that these may be captured C-class asteroids from a postaccretional dispersion. The proposed testing technique could be applied to an Earth-based survey of asteroids or a space-probe study of an asteroid's surface characteristic in order to identify a potential water source.

Description: The infrared transmission spectra and photochemical behavior of various organic compounds isolated in solid N2 ices, appropriate for applications to Triton ad Pluto, are presented. It is shown that excess absorption in the surface spectra of Triton and Pluto, i.e., absorption not explained by present models incorporating molecules already identified on these bodies (N2, CH4, CO, and CO2), that starts near 4450/cm (2.25 microns) and extends to lower frequencies, may be due to alkanes (C(n)H(2n+2)) and related molecules frozen in the nitrogen. Branched and linear alkanes may be responsible. Experiments in which the photochemstry of N2: CH4 and N2: CH4: CO ices was explored demonsrtrate that the surface ices of Triton and Pluto may contain a wide variety of additional species containing H, C, O, and N. Of these, the reactive molecule diazomethane, CH2N2, is particularly important since it may be largely responsible for the synthesis of larger alkanes from CH4 and other small alkanes. Diazomethane would also be expected to drive chemical reactions involving organics in the surface ices of Triton and Pluto toward saturation, i.e., to reduce multiple CC bonds. The positions and intrinsic strengths (A values) of many of the infrared absorption bands of N2 matrix-isolated molecules of relevance to Triton and Pluto have also been determined. These can be used to aid in their search and to place constraints on their abundances.

Description: This report presents the most recent spherical harmonic topography model of Venus developed at Jet Propulsion Laboratory. It was produced by a spherical harmonic analysis of the most complete set of Magellan altimetry data, augmented by Pioneer Venus and Venera data. The harmonic coefficients of the topography were computed to degree and order 360. Compared to previous topography models, this one has the highest correlation with the gravity field of Venus.

Description: The 500-Myr average crater retention age for Venus has raised questions about the present-day level of tectonic activity. In this study we examine the relationship between the gravity and topography of four large volcanic swells, Beta, Atla, Bell, and Western Eistla Regiones, for clues about their stage evolution. The Magellan line-of-sight gravity data are inverted using a point mass model of the anomalous mass to solve for the local vertical gravity field. Spectral admittance calculated from both the local gravity inversions and a spherical harmonic model is compared to three models of compensation: local compensation, a 'flexural' model with local and regional compensation of surface and subsurface loads, and a 'hotspot' model of compensation that includes top loading by volcanoes and subsurface loading due to a deep, low density mass anomaly. The coherence is also calculated in each region, but yields an elastic thickness estimate only at Bell Regio. In all models, the long wavelengths are compensated locally. Our results may indicate a relatively old, possibly inactive plume.

Description: Recent mapping studies west of Elysium Mons, Mars, have pinpointed subice features that suggest the existence of a frozen paleolake in Utopia Planitia as recently as 1.8 billion years ago. The subice features are interpreted to be hyaloclastic ridges and hills, table moutains, associated joekulhalaup deposits, and fluvial channels. Photoclinometric studies of these features and of a basal scarp around the northwest flank of Elysium Mons interpreted to have been an ice-sheet boundary indicate that the maximum thickness of ice within the basin may have been about 180 m. This thickness of ice during a relatively late stage of Martian geologic history would have important implications concerning the atmospheric, the climatic, and possibly the exobiologic history of the planet.

Description: The ion and electron momentum equations, along with Ampere's law, are solved for the ion and electron drift velocities and the electric field in the subsolar Venus ionosphere, assuming a partially ionized gas and a single ion species having the ion mean mass. All collision terms among the ions, electrons and neutral particles are retained in the equations. A general expression for the evolution of the magnetic field is derived and compared with earlier expressions. Subsolar region data in the altitude range 150-300 km from the Pioneer Venus Orbiter are used to calculate altitude profiles of the components of the current due to the electric field, gradients of pressure, and gravity. Altitude profiles of the ion and electron velocities as well as the electric field, electrodynamic heating, and the energy density are determined. Only orbits having a complete set of measured plasma temperatures and densities, neutral densities, and magnetic field were considered for analysis; the results are shown only for orbit 202. The vertical velocity at altitudes above 220 km is upgoing for orbit 202. This result is consistent with observations of molecular ions at high altitudes and of plasma flow to the nightside, both of which require upward velocity of ions from the dayside ionosphere. Above about 230 km the momentum equations are extremely sensitive to the altitude profiles of density, temperature, and magnetic field.

Description: We report here analyses of olivines and pyroxenes, and petrofabrics of 27 chondritic interplanetary dust particles (IDPs), comparing those from anhydrous and hydrous types. Approximately 40% of the hydrous particles contain diopside, a probable indicator of parent body thermal metamorphism, while this mineral is rarely present in the anhydrous particles. Based on this evidence, we find that hydrous and anhydrous IDPs are, in general, not directly related, and we conclude that olivine and pyroxene major-element compositions can be used to help discriminate between IDPs that are (1) predominantly nebular condensates, and lately resided in anhydrous or icy (no liquids) primitive parent bodies, and (2) those originating from more geochemically active parent bodies (probably hydrous and anhydrous asteroids).

Description: This paper reports on a magnetic field phenomenon, hereafter referred to as null fields, which were discovered during the inbound pass of the recent flyby of Jupiter by the Ulysses spacecraft. These null fields which were observed in the outer dayside magnetosphere are characterised by brief but sharp decreases of the field magnitude to values less than 1 nT. The nulls are distinguished from the current sheet signatures characteristic of the middle magnetosphere by the fact that the field does not reverse across the event. A field configuration is suggested that accounts for the observed features of the events.

Description: The analysis of balloon envelopes by the finite element (FE) method is plagued by convergence problems. A pratical FE analysis approach is based on the fact that in thin shells with non-zero Gaussian curvature the membrane solution component is essentially decoupled from the bending solution component. A proxy-problem is solved by using a small artificial bending stiffness that assures convergence without significantly affecting the membrane solution component. This approach has been previously validated on slightly overpressurized balloon envelopes. Extensions of this approach to more difficult problems in the structural analysis of balloon envelopes are presented. The convergence forcing modelling measures are discussed. Implications of the findings of the analysis results to future balloon designs are also discussed.

Description: Geochemical profiles of surface units, impact, and volcanic features are studied in detail to determine the underlying structure in an area of extensive mare/highland interface, Sinus Amoris. This study region includes and surrounds the northeastern embayment of Mare Tranquillitatis. The concentrations of two major rock-forming elements (Mg and Al), which were derived from the Apollo 15 orbital geochemical measurements, were used in this study. Mapped units and deposits associated with craters in the northwestern part of the region tend to have correlated low Mg and Al concentrations, indicating the presence of Potassium (K)-Rare Earth Elements (REE)-Phosphorus (P) (KREEP)-enriched basalt. Found along the northeastern rim of Tranquillitatis were areas with correlated high Mg and Al concentration, indicating the presence of troctolite. Distinctive west/east and north/south trends were observed in the concentrations of Mg and Al, and, by implication, in the distribution of major rock components on the surface. Evidence for a systematic geochemical transition in highland or basin-forming units may be observed here in the form of distinctive differences in chemistry in otherwise similar units in the western and eastern portions of the study region.

Description: The aerobraking orbital activities of Magelland during the gravity mapping of Venus are discussed. The goal of aerobraking was to circularize Magellan's orbit. By aerobraking the spacecraft into a nearly circula orbit, the Magellan team was able to provide scientists with a different data set to deepen their understanding of what is going on beneath Venus' surface. Before undertaking its gravity-mapping mission, Magellan completed three cycles of radar mapping. This repeated coverage allowed the spacecraft to see some of Venus' geologic features from different viewing angles. Various aspects of the mission are discussed, and maps of Venus are presented.

Description: Spatial correlation among densely packed particles can substantially change their single-scattering properties, thus making questionable the applicability of the independent scattering approximation in calculations of light scattering by planetary regoliths. The same problem arises in geophysics in light scattering computations for snow, frosts, and bare soil. In this paper, we use a dense-medium light-scattering theory based on the introduction of the static structure factor to calculate asymmetry parameters of the phase function for densely packed particles with real refractive indices 1.31 and 1.66, approximating water ice and soil particles, respectively, and imaginary refractive indices 0, 0.01, and 0.3. For sparsely distributed, independently scattering grains, the calculated asymmetry parameters are always positive and always larger than those for densely packed particles. For densely packed grains, the asymmetry parameters may be negative but only for radius-to-wavelength ratios from about 0.1 to about 0.4. With decreasing particle size, the calculated asymmetry parameters tend to zero independently of the compaction state. In the geometrical optics regime, the asymmetry parameters for densely packed scatterers are positive and very close to those for independently scattering grains. These results may have important implications for remote sensing of the Earth and solid planetary surfaces. In particular, it is demonstrated that negative asymmetry parameters derived with some approximate multiple-scattering theories may be physically irrelevant and can be the result of using an inaccurate bidirectional reflection function combined with the ill-conditionally of the inverse scattering problem.

Description: The first comprehensive discussion of the south seasonal polar cap spectra obtained by the Mariner 7 infrared spectrometer in the short-wavelength region (2-4 microns) is presented. The infrared spectra is correlated with images acquired by the wide-angle camera. Significant spectral variation is noted in the cap interior and regions of varying water frost abundance, CO2 ice/frost cover, and CO2-ice path length can be distinguished. Many of these spectral variations correlate with heterogeneity noted in the camera images, but certain significant infrared spectral variations are not discernible in the visible. Simple reflectance models are used to classify the observed spectral variations into four regions. Region I is at the cap edge, where there is enhanced absorption beyond 3 microns inferred to be caused by an increased abundance of water frost. The increase in water abundance over that in the interior is on the level of a few parts per thousand or less. Region II is the typical cap interior characterized by spectral features of CO2 ice at grain sizes of several millimeters to centimeters. These spectra also indicate the presence of water frost at the parts per thousand level. A third, unusual region (III), is defined by three spectra in which weak CO2 absorption features are as much as twice as strong as in the average cap spectra and are assumed to be caused by an increased path length in the CO2. Such large paths are inconsistent with the high reflectance in the visible and at 2.2 microns and suggest layered structures or deposition conditions that are not accounted for in current reflectance models. The final region (IV) is an area of thinning frost coverage or transparent ice well in the interior of the seasonal cap. These spectra are a combination of CO2 and ground signatures.

Description: We report new measurements of the sodium emission intensity seen in a line of sight just above the surface of the Moon. These data show a strong dependence on lunar phase. The emission intensity decreases from a maximum around first quarter (phase angle 90 deg) to very small values near full Moon (phase angle 0 deg). This suggests that the rate of sodium vapor production from the lunar surface is largest at the subsolar point and becomes small near the terminator. However, the sodium emission near full Moon falls below that which would be expected for solar photon-driven processes. Since the solar wind flux decreases substantially when the Moon enters the Earth's magnetotail near full Moon, while the global solar photon flux is undiminished, we suggest that solar wind sputtering is the dominant process for sodium production.

Description: Infrared diffuse reflectance spectra (2.53-25 microns) of some carbonaceous (C) chondrites were measured. The integrated intensity of the absorption bands near 3 microns caused by hydrous minerals were compared with the modal content of hydrous minerals for the meteorites. The CM and CI chondrites show larger values of the intergated intensity than those of the unique C chondrites Y82162, Y86720 and B7904, suggesting that the amount of hydrous minerals in the CM and CI chondrites is larger, which supports the contention that hydrous minerals were dehydrated by thermal metamorphism in the unique chondrites. Orgueil (CI) has the largest value of the integrated intensity among the C chondrites we measured and shows a sharp absorption band at 3685/cm (2.71 microns) that is not seen in the spectra of the CM chondrites. There is an excellent correlation between the observed hydrogen content in C chondrites and the integrated intensity. The CM chondrites show a wide variation in the strength of absorption bands at 1470/cm (6.8 microns), despite the similarity in absorption features near 3 micron for all CM chondites. The 1470/cm band could be due to the presence of some hydrocarbons but may also be a result of terrestrial alteration processes.

Description: Oxygen production from a lunar rock has been experimentally demonstrated for the first time. A 10 g sample of high-Ti basalt 70035 was reduced with hydrogen in seven experiments at temperatures of 900-1050 C and pressures of 14.7-150 psia. In all experiments, water evolution began almost immediately and was essentially complete in tens of minutes. Oxygen yields ranged from 2.93 to 4.61% of the starting sample weight, and showed weak dependence on temperature and pressure. Analysis of the solid samples demonstrated total reduction of Fe(2+) in ilmenite and small degrees of reduction in olivine and pyroxene. Ti O2 was also partially reduced to one or more suboxides. Data from these experiments provide a basis for predicting the yield of oxygen from lunar basalt as well as new constraints on natural reduction in the lunar regolith.

Description: The ORBSIM program was developed for the accurate extraction of geophysical model parameters from Doppler radio tracking data acquired from orbiting planetary spacecraft. The model of the proposed planetary structure is used in a numerical integration of the spacecraft along simulated trajectories around the primary body. Using line of sight (LOS) Doppler residuals, ORBSIM applies fast and efficient modelling and optimization procedures which avoid the traditional complex dynamic reduction of data. ORBSIM produces quantitative geophysical results such as size, depth, and mass. ORBSIM has been used extensively to investigate topographic features on the Moon, Mars, and Venus. The program has proven particulary suitable for modelling gravitational anomalies and mascons. The basic observable for spacecraft-based gravity data is the Doppler frequency shift of a transponded radio signal. The time derivative of this signal carries information regarding the gravity field acting on the spacecraft in the LOS direction (the LOS direction being the path between the spacecraft and the receiving station, either Earth or another satellite). There are many dynamic factors taken into account: earth rotation, solar radiation, acceleration from planetary bodies, tracking station time and location adjustments, etc. The actual trajectories of the spacecraft are simulated using least squares fitted to conic motion. The theoretical Doppler readings from the simulated orbits are compared to actual Doppler observations and another least squares adjustment is made. ORBSIM has three modes of operation: trajectory simulation, optimization, and gravity modelling. In all cases, an initial gravity model of curved and/or flat disks, harmonics, and/or a force table are required input. ORBSIM is written in FORTRAN 77 for batch execution and has been implemented on a DEC VAX 11/780 computer operating under VMS. This program was released in 1985.

Description: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Description: The Composite Plate Buckling Analysis Program (COMPPAP) was written to help engineers determine buckling loads of orthotropic (or isotropic) irregularly shaped plates without requiring hand calculations from design curves or extensive finite element modeling. COMPPAP is a one element finite element program that utilizes high-order displacement functions. The high order of the displacement functions enables the user to produce results more accurate than traditional h-finite elements. This program uses these high-order displacement functions to perform a plane stress analysis of a general plate followed by a buckling calculation based on the stresses found in the plane stress solution. The current version assumes a flat plate (constant thickness) subject to a constant edge load (normal or shear) on one or more edges. COMPPAP uses the power method to find the eigenvalues of the buckling problem. The power method provides an efficient solution when only one eigenvalue is desired. Once the eigenvalue is found, the eigenvector, which corresponds to the plate buckling mode shape, results as a by-product. A positive feature of the power method is that the dominant eigenvalue is the first found, which is this case is the plate buckling load. The reported eigenvalue expresses a load factor to induce plate buckling. COMPPAP is written in ANSI FORTRAN 77. Two machine versions are available from COSMIC: a PC version (MSC-22428), which is for IBM PC 386 series and higher computers and compatibles running MS-DOS; and a UNIX version (MSC-22286). The distribution medium for both machine versions includes source code for both single and double precision versions of COMPPAP. The PC version includes source code which has been optimized for implementation within DOS memory constraints as well as sample executables for both the single and double precision versions of COMPPAP. The double precision versions of COMPPAP have been successfully implemented on an IBM PC 386 compatible running MS-DOS, a Sun4 series computer running SunOS, an HP-9000 series computer running HP-UX, and a CRAY X-MP series computer running UNICOS. COMPPAP requires 1Mb of RAM and the BLAS and LINPACK math libraries, which are included on the distribution medium. The COMPPAP documentation provides instructions for using the commercial post-processing package PATRAN for graphical interpretation of COMPPAP output. The UNIX version includes two electronic versions of the documentation: one in LaTex format and one in PostScript format. The standard distribution medium for the PC version (MSC-22428) is a 5.25 inch 1.2Mb MS-DOS format diskette. The standard distribution medium for the UNIX version (MSC-22286) is a .25 inch streaming magnetic tape cartridge (Sun QIC-24) in UNIX tar format. For the UNIX version, alternate distribution media and formats are available upon request. COMPPAP was developed in 1992.

Description: Structural flaws and cracks may grow under fatigue inducing loads and, upon reaching a critical size, cause structural failure to occur. The growth of these flaws and cracks may occur at load levels well below the ultimate load bearing capability of the structure. The Fatigue Crack Growth Computer Program, NASA/FLAGRO, was developed as an aid in predicting the growth of pre-existing flaws and cracks in structural components of space systems. The earlier version of the program, FLAGRO4, was the primary analysis tool used by Rockwell International and the Shuttle subcontractors for fracture control analysis on the Space Shuttle. NASA/FLAGRO is an enhanced version of the program and incorporates state-of-the-art improvements in both fracture mechanics and computer technology. NASA/FLAGRO provides the fracture mechanics analyst with a computerized method of evaluating the "safe crack growth life" capabilities of structural components. NASA/FLAGRO could also be used to evaluate the damage tolerance aspects of a given structural design. The propagation of an existing crack is governed by the stress field in the vicinity of the crack tip. The stress intensity factor is defined in terms of the relationship between the stress field magnitude and the crack size. The propagation of the crack becomes catastrophic when the local stress intensity factor reaches the fracture toughness of the material. NASA/FLAGRO predicts crack growth using a two-dimensional model which predicts growth independently in two directions based on the calculation of stress intensity factors. The analyst can choose to use either a crack growth rate equation or a nonlinear interpolation routine based on tabular data. The growth rate equation is a modified Forman equation which can be converted to a Paris or Walker equation by substituting different values into the exponent. This equation provides accuracy and versatility and can be fit to data using standard least squares methods. Stress-intensity factor numerical values can be computed for making comparisons or checks of solutions. NASA/FLAGRO can check for failure of a part-through crack in the mode of a through crack when net ligament yielding occurs. NASA/FLAGRO has a number of special subroutines and files which provide enhanced capabilities and easy entry of data. These include crack case solutions, cyclic load spectrums, nondestructive examination initial flaw sizes, table interpolation, and material properties. The materials properties files are divided into two types, a user defined file and a fixed file. Data is entered and stored in the user defined file during program execution, while the fixed file contains already coded-in property value data for many different materials. Prompted input from CRT terminals consists of initial crack definition (which can be defined automatically), rate solution type, flaw type and geometry, material properties (if they are not in the built-in tables of material data), load spectrum data (if not included in the loads spectrum file), and design limit stress levels. NASA/FLAGRO output includes an echo of the input with any error or warning messages, the final crack size, whether or not critical crack size has been reached for the specified stress level, and a life history profile of the crack propagation. NASA/FLAGRO is modularly designed to facilitate revisions and operation on minicomputers. The program was implemented on a DEC VAX 11/780 with the VMS operating system. NASA/FLAGRO is written in FORTRAN77 and has a memory requirement of 1.4 MB. The program was developed in 1986.

Description: Delaminations near the outer surface of a laminate are susceptible to local buckling and buckling-induced delamination propagation when the laminate is subjected to transverse impact loading. This results in a loss of stiffness and strength. TRBUCKL is an unique dynamic delamination buckling and delamination propagation analysis capability that can be incorporated into the structural analysis program, NASTRAN. This capability will aid engineers in the design of structures incorporating composite laminates. The capability consists of: (1) a modification of the direct time integration solution sequence which provides a new analysis algorithm that can be used to predict delamination buckling in a laminate subjected to dynamic loading; and (2) a new method of modeling the composite laminate using plate bending elements and multipoint constraints. The capability now exists to predict the time at which the onset of dynamic delamination buckling occurs, the dynamic buckling mode shape, and the dynamic delamination strain energy release rate. A procedure file for NASTRAN, TRBUCKL predicts both impact induced buckling in composite laminates with initial delaminations and the strain energy release rate due to extension of the delamination. In addition, the file is useful in calculating the dynamic delamination strain energy release rate for a composite laminate under impact loading. This procedure simplifies the simulation of progressive crack extension. TRBUCKL has been incorporated into COSMIC NASTRAN. TRBUCKL is a DMAP Alter for NASTRAN. It is intended for use only with the COSMIC NASTRAN Direct Transient Analysis (RF 9) solution sequence. The program is available as a listing only. TRBUCKL was developed in 1987.

Description: Predictions of fatigue crack growth behavior can be made with the Fatigue Crack Growth Structural Analysis (FASTRAN II) computer program. As cyclic loads are applied to a selected crack configuration with an initial crack size, FASTRAN II predicts crack growth as a function of cyclic load history until either a desired crack size is reached or failure occurs. FASTRAN II is based on plasticity-induced crack-closure behavior of cracks in metallic materials and accounts for load-interaction effects, such as retardation and acceleration, under variable-amplitude loading. The closure model is based on the Dugdale model with modifications to allow plastically deformed material to be left along the crack surfaces as the crack grows. Plane stress and plane strain conditions, as well as conditions between these two, can be simulated in FASTRAN II by using a constraint factor on tensile yielding at the crack front to approximately account for three-dimensional stress states. FASTRAN II contains seventeen predefined crack configurations (standard laboratory fatigue crack growth rate specimens and many common crack configurations found in structures); and the user can define one additional crack configuration. The baseline crack growth rate properties (effective stress-intensity factor against crack growth rate) may be given in either equation or tabular form. For three-dimensional crack configurations, such as surface cracks or corner cracks at holes or notches, the fatigue crack growth rate properties may be different in the crack depth and crack length directions. Final failure of the cracked structure can be modelled with fracture toughness properties using either linear-elastic fracture mechanics (brittle materials), a two-parameter fracture criterion (brittle to ductile materials), or plastic collapse (extremely ductile materials). The crack configurations in FASTRAN II can be subjected to either constant-amplitude, variable-amplitude or spectrum loading. The applied loads may be either tensile or compressive. Several standardized aircraft flight-load histories, such as TWIST, Mini-TWIST, FALSTAFF, Inverted FALSTAFF, Felix and Gaussian, are included as options. FASTRAN II also includes two other methods that will help the user input spectrum load histories. The two methods are: (1) a list of stress points, and (2) a flight-by-flight history of stress points. Examples are provided in the user manual. Developed as a research program, FASTRAN II has successfully predicted crack growth in many metallic materials under various aircraft spectrum loading. A computer program DKEFF which is a part of the FASTRAN II package was also developed to analyze crack growth rate data from laboratory specimens to obtain the effective stress-intensity factor against crack growth rate relations used in FASTRAN II. FASTRAN II is written in standard FORTRAN 77. It has been successfully compiled and implemented on Sun4 series computers running SunOS and on IBM PC compatibles running MS-DOS using the Lahey F77L FORTRAN compiler. Sample input and output data are included with the FASTRAN II package. The UNIX version requires 660K of RAM for execution. The standard distribution medium for the UNIX version (LAR-14865) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. The standard distribution medium for the MS-DOS version (LAR-14944) is a 5.25 inch 360K MS-DOS format diskette. The contents of the diskette are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. The program was developed in 1984 and revised in 1992. Sun4 and SunOS are trademarks of Sun Microsystems, Inc. IBM PC is a trademark of International Business Machines Corp. MS-DOS is a trademark of Microsoft, Inc. F77L is a trademark of the Lahey Computer Systems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. PKWARE and PKUNZIP are trademarks of PKWare, Inc.

Description: The Panel Analysis and Sizing Code (PASCO) was developed for the buckling and vibration analysis and sizing of prismatic structures having an arbitrary cross section. PASCO is primarily intended for analyzing and sizing stiffened panels made of laminated orthotropic materials and is of particular value in analyzing and sizing filamentary composite structures. When used in the analysis mode, PASCO calculates laminate stiffnesses, lamina stress and strains (including the effects of temperature and panel bending), buckling loads, vibration frequencies, and overall panel stiffness. When used in the sizing mode, PASCO adjusts sizing variables to provide a low-mass panel design that carries a set of specified loadings without exceeding buckling or material strength allowables and that meets other design requirements such as upper and lower bounds on sizing variables, upper and lower bounds on overall bending, extensional and shear stiffnesses, and lower bounds on vibration frequencies. Although emphasis in PASCO is placed on flat panels having several identical bays, the only restriction on configuration modeling is that the structure is assumed to be prismatic. In addition, it is assumed that loads and temperatures do not vary along the length of a panel. Because of their wide application in aerospace structures, stiffened panels are readily handled by PASCO. The panel cross section may be composed of an arbitrary assemblage of thin, flat, rectangular plate elements that are connected together along their longitudinal edges. Each plate element consists of a balanced symmetric laminate of any number of layers of orthotropic material. Any group of element widths, layer thicknesses, and layer orientation angles may be selected as sizing variables. Substructuring is available to increase the efficiency of the analysis and to simplify the modeling of complicated structures. The Macintosh version of PASCO includes an interactive, graphic preprocessor called MacPASCO. The main objective of MacPASCO is to make the use of PASCO faster, simpler, and less error-prone. By using a graphical user interface (GUI), MacPASCO simplifies the specification of panel geometry and reduces user input errors, thus making the modeling and analysis of panel designs more efficient. The user draws the initial structural geometry on the computer screen, then uses a combination of graphic and text inputs to: refine the structural geometry, specify information required for analysis such as panel load conditions, and define design variables and constraints for minimum-mass optimization. Composite panel design is an ideal application because the graphical user interface can: serve as a visual aid, eliminate the tedious aspects of text-based input, and eliminate many sources of input errors. The current version of MacPASCO does not implement all the modeling features of PASCO, but has been found to be sufficient for many users. Many difficulties common to text-based inputs are avoided because MacPASCO uses a GUI. First, the graphic displays eliminate syntax errors, like misplaced commas and incorrect command names, because there is no text-based syntax. Second, graphic displays allow the user to see the geometry as it is created and immediately detect and correct any errors. Third, MacPASCO's drawing tools have been designed to avoid modeling errors. Fourth, the graphic displays make revisions to existing structural designs much easier and less error-prone by eliminating the need for the user to conceptualize the text input as geometry. The user can work directly with the geometry displayed on the screen. Finally, MacPASCO automatically generates the correct PASCO input from the geometry displayed on the screen. This input file can be used with any machine version of PASCO to actually perform the analysis and sizing and to output results. The DEC VAX version of PASCO is written in FORTRAN IV for batch execution and has been implemented on a DEC VAX series computer. The Macintosh version of PASCO was developed for Macintosh II series computers with at least 2Mb of RAM running MPW Pascal 3.0 and Language Systems FORTRAN 2.0 under the MPW programming environment. It includes MPW compatible makefiles for compiling the source code. The Macintosh version uses input files compatible with versions of PASCO running on different platforms. MacPASCO is written in Macintosh Programmers Workbench 3.0, MPW Pascal 3.0, and MacAPP 2.0. The Pascal source code is included on the distribution diskette. MacAPP is a development library which is not included. MacPASCO requires a Mac Plus, SE/30, or MacII, IIx, IIcx, IIci, or IIfx running System 6.0 or greater. MacPASCO is System 7.0 compatible. A minimum of 2Mb of RAM is required for execution. The Macintosh version of PASCO is distributed on four 3.5 inch 800K Macintosh format diskettes. The DEC VAX version is distributed on a 9-track 1600 BPI magnetic tape. The PASCO program was developed in 1981, adapted to the DEC VAX in 1983 and to the Macintosh in 1991. MacPASCO was released in 1992.