ALBERT

Description: A FORTRAN program for calculating the stresses of n cracks embedded in an isotropic plate is presented. Formulas are given for calculating the stresses of one crack, two cracks, and n cracks in an isotropic plate. Then the program code that accomplishes this is provided.

Description: Recently, extremely interesting high-altitude cloud-ionosphere electrical discharges, like lightning above thunderstorms, have been observed from NASA's space shuttle missions and during airborne and ground-based experiments. To understand these discharges, a new experiment was conceived to simulate a thundercloud in a vacuum chamber using a dielectric in particulate form into which electrodes were inserted to create charge centers analogous to those in an electrified cloud. To represent the ionosphere, a conducting medium (metallic plate) was introduced at the top of the chamber. It was found that for different pressures between approximately 1 and 300 mb, corresponding to various upper atmospheric altitudes, different discharges occurred above the simulated thundercloud, and these bore a remarkable similarity to the observed atmospheric phenomena. At pressures greater than 300 mb, these discharges were rare and only discharges within the simulated thundercloud were observed. Use of a particulate dielectric was critical for the successful simulation of the high-altitude lightning.

Description: The objects of the first, exploratory, stage of the project were listed as: (1) to make a detailed and critical review of the Boundary Element method as already published and with regard to elastic-plastic fracture mechanics, to assess its potential for handling present concepts in two-dimensional and three-dimensional cases. To this was subsequently added the Finite Volume method and certain aspects of the Finite Element method for comparative purposes; (2) to assess the further steps needed to apply the methods so far developed to the general field, covering a practical range of geometries, work hardening materials, and composites: to consider their application under higher temperature conditions; (3) to re-assess the present stage of development of the energy dissipation rate, crack tip opening angle and J-integral models in relation to the possibilities of producing a unified technology with the previous two items; and (4) to report on the feasibility and promise of this combined approach and, if appropriate, make recommendations for the second stage aimed at developing a generalized crack growth technology for its application to real-life problems.

Description: An experimental program was conducted to study the damaging effects of tensile and compressive prestrains on the fatigue life of nickel-base, Inconel 718 superalloy at room temperature. To establish baseline fatigue behavior, virgin specimens with a solid uniform gage section were fatigued to failure under fully-reversed strain-control. Additional specimens were prestrained to 2 percent, 5 percent, and 10 percent (engineering strains) in the tensile direction and to 2 percent (engineering strain) in the compressive direction under stroke-control, and were subsequently fatigued to failure under fully-reversed strain-control. Experimental results are compared with estimates of remaining fatigue lives (after prestraining) using three life prediction approaches: (1) the Linear Damage Rule; (2) the Linear Strain and Life Fraction Rule; and (3) the nonlinear Damage Curve Approach. The Smith-Watson-Topper parameter was used to estimate fatigue lives in the presence of mean stresses. Among the cumulative damage rules investigated, best remaining fatigue life predictions were obtained with the nonlinear Damage Curve Approach.

Description: The latest version of the Global Reference Atmospheric Model (GRAM-95) is presented and discussed. GRAM-95 uses the new Global Upper Air Climatic Atlas (GUACA) CD-ROM data set, for 0- to 27-km altitudes. As with earlier versions, GRAM-95 provides complete geographical and altitude coverage for each month of the year. Individual years 1985 to 1991 and a period-of-record (1980 to 1991) can be simulated for the GUACA height range. GRAM-95 uses a specially developed data set, based on Middle Atmosphere Program (MAP) data, for the 20- to 120-km height range, and the NASA Marshall Engineering Thermosphere (MET) model for heights above 90 km. Fairing techniques assure a smooth transition in the overlap height ranges (20 to 27 km and 90 to 120 km). In addition to the traditional GRAM variables of pressure, density, temperature and wind components, GRAM-95 now includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He, and H). A new, variable-scale perturbation model provides both large-scale and small-scale deviations from mean values for the thermodynamic variables and horizontal and vertical wind components. The perturbation model includes new features that simulate intermittency (patchiness) in turbulence and small-scale perturbation fields. The density perturbations and density gradients (density shears) computed by the new model compare favorably in their statistical characteristics with observed density perturbations and density shears from 32 space shuttle reentry profiles. GRAM-95 provides considerable improvement in wind estimates from the new GUACA data set, compared to winds calculated from the geostrophic wind relations previously used in the 0- to 25-km height range. The GRAM-95 code has been put into a more modular form, easier to incorporate as subroutines in other programs (e.g., trajectory codes). A complete user's guide for running the program, plus sample input and output, is provided.

Description: Measurements by single-frequency satellite altimeter (Geosat, ERS-1) require a ionospheric correction to account for the signal time delay in the ionosphere. We propose using the International Reference Ionosphere (IRI) for the determination of this time delay. To investigate the effectiveness of an IRI correction, we have compared the IRI values with ionospheric corrections deduce from measurements by the dual-frequency Topex altimeter. By measuring at two frequencies, the Topex instrument can record (and thus eliminate) the ionospheric influence. We find that IRI agrees with the Topex data much better than the model that is currently used in Geosat data analysis. In particular the earlier model does not represent the equatorial double-peak (equator anomaly) clearly seen in the Topex data. Overall, the use of IRI results in a 30% improvement (over the older model) in the accuracy of ionospheric corrections computed for the first year of the Topex mission.

Description: Papers from the conference are presented and cover the following topics: ion and electron beams; ionospheric modification; spacecraft interactions; chemical releases; and plasma waves. Auroras and plasma emissions are reported from electron beam injection experiments on the EXCEDE 3 rocket and APEX satellite respectively. The important parameters affecting the charging of spacecraft during the operation of electron guns is covered. The Active Magnetospheric Particle Acceleration Satellite (AMPAS) mission utilizing dual-payload tethered satellites and both up and downward directed electron beams is proposed to study the magnetosphere. Recent results and associated theories from the Sura, Arecibo and Troms ionospheric heating facitlites are presented. The effects of neutral gases on spacecraft charging are examined in a series of rocket flights. Many results from the Combined Release and Radiation Effects Satellite chemical release experiments are presented.

Description: The main subject of the symposium was 'off-median phenomena'. The title denotes a range of problems that are rarely considered in 'pure science' studies of the ionosphere. The appearance of regular ionospheric variations is well known. Most of these depend on evident solar-geophysical influences like day and night, solar zenith angle, the seasons, geomagnetic control, solar activity, etc. Applicants and theoreticians as well used to work with monthly medians so that the in fact existing and quite important day-to-day variability is systematically overlooked. It is evident that a descriptive model like the International Reference Ionosphere (IRI) would be inadequate if this variability were denied. Interesting contributions from the symposium on 'off-median phenomena' and the IRI are presented.

Description: This serves as a final report entitled Energization and Transport of Ions of Ionospheric Origin in the Terrestrial Magnetosphere. The work has been predominantly focused on ion outflows identified in two data sets: (1) Prognoz 7; and (2) Dynamics Explorer. The study analyzed ion densities, temperatures, and flow velocities in the magnetotail. The work performed under this contract consisted of developing a program to load the raw data, compute the background subtraction of a strong sun pulse, and use the net counts to calculate the low order moments of the distribution function. The study confirms the results of ISEE that the the cusp is a major source of plasmasheet plasma and goes beyond this to discuss the use of ion velocities as a way to examine the motions of the magnetotail.

Description: Design loads are presented for the General Electric MOD-SA wind turbine. The MOD-SA system consists of a 400 ft. diameter, upwind, two-bladed, teetered rotor connected to a 7.3 mW variable-speed generator. Fatigue loads are specified in the form of histograms for the 30 year life of the machine, while limit (or maximum) loads have been derived from transient dynamic analysis at critical operating conditions. Loads prediction was accomplished using state of the art aeroelastic analyses developed at General Electric. Features of the primary predictive tool - the Transient Rotor Analysis Code (TRAC) are described in the paper. Key to the load predictions are the following wind models: (1) yearly mean wind distribution; (2) mean wind variations during operation; (3) number of start/shutdown cycles; (4) spatially large gusts; and (5) spatially small gusts (local turbulence). The methods used to develop statistical distributions from load calculations represent an extension of procedures used in past wind programs and are believed to be a significant contribution to Wind Turbine Generator analysis. Test/theory correlations are presented to demonstrate code load predictive capability and to support the wind models used in the analysis. In addition MOD-5A loads are compared with those of existing machines. The MOD-5A design was performed by the General Electric Company, Advanced Energy Program Department, under Contract DEN3-153 with NASA Lewis Research Center and sponsored by the Department of Energy.

Description: The Compton Observatory commonly encounters fluxes of energetic electrons which have been scattered from the inner radiation belt to the path of the satellite by resonant interactions with VLF waves from powerful man-made transmitters. The present investigation was motivated by the fact that in the Fall of 1993, the Gamma Ray Observatory was boosted from a 650 km altitude circular orbit to a 750 km altitude circular orbit. This was an opportunity, for the first time, to make observations at two different altitudes using the same instrument. We have examined DISCLA data from the Burst & Transient Source Experiment (BATSE) experiment from 1 Sep. 1993 to 29 Jan. 1994. During the period of study we identified 48 instances of the satellite encountering a cloud of energetic electrons which had been scattered by VLF transmitters. We find that boosting the altitude of the circular orbit from 650 km to 750 km increased the intensity of cyclotron resonance scattered electrons by a factor of two. To search for long term changes in the cyclotron resonance precipitation, we have compared the approx. 750 km altitude data from 106 days at the end of 1993 with data at the same altitudes and time of year in 1991. The cyclotron resonance events in 1991 were three times more frequent and 25% of those cases were more intense than any seen in the 1993 data. We attribute this difference to increased level of geomagnetic activity in 1991 near the Solar Maximum.

Description: Mesoscale model simulations provide insight into the complex jet streak adjustments on 11-12 July 1981 that preceded the first of two significant gravity wave events to have been generated over the Rocky Mountains in Montana. Simulations employing a variety of terrain treatments indicate that prior to wave formation, geostrophic adjustment processes modified the structure of the mid-upper tropospheric jet streak by creating secondary jetlets to the southeast of the polar jet streak in proximity to the gravity wave generation region. This simulated restructuring of the mid-upper tropospheric jet streak is the result of a four stage process. During stage 1, the wind adjusts to the mass field as the jet streak exit region propagates into the inflection point between the upstream trough and downstream ridge in the height field. Stage 2 is initiated as the mass field is forced to adjust to the new ageostrophic wind field created during stage 1. Stage 3 is defined by a second geostrophic adjustment process occurring in a similar manner but to the south and east of the adjustment which occurs during stage 1. A low-level mesoscale jetlet is formed during stage 4 in response to the low-level pressure falls that are established during stage 3. The perturbation of this jetlet, caused by orographically-induced adiabatic and diabatic physical processes, is the likely mechanism responsible for the generation of the first and second episode of observed gravity waves. The dynamics responsible for this wave episode are discussed as differential surface sensible heating inducing an orographically-forced mountain-plains solenoid, resulting in the formation of additional mesoscale jetlets and internal gravity waves. Also discussed is how convective latent heating modifies the numerically simulated terrain-induced internal gravity waves, especially their amplitude and phase velocities, which provide better agreement with those wave characteristics observed in nature. Finally, the three-dimensional linear response of a zonally uniform barotropic flow in a vertically unbounded, continuously stratified, Boussinesq atmosphere which is perturbed from geostrophic equilibrium is investigated.

Description: The implementation of an interface element capability within the COMET-AR software system is described. The report is intended for use by both users of currently implemented interface elements and developers of new interface element formulations. Guidance on the use of COMET-AR is given. A glossary is provided as an Appendix to this report for readers unfamiliar with the jargon of COMET-AR. A summary of the currently implemented interface element formulation is presented in Section 7.3 of this report.

Description: The retarding ion mass spectrometer (RIMS) experiment onboard the Dynamics Explorer 1 (DE 1) satellite was designed to perform energy and mass-per-charge analysis on low-energy ions (less than 50 eV) with mass/charge ratios ranging from 1 to 40 amu/Z. the DE 1 satellite, carrying the RIMS experiment, was launched into an elliptical polar orbit on August 3, 1981. The approximately 7.5 hour orbit has perigee of 675 km altitude and apogee of 24,875 km altitude. This document and those that follow in this series, contain summary RIMS data spectrograms for each orbit for which RIMS data are available. The RIMS instrument began returning science data on day 280 of 1981 and continued to return usable data until the end of the DE mission in March 1991. It should be noted that studies of the RIMS data set should be conducted only with a thorough awareness of the material described in the introduction section presented here, or in collaboration with a scientist familiar with RIMS data analysis.

Description: Impact dampers belong to the category of passive vibration devices used to attenuate the vibration of discrete and continuous systems. An impact damper generally consists of a mass which is allowed to travel freely between two defined stops. Under the right conditions, the vibration of the structure to which the impact damper is attached will cause the mass of the impact damper to strike the structure. Previous analytical and experimental research work on the effect of impact dampers in attenuating the vibration of discrete and continuous systems have demonstrated their effectiveness. It has been shown in this study that impact dampers can increase the intrinsic damping of a lightly-damped flexible structure. The test structure consists of a slender flexible beam supported by a pin-type support at one end and supported by a linear helical flexible spring at another location. Sinusoidal excitation spanning the first three natural frequencies was applied in the horizontal plane. The orientation of the excitation and the test structure in the horizontal plane minimizes the effect of gravity on the behavior of the test structure. The excitation was applied using a linear sine sweep technique. The span of the test structure, the mass of the impact damper, the distance of travel, and the location of the impact damper along the span of the test structure were varied. The damping ratio are estimated for sixty test configurations. The results show that the impact damper significantly increases the damping ratio of the test structure. Statistical analysis of the results using the method of multiple linear regression indicates that a reasonable fit has been accomplished. It is concluded that additional experimental analysis of flexible structures in microgravity environment is needed in order to achieve a better understanding of the behavior of impact damper under conditions of microgravity. Numerical solution of the behavior of flexible structures equipped with impact dampers is also needed to predict stresses and deformations under operating conditions of microgravity in space applications.

Description: The FPCAS3D computer code has been developed for aeroelastic stability analysis of bladed disks such as those in fans, compressors, turbines, propellers, or propfans. The aerodynamic analysis used in this code is based on the unsteady three-dimensional full potential equation which is solved for a blade row. The structural analysis is based on a finite-element model for each blade. Detailed explanations of the aerodynamic analysis, the numerical algorithms, and the aeroelastic analysis are not given in this report. This guide can be used to assist in the preparation of the input data required by the FPCAS3D code. A complete description of the input data is provided in this report. In addition, six examples, including inputs and outputs, are provided.

Description: A probabilistic design method is described and demonstrated using a smart composite wing. Probabilistic structural design incorporates naturally occurring uncertainties including those in constituent (fiber/matrix) material properties, fabrication variables, structure geometry and control-related parameters. Probabilistic sensitivity factors are computed to identify those parameters that have a great influence on a specific structural reliability. Two performance criteria are used to demonstrate this design methodology. The first criterion requires that the actuated angle at the wing tip be bounded by upper and lower limits at a specified reliability. The second criterion requires that the probability of ply damage due to random impact load be smaller than an assigned value. When the relationship between reliability improvement and the sensitivity factors is assessed, the results show that a reduction in the scatter of the random variable with the largest sensitivity factor (absolute value) provides the lowest failure probability. An increase in the mean of the random variable with a negative sensitivity factor will reduce the failure probability. Therefore, the design can be improved by controlling or selecting distribution parameters associated with random variables. This can be implemented during the manufacturing process to obtain maximum benefit with minimum alterations.

Description: If current plans are realized, within the next few years, an extraordinary set of coordinated research efforts focusing on energy flows in the Arctic will be implemented. All are motivated by the prospect of global climate change. SHEBA (Surface Energy Budget of the Arctic Ocean), led by the National Science Foundation (NSF) and the Office of Naval Research (ONR), involves instrumenting an ice camp in the perennial Arctic ice pack, and taking data for 12-18 months. The ARM (Atmospheric Radiation Measurement) North Slope of Alaska and Adjacent Arctic Ocean (NSA/AAO) Cloud and Radiation Testbed (CART) focuses on atmospheric radiative transport, especially in the presence of clouds. The NSA/AAO CART involves instrumenting a sizeable area on the North Slope of Alaska and adjacent waters in the vicinity of Barrow, and acquiring data over a period of about 10 years. FIRE (First ISCCP (International Satellite Cloud Climatology Program) Regional Experiment) Phase 3 is a program led by the National Aeronautics and Space Administration (NASA) which focuses on Arctic clouds, and which is coordinated with SHEBA and ARM. FIRE has historically emphasized data from airborne and satellite platforms. All three program anticipate initiating Arctic data acquisition during spring, 1997. In light of his historic opportunity, the authors discuss a strawman atmospheric radiative transfer experimental plan that identifies which features of the radiative transport models they think should be tested, what experimental data are required for each type of test, the platforms and instrumentation necessary to acquire those data, and in general terms, how the experiments could be conducted. Aspects of the plan are applicable to all three programs.

Description: A survey of bulk parameters of analyzable O(+) outward streams in the mid-altitude (3-4.7 R(sub E) geocentric distance) polar cap magnetosphere is obtained from measurements by the Retarding Ion Mass Spectrometer (RIMS) aboard the Dynamics Explorer-1 (DE-1) spacecraft. There is wide scatter in the obtained densities, but they do display discernible trends: the average O(+) density in these streams decreases from about 60 ions/cc at 3.5 R(sub E) to about 1 ion/cc at 4.6 R(sub E). The streaming velocities are somewhat more defined, and their average increases from about 8 km/s at 3.5 R(sub E) to about 12 km/s at 4.6 R(sub E). The densities and bulk velocities are inversely correlated. We have further compared these observational trends with model profiles for the centrifugally-accelerated polar wind as recently described by Horwitz et al. (1994). The large outflow velocities observed can be understood in part as centrifugally-driven by convection with ionospheric electric field magnitudes of the order 50-70 mV/m, perhaps including plasma expansion effects.

Description: The Generalized Method of Cells (GMC), a micromechanics based constitutive model, is implemented into the finite element code MARC using the user subroutine HYPELA. Comparisons in terms of transverse deformation response, micro stress and strain distributions, and required CPU time are presented for GMC and finite element models of fiber/matrix unit cell. GMC is shown to provide comparable predictions of the composite behavior and requires significantly less CPU time as compared to a finite element analysis of the unit cell. Details as to the organization of the HYPELA code are provided with the actual HYPELA code included in the appendix.

Description: Because the 2195 aluminum-lithium material of the super lightweight external tank (SLWT ET) has a lower toughness than the 2219 aluminum used in previous ET's, careful attention must be paid to stress concentrations. This report details the analysis performed on some of the stress concentrations in the orthogrid panels of the liquid hydrogen tank.

Description: Passive isolator, active vibration absorber, and an integrated passive/active (hybrid) control are studied for their effectiveness in reducing structural vibration under seismic excitations. For the passive isolator, a laminated rubber bearing base isolator which has been studied and used extensively by researchers and seismic designers is considered. An active vibration absorber concept, which can provide guaranteed closed-loop stability with minimum knowledge of the controlled system, is used to reduce the passive isolator displacement and to suppress the top floor vibration. A three-story building model is used for the numerical simulation. The performance of an active vibration absorber and a hybrid vibration controller in reducing peak structural responses is compared with the passively isolated structural response and with absence of vibration control systems under the N00W component of El Centro 1940 and N90W component of the Mexico City earthquake excitation records. The results show that the integrated passive/active vibration control system is most effective in suppressing the peak structural acceleration for the El Centro 1940 earthquake when compared with the passive or active vibration absorber alone. The active vibration absorber, however, is the only system that suppresses the peak acceleration of the structure for the Mexico City 1985 earthquake.

Description: A high-temperature optical window is essential to the optical diagnostics of high-temperature combustion rigs. Laser Doppler velocimetry, schlieren photography, light sheet visualization, and laser-induced fluorescence spectroscopy are a few of the tests that require optically clear access to the combustor flow stream. A design was developed for a high-temperature window that could withstand the severe environment of the NASA Lewis 3200 F Lean Premixed Prevaporized (LPP) Flame Tube Test Rig. The development of this design was both time consuming and costly. This report documents the design process and the lessons learned, in an effort to reduce the cost of developing future designs for high-temperature optical windows.

Description: The Space Power Institute (SPI) at Auburn University has conducted experiments on the effects of impact angle on crater morphology and impactor residue retention for hypervelocity impacts. Copper target plates were set at angles of 30 deg, 45 deg, 60 deg, and 75 deg from the particle flight path. For the 30 deg and 45 deg impacts, in the velocity regime greater than 8 km s(exp -1) the resultant craters are almost identical to normal incidence impacts. The only difference found was in the apparent distribution of particle residue within the crater, and further research is needed to verify this. The 60 deg and 75 deg impacts showed marked differences in crater symmetry, crater lip shape, and particle residue distribution in the same velocity regime. Impactor residue shock fractionation effects have been quantified in first-order. It is concluded that a combination of analysis techniques can yield further information on impact velocity, direction, and angle of incidence.

Description: This report describes the analysis component of the Goddard Earth Observing System, Data Assimilation System, Version 1 (GEOS-1 DAS). The general features of the data assimilation system are outlined, followed by a thorough description of the statistical interpolation algorithm, including specification of error covariances and quality control of observations. We conclude with a discussion of the current status of development of the GEOS data assimilation system. The main components of GEOS-1 DAS are an atmospheric general circulation model and an Optimal Interpolation algorithm. The system is cycled using the Incremental Analysis Update (IAU) technique in which analysis increments are introduced as time independent forcing terms in a forecast model integration. The system is capable of producing dynamically balanced states without the explicit use of initialization, as well as a time-continuous representation of non- observables such as precipitation and radiational fluxes. This version of the data assimilation system was used in the five-year reanalysis project completed in April 1994 by Goddard's Data Assimilation Office (DAO) Data from this reanalysis are available from the Goddard Distributed Active Center (DAAC), which is part of NASA's Earth Observing System Data and Information System (EOSDIS). For information on how to obtain these data sets, contact the Goddard DAAC at (301) 286-3209, EMAIL daac@gsfc.nasa.gov.

Description: The Center for Computational Structures Technology (CST) is intended to serve as a focal point for the diverse CST research activities. The CST activities include the use of numerical simulation and artificial intelligence methods in modeling, analysis, sensitivity studies, and optimization of flight-vehicle structures. The Center is located at NASA Langley and is an integral part of the School of Engineering and Applied Science of the University of Virginia. The key elements of the Center are: (1) conducting innovative research on advanced topics of CST; (2) acting as pathfinder by demonstrating to the research community what can be done (high-potential, high-risk research); (3) strong collaboration with NASA scientists and researchers from universities and other government laboratories; and (4) rapid dissemination of CST to industry, through integration of industrial personnel into the ongoing research efforts.

Description: Consensus on the timing and mapping of substorm features has permitted a synthesis of substorm models. Within the synthesis model the mechanism for onset of substorm expansion is still unknown. Possible mechanisms are: growth of an ion tearing mode, current disruption by a cross-field current instability, and magnetosphere-ionosphere coupling. While the synthesis model is consistent with overall substorm morphology, including near-Earth onset, none of the onset theories, taken individually, appear to account for substorm expansion onset. A grand synthesis with unification of the underlying onset theories appears necessary.

Description: This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present first part of the report, we focus on the theoretical developments, and discussions of the results of numerical-performance studies using the integration schemes for GVIPS and NAV models.

Description: In this report, we address the intercomparison of precipitation (P), evaporation (E), and surface hydrologic forcing (P-E) for 23 Atmospheric Model Intercomparison Project (AMIP) general circulation models (GCM's) including relevant observations, over a variety of spatial and temporal scales. The intercomparison includes global and hemispheric means, latitudinal profiles, selected area means for the tropics and extratropics, ocean and land, respectively. In addition, we have computed anomaly pattern correlations among models and observations for different seasons, harmonic analysis for annual and semiannual cycles, and rain-rate frequency distribution. We also compare the joint influence of temperature and precipitation on local climate using the Koeppen climate classification scheme.

Description: This two-part report is concerned with the development of a general framework for the implicit time-stepping integrators for the flow and evolution equations in generalized viscoplastic models. The primary goal is to present a complete theoretical formulation, and to address in detail the algorithmic and numerical analysis aspects involved in its finite element implementation, as well as to critically assess the numerical performance of the developed schemes in a comprehensive set of test cases. On the theoretical side, the general framework is developed on the basis of the unconditionally-stable, backward-Euler difference scheme as a starting point. Its mathematical structure is of sufficient generality to allow a unified treatment of different classes of viscoplastic models with internal variables. In particular, two specific models of this type, which are representative of the present start-of-art in metal viscoplasticity, are considered in applications reported here; i.e., fully associative (GVIPS) and non-associative (NAV) models. The matrix forms developed for both these models are directly applicable for both initially isotropic and anisotropic materials, in general (three-dimensional) situations as well as subspace applications (i.e., plane stress/strain, axisymmetric, generalized plane stress in shells). On the computational side, issues related to efficiency and robustness are emphasized in developing the (local) interative algorithm. In particular, closed-form expressions for residual vectors and (consistent) material tangent stiffness arrays are given explicitly for both GVIPS and NAV models, with their maximum sizes 'optimized' to depend only on the number of independent stress components (but independent of the number of viscoplastic internal state parameters). Significant robustness of the local iterative solution is provided by complementing the basic Newton-Raphson scheme with a line-search strategy for convergence. In the present second part of the report, we focus on the specific details of the numerical schemes, and associated computer algorithms, for the finite-element implementation of GVIPS and NAV models.

Description: A new method for commanding machines to move with increased dynamic performance was developed. This method is an enhanced version of input shaping, a patented vibration suppression algorithm. This technique intercepts a command input to a system command that moves the mechanical system with increased performance and reduced residual vibration. This document describes many advanced methods for generating highly optimized shaping sequences which are tuned to particular systems. The shaping sequence is important because it determines the trade off between move/settle time of the system and the insensitivity of the input shaping algorithm to variations or uncertainties in the machine which can be controlled. For example, a system with a 5 Hz resonance that takes 1 second to settle can be improved to settle instantaneously using a 0.2 shaping sequence (thus improving settle time by a factor of 5). This system could vary by plus or minus 15% in its natural frequency and still have no apparent vibration. However, the same system shaped with a 0.3 second shaping sequence could tolerate plus or minus 40% or more variation in natural frequency. This document describes how to generate sequences that maximize performance, sequences that maximize insensitivity, and sequences that trade off between the two. Several software tools are documented and included.

Description: Solutions to increasingly larger structural optimization problems are desired. However, computational resources are strained to meet this need. New methods will be required to solve increasingly larger problems. The present approaches to solving large-scale problems involve approximations for the constraints of structural optimization problems and/or decomposition of the problem into multiple subproblems that can be solved in parallel. An area of game theory, equilibrium programming (also known as noncooperative game theory), can be used to unify these existing approaches from a theoretical point of view (considering the existence and optimality of solutions), and be used as a framework for the development of new methods for solving large-scale optimization problems. Equilibrium programming theory is described, and existing design techniques such as fully stressed design and constraint approximations are shown to fit within its framework. Two new structural design formulations are also derived. The first new formulation is another approximation technique which is a general updating scheme for the sensitivity derivatives of design constraints. The second new formulation uses a substructure-based decomposition of the structure for analysis and sensitivity calculations. Significant computational benefits of the new formulations compared with a conventional method are demonstrated.

Description: Because the 2195 aluminum-lithium of the super lightweight external tank (SLWT ET) has a lower toughness than the 2219 aluminum used in previous ET's, careful attention must be paid to stress concentration in the SLWT ET. This report details the initial analysis performed by NASA to determine the material properties required to ensure structural integrity in these critical areas.

Description: This guide describes the input data required for using ECAP2D (Euler Cascade Aeroelastic Program-Two Dimensional). ECAP2D can be used for steady or unsteady aerodynamic and aeroelastic analysis of two dimensional cascades. Euler equations are used to obtain aerodynamic forces. The structural dynamic equations are written for a rigid typical section undergoing pitching (torsion) and plunging (bending) motion. The solution methods include harmonic oscillation method, influence coefficient method, pulse response method, and time integration method. For harmonic oscillation method, example inputs and outputs are provided for pitching motion and plunging motion. For the rest of the methods, input and output for pitching motion only are given.

Description: A recently developed micromechanical theory for the thermoelastic response of functionally graded composites with nonuniform fiber spacing in the through-thickness direction is further extended to enable analysis of material architectures characterized by arbitrarily nonuniform fiber spacing in two directions. In contrast to currently employed micromechanical approaches applied to functionally graded materials, which decouple the local and global effects by assuming the existence of a representative volume element at every point within the composite, the new theory explicitly couples the local and global effects. The analytical development is based on volumetric averaging of the various field quantities, together with imposition of boundary and interfacial conditions in an average sense. Results are presented that illustrate the capability of the derived theory to capture local stress gradients at the free edge of a laminated composite plate due to the application of a uniform temperature change. It is further shown that it is possible to reduce the magnitude of these stress concentrations by a proper management of the microstructure of the composite plies near the free edge. Thus by an appropriate tailoring of the microstructure it is possible to reduce or prevent the likelihood of delamination at free edges of standard composite laminates.

Description: This paper describes the work that is done at the Lehrstuhl fur Raumfahrttechnik (lrt) at the Technische Universitat Munchen to examine particle impacts into germanium surfaces which were flown on board the LDEF satellite. Besides the description of the processing of the samples, a brief overview of the particle launchers at our institute is given together with descriptions of impact morphology of high- and hypervelocity particles into germanium. Since germanium is a brittle, almost glass-like material, the impact morphology may also be interesting for anyone dealing with materials such as optics and solar cells. The main focus of our investigations is to learn about the impacting particle's properties, for example mass, velocity and direction. This is done by examining the morphology, various geometry parameters, crater obliqueness and crater volume.

Description: The investigation of the vibrational disturbances of the Hubble Space Telescope that were discovered soon after deployment in orbit is described in detail. It was found that the disturbances were particularly evident during orbital day-night crossings, and that the magnitudes of the disturbances were considerably larger than the design jitter requirements. This paper describes the process by which the vibrations were characterized and isolated to a particular mechanism. The analysis of the flight data and comparisons with computer simulation results showed that the source of the disturbances was the thermally driven deformation of the solar arrays in conjunction with frictional effects in the array mechanisms. The control system was successfully modified to attenuate the disturbances to tolerable levels pending mechanical and thermal redesign of the solar arrays. The new arrays were installed during the first space telescope servicing mission and, in combination with the enhanced control system algorithm, reduced the disturbances to satisfactory levels.

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: A detailed description of the numerical formulation of Version 2 of the ARIES/GEOS 'dynamical core' is presented. This code is a nearly 'plug-compatible' dynamics for use in atmospheric general circulation models (GCMs). It is a finite difference model on a staggered latitude-longitude C-grid. It uses second-order differences for all terms except the advection of vorticity by the rotation part of the flow, which is done at fourth-order accuracy. This dynamical core is currently being used in the climate (ARIES) and data assimilation (GEOS) GCMs at Goddard.

Description: The retarding ion mass spectrometer (RIMS) experiment onboard the Dynamics Explorer 1 (DE 1) satellite was designed to perform energy and mass-per-charge analysis on low-energy ions (less than 50 eV) with mass/charge ratios ranging from 1 to 40 amu/Z. The DE 1 satellite, carrying the RIMS experiment, was launched into an elliptical polar orbit on August 3, 1981. The approximately 7.5 hour orbit has perigee of 675 km altitude and apogee of 24,875 km altitude. this document and those that following in this series, contains summary RIMS data spectrograms for each orbit for which RIMS data are available. The RIMS instrument began returning science data on day 280 of 1981 and continued to return usable data until the end of the DE mission in March 1991. It should be noted that studies of the RIMS data set should be conducted only with a thorough awareness of the material described in the introduction section presented here, or in collaboration with a scientist familiar with RIMS data analysis.

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: Although many properties of the Earth's magnetosphere have been measured and quantified in the past 30 years since it was discovered, one fundamental measurement (for zeroth order MHD equilibrium) has been made infrequently and with poor spatial coverage - the global electric field. This oversight is due in part to the neglect of theorists. However, there is renewed interest in the convection electric field because it is now realized to be central to many magnetospheric processes, including the global MHD equilibrium, reconnection rates, Region 2 Birkeland currents, magnetosphere ionosphere coupling, ring current and radiation belt transport, substorm injections, and several acceleration mechanisms. Unfortunately the standard experimental methods have not been able to synthesize a global field (excepting the pioneering work of McIlwain's geostationary models) and we are left with an overly simplistic theoretical field, the Volland-Stern electric field model. Single point measurements of the plasmapause were used to infer the appropriate amplitudes of this model, parameterized by K(sub p). Although this result was never intended to be the definitive electric field model, it has gone nearly unchanged for 20 years. The analysis of current data sets requires a great deal more accuracy than can be provided by the Volland-Stern model. The variability of electric field shielding has not been properly addressed although effects of penetrating magnetospheric electric fields has been seen in mid-and low-latitude ionospheric data sets. The growing interest in substorm dynamics also requires a much better assessment of the electric fields responsible for particle injections. Thus we proposed and developed algorithms for extracting electric fields from particle data taken in the Earth's magnetosphere. As a test of the effectiveness of these new techniques, we analyzed data taken by the AMPTE/CCE spacecraft in equatorial orbit from 1984 to 1989.

Description: The Data Assimilation Office (DAO) at Goddard Space Flight Center has produced a multiyear global assimilated data set with version 1 of the Goddard Earth Observing System Data Assimilation System (GEOS-1 DAS). One of the main goals of this project, in addition to benchmarking the GEOS-1 system, was to produce a research quality data set suitable for the study of short-term climate variability. The output, which is global and gridded, includes all prognostic fields and a large number of diagnostic quantities such as precipitation, latent heating, and surface fluxes. Output is provided four times daily with selected quantities available eight times per day. Information about the observations input to the GEOS-1 DAS is provided in terms of maps of spatial coverage, bar graphs of data counts, and tables of all time periods with significant data gaps. The purpose of this document is to serve as a users' guide to NASA's first multiyear assimilated data set and to provide an early look at the quality of the output. Documentation is provided on all the data archives, including sample read programs and methods of data access. Extensive comparisons are made with the corresponding operational European Center for Medium-Range Weather Forecasts analyses, as well as various in situ and satellite observations. This document is also intended to alert users of the data about potential limitations of assimilated data, in general, and the GEOS-1 data, in particular. Results are presented for the period March 1985-February 1990.

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: Evidence for the probable existence of magnetospheric boundary layers was first presented by Hones, et al. (1972), based on VELA satellite plasma observations (no magnetic field measurements were obtained). This magnetotail boundary layer is now known to be the tailward extension of the high-latitude boundary layer or plasma mantle (first uniquely identified using HEOS 2 plasma and field observations by Rosenbauer et al., 1975) and the low-latitude boundary layer (first uniquely identified using IMP 6 plasma and field observations by Eastman et al., 1976). The magnetospheric boundary layer is the region of magnetosheath-like plasma located Earthward of, but generally contiguous with the magnetopause. This boundary layer is typically identified by comparing low-energy (less than 10 keV) ion spectra across the magnetopause. Low-energy electron measurements are also useful for identifying the boundary layer because the shocked solar wind or magnetosheath has a characteristic spectral signature for electrons as well. However, there are magnetopause crossings where low-energy electrons might suggest a depletion layer outside the magnetopause even though the traditional field-rotation signature indicates that this same region is a boundary layer Earthward of the current layer. Our analyses avoided crossings which exhibit such ambiguities. Pristine magnetopause crossings are magnetopause crossings for which the current layer is well defined and for which there is no adjoining magnetospheric boundary layer as defined above. Although most magnetopause models to date apply to such crossings, few comparisons between such theory and observations of pristine magnetopause crossings have been made because most crossings have an associated magnetospheric boundary layer which significantly affects the applicable boundary conditions for the magnetopause current layer. Furthermore, almost no observational studies of magnetopause microstructure have been done even though key theoretical issues have been discussed for over two decades. This is because plasma instruments deployed prior to the ISEE and AMPTE missions did not have the required time resolution and most ISEE investigations to-date have focused on tests of MHD plasma models, especially reconnection. More recently, many phenomenological and theoretical models have been developed to explain the existence and characteristics of the magnetospheric boundary layers with only limited success to date. The cases with no boundary layer treated in this study provide a contrary set of conditions to those observed with a boundary layer. For the measured parameters of such cases, a successful boundary layer model should predict no plasma penetration across the magnetopause. Thus, this research project provides the first direct observational tests of magnetopause models using pristine magnetopause crossings and provides important new results on magnetopause microstructure and associated kinetic processes.

Description: The work of this grant has been predominantly focused on ion outflows from two data sets: Prognoz 7 and Dynamics Explorer. The Prognoz analysis studied ion densities, temperatures, and flow velocities in the magnetotail. The work performed under this contract consisted of developing a program to load the raw data, computing the background subtraction of a strong sun pulse, and using the net count to calculate the low order moments of the distribution function. The study confirms the results of ISEE with regard to the supply of plasma from the cusp as a major source of plasmasheet plasma and goes beyond this to discuss the use of ion velocities as a way to examine the motions of the magnetotail. The abstract of the work to be reported is included as an appendix. The work on the DE/Retarding Ion Mass Spectrometer is separated into two categories: (1) classification of low-energy ion flows from high-latitudes, and (2) studies of the polar wind. Major publications resulting from this work are also included as an appendix to this report. The polar wind is in a category by itself as a result of the thermal escape of hydrogen and helium because of ambipolar diffusion through the heavier, oxygen-dominated topside ionosphere. The analysis of the polar wind reports the flux variability as a function of season, magnetic activity, etc. Much effort has been expended under this grant to complete a follow on study of the thermal structure of the polar wind. Extensive display tools and analysis software have been developed and used in an attempt to carry out this thermal analysis. The present work uses a constrained fit scheme that combines the ion densities and flow velocities derived from Chandler et al. and a spacecraft potential derived from an empirical relation to the total ion density to determine the remaining fit parameter, the ion temperature, via a least squares fit to the RIMS data.

Description: This report covers the time period 1 January 1994 to 31 December 1994. During this reporting period we had our fourth Upper Atmosphere Research Satellite (UARS) correlative balloon flight; the data from this flight have been reduced and submitted to the UARS CDHF. We have spent most of the past year analyzing data from this and past flights. For example, using data from our September 1989 balloon flight we have demonstrated for the first time ever that the rates of production and loss of ozone are in balance in the upper stratosphere. As part of this analysis, we have completed the most detailed study to date of radical partitioning throughout the stratosphere. We have also produced the first measurement of HBr and HOBr mixing ratio profiles over a full diurnal cycle.

Description: Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Description: The 'radial' transport of energy by internal ULF waves, stimulated by dayside magnetospheric boundary oscillations, is analyzed in the framework of one-fluid magnetohydrodynamics. (the term radial is used here to denote the direction orthogonal to geomagnetic flux surfaces.) The model for the inhomogeneous magnetospheric plasma and background magnetic field is axisymmetric and includes radial and parallel variations in the magnetic field, magnetic curvature, plasma density, and low but finite plasma pressure. The radial mode structure of the coupled fast and intermediate MHD waves is determined by numerical solution of the inhomogeneous wave equation; the parallel mode structure is characterized by a Wentzel-Kramer-Brillouin (WKB) approximation. Ionospheric dissipation is modeled by allowing the parallel wave number to be complex. For boudnary oscillations with frequencies in the range from 10 to 48 mHz, and using a dipole model for the background magnetic field, the combined effects of magnetic curvature and finite plasma pressure are shown to (1) enhance the amplitude of field line resonances by as much as a factor of 2 relative to values obtained in a cold plasma or box-model approximation for the dayside magnetosphere; (2) increase the energy flux delivered to a given resonance by a factor of 2-4; and (3) broaden the spectral width of the resonance by a factor of 2-3. The effects are attributed to the existence of an 'Alfven buoyancy oscillation,' which approaches the usual shear mode Alfven wave at resonance, but unlike the shear Alfven mode, it is dispersive at short perpendicular wavelengths. The form of dispersion is analogous to that of an internal atmospheric gravity wave, with the magnetic tension of the curved background field providing the restoring force and allowing radial propagation of the mode. For nominal dayside parameters, the propagation band of the Alfven buoyancy wave occurs between the location of its (field line) resonance and that of the fast mode cutoff that exists at larger radial distances.

Description: A statistical study has been made of the high-latitude impulsive events that were observed during the 1985-1986 South Pole Balloon Campaign. The events were selected by searching for unipolar pulses greater than or equal to 10 nT above background in the vertical component of the magnetic field on the ground and/or pedestal or 'W' shaped horizontal electric field perturbations greater than or equal to 10 mV/m in amplitude and accompanied by perturbations in the vertical electric field at balloon altitude. A main event list comprising 112 events was compiled from the 468 hours of data available. Three aspects of the events were examined: the solar wind conditions prior to the event, local time of observation, and intrinsic properties of the events. The local time distribution was obtained from the 112 entry main event list and was found to be nearly uniform across the dayside, with no midday gap. The event rate found using our low-amplitude selection criteria was 0.7 event/hr, comparable to expectations based on in situ studies of the magnetopause. A total of 42 events were found for which data were available from Interplanetary Monitoring Platform (IMP) 8. Of these events, 12 occurred when the Z(sub GSM) component (B(sub Z)) of the interplanetary magnetic field (IMF) was northward and 30 occurred when B(sub Z) was southward or fluctuating. Only three of the B(sub Z) northward cases and only five of the B(sub Z) southward cases were preceded by pressure pulses greater than 0.4 nPa in amplitude. Ten of the events were studied in detail by means of a model-fitting method discussed elsewhere. This method infers values of several parameters, including the total current flowing in a coaxial or monopole system and a two-dimensional dipole system. The intrinsic properties of the events showed that only approximately 10% of the total current contributed to momentum transfer to the high-latitude ionosphere, that the direction of the motion depended more on local time of observation than IMF B(sub y), and that events were usually several hundred kilometers in size. The observed B(sub z) control found in the 42 event list and the prevalence of coaxial current dominated events are inconsistent with the predictions of the pressure pulse model.

Description: The peak vertical velocities predicted by three realistic, but contrasting, present-day scenarios of Antarctic ice sheet mass balance are found to be of the order of several mm/a. One scenario predicts local uplift rates in excess of 5 mm/a. These rates are small compared to the peak Antarctic vertical velocities of the ICE-3G glacial rebound model, which are in excess of 20 mm/a. If the Holocene Antarctic deglaciation history protrayed in ICE-3G is realistic, and if regional upper mantle viscosity is not an order of magnitude below 10(exp 21) Pa(dot)s, then a vast geographical region in West Antarctica is uplifting at a rate that could be detected by a future Global Positioning System (GPS) campaign. While present-day scenarios predict small vertical crustal velocities, their overall continent-ocean mass exchange is large enough to account for a substantial portion of the observed secular polar motion (omega m(arrow dot)) and time-varying zonal gravity field.

Description: We demonstrate connections between decadal and secular global climatic variations, and historical variations in the volume of the Great Salt Lake. The decadal variations correspond to a low-frequency shifting of storm tracks which influence winter precipitation and explain nearly 18% of the interannual and longer-term variance in the record of monthly volume change. The secular trend accounts for a more modest approximately 1.5% of the variance.

Description: During the period July-August 1991, observations were made of Polar Mesospheric Summer Echoes (PMSE) at 46.9 MHz and 224 MHz by the CUPRI and EISCAT radars, respectively, at two sites in northern Scandinavia. Those observations are compared here with observations of noctilucent clouds, nergetic particle precipitation and magnetic disturbances. The appearance and morphology of PMSE are found to be closely correlated at the two frequencies and the two sites, 200 km apart. No correlation is found between PMSE and noctilucent clouds or magnetic disturbance. No correlation is found between energetic particle precipitation and the appearance of PMSE at 46.9 MHz for the whole time period. At 224 MHz, there is no evidence for a correlation before the beginning of August and only one event suggesting a possible correlation after the beginning of August. A minimum in occurrence for PMSE is found between 16 and 21 UT (17-22 LST) which may be related to an expected minimum in background wind strength in that time interval.

Description: A new analytical radiative transfer model of a leaf canopy is developed that approximates multiple-scattering radiance by a four-stream formulation. The canopy model is coupled to a homogeneous atmospheric model as well as a non-Lambertian lower boundary soil surface. The same four-stream formulation is also used for the calculation of multiple scattering in the atmosphere. Comparisons of radiance derived from the four-stream model with those calculated by an iterative numerical solution of the radiative transfer equation show that the analytic model has a very high accuracy, even with a turbid atmosphere and a very dense canopy in which multiple scattering dominates. Because the coupling of canopy and atmospheric models fully accommodates anisotropic surface reflectance and atmospheric scattering and its effect on directional radiance, the model is especially useful for application to directional radiance and measurements obtained by remote sensing. Retrieval of biophysical parameters using this model is under investigation.

Description: As part of the Cooperative High Altitude Rocket Gun Experiment (CHARGE-2B) rocket mission, an electron beam was injected into the ionosphere with a modulated beam current in an effort to generate very low frequency (VLF) waves. The propagation of the beam-driven VLF waves through the ionosphere is examined here to determine whether it is possible to detect these wave emissions with ground receivers. The paths of the VLF waves from where they were generated near the rocket were followed to the bottom of the ionosphere and the decrease in wave amplitude due to wave-particle resonance and collisional damping was calculated. It was found that due to collisional damping, which for these VLF waves becomes large at altitudes below about 150 km, wave amplitudes were decreased below the background atmospheric noise level. A number of different beam injection events have been examined and in all of these cases studied the waves were sufficiently damped such that detection on the ground would not be possible. This is in agreement with observations on the ground in which no wave emissions were observed during the CHARGE-2B mission. Control parameters that would be more favorable for beam-generated VLF propagation to the ground are discussed for future experiments of this type.

Description: The rapid formation of a new proton radiation belt at L approximately = 2.5 following the March 24, 1991 Storm Sudden Commencement (SSC) observed at the Combined Release and Radiation Effects Satellite (CRRES) satellite is modeled using a relativistic guiding center test particle code. The SSC is modeled by a bipolar electric field and associated compression and relaxation in the magnetic field, superimposed on a dipole magnetic field. The source population consists of both solar and trapped inner zone protons. The simulations show that while both populations contribute to drift echoes in the 20-80 MeV range, primary conditions is from the solar protons. Proton acceleration by the SSC differs from relativistic electron acceleration in that different source populations contribute and nonrelativistic conservation of the first adiabatic invariation leads to greater energization of protons for a given decrease in L. Model drift echoes and flux distribution in L at the time of injection compare well with CRRES observations.

Description: The properties of the 1-micrometer volume extinction coefficient of two geographically different high-altitude cloud systems have been examined for the posteruption period (1985-1990) of the April 1982 El Chichon volcanic event with emphasis on the effect of volcanic aerosols on clouds. These two high-altitude cloud systems are the tropical clouds in the tropopause region observed by the Stratospheric Aerosol and Gas Experiment (SAGE) 2 and the polar stratospheric clouds (PSCs) sighted by the Stratospheric Aerosol Measurement (SAM) 2. The results indicate that volcanic aerosols alter the frequency distributions of these high-altitude clouds in such a manner that the occurrence of clouds having high extinction coefficients (6 x 10(exp -3) - 2 x 10(exp -2)/km) is suppressed, while that of clouds having low extinction coefficients (2 x 10(exp -3) - 6 x 10(exp -2)/km) is enhanced. This influence of the volcanic aerosols appears to be opposite to the increase in the extinction coefficient of optically thick clouds observed by the Earth Radiation Budget Experiment (ERBE) during the initial posteruption period of the June 1991 Pinatubo eruption. A plausible explanation of this difference, based on the Mie theory, is presented. As a consequence of the Mie theory, the effective radius of most, if not all, of the high-altitude clouds, measured by the SAGE series of satellite instruments must be less than about 0.8 micrometers. This mean cloud particle size implied by the satellite extinction-coefficient data at a single wavelength (1 micrometer) is further substantiated by the particle size analysis based on cloud extinction coefficient at two wavelengths (0.525 and 1.02 micrometers) obtained by the SAGE 2 observations. Most of the radiation measured by ERBE is reflected by cloud systems comprised of particles having effective radii much greater than 1 micrometer. A reduction in the effective radius of these clouds due to volcanic aerosols is expected to increase their extinction-coefficient values, opposite the effect observed by SAGE 2 and SAM 2. This work further illustrates the capability of the solar occultation satellite sensor to provide particulate extinction-coefficient measurements important to the study of the aerosol-cloud interactions. It is important to examine the variations of the extinction coefficient of these two high-altitude cloud systems for the posteruption years of the Pinatubo volcanic event for further evidence of the impact of volcanic aerosols on high-altitude clouds.

Description: We use global, multiyear observations of the properties of clouds, the atmosphere, and the surface to calculate global shortwave (SW) and longwave (LW) fluxes at the top of the atmosphere and at the surface at a resolution of 280 km and 3 hours for every third month from April 1985 to January 1989. Our validation studies suggest that the specification of cloud effects is no longer the dominant uncertainty in reconstructing the radiative fluxes at the top of atmosphere and at the surface. Rather cloud property uncertainties are now roughly equal contributors to the flux uncertainty, along with surface and atmospheric properties. The resulting SW and LW flux data sets suggest the following conclusions: (1) The net SW heating of Earth appears predominantly at the surface, whereas the net LW cooling arises predominantly from the atmosphere. The net cooling effect of clouds on top of atmospheric radiation appears primarily at the surface rather than in the atmosphere. (2) Clouds have almost no net effect on the global mean radiation balance of the atmosphere, but they enhance the latitudinal gradient in the LW cooling and reinforce the radiative forcing for the mean atmospheric circulation. Clouds act to mute seasonal contrasts however. (3) Clouds enhance the land-ocean contrasts of the atmospheric cooling, reinforcing the growth of standing eddy motions; but reduce land-ocean contrasts of the surface heating.

Description: The largest uncertainty in upwelling shortwave (SW) fluxes (approximately equal 10-15 W/m(exp 2), regional daily mean) is caused by uncertainties in land surface albedo, whereas the largest uncertainty in downwelling SW at the surface (approximately equal 5-10 W/m(exp 2), regional daily mean) is related to cloud detection errors. The uncertainty of upwelling longwave (LW) fluxes (approximately 10-20 W/m(exp 2), regional daily mean) depends on the accuracy of the surface temperature for the surface LW fluxes and the atmospheric temperature for the top of atmosphere LW fluxes. The dominant source of uncertainty is downwelling LW fluxes at the surface (approximately equal 10-15 W/m(exp 2)) is uncertainty in atmospheric temperature and, secondarily, atmospheric humidity; clouds play little role except in the polar regions. The uncertainties of the individual flux components and the total net fluxes are largest over land (15-20 W/m(exp 2)) because of uncertainties in surface albedo (especially its spectral dependence) and surface temperature and emissivity (including its spectral dependence). Clouds are the most important modulator of the SW fluxes, but over land areas, uncertainties in net SW at the surface depend almost as much on uncertainties in surface albedo. Although atmospheric and surface temperature variations cause larger LW flux variations, the most notable feature of the net LW fluxes is the changing relative importance of clouds and water vapor with latitude. Uncertainty in individual flux values is dominated by sampling effects because of large natrual variations, but uncertainty in monthly mean fluxes is dominated by bias errors in the input quantities.

Description: A 3-km-long gap in the dextral surficial rupture of the 1992 M(sub W) = 7.3 Landers earthquake occurs at the north end of a major fault stepover between the Johnson Valley and Homestead Valley faults. This gap is situated along a segment of the Landers rupture that has been modeled geophysically as having a deficit in average slip at depth. To better evaluate the nature of the slip gap, we document in detail the character and distribution of surficial rupture within it. Along the gap, is a northwest trending thrust fault rupture with an average of less than 1 m of northeast directed reverse-slip and nearly no oblique right slip. We interpret this rupture to be limited to the shallow crust of the northern end of the stepover and to have been the secondary result of dextral shear, rather than a mechanism of rigid-block slip-transfer from the Landers-Kickapoo fault. A zone of en echelon extensional ruptures also occurs along the slip gap, which we interpret as the secondary result of diffuse dextral shear that accommodated less than 0.5 m of west-northwest extension. These secondary ruptures represent a discontinuity in the surficial dextral ruptre of the Landers earthquake, which we propose resulted from the lack of a mature fault connection between the Johnson Valley and Homestead Valley faults. The rupture pattern of the slip gap implies a significant deficit in net surficial slip, which compares favorably with some geophysical models. Aspects of this rupture pattern also suggest a temporal sequence of rupture that compares favorably with geophysical interpretations of the dynamic rupture propagation.

Description: This report covers the activities performed under NAS5-32572. The results of those activities are included in this Final Report. TIMED Science Objectives: (1) To determine the temperature, density, and wind structure of the MLTI (mixed layer thermal inertia), including the seasonal and latitudinal variations; and (2) To determine the relative importance of the various radiative, chemical, electrodynamical, and dynamical sources and sinks of energy for the thermal structure of the MLTI. GUVI Science Goals: (1) Determine the spatial and temporal variations of temperature and constituent densities in the lower thermosphere; and (2) Determine the importance of auroral energy sources and solar EUV (extreme ultraviolet) to the energy balance of the region.

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: 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: This report describes progress on the first year of a research program on the infrared radiation of air plasmas conducted in the High Temperature Gasdynamics Laboratory at Stanford University. This program is intended to investigate the masking of infrared signatures by the air plasma formed behind the bow shock of high velocity missiles. To this date, the radiative emission of air plasmas in the infrared has been the object of few experimental investigations, and although several infrared systems are already modeled in radiation codes such as NEQAIR, measurements are required to validate numerical predictions and indicate whether all transitions of importance are accounted for. The present program is motivated by the fact that 9 excited states (A, B, C, D, B', F, H, and H') of NO radiate in the infrared, especially between 1 and 1.5 microns where at least 9 transitions involving can be observed. Because these IR transitions are relatively well separated from each other, excited NO states concentrations can be easily measured, thus providing essential information on excited-state chemistry for use in optical diagnostics or in electronic excitation model validation. Developing accurate collisional-radiative models for these excited NO states is of importance as the UV-VUV transitions of NO (beta, gamma, epsilon, beta prime, gamma prime) produce a major, if not dominant, fraction of the radiation emitted by air plasmas. During the first year of the program, research has focused on the spectral range 1.0 to 1.5 microns, as detailed in Section 2 of this report. The measurements, conducted in a 50 kW radio-frequency inductively coupled plasma torch operating on air at atmospheric pressure, extend previous shock tube investigations by Wray to a wider spectral range (1.0 to 1.5 microns vs 0.9 to 1.2 microns) and higher temperatures (7600 K in the plasma torch versus 6700 K in the shock-tube). These higher temperatures in the present experiment have made it possible to observe high-lying excited NO states that were previously undetectable. These measurements are currently being extended up to 5 microns, with particular attention paid to the rovibronic bands of ground state NO, molecular continua, CO transitions, and other systems of importance. Publications and presentations resulting from or related to this work are cited in Section 3, and Section 4 lists the personnel who contributed to this report.

Description: Work at SRI involved modeling the exchange of electromagnetic energy between the ionosphere and magnetosphere to help interpret the DE-B Poynting flux observations. To describe the electrical properties of the high-latitude ionosphere, we constructed a numerical model, from the framework provided by the Vector Spherical Harmonic (VSH) model, that determines the ionospheric currents, conductivities, and electric fields including both magnetospheric inputs and neutral wind dynamo effects. This model development grew from the earlier question of whether an electrical energy source in the ionosphere was capable of providing an upward Poynting flux. The model solves the steady-state neutral wind dynamo equations and the Poynting flux equation to provide insight into the electrodynamic role of the neutral winds. The modeling effort to determine the high-latitude energy flux has been able to reproduce many of the large-scale features observed in the Poynting flux measurements made by DE-2. Because the Poynting flux measurement is an integrated result of energy flux into or out of the ionosphere, we investigated the ionospheric properties that may contribute to the observed flux of energy measured by the spacecraft. During steady state the electromagnetic energy flux, or DC Poynting flux, is equal to the Joule heating rate and the mechanical energy transfer rate in the high-latitude ionosphere. Although the Joule heating rate acts as an energy sink, transforming electromagnetic energy into thermal or internal energy of the gas, the mechanical energy transfer rate may be either a sink or source of electromagnetic energy. In the steady state, it is only the mechanical energy transfer rate that can generate electromagnetic energy and result in a DC Poynating flux that is directed out of the ionosphere.

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 objective of this project was to investigate the utility of satellite synthetic aperture radar (SAR) imagery for measurement of geophysical parameters on Alaskan glaciers relevant to their mass balance and dynamics, including: (1) the positions of firn lines (late-summer snow lines); (2) surface velocities on fast-flowing (surging) glaciers, and also on slower steady-flow glaciers; and (3) the positions and changes in the positions of glacier termini. Preliminary studies of topography and glacier surface velocity with SAR interferometry have also been carried out. This project was motivated by the relationships of multi-year to decadal changes in glacier geometry to changing climate, and the probable significant contribution of Alaskan glaciers to rising sea level.

Description: The poleward arc system of a double oval distribution is shown to activate at the end of the optical expansion phase signifying the beginning of substorm recovery. The velocity dispersed ion signature (VDIS) can exist coincident with this discrete aurora developing on the most poleward oval. Although the VDIS is usually associated with ion beams in the plasma sheet boundary layer, it is demonstrated that the ionospheric signature is not beamlike but distributed in pitch angle. At the time when the double oval begins to form, the magnetic field in the magnetotail lobe becomes less flared and can show Pc 5 period oscillations. Similar pulsations also exist in the ionosphere associated with the most poleward oval and with stationary surge formation. Theoretical considerations link this phenomenon with a wave source tailward of x(sub GSE) = -30R(sub E) and fast mode evanescent waves propagating earthward in the tail lobe region. In this case the magnetotail appears to act like a waveguide and the plasma sheet boundary layer as a resonance region. This implies that the coupling of this fast mode waves is with the plasma sheet boundary layer and not with dipolar like field lines. The implications of this for the reconnection model of substorms are discussed.

Description: During the later stages of the auroral substorm the luminosity distribution frequently resembles a double oval, one oval lying poleward of the normal or main UV auroral oval. We interpret the double oval morphology as being due to the plasma sheet boundary layer becoming active in the later stages of the substorm process. If the disturbance engulfs the nightside low-latitude boundary layers, then the double oval configuration extends into the dayside ionospheric region. The main UV oval is associated with the inner portion of the central plasma sheet and can rapidly change its auroral character from being diffuse to discrete. This transition is associated with the substorm process and is fundamental to understanding the near-Earth character of substorm onset. On the other hand, the poleward arc system in the nightside ionosphere occurs adjacent to or near the open-closed field line boundary. This system activates at the end of the optical expansion phase and is a part of the recovery phase configuration in substorms where it occurs. These two source regions for nightside discrete auroral arcs are important in resolving the controversy concerning the mapping of arcs to the magnetosphere. The dayside extension of this double oval configuration is also investigated and shows particle signatures which differ considerably from those on the nightside giving clues to the magnetospheric source regions of the aurora in the two local time sectors. Near-Earth substorm onsets are shown to be coupled to processes occurring much further tailward and indicate the importance of understanding the temporal development of features within the double oval. Using 'variance images,' a new technqiue for the investigation of these dynamics is outlined.

Description: Sources such as atmospheric or buried explosions and shallow earthquakes producing strong vertical ground displacements produce pressure waves that propagate at infrasonic speeds in the atmosphere. At ionospheric altitudes low frequency acoustic waves are coupled to ionispheric gravity waves and induce variations in the ionoispheric electron density. Global Positioning System (GPS) data recorded in Southern California were used to compute ionospheric electron content time series for several days preceding and following the January 17, 1994, M(sub w) = 6.7 Northridge earthquake. An anomalous signal beginning several minutes after the earthquake with time delays that increase with distance from the epicenter was observed. The signal frequency and phase velocity are consistent with results from numerical models of atmospheric-ionospheric acoustic-gravity waves excited by seismic sources as well as previous electromagnetic sounding results. It is believed that these perturbations are caused by the ionospheric response to the strong ground displacement associated with the Northridge earthquake.

Description: Observations of waves at frequencies below approximately 200 Hz obtained near the magnetopause are presented. For one case identified in the ISEE 1 data as a period when steady state reconnection was occurring, there were waves below the lower hybrid frequency with amplitudes up to approximately 7 mV/m. Intense low-frequency waves with amplitudes up to approximately 20 mV/m at the subsolar magnetopause have also been observed by the Geotail electric field instrument. In some cases, large spiky fields were embedded in the waves. The waves observed by ISEE 1 and Geotail were large enough to provide the dissipation required for reconnection to occur.

Description: Pyranometers have been used for many years to measure broadband surface incoming solar irradiance, data that is necessary for surface energy budget, cloud forcing, and satellite validation research. Because such measurements are made at a specific location, it is unclear how representative they may be of a larger area. This study attempts to determine a reasonable spacing between measurement sites for such research by computing the correlation, and standard deviation from perfect correlation, between simultaneous measurements of incoming solar irradiance for a network of surface measurement sites covering a 75 km x 75 km area. Using 1-min data collected from this network of 11 sites during the NASA First ISSCP Radiation Experiment/Surface Radiation Budget (FIRE/SRB) Project temporal averages were calculated. The correlation between any two of these sites was determined by comparing simultaneous measurement averages for the 55 possible combinations of site pairs, along with the distances between them. In an attempt to remove the effect of the diurnal cycle, thus leaving clouds as the primary influence on correlation of the radiation field, model results for a clear day were used to normalize measured irradiances and correlations were again calculated.

Description: We examine the consequences of magnetic reconnection at the high-latitude magnetopause using a three-dimensional global magnetohydrodynamic simulation of the solar wind interaction with the Earth's magnetosphere. Magnetic field lines from the simulation reveal the formation of magnetic flux ropes during periods with northward interplanetary magnetic field. These flux ropes result from multiple reconnection processes between the lobes field lines and draped magnetosheath field lines that are convected around the flank of the magnetosphere. The flux ropes identified in the simulation are consistent with features observed in the magnetic field measured by Hawkeye-1 during some high-latitude magnetopause crossings.

Description: Analysis of partical data from the CRIT 2 experiment, studying Alfven's critical ionization velocity (CIV) effect, shows that the density of newly created ions (presumably Ba(+) from the shaped-charge beam) is consistent with the increase in total plasma density measured by the independent RF plasma probe on board (Swenson et al., 1990) at the most active time period. We model this ion production using the measured electron flux data and the neutral barium model of Stenbaek-Nielsen et al. (1990a). To identify the main source mechanisms which may contribute most to the barium ionization, a simple model for the barium ion density at the payload location is developed based on Liouvilles theorem. We estimate that the electron impact ionization is responsible for 90% of the barium ion production observed by CRIT 2 in the first release and up to 45% in the second release. By employing a two-state approximation calculation (Rapp and Francis, 1962), the Ba-O(+) charge exchange cross section is found to range from about 2.0 X 10(exp -17) sq cm at a velocity of 4 km/s to 2.0 X 10(exp -15) sq cm at a velocity of 20 km/s. This result suggests that the Ba-O(+) charge exchange is probably dominant among all the non-CIV ionization processes. By considering the charge exchange process in our density model, the barrium ion densities are calculated for the two releases on CRIT II. The comparison between the model results and the observed data is found to be resonably consistent if the cross sections, as calculated above, are multiplied by 0.3 for the first release and 1.0 for the second release. Our result suggests that the charge exchange process could be the most important non-CIV ionization mechanism in the CRIT II experiment and it should be considered carefully case by case in CIV experiments.

Description: A selected set of far ultra violet images of Earth have been analyzed quanitatively to establish their validity for studying thermospheric weather. The set of images chosen for the study was restricted to mostly geomagnetically quiet conditions in order to obtain a baseline understanding of the relationship between the observations and thermospheric phenomenology. The images included low to modrerate solar activity levels. A new model was developed to generate global dayglow images using first principles methods. The mass Spectrometer/incoherent scatter (MSIS-86) model was used to predict the thermospheric concentrations. The analyses of thermospheric images observed in the 123 to 160-nm nominal passband show that the spectral composition for observations on the projected earth disk is dominated by O I 130.4-nm radiation (85-90%), with concentrations from O I 135.6-nm and N2 Lyman-Birge-Hopefield (LBH) bands of about 5-8% each. The synthetic images reproduce the global features of the observed images rather well. Differences between the model and the data are attributed to real atmospheric effects, such as atomic oxygen depletions which are not well reproduced by the MSIS model when geomagnetic activity is elevated. The absoute values recorded were 38-54% higher than predicted. We attribute this discrepancy to low values of the solar extreme ultraviolet irradiances used in the model. Images obtained in the 136 to 165-nm nominal passband are a factor of 2.7 greater than the model. The excess signal observed is most likely due to a long wavelength tail in the instrument sensitivity which allowed Rayleigh scattered sunlight between 180and 250nm to be detected. The understanding of the Dynamics Explorer (DE 1) images gained by this study provides the basis for future work on the global response of the theremosphere to geomagnetic forcing.

Description: There are four low-frequency modes which may propagate in a high-beta nearly bi-Maxwellian plasma. These are the magnetosonic, Alfven, ion acoustic, and mirror modes. This manuscript defines a procedure based on linear Vlasov theory for the unique identification of these modes by use of transport ratios, dimensionless ratios of the fluctuating field and plasma quantities. A single parameter, the mode deviation is calculated using the plasma and magnetic field data gathered by the Active Magnetospheric Particle Tracer Explorers/Ion Release Module (AMPTE/IRM) spacecraft to identify the modes observed in the terrestial magnetosheath near the magnetopause. As well as determining the mode which best describes the observed fluctuations, it gives us a measure of whether or not the resulting identification is unique. Using 17 time periods temporally close to a magnetopause crossing, and confining our study to the frequency range from 0.01 to 0.04 Hz, we find that the only clearly identified mode in this frequency range is the mirror mode. Most commonly, the quasi-perpendicular mirror mode (with wave vector k roughly perpendicular to the background magnetic field B(sub zero) is observed. In two events the quasi-parallel mirror mode k parallel B(sub zero) was identified.

Description: Low-altitude spacecraft on magnetotail field lines often detect a distinctive signature in the precipitating ion flux. A velocity-dispersed ion structure is often observed near the poleward boundary of the auroral oval. At the low-latitude edge of this structure an absence of precipitating ions is seen, previously referred to as 'the gap,' separating the velocity -- dispersed ions at the higher latitudes from the more diffuse, plasma sheet-like ions at lower latitudes. We present a model of low-altitude particle precipitation that reproduces these observed features in the ion spectra and provides a quantitative estimate of the downtail plasma sheet properties. The model calculations are compared with observations from the Akebono spacecraft. In this model, the dispersed ion velocity signature maps to a region in the distant plasma sheet where the plasma has a field-aligned bulk flow. The gap maps to a region in the distant magnetotail where the ion fluxes are below the detection threshold of the instrument, due to the low plasma sheet density and temperature in that region.

Description: Earthquake recurrence data from the Pallett Creek and Wrightwood paleoseismic sites on the San Andreas fault appear to show temporal variations in repeat interval. We investigate the interaction between strike-slip faults and auxiliary reverse and normal faults as a physical mechanism capable of producing such variations. Under the assumption that fault strength is a function of fault-normal stress (e.g. Byerlee's Law), failure of an auxiliary fault modifies the strength of the strike-slip fault, thereby modulating the recurrence interval for earthquakes. In our finite element model, auxiliary faults are driven by stress accumulation near restraining and releasing bends of a strike-slip fault. Earthquakes occur when fault strength is exceeded and are incorporated as a stress drop which is dependent on fault-normal stress. The model is driven by a velocity boundary condition over many earthquake cycles. Resulting synthetic strike-slip earthquake recurrence data display temporal variations similar to observed paleoseismic data within time windows surrounding auxiliary fault failures. Our simple model supports the idea that interaction between a strike-slip fault and auxiliary reverse or normal faults can modulate the recurrence interval of events on the strike-slip fault, possibly producing short term variations in earthquake recurrence interval.

Description: Global Positioning System (GPS) measurements spanning approximately 3 years have been used to determine velocities for 7 sites on the Australian, Pacific and Antarctic plates. The site velocities agree with both plate model predictions and other space geodetic techniques. We find no evidence for internal deformation of the interior of the Australian plate. Wellington, New Zealand, located in the Australian-Pacific plate boundary zone, moves 20 +/- 5 mm/yr west-southwest relative to the Australian plate. Its velocity lies midway between the predicted velocities of the two plates. Relative Euler vectors for the Australia-Antarctica and Pacific-Antarctica plates agree within one standard deviation with the NUVEL-1A predictions.

Description: Coseismic surface deformation associated with the M(sub w) 6.1, April 23, 1992, Joshua Tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from Global Positioning System and trilateration measurements. The rms signal to noise ratio for these inferred displacements is 1.8 with near-fault displacement estimates exceeding 40 mm. In order to determine the long-wavelength distribution of slip over the plane of rupture, a Tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right-lateral slip and zero surface slip constraints. The resulting slip distribution yields a geodetic moment estimate of 1.7 x 10(exp 18) N m with corresponding maximum slip around 0.8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations. From empirical Green's functions analyses, a rupture duration of 5 s is obtained which implies a rupture radius of 6-8 km. Most of the inferred slip lies to the north of the hypocenter, consistent with northward rupture propagation. Stress drop estimates are in the range of 2-4 MPa. In addition, predicted Coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters; most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0.1 MPa. In contrast, predicted stress changes are near zero at the hypocenter of the M(sub w) 7.3, June 28, 1992, Landers earthquake which nucleated about 20 km beyond the northernmost edge of the Joshua Tree rupture. Based on aftershock migrations and the predicted static stress field, we speculate that redistribution of Joshua Tree-induced stress perturbations played a role in the spatio-temporal development of the earth sequence culminating in the Landers event.

Description: Ions that are energized at quasi-parallel collisionless shocks and move back upstream generate low-frequency waves, largely on the fast/magnetosonic branch. At sufficient Mach number, the waves are convected back into the shock, lead to shock re-formation, and are mode-converted into downstream (magnetosheath) Alfvenic turbulence. Other waves are generated more locally at the interface of the incoming solar wind and the partially thermalized plasma. This paper reviews how recent simulation studies of collisionless shocks in conjunction with linear kinetic theory and proper wave diagnostics have aided in our understanding of the upstream and magnetosheath waves.

Description: A variety of suprathermal and energetic ion distributions are found upstream from shocks. Some distributions, such as field-aligned beams, are generated directly at the shock either through reflection processes or through leakage from the hotter downstream region. Other distributions, such as intermediate distributions, evolve from these parent distributions through wave-particle interactions. This paper reviews our current understanding of the creation and evolution of suprathermal distributions at shocks. Examples of suprathermal ion distributions are taken from observations at the Earth's bow shock. Particular emphasis is placed on the creation of field-aligned beams and specularly reflected ion distributions and on the evolution of these distributions in the Earth's ion foreshock. However, the results from this heavily studied region are applicable to interplanetary shocks, bow shocks at other planets, and comets.

Description: In this study, Atmosphere Explorer data and model results for the ion and electron temperature and the density of N(+), O(+), H(+), and He(+) between 120 and 1400 km altitude are compared for two midlatitude ranges (L=2 and L=4), noon and midnight local time, winter and summer, at solar minimum. The data for the heavy atomic ions (O(+) and N(+)) show that their densities are greater at noon than at midnight for a given season and greater in summer than winter for a given local time. There is only a weak latitudinal variation in the density of these ions. The data show that the light ion (H(+) and He(+)) densities are greater at midnight than at noon and are generally greater in winter than summer. There is a strong latitudinal variation of the light ion densities, with the densities decreasing with increasing latitude. The model densities are in good agreement with the AE densities for N(+), O(+), and H(+). Model He(+) densities are lower, by a factor of 2 or more, than the measured densities. Model ion and electron temperatures agree well with the measured temperatures with only a modest increase in plasmapheric heating.

Description: From a study of the 21 largest geomagnetic storms during solar cycle 21, a strong correlation is established between the ring current index Dst and the time-weighted accumulation of the 1-hour auroral electrojets indices, AE and AL. The time-weighted accumulation corresponds to convolution of the auroral electrojet indices with an exponential weighting function with an e-folding time of 9.4 hours. The weighted indices AE(sub w) and AL(sub w) have correltation coefficients against Dst ranging between 0.8 and 0.95 for 20 of the 21 storms. Correlation over the entire solar cycle 21 database is also strong but not as strong as for an individual storm. A set of simple Dst prediction functions provide a first approximation of the inferred dependence, but the specific functional relationship of Dst (AL(sub w)) or Dst (AL(sub w)) varies from one storm to the next in a systematic way. This variation reveals a missing parametric dependence in the transfer function. However, our results indicate that auroral electroject indices are potentially useful for predicting storm time enhancements of ring current intensity with a few hours lead time.

Description: The poloidal mode field line resonance in the Earth's dipole magnetic field is investigated using cold plasma ideal MHD simulations in dipole geometry. In order to excite the poloidal mode resonance, we use either an initial or a continuous velocity perturbation to drive the system. The perturbation is localized at magnetic shell L = 7 with plasma flow in the radial direction (electric field component in the azimuthal direction). It is found that with the initial perturbation alone, no polodial mode resonance can be obtained and the initially localized perturbation spreads out across all magnetic L shells. With the continuous perturbation, oscillating near the poloidal resonance frequency, a global-scale poloidal cavity mode can be obtained. For the first time, a localized guided poloidal mode resonance is obtained when a radial component of electric field is added to the initial perturbation such that the curl of the electric field is everywhere perpendicular to the background dipole magnetic field. During the localized poloidal resonance, plasma vortices parallel/antiparallel to the background dipole magnetic field B(sub 0). This circular flow, elongated radially, results in twisting of magnetic field flux tubes, which, in turn, leads to the slowdown of the circular plasma flow and reversal of the plasma vortices. The energy associated with the localized poloidal resonance is conserved as it shifts back and forth between the oscillating plasma vortices and the alternately twisted magnetic flux tubes. In the simulations the eigenfunctions associated with the localized poloidal resonance are grid-scale singular functions. This result indicates that ideal MHD is inadequate to describe the underlying problem and nonideal MHD effects are needed for mode broadening.

Description: The long-term time series of global ozone from the Nimbus-7 Solar Backscatter Ultraviolet instrument (SBUV) (Nov. 1978-June 1990) are extended through June 1994 by using measurements from the NOAA-11 SBUV/2. The data sets are merged based upon comparisons during the 18-month overlap period in which both instruments were operational. During this period, the average offset between the two time series is less than 2% in total ozone, and less than 6% in Umkehr layers 1-10. A linear-regression trend model is applied to the extended time series to calculate updated trends as a function of latitude and altitude. Trends through June 1994 are 1.5-2% per decade less negative than through June 1990 in the tropical middle stratosphere (35-40 km) and in the upper stratosphere (45-50 km) at mid-latitudes. In the lower stratosphere, the trends are nearly 1.5% per decade more negative in the southern hemisphere tropical regions to 25 deg S, but remain relatively unchanged elsewhere. The seasonal structure of the total ozone trends is similar to past trend study results, but the magnitude of the seasonal trend can vary by 2% per decade depending on the length of the time series. Both Total Ozone Mapping Spectrometer (TOMS) (through April 1993) and SBUV total ozone time series show small negative trends in the equatorial region, though they are not statistically at the 2-sigma level.

Description: A new technique for estimating electron density perturbation amplitudes of traveling ionospheric disturbances (TIDs), using HF radar data, is presented. TIDs are observed in HF radar data as enhancements of the ground-scattered power which propagate through the radar's field of view. These TIDs are the ionospheric manifestation of atmospheric acoustic-gravity waves. TID electron density perturbation amplitudes were estimated by simulating the radar returns, using HF ray tracing through a model ionosphere perturbed by a model gravity wave. The simulation determined the return power in the ground-scattered portion of the signal as a function of range, and this was compared to HF radar data from the Goose Bay HF radar at a time when evidence of gravity waves was present in the data. By varying the amplitude of the electron density perturbation in the model it was possible to estimate the perturbation of the actual wave. It was found that the perturbations that are observed by the Goose Bay HF radar are of the order of 20% to 35%. It was also found that the number of observable power enhancements, and the relative amplitudes of these enhancements, depended on the vertical thickness of the gravity wave's source region. From the simulations and observations it was estimated that the source region for the case presented here was approximately 20 km thick. In addition, the energy in the wave packet was calculated and compared to an estimate of the available energy in the source region. It was found that the wave energy was about 0.2% of the estimated available source region energy.

Description: Distribution functions of ions heated in quasi-perpendicular bow shocks have a large perpendicular temperature anisotropy that provides free energy for the growth of Alfven ion cyclotron (AIC) waves and mirror waves. Both types of waves have been observed in the Earth's magnetosheath downstream of quasi-perpendicular shocks. We use a two-dimensional hybrid simulations to give a self-consistent description of the evolution of the wave spectra downstream of quasi-perpendicular shocks. Both mirror and AIC waves are identified in the simulated magnetosheath. They are generated at or near the shock front and convected away from it by the sheath plasma. Near the shock, the waves have a broad spectrum, but downstream of the shock, shorter-wavelength modes are heavily damped and only longer-wavelength modes persist. The characteristics of these surviving modes can be predicted with reasonable accuracy by linear kinetic theory appropriate for downstream conditions. We also follow the evolution of the ion distribution function. The shocked ions that provide the free energy for wave growth have a two-component distribution function. The halo is initially gyrophase-bunched and extremely anisotropic. Within a relatively short distance downstream of the shock (of the order of 10 ion inertial lengths), wave-particle interactions remove these features from the halo and reduce the anisotropy of the distribution to near-threshold levels for the mirror and AIC instabilities. A similar evolution has been observed for ions at the Earth's bow shock.

Description: The following areas of interest in this progress report are: (1) the continuation of software development in the examination of F-region gravity-wave power using in-situ data from the Atmosphere Explorer (AE-E); (2) the inquiry into the use of the San Marco data for the study of the initiation and growth of bubbles, particularly when the satellite passes through the early evening hours at relatively high altitudes, and the development of bubbles using not only the San Marco data but includes the use of airglow observations made in Hawaii; and (3) the promising development in the observation of distinct well formed waves at about 400 km altitude in the equatorial region. These waves look very much like waves seen over the polar cap that are attributed to internal gravity waves in the neutral atmosphere driving ionization up and down the magnetic field lines. These equatorial waves show no modulation of the total ion concentration.

Description: This final report was concerned with the ideas that: (1) magnetospheric boundary layers link disparate regions of the magnetosphere-solar wind system together; and (2) global behavior of the magnetosphere can be understood only by understanding its internal linking mechanisms and those with the solar wind. The research project involved simultaneous research on the global-, meso-, and micro-scale physics of the magnetosphere and its boundary layers, which included the bow shock, the magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical, and simulation projects were performed on these subjects, as well as comparisons of theoretical results with observational data. Other related activity included in the research included: (1) prediction of geomagnetic activity; (2) global MHD (magnetohydrodynamic) simulations; (3) Alfven resonance heating; and (4) Critical Ionization Velocity (CIV) effect. In the appendixes are list of personnel involved, list of papers published; and reprints or photocopies of papers produced for this report.

Description: A brief description of the major activities pursued during the last year (March 1994 - February 1995) of this grant are: (1) the development of a 200 km to 1 Re, O(+) H(+) Model; (2) the extension of the E x B convection heating study to include centrifugal effects; (3) the study of electron precipitation effects; (4) the study of wave heating of O(+); and (5) the polar wind acceleration study. A list of both papers published and papers submitted, along with a proposal for next year's study and a copy of the published paper is included.

Description: The Langmuir probe flown as part of the Solar Array Module Plasma Interactions Experiment (SAMPIE) package aboard the space shuttle flight STS-62 was used to determine plasma potential fluctuations in the vicinity of the shuttle. The broadband noise was observed at frequencies 250 - 20,000 Hz. Measurements were performed in ram conditions; thus, it seems reasonable to believe that the influence of spacecraft operations on plasma parameters was absolutely negligible. The average spectrum of fluctuations is in agreement with theoretical predictions. The influence on the observed spectra of arcing generated by high negative bias voltages applied to solar cell samples is briefly discussed.

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: The SBUV/2 instrument onboard the NOAA-11 satellite made daily solar spectral irradiance measurements in the wavelength region 160405 nm at 1.1 nm resolution between January 1989 and October 1994. These observations continued the uninterrupted series of solar measurements begun by the Nimbus-7 SBUV in 1978 and continued by NOAA-9 SBUV/2. While the measurements made by the SBUV-series instruments furnish an excellent data base for studies of solar UV variability, these instruments do not have an internal mew to evaluate and correct for long-term instrument sensitivity degradation, needed to evaluate solar cycle timescale irradiance change. During yearly Shuttle flights the Shuttle SBUV (SSBUV) also performs solar spectral irradiance measurements in the wavelength region 200 to 400 nm with an instrument that is calibrated preflight, inflight, and postflight. Comparisons between the simultaneous NOAA-11 SBUV/2 and SSBUV solar measurements are used to identify and correct long term sensitivity changes in the satellite instrument. The NOAA-11 data will then be used to evaluate long-term solar change. We present a progress report on the above process. At this preliminary stage uncertainties in the calibration transfer between SSBUV and NOAA-11 SBUV/2 are too large to accurately evaluate long-term solar change near the A1 edge, but solar rotational activity variations can be evaluated. We find that rotational activity declined from roughly 6% peak-to-peak (p-p) near the maximum of solar cycle 22 in 1989-1991 to approximately 3% p-p in mid 1992 and 2% p-p by mid 1994. Emphasizing rotational variations, comparisons between the 200 nm data and the NOAA-11 Mg II proxy index are presented.

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: A far-wing line shape theory based on the binary collision and quasistatic approximations that is applicable for both the low- and high-frequency wings of the vibration-rotational bands has been developed. This theory has been applied in order to calculate the frequency and temperature dependence of the continuous absorption coefficient for frequencies up to 10,000 cm(exp -1) for pure H2O and for H2O-N2 mixtures. The calculations were made assuming an interaction potential consisting of an isotropic Lennard-Jones part with two parameters that are consistent with values obtained from other data, and the leading long-range anisotropic part, together with the measured line strengths and transition frequencies. The results, obtained without the introduction of adjustable parameters, compare well with the existing laboratory data, both in magnitude and in temperature dependence. This leads us to the conclusion that the water continuum can be explained in terms of far-wing absorption. Current work in progress to extend the theory and to validate the theoretically calculated continuum will be discussed briefly.

Description: The Argon Release for Controlled Studies (ARCS) 4 sounding rocket was launched northward into high altitude from Poker Flat Research Range on February 23, 1990. The vehicle crossed geomagnetic field lines containing discrete auroral activity. An instrumented subpayload released 100-eV and 200-eV Ar(+) ion beams sequentially, in a direction largely perpendicular to both the local geomagnetic field and the subpayload spin axis. The instrumented main payload was separated along field lines from the beam emitting subpayload by a distance which increased at a steady rate of approximately 2.4 m/s. Three dimensional mass spectrometric ion observations of ambient H(+) and O(+) ions, obtained on board the main payload, are presented. Main payload electric field observations in the frequency range 0-16 kHz, are also presented. These observations are presented to demonstrate the operation of transverse ion acceleration, which was differential with respect to ion mass, primarily during 100-eV beam operations. The preferential transverse acceleration of ambient H(+) ions, as compared with ambient O(+) ions, during the second, third, fourth, and fifth 100-eV beam operations, is attributed to a resonance among the injected Ar(+) ions, beam-generated lower hybrid waves, and H(+) ions in the tail of the ambient thermal distribution. This work provides experimental support of processes predicted by previously published theory and simulations.

Description: Spectral and spatial images obtained with the Imaging Spectrometric Observatory on the ATLAS 1 and Spacelab 1 missions are used to study the ultraviolet emissions of nitric oxide in the thermosphere. By synthetically fitting the measured NO gamma bands, intensities are derived as a function of altitude and latitude. We find that the NO concentrations inferred from the ATLAS 1 measurements are higher than predicted by our thermospheric airglow model and tend to lie to the high side of a number of earlier measurements. By comparison with synthetic spectral fits, the shape of the NO gamma bands is used to derive temperature as a function of altitude. Using the simultaneous spectral and spatial imaging capability of the instrument, we present the first simultaneously acquired altitude images of NO gamma band temperature and intensity in the thermosphere. The lower thermospheric temperature images show structure as a function of altitude. The spatial imaging technique appears to be a viable means of obtaining temperatures in the middle and lower thermosphere, provided that good information is also obtained at the higher altitudes, as the contribution of the overlying, hotter NO is nonnegligible. By fitting both self-absorbed and nonabsorbed bands of the NO gamma system, we show that the self absorption effects are observable up to 200 km, although small above 150 km. The spectral resolution of the instrument (1.6 A) allows separation of the N(+)(S-5) doublet, and we show the contribution of this feature to the combination of the NO gamma (1, 0) band and the N(+)(S-5) doublet as a function of altitude (less than 10% below 200 km). Spectral images including the NO delta bands support previous findings that the fluorescence efficiency is much higher than that determined from laboratory measurements. The Spacelab 1 data indicate the presence of a significant population of hot NO in the vehicle environment of that early shuttle mission.