ALBERT

Description: The thinning and intensification of the cross tail current sheet during the substorm growth phase are analyzed during the CDAW 6 substorm (22 Mar. 1979) using two complementary methods. The magnetic field and current sheet development are determined using data from two spacecraft and a global magnetic field model with several free parameters. These results are compared with the local calculation of the current sheet location and structure previously done by McPherron et al. Both methods lead to the conclusion that an extremely thin current sheet existed prior to the substorm onset, and the thicknesses estimated by the two methods at substorm onset agree relatively well. The plasma data from the ISEE 1 spacecraft at 13 R(sub E) show an anisotropy in the low energy electrons during the growth phase which disappears just before the substorm onset. The global magnetic model results suggest that the field is sufficiently stretched to scatter such low energy electrons. The strong stretching may improve the conditions for the growth of the ion tearing instability in the near Earth tail at substorm onset.

Description: An approach to the study of the solar wind-magnetosphere interaction by signal type, that is, by examining the effect in the magnetosphere of well defined interplanetary structures, is presented. Focus is on the response of the magnetosphere to interplanetary magnetic clouds. Among their properties are: the slow and smooth variation of the magnetic field vector, with fluctuation level well below common interplanetary values; the similarly well behaved bulk flow; the wide range of field and flow parameters; and the longevity of passage (1 to 2 days). If the magnetic cloud is oriented such that a long period of uninterruptedly northward pointing field is followed by a long interval of continuously southward pointing field, then the transition of the magnetosphere from a quiescent state (the 'ground state') to a very active state can be studied, the latter being sustained by continued forcing from the magnetic cloud. A synopsis of the main findings of a recent study in such an interaction is given, concentrating on the substorm activity attending the second part of cloud passage.

Description: The stability of the geomagnetic tail is investigated on the basis of three dimensional resistive magnetohydrodynamic simulations, using different dynamic constraints and different initial equilibria. Different forms of the energy equation for isotropic pressure are found to have no significant effect on the dynamic growth of a resistive tearing instability, which is responsible for near Earth reconnection, plasmoid formation and ejection, and the generation of fast plasma flows. The constraints of a modified double adiabatic approach, however, can quench the tearing instability through the development of large, mirror type, anisotropies in the boundary regions of the plasma sheet, unless isotropization occurs on fast, nearly Alfvenic, time scales. The presence of a net cross tail magnetic field component B(sub yN) can reduce the growth of the instability without complete stabilization. An increase of B(sub z) from midnight toward the tail flanks, however, by more than a factor of about 3, apparently completely stabilizes the tearing mode. Stabilization and destabilization thus may depend on properties and constraints (and their release) in regions other than the neutral sheet where reconnection is initiated.

Description: Issues concerning the 'driven' versus 'unloading' nature of substorms are presented. The original concepts attendant to this debate are presented and substorms are concluded to inextricably combine aspects that are driven with aspects that represent a loading-unloading system. For isolated substorms, the magnetosphere-ionosphere system is shown to exhibit a bimodal response to solar wind changes. A 20 min response characteristic is associated with the driven aspect of substorms, while a 1 hr response time is associated with unloading. It is found that for strong solar wind input conditions, the magnetospheric response becomes more nearly unimodal. This is interpreted in terms of a nonlinear dynamical evolution of the system. Simple analog models are described which capture the essence of the nonlinear magnetospheric behavior. These models exhibit chaotic transitions for strong driving conditions: this may explain the observed behavior of the magnetosphere during strong geomagnetic activity.

Description: A pseudobreakup is a phenomenon similar to the substorm expansive phase onset, including an activation of an auroral arc, a burst of Pi2 micropulsations, and enhancement of the westward electrojet. However, these effects are weak and a pseudobreak is generally assumed to be very localized. The pseudobreakups are discussed based on simultaneous observations made in space and on the ground during the substorm growth phase. In the events studied the main features listed above are found, but the significance of the localization is unclear. The optical pseudobreakup, with associated magnetic perturbations, is highly localized, but simultaneously a wide local time sector of the auroral oval may be activated. The major differences between pseudobreakups and substorm expansive phase onsets are concluded to be the intensity and the development that follows. Careful study of pseudobreakups may help to determine phase initiation, and the role of the ionosphere-magnetosphere coupling in the substorm process.

Description: A 54.95-MHz coherent backscatter radar, an ionosonde and the magnetometer located at Trivandrum in India (8.5 deg N, 77 deg E, 0.5 deg N dip angle) recorded large-amplitude ionospheric fluctuations and magnetic field fluctuations associated with a Pc5 micropulsation event, which occurred during an intense magnetic storm on 24 March 1991 (A(sub p) = 161). Simultaneous 100-n T-level fluctuations are also observed in the H-component at Brorfelde, Denmark (55.6 deg N gm) and at Narsarsuaq, Greenland (70.6 deg N gm). Our study of the above observations shows that the E-W electric field fluctuations in the E- and F-regions and the magnetic field fluctuations at Thumba are dominated by a near-sinusoidal oscillation of 10 min during 1730-1900 IST (1200-1330 UT), the amplitude of the electric field oscillation in the equatorial electrojet (EEJ) is 0.1-0.25 mV/m and it increases with height, while it is about 1.0 mV/m in the F-region, the ground-level H-component oscillation can be accounted for by the ionospheric current oscillation generated by the observed electric field oscillation in the EEJ and the H-component oscillations at Trivandrum and Brofelde are in phase with each other. The observations are interpreted in terms of a compressional cavity mode resonance in the inner magnetosphere and the assoicated ionospheric electric field penetrating from high latitudes to the magnetic equator.

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

Description: The Coast Ranges of the Cascadia margin are overriding the subducted Juan de Fuca/Gorda plate. We investigate the extent to which the latitudinal change in attributes related to the subduction process. These attributes include the varibale age of the subducted slab that underlies the Coast Ranges and average vertical crustal velocities of the western margin of the Coast Rnages for two markedly different time periods, the last 45 years and the last 100 kyr. These vertical crustal velocities are computed from the resurveying of highway bech marks and from the present elevation of shore platforms that have been uplifted in the late Quaternary, respectively. Topogarphy of the Coast Ranges is in part a function of the age and bouyancy of the underlying subducted plate. This is evident in the fact that the two highest topographic elements of the Coast Rnages, the Klamath Mountains and the Olympic Mountains, are underlain by youngest subducted oceanic crust. The subducted Blanco Fracture Zone in southernmost Oregon is responsible for an age discontinuity of subducted crust under the Klamath Mountains. The norhtern terminus of hte topographically higher Klamaths is offset to the north relative to the position of the underlying Blanco Fracture Zone, teh offset being in the direction of migration of the farcture zone, as dictated by relative plate motions. Vertical crustal velocities at the coast, derived from becnh mark surveys, are as much as an order of magnitude greater than vertical crustal velocities derived from uplifted shore platforms. This uplift rate discrepancy indicates that strain is accumulating on the plate margin, to be released during the next interplate earthquake. In a latitudinal sense, average Coast Rnage topography is relatively high where bench mark-derived, short-term vertical crustal velocities are highest. Becuase the shore platform vertical crustal velocities reflect longer-term, premanent uplift, we infer that a small percentage of the interseismic strain that accumulates as rapid short-term uplift is not recovered by subduction earthquakes but rather contributes to rock uplift of the Coast Ranges. The conjecture that permanent rock uplift is related to interseismic uplift is consistent with the observation that those segments of the subduction zone subject to greater interseismic uplift rates are at approximately the same latitudes as those segments of the Coast Ranges that have higher magnitudes of rock uplift over the long term.

Description: A comprehensive review is presented of the mathematical models used to represent magnetic fields in the Earth's magnetosphere, of the way existing data-based models use these methods and of the associated problems and concepts. The magnetic field has five main components: the internal field, the magnetopause, the ring current, the tail and Birkeland currents. Methods of representing separately each of these are discussed, as is the deformation of magnetic fields; Appendix B traces the connection between deformations and the Cauchy integral. A summary section lists the uses of data-based models and their likely future evolution, and Appendix A supplements the text with a set of problems.

Description: We present a map of the coseimic displacement field resulting from the Landers, California, June 28, 1992, earthquake derived using data acquired from an orbiting high-resolution radar system. We achieve results more accurate than previous space studies and similar in accuracy to those obtained by conventional field survey techniques. Data from the ERS 1 synthetic aperture radar instrument acquired in April, July, and August 1992 are used to generate a high-resolution, wide area map of the displacements. The data represent the motion in the direction of the radar line of sight to centimeter level precision of each 30-m resolution element in a 113 km by 90 km image. Our coseismic displacement contour map gives a lobed pattern consistent with theoretical models of the displacement field from the earthquake. Fine structure observed as displacement tiling in regions several kilometers from the fault appears to be the result of local surface fracturing. Comparison of these data with Global Positioning System and electronic distance measurement survey data yield a correlation of 0.96; thus the radar measurements are a means to extend the point measurements acquired by traditional techniques to an area map format. The technique we use is (1) more automatic, (2) more precise, and (3) better validated than previous similar applications of differential radar interferometry. Since we require only remotely sensed satellite data with no additioanl requirements for ancillary information. the technique is well suited for global seismic monitoring and analysis.

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

Description: A search was conducted for the signatures of Birkeland currents in the Earth's magnetic tail, using observed values of B(sub x) and B(sub y) from large sets of spacecraft data. The data were binned by x and y for -10 greater than x(sub GSM) greater than -35 and absolute value of y(sub GSM) less than or equal to 20 R(sub E) (less than or equal to 30 R(sub E) for x(sub GSM) less than or equal to -25 R(sub E)) and in each bin their distribution in the (B(sub x), B(sub y)) plane was fitted by least squares to a piecewise linear function. That gave average x-y distributions of the flaring angle between B(sub xy) and the x direction, as well as that angle's variation across the thickness of the plasma sheet. Angles obtained in the central plasma sheet differed from those derived near the lobe boundary. That is the expected signature if earthward or tailward Birkeland current sheets are embedded in the plasma sheet, and from this dfiference we derived the dawn-dusk profiles of the tail Birkeland currents for several x(sub GSM) intervals. It was found that (1) the Birkeland currents have the sense of region 1 currents, when mapped to the ionosphere; (2) both the linear current density (kiloamperes/R(sub E)) and the net magnitude of the field-aligned currents decrease rapidly down the tail; (3) the total Birkeland current at x approximately equals -10 R(sub E) equals approximately equals 500-700 kA, which is approx. 30% of the net region 1 current observed at ionospheric altitudes, in agreement with model mapping results; and (4) the B(sub z) and B(sub y) components of the interplanetary magnetic field influence the distribution of Birkeland currents in the tail.

Description: The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Description: The effect of the porous leading edge of an airfoil on the blade-vortex interaction noise, which dominates the far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near field are calculated directly from the computation. The porosity effect on the surface is modeled in two ways: (1) imposition of prescribed transpiration velocity distribution and (2) calculation of transpiration velocity distribution by Darcy's law. Results show leading-edge transpiration can suppress pressure fluctuations at the leading edge during blade-vortex interaction and consequently reduce the amplitude of propagating noise by 30% at a maximum in the near field.

Description: A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered-grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the noniterative solution of the flowfield and satisfies the continuity equation to machine zero at each time step. The method is evaluated in terms of its stability to predict two-dimensional flow about an airfoil at -90-deg incidence for varying Reynolds number and laminar/turbulent models. The variations of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisom indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

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

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

Description: New line parameters for two heavy odd nitrogen molecules HNO3 in the nu(sub 5)/2nu(sub 9) region, and ClONO2 in the nu(sub 4) region are incorporated in the analysis of high resolution i.r. atmospheric spectra. The line parameters are tested and renormalized vs laboratory spectra, and then applied to retrievals from balloon-borne and ground-based solar absorption spectra.

Description: The spectroscopic identification for the HNO3 3 nu(sub 9) - nu(sub 9) band Q branch at 830.4/cm is reported based on 0.01/cm resolution solar occultation spectra of the lower stratosphere recorded by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer and a recent analysis of this band. Least-squares fits to 0.0025/cm resolution laboratory spectra in the Q branch region indicate an integrated intensity of 0.529 x 10(exp -18)/cm/mol/sq cm at 296 K for this weak band. Stratospheric HNO3 retrievals derived from the ATMOS data are consistent with this value within its estimated uncertainty of about +/- 30%. A set of spectroscopic line parameters suitable for atmospheric studies has been generated.

Description: About 200 i.r. solar spectra recorded at 0.01/cm resolution on 71 days between November 1991 and July 1993 at the Network for the Detection of Stratospheric Change (NDSC) station at Mauna Loa, Hawaii (latitude 19.53 deg N, longitude 155.58 deg W, elevation 3.459 km) have been analyzed with a nonlinear least-squares spectral fitting technique to study temporal variations in the total column of atmospheric ethane (C2H6) above the site. The results were derived from the analysis of the unresolved nu(sub 7) band (P)Q(sub 3) subbranch at 2976.8/cm. A distinct seasonal cycle is observed with a factor of 2 variation, a maximum total column of 1.16 x 10(exp 16) mol/sq cm at the end of winter, and a minimum total column of 0.53 x 10(exp 16) mol/sq cm at the end of summer. Our measurements are compared with previous observations and model predictions.

Description: Quantitative measurements of the wavelength dependence of aerosol extinction in the 750-3400/cm spectral region have been derived from 0.01/cm resolution stratospheric solar occultation spectra recorded by the ATMOS (Atmospheric Trace Molecule Spectroscopy) Fourier transform spectrometer about 9 1/2 months after the Mt Pinatubo volcanic eruption. Strong, broad aerosol features have been identified near 900, 1060, 1190, 1720, and 2900/cm below a tangent height of approximately 30 km. Aerosol extinction measurements derived from approximately 0.05/cm wide microwindows nearly free of telluric line absorption in the ATMOS spectra are compared with transmission calculations derived from aerosol size distribution profiles retrieved from correlative SAGE (Stratospheric Aerosol and Gas Experiment) II visible and near i.r. extinction measurements, seasonal and zonally averaged H2SO4 aerosol weight percentage profiles, and published sulfuric acid optical constants derived from room temperature laboratory measurements. The calculated shapes and positions of the aerosol features are generally consistent with the observations, thereby confirming that the aerosols are predominantly concentrated H2SO4-H2O droplets, but there are significant differences between the measured and calculated wavelength dependences of the aerosol extinction. We attribute these differences as primarily the result of errors in the calculated low temperature H2SO4-H2O optical constants. Errors in both the published room temperature optical constants and the limitations of the Lorentz-Lorenz relation are likely to be important.

Description: We present vertical column abundances of H2O, N2O, HNO3, NO2, O3, HF, HCl, and ClNO3, determined from solar absorption spectra measured by the JPL MkIV interferometer from the NASA DC-8 aircraft. These observations, taken in 1987 and 1992, covered latitudes ranging from 85 deg S to 85 deg N. Although most gases display latitude symmetry, large asymmetries in H2O, HNO3, and O3 are apparent, which can be ascribed to processes enhanced by the colder Antarctic winter temperatures.

Description: Simultaneous in situ measurements of NO2, NO, O3, ClO, pressure and temperature have been made for the first time, presenting a unique opportunity to test our current understanding of the photochemistry of the lower stratospere. Data were collected from several flights of the ER-2 aircraft at mid-latitudes in May 1993 during NASA's Stratospheric Photochemistry, Aerosols and Dynamics Expedition (SPADE). The daytime ratio of NO2/NO remains fairly constant at 19 km with a typical value of 0.68 and standard deviation of +/- 17. The ratio observations are compared with simple steady-state calculations based on laboratory-measured reaction rates and modeled NO2 photolysis rates. At each measurement point the daytime NO2/NO with its measurements uncertainty overlap the results of steady-state caculations and associated uncertainty. Possible sources of error are examined in both model and measurements. It is shown that more accurate laboratory determinations of the NO + 03 reaction rate and of the NO2 cross-sections in the 200-220 K temperature range characteristic of the lower stratosphere would allow for a more robust test of our knowledge of NO(X) phtochemistry by reducing significant sources if uncertainties in the interpretation of statospheric measurements.

Description: The ability to predict short-term variations in the Earth's rotation has gained importance in recent years owing to more precise spacecraft tracking requirements. Universal time (UT1), that component of the Earth's orientation corresponding to the rotation angle, can be measured by number of high-precision space geodetic techniques. A Kalman filter developed at the Jet Propulsion Laboratory (JPL) optimally combines these different data sets and generates a smoothed times series and a set of predictions for UT1, as well as for additional Earth orientation components. These UT1 predictions utilize an empirically derived random walk stochastic model for the length of the day (LOD) and require frequent and up-to-date measurements of either UT1 or LOD to keep errors from quickly accumulating. Recent studies have shown that LOD variations are correlated with changes in the Earth's axial atmospheric angular momentum (AAM) over timescales of several years down to as little as 8 days. AAM estimates and forecasts out to 10 days are routinely available from meteorological analysis centers; these data can supplement geodetic measurements to improve the short-term prediction of LOD and have therefore been incorporated as independent data types in the JPL Kalman filter. We find that AAM and, to a lesser extent, AAM forecast data are extremely helpful in generating accurate near-real-time estimates of UT1 and LOD and in improving short-term predictions of these quantities out to about 10 days.

Description: Rotor noise prediction codes predict the thickness and loading noise produced by a helicopter rotor, given the blade motion, rotor operating conditions, and fluctuating force distribution over the blade surface. However, the criticality of these various inputs, and their respective effects on the predicted acoustic field, have never been fully addressed. This paper examines the importance of these inputs, and the sensitivity of the acoustic predicitions to a variation of each parameter. The effects of collective and cyclic pitch, as well as coning and cyclic flapping, are presented. Blade loading inputs are examined to determine the necessary spatial and temporal resolution, as well as the importance of the chordwise distribution. The acoustic predictions show regions in the acoustic field where significant errors occur when simplified blade motions or blade loadings are used. An assessment of the variation in the predicted acoustic field is balanced by a consideration of Central Processing Unit (CPU) time necessary for the various approximations.

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

Description: During the final, low solar activity phase of the Pioneer Venus (PV) mission, the Orbiter Ion Mass Spectrometer (OIMS) measurements found all ion species, in the midnight-dusk sector, reduced in concentration relative to that observed at solar maximum. Molecular ion species comprised a greater part of the total ion concentration as O(+) and H(+) had the greatest depletions. The nightside ionospheric states were strikingly similar to the isolated solar maximum 'disappearing' ionospheres. Both are very dynamic states characterized by a rapidly drifting plasma and 30-100 eV superthermal O(+) ions.

Description: In October, 1992 the Pioneer Venus Orbiter entered the atmosphere of Venus, ending nearly 14 years of observations at Venus. Prior to the entry into the atmosphere and subsequent loss of the spacecraft careful management of spacecraft resources had allowed the acquisition of much low altitude data over the nightside of the planet. The long duration of the Pioneer Venus mission has enabled us to study the ionosphere and atmosphere of Venus under different levels of solar activity.

Description: Atmospheric drag measurements from the orbital decay of the Pioneer Venus Orbiter and Magellan spacecraft have recently been obtained of the Venus dayside and nightside atmosphere between 130 and 210 km during a period of low solar activity. These new measurements, combined with the earlier Pioneer Venus drag measurements (1978-80) obtained near the maximum of the 11-year solar cycle, have allowed the detection of the detailed response of temperature, atomic oxygen and carbon dioxide to solar variations. We have found a weak but detectable temperature response on the dayside which is in accord with the response predicted by Keating and Bougher when they assumed very strong CO2 radiative cooling resulting from atomic oxygen exciting CO2 into 15 micron emission. This same radiative process may cause strong cooling in the Earth's upper atmosphere with the doubling of CO2 in the future. With decreasing solar activity, the O/CO2 ratio in the lower thermosphere is found to decrease, apparently due to decreased photodissociation of CO2 and lower temperatures. The percent decrease in atomic oxygen with decreasing solar activity on the dayside is found to be approximately the same as the percent decreases of atomic oxygen transported to the nightside. A very weak response of nightside temperatures to solar activity variations has also been detected.

Description: Disk-shaped current distributions are useful tools in modeling the magnetospheres of Earth and other planets and have also been adapted for modeling the geotail. Such models usually start with an axisymmetric vector potential but may be modified to account for observed asymmetries, variable thickness, and warping in response to an inclined orientation of the planetary dipole axis. Models of this type until now either have lacked the ability to simulate a sharp inner edge of the current and to control accurately its falloff with distance or did not allow a simple analytical representation. Here existing methods will be reviewed, after which a new class of models which overcomes the above deficiencies and also allows the modeling of current disks of finite thickness flanked by current-free regions will be presented.

Description: An improved model of Earth's gravitational field, Goddard Earth Model T-3 (GEM-T3), has been developed from a combination of satellite tracking, satellite altimeter, and surface gravimetric data. GEM-T3 provides a significant improvement in the modeling of the gravity field at half wavelengths of 400 km and longer. This model, complete to degree and order 50, yields more accurate satellite orbits and an improved geoid representation than previous Goddard Earth Models. GEM-T3 uses altimeter data from GEOS 3 (1975-1976), Seasat (1978) and Geosat (1986-1987). Tracking information used in the solution includes more than 1300 arcs of data encompassing 31 different satellites. The recovery of the long-wavelength components of the solution relies mostly on highly precise satellite laser ranging (SLR) data, but also includes Tracking Network (TRANET) Doppler, optical, and satellite-to-satellite tracking acquired between the ATS 6 and GEOS 3 satellites. The main advances over GEM-T2 (beyond the inclusion of altimeter and surface gravity information which is essential for the resolution of the shorter wavelength geoid) are some improved tracking data analysis approaches and additional SLR data. Although the use of altimeter data has greatly enhanced the modeling of the ocean geoid between 65 deg N and 60 deg S latitudes in GEM-T3, the lack of accurate detailed surface gravimetry leaves poor geoid resolution over many continental regions of great tectonic interest (e.g., Himalayas, Andes). Estimates of polar motion, tracking station coordinates, and long-wavelength ocean tidal terms were also made (accounting for 6330 parameters). GEM-T3 has undergone error calibration using a technique based on subset solutions to produce reliable error estimates. The calibration is based on the condition that the expected mean square deviation of a subset gravity solution from the full set values is predicted by the solutions' error covariances. Data weights are iteratively adjusted until this condition for the error calibration is satisfied. In addition, gravity field tests were performed on strong satellite data sets withheld from the solution (thereby ensuring their independence). In these tests, the performance of the subset models on the withheld observations is compared to error projections based on their calibrated error covariances. These results demonstrate that orbit accuracy projections are reliable for new satellites which were not included in GEM-T3.

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

Description: We present predictions of the signatures of magnetosheath particle precipitation (in the regions classified as open low-latitude boundary layer, cusp, mantle and polar cap) for periods when the interplanetary magnetic field has a southward component. These are made using the 'pulsating cusp' model of the effects of time-varying magnetic reconnection at the dayside magnetopause. Predictions are made for both low-altitude satellites in the topside ionosphere and for midaltitude spacecraft in the magnetosphere. Low-altitude cusp signatures, which show a continuous ion dispersion signature, reveal 'quasi-steady reconnection' (one limit of the pulsating cusp model), which persists for a period of at least 10 min. We estimate that 'quasi-steady' in this context corresponds to fluctuations in the reconnection rate of a factor of 2 or less. The other limit of the pulsating cusp model explains the instantaneous jumps in the precipitating ion spectrum that have been observed at low altitudes. Such jumps are produced by isolated pulses of reconnection: that is, they are separated by intervals when the reconnection rate is zero. These also generate convecting patches on the magnetopause in which the field lines thread the boundary via a rotational discontinuity separated by more extensive regions of tangential discontinuity. Predictions of the corresponding ion precipitation signatures seen by midaltitude spacecraft are presented. We resolve the apparent contradiction between estimates of the width of the injection region from midaltitude data and the concept of continuous entry of solar wind plasma along open field lines. In addition, we reevaluate the use of pitch angle-energy dispersion to estimate the injection distance.

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

Description: Simultaneous profiles of aerosol backscatter ratio were measured over Lauder, New Zealand (45 deg S, 170 deg E) on the night of November 24, 1992. Instrumentation comprised two complementary lidar systems and a backscattersonde, to give measurements at wavelengths 351, 490, 532, and 940 nm. The data from the lidars and the backscattersonde were self-consistent, enabling the wavelength dependence of aerosol backscatter to be determined as a function of altitude. This wavelength-dependence is a useful parameter in radiative transfer calculations. In the stratosphere, the average wavelength exponent between 351 and 940 nm was -1.23 +/- 0.1, which was in good agreement with values derived from measured physical properties of aerosols.

Description: High-energy electrons have been measured systematically in a low-altitude (520 x 675 km), nearly polar (inclination = 82 deg) orbit by sensitive instruments onboard the Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX). Count rate channels with electron energy thresholds ranging from 0.4 MeV to 3.5 MeV in three different instruments have been used to examine relativistic electron variations as a function of L-shell parameter and time. A long run of essentially continuous data (July 1992 - July 1993) shows substantial acceleration of energetic electrons throughout much of the magnetosphere on rapid time scales. This acceleration appears to be due to solar wind velocity enhancements and is surprisingly large in that the radiation belt 'slot' region often is filled temporarily and electron fluxes are strongly enhanced even at very low L-values (L aprroximately 2). A superposed epoch analysis shows that electron fluxes rise rapidly for 2.5 is approximately less than L is approximately less than 5. These increases occur on a time scale of order 1-2 days and are most abrupt for L-values near 3. The temporal decay rate of the fluxes is dependent on energy and L-value and may be described by J = Ke-t/to with t(sub o) approximately equals 5-10 days. Thus, these results suggest that the Earth's magnetosphere is a cosmic electron accelerator of substantial strength and efficiency.

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

Description: The U.S. National Aeronautics and Space Administration (NASA) Balloon Program has been highly successful since recovering from the catastrophic balloon failure problems of the early to mid 1980s. Balloons have continued to perform at unprecedented success rates. The comprehensive research and development (R&D) effort has continued with advances being made across the spectrum of balloon related disciplines. The long duration balloon project will be transitioning from a development effort to an operational capability this year. Recently, emphasis has been placed on the development and implementation of new support systems and facilities. A new permanent launch facility at Fort Sumner, New Mexico has been established. New ground station support equipment is being implemented, and a new heavy load launch vehicle is scheduled to be implemented in 1992. The progress, status and future plans for these and other aspects of the NASA program will be presented.

Description: The catastrophic balloon failure during the first half of the 1980's identified the need for a comprehensive and continuing balloon research and development (R&D) commitment by NASA. Technical understanding was lacking in many of the disciplines and processes associated with scientific ballooning. A comprehensive balloon R&D plan was developed in 1986 and implemented in 1987. The objectives were to develop the understanding of balloon system performance, limitations, and failure mechanisms. The program consisted of five major technical areas: structures, performance and analysis, materials, chemistry and processing, and quality control. Research activitites have been conducted at NASA/Goddard Space Flight Center (GSFC)-Wallops Flight Facility (WFF), other NASA centers and government facilities, universities, and the balloon manufacturers. Several new and increased capabilities and resources have resulted from this activity. The findings, capabilities, and plan of the balloon R&D program are presented.

Description: Caps have been used to structurally reinforce scientific research balloons since the late 1950's. The scientific research balloons used by the National Aeronautics and Space Administration (NASA) use internal caps. A NASA cap placement specification does not exist since no empirical information exisits concerning cap placement. To develop a cap placement specification, NASA has completed two in-hangar inflation tests comparing the structural contributions of internal caps and external caps. The tests used small scale test balloons designed to develop the highest possible stresses within the constraints of the hangar and balloon materials. An externally capped test balloon and an internally capped test balloon were designed, built, inflated and simulated to determine the structural contributions and benefits of each. The results of the tests and simulations are presented.

Description: Rocket data have been used to evaluate the characteristics of precipitating relativistic electrons and their effects on the electrodynamic structure of the middle atmosphere. These data were obtained at Poker Flat, Alaska, on May 13 and 14, 1990, during a midday, highly relativistic electron (HRE) precipitation event. Solid state detectors were used to measure the electron fluxes and their energy spectra. An X ray scintillator was included on each flight to measure bremsstrahlung X rays produced by energetic electrons impacting on the upper atmosphere. However, these were found the be of negligible importance for this particular event. The energy deposition by the electrons has been determined from the flux measurements and compared with in situ measurements of the atmospheric electrical response. The electrodynamic measurements were obtained by the same rockets and additionally on May 13, with an accompanying rocket. The impact flux was highly irregular, containing short-lived bursts of relativistic electrons, mainly with energies below 0.5 MeV and with fluxes most enhanced between pitch angles of 0 deg - 20 deg. Although the geostationary counterpart of this measured event was considered to be of relatively low intensity and hardness, energy deposition peaked near 75 km with fluxes approaching an ion pair production rate in excess of 100/cu cm s. This exceeds peak fluxes in relativistic electron precipitation (REP) events as observed by us in numerous rocket soundings since 1976. Conductivity measurements from a blunt probe showed that negative electrical conductivities exceeded positive conductivities down to 50 km or lower, consistent with steady ionization by precipitating electrons above 1 MeV. These findings imply that the electrons from the outer radiation zone can modulate the electrical properties of the middle atmosphere to altitudes below 50 km. During the decline and activity minimum of the current solar cycle, we anticipate the occurence of similar events but with fluxes 1-2 orders of magnitude above that reported here, based on studies of earlier solar cycles (e.g., Baker et al., 1993).

Description: Nitrous oxide (N2O) measured on board the ER-2 aircraft during the Airborne Arctic Stratospheric Expedition 2 (AASE 2) has been used to monitor descent of air inside the Arctic vortex between October 1991 and March 1992. Monthly mean N2O fields are calculated from the flight data and then compared with mean fields calculated from the high-resolution Geophysical Fluid Dynamics Laboratory general circulation model SKYHI in order to evaluate the model's simulation of the polar vortex. From late fall through winter the model vortex evolves in much the same way as the 1991-1992 vortex, with N2O gradients at the edge becoming progressively steeper. The October to March trends in N2O profiles inside the vortex are used to verify daily net heating rates in the vortex that were computed from clear sky radiative heating rates and National Meteorological Center temperature observations. The computed heating rates successfully estimate the descent of vortex air from December through February but suggest that before December, air at high latitudes may not be isolated from the midlatitudes. SKYHI heating rates are in good agreement with the computed rates but tend to be slightly higher (i.e., less cooling) due to meteorological differences between SKYHI and the 1991-1992 winter. Three ER-2 flights measured N2O just north of the subtropical jet. These low-midlatitude profiles show only slight differences from the high-midlatitude profiles (45 deg - 60 deg N), indicating strong meridional mixing in the midlatitude 'surf zone.' Mean midwinter N2O profiles inside and outside the vortex calculated from AASE 2 data are shown to be nearly identical to 1989 AASE profiles, pointing to the N2O/potential temperature relationship as an excellent marker for vortex air.

Description: The combined Nimbus 7 solar backscattered ultraviolet (SBUV) and NOAA 11 SBUV/2 ozone data, covering a period of more than a solar cycle (about 15 years), are used to study the UV response of ozone in the stratosphere. The study shows that about 2% change in total column ozone and about 5-7% change in ozone mixing ratio in the upper stratosphere (0.7 to 2 hPa) may be attributed to the change in the solar UV flux over a solar cycle. In the upper stratosphere, where photochemical processes are expected to play a major role, the measured solar cycle variation of ozone is significantly larger than inferred either from the photochemical models or from the ozone response to the 27-day solar UV modulation. For example, the observed solar cycle related change in ozone mixing ratio at 2 hPa is about 1% for 1% change in the solar UV flux near 200 nm. The inferred change in ozone from either the photochemical models or from the 27-day ozone-UV response is about a factor of 2-3 lower than this value.

Description: Vertical profiles of N2O and NO(y) taken by the ER-2 outside the vortex are used to construct average vertical profiles of F(NO(y)) = NO(y)/(A-N2O), where A is the tropospheric content of N2O three years prior to the measurements. The southern hemisphere had less nitrous oxide in the range 400 less than Theta less than 470 K, by up to 25% relative to the northern hemisphere. F(NO(y)) is the ratio of NOy produced to N2O lost in a stratospheric air mass since entry from the troposphere. The profiles of F(NO(y)) have the following characteristics: (1) Relative to 1991-1992, a year without denitrification inside or outside the vortex, the northern hemisphere in 1988-1989 showed denitrification outside the vortex ranging up to 25% and averaging 17% above Theta = 425 K. (2) Relative to the northern hemisphere in 1991-1992, the southern hemisphere in 1987 showed denitrification outside the vortex ranging up to 32% and averaging 20% above Theta = 400 K. (3) Below Theta = 400 K the southern hemisphere showed enhancements of F(NO(y)) relative to the northern hemisphere in 1991-1992 ranging up to 200% at Theta = 375 K, outside the vortex. Corresponding profiles of residual water, R(H2O) = H2O - 2(1.6 - CH4), are considered and shown to be consistent with those of F(NO(y)) in the sense that they show deficits outside the Antarctic vortex, which was both dehydrated and denitrified, but not outside the 1988-1989 Arctic vortex, which was denitrified but not dehydrated. R(H2O) is the water content of stratospheric air with the contribution from methane oxidation subtracted. Comparison of F(NO(y)) and R(H2O) below 400 K outside the Antarctic vortex leads to the suggetion that dehydration in the Antarctic vortex occurs by the sedimentation of ice crystals large enough to fall out of the stratosphere, whereas denitrification occurs mainly on mixed nitric acid-water crystals which evaporate below the base of the vortex at Theta = 400 K but above the tropopause.

Description: During January and August 1985, the scanning radiometers of the Earth Radiation Budget Experiment(ERBE) aboard the Earth Radiation Budget Satellite (ERBS) and the NOAA-9 satellite were operated in along-track scanning modes. Along-track scanning permits the study of many measurement problems. It provides the data for developing a limb-darkening model for a single site over a short period of time and also permits the indentification of the scene from data taken at smaller nadir angles. The earth-emitted radiation measured by the scanners has been analyzed to produce limb-darkening models for a variety of scene types. Limb-darkening models relate the radiance in any given direction to the radiant flux. The scene types were computed using measurements within 10 deg of zenith. The models have values near zenith of 1.02-1.09. The typical zenith values of the model are 1.06 for both day and night for ERBS, and for NOAA-9, 1.06 for day and 1.05 for night. Mean models are formed for the ERBS and NOAA-9 results and are found to differ less than 1%, the ERBS results being the higher. The models vary about 1% with latitude near zenith and agree with earlier models that were used to analyze ERBE data typically to 2%.

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

Description: One of the oldest mysteries in geomagnetism is the linkage between solar and geomagnetic activity. In investigating the causes of geomagnetic storms occurring during solar maximum, the following topics are discussed: solar phenomena; types of solar wind; magnetic reconnection and magnetic storms; an interplanetary example; and future space physics missions.

Description: Atmospheric mass loading produces a primarily vertical displacement of the Earth's crust. This displacement is correlated with surface pressure and is large enough to be detected by very long baseline interferometry (VLBI) measurements. Using the measured surface pressure at VLBI stations, we have estimated the atmospheric loading term for each station location directly from VLBI data acquired from 1979 to 1992. Our estimates of the vertical sensitivity to change in pressure range from 0 to -0.6 mm/mbar depending on the station. These estimates agree with inverted barometer model calculations (Manabe et al., 1991; vanDam and Herring, 1994) of the vertical displacement sensitivity computed by convolving actual pressure distributions with loading Green's functions. The pressure sensitivity tends to be smaller for stations near the coast, which is consistent with the inverted barometer hypothesis. Applying this estimated pressure loading correction in standard VLBI geodetic analysis improves the repeatability of estimated lengths of 25 out of 37 baselines that were measured at least 50 times. In a root-sum-square (rss) sense, the improvement generally increases with baseline length at a rate of about 0.3 to 0.6 ppb depending on whether the baseline stations are close to the coast. For the 5998-km baseline from Westford, Massachusetts, to Wettzell, Germany, the rss improvement is about 3.6 mm out of 11.0 mm. The average rss reduction of the vertical scatter for inland stations ranges from 2.7 to 5.4 mm.

Description: A theoretical model is used to describe and investigate the effects of simultaneous crystallization, radiation loss, and entrainment of cooler material on the temperature of a well-mixed core of an active aa lava flow. Entrainment of crust, levee debris, and base material into the interior of active flows has been observed, but the degree of assimilation and the thermal consequences are difficult to quantify. The rate of entrainment can be constrained by supplementing the theoretical model with information on the crystallization along the path of the flow and estimation of the radiative loss from the flow interior. Application of the model is demonstrated with the 1984 Mauna Loa flow, which was erupted about 30 C undercooled. Without any entrainment of cooler material, the high crystallization rates would have driven temperatures in the core wall above temperatures measured by thermocouple and estimated from glass geothermometry. One plausible scenario for this flow, which agrees with available temperature and crystallinity measurements, has a high initial rate of entrainment during the first 8 hours of travel (a mass ratio of entrained material to fluid core of about 15% if the average temperature of the entrained material was 600 C), which counterbalances the latent heat from approximately 40% crystallization. In this scenario, the model suggests an additional 5% crystallization and a 5% entrainment mass ratio over the subsequent 16-hour period. Measurements of crystallization, radiative losses, and entrainment factors are necessary for understanding the detailed thermal histories of active lava flows.

Description: The purpose of this Note is to present results from an analytic/experimental study that investigated the potential for passively changing blade twist through the use of extension-twist coupling. A set of composite model rotor blades was manufactured from existing blade molds for a low-twist metal helicopter rotor blade, with a view toward establishing a preliminary proof concept for extension-twist-coupled rotor blades. Data were obtained in hover for both a ballasted and unballasted blade configuration in sea-level atmospheric conditions. Test data were compared with results obtained from a geometrically nonlinear analysis of a detailed finite element model of the rotor blade developed in MSC/NASTRAN.

Description: The data from the infrared telescope (IRT), which was flown on space shuttle Challenger Spacelab 2 mission (July 1985), were originally reported by Koch et al. (1987) as originating from near orbital emissions, primarily H2O. In this study, analysis of this data was extended to determine the collisional cross sections for the excitation of the low lying vibrational levels of H2O, present in the orbiter cloud, by atmospheric O(3P). The evaluation of the contribution to the measured signal from solar excitation and ram O excitation of outgassing H2O permits the determination of the H2O column density and the excitation cross section of the (101) level at an O(3P) velocity of approximately 7.75 km/s. Contributions to the radiation in the 1.7-3.0 micron band by transitions from the (100), (001), and multiquantum excited levels are discussed. The findings of the study are (1) the IRT data for the 4.5-9.5 micron and the nighttime data for the 1.7-3.0 micron sensors are consistent with being explained by collision excitation of H2O by O(3P), (2) diurnal variations of 4.5-9.5 micron intensities follow the model predicted O density for a full orbit, (3) daytime increases in the H2O cloud density were not evident, (4) the cross sections for the collisional excitation process are derived and compared to values computated by Johnson (1986) and Redmon et al. (1986), (5) theoretical investigation suggests greater than 60% of the radiation from H2O is a result of multiphoton emission resulting from collisional multiquanta excitation, and (6) the large daytime increase in the 1.7-3.0 micron intensity data suggests that O(+) may likely be instrumental in producing excited H2O(+) through charge exchange.

Description: The paper considers the compressible Rayleigh equation as a model for the Mach wave emission mechanism associated with high-temperature supersonic jets. Solutions to the compressible Rayleigh equation reveal the existence of several families of supersonically convecting instability waves. These waves directly radiate noise to the jet far field. The predicted noise characteristics are compared to previously acquired experimental data for an axisymmetric Mach 2 fully pressure balanced jet operating over a range of jet total temperatures from ambient to 1370 K. The results of this comparison show that the first-order supersonic instability wave and the Kelvin-Hemlhlotz first-, second-, and third-order modes have directional radiation characteristics that are in agreement with observed data. The assumption of equal initial amplitudes for all of the waves leads to the conclusion that the flapping mode of instability dominates the noise radiatio process of supersonic jets. At a jet temperature of 1370 K, supersonic instability waves are predicted to dominate the noise radiated at high frequency at narrow angles to the jet axis.

Description: The objective of the present work is to study the mixing characteristics of a linear array of supersonic rectangular jets under conditions of screech synchronization. The screech synchronization at a fully expanded jet Mach number of 1.61 is achieved by a precise adjustment of the internozzle spacing. To our knowledge, such an experiment on the resonant mixing of screech synchronized multiple rectangular jets has not been reported before. The results are compared with the case where the screech was suppressed in the multijet configuration.

Description: The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

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

Description: We have combined the most recent Pioneer Venus Orbiter (PVO) and Magellan (MGN) data with the earlier 1978-1982 PVO data set to obtain a new 60th degree and order spherical harmonic gravity model and a 120th degree and order spherical harmonic topography model. Free-air gravity maps are shown over regions where the most marked improvement has been obtained (Ishtar-Terra, Alpha, Bell and Artemis). Gravity versus topography relationships are presented as correlations per degree and axes orientation.

Description: Observations of stratospheric column amounts of nitrogen dioxide (NO2), nitric oxide (NO) and nitric acid (HNO3) have been made following major eruptions of the El Chichon and Mt. Pintatubo volcanoes. Midlatitude abundances of NO2 and NO were reduced by as much as 70% in the months following the appearance of the volcanic aerosols as compared to volcanically quite periods. There are heterogeneous reactions which could occur on the volcanic aerosols to convert NO2 into HNO3 but no commensurate increase in HNO3 column amounts was observed at the times of NO2 decrease.

Description: The benefits of using a hypersonic waverider for spacecraft trajectory modification are presented. A waverider is a hypersonic vehicle specifically designed so that the undersurface bow shock is attached to the leading edge, which provides for the highest known lift-to-drag ratios achievable at high Mach number flight. Several viable space missions are suggested which could use such configurations for low-drag aero-assisted maneuvers in planetary atmospheres. It is shown that large changes in the spacecraft velocity vector can be accomplished with acceptably small losses in energy due to drag using a waverider aeroshell. The primary advantage of an aero-assist maneuver is suggested by comparison to a traditional gravity-assist trajectory. Some scaling laws are presented for comparing waveriders designed for different planetary atmospheres, and it is shown that the compositional differences between the terrestrial planets has a minimal impact on waverider design.

Description: The devolatilization of calcium sulfate, which is present in the target rock of the Chicxulub, Mexico impact structure, and dispersal in the stratosphere of the resultant sulfuric acid aerosol have been suggested as a possible mechanism for the Cretaceous-Tertiary extinctions. We measured the amount of SO2 produced from two shock-induced devolatilization reactions of calcium sulfate up to 42 GPa in the laboratory. We found both to proceed to a much lower extent than calculated by equilibrium thermodynamic calculations. Reaction products are found to be approx. 10(exp -2) times those calculated for equilibrium. Upon modeling the quantity of sulfur oxides degassed into the atmosphere from shock devolatilization of CaSO4 in the Chicxulub lithographic section, the resulting 9 x 10(exp 16) to 6 x 10(exp 17) g (in sulfur mass) is lower by a factor of 10-100 than previous upper limit estimates, the related environmental stress arising from the resultant global cooling and fallout of acid rain is insufficient to explain the widespread K-T extinctions.

Description: Energetic atomic (O(+1) and N(+1)) and molecular (O2(+1), NO(+1), and N2(+1)) ions of ionospheric origin were observed in Earth's magnetotail at X approximately -146 R(sub E) during two plasma sheet sunward/tailward flow-reversal events measured by instruments on the GEOTAIL spacecraft. These events were associated with concurrent ground-measured geomagnetic disturbance intensification at auroral-and mid-latitudes (Kp = 7(-)). Energetic ions in the sunward-component and tailward flows were from both the solar wind and ionosphere. Plasma and energetic ions participated in the flows. During tailward flow, ionospheric origin ion abundance ratios at approximately 200-900 km/s in the rest frame were N(+1)/O(+1) = approximately 25-30% and ((O2(+1), NO(+1), and N2(+1))/O(+1) = approximately 1-2%. We argue that tailward flow most likely initiated approximately 80-100 R(sub E) tailward of Earth and molecular ions were in the plasma sheet prior to geomagnetic intensification onset.

Description: Wave forms of BEN (Broadband Electrostatic Noise) in the geomagnetic tail were first detected by the Wave Form Capture reciever on the GEOTAIL spacecraft. The results show that most of the BEN in the plasma sheet boundary layer (PSBL) are not continuous broadband noise but are composed of a series of solitary pulses having a special form which we term 'Electrostatic Solitary Waves (ESW)'. A nonlinear BGK potential model is proposed as the generation mechanism for the ESW based upon a simple particle simulation which considers the highly nonlinear evolution of the electron beam instability. The wave forms produced by this simulation are very similar to those observed by GEOTAIL and suggest that the nonlinear dynamics of the electron beam play an essential role in the generation of ESW.

Description: Radiometric (K-Ar and Ar-40/Ar-39) age determinations of volcanic and volcaniclastic rocks, combined with structural, gravity, and seismic reflection data, are used to constrain the age of sedimentary strata contained within the seismically and volcanically active northern Malawi (Nyasa) rift and to characterize changes in basin and flank morphologies with time. Faulting and volcanism within the Tukuyu-Karonga basin began at approximately 8.6 Ma, when sediments were deposited in abroad, initially asymmetric lake basin bounded on its northeastern side by a border fault system with minor topographic relief. Extensions, primarily by a slip along the border fault, and subsequent regional isostatic compensation led to the development of a 5-km-deep basin bounded by broad uplifted flanks. Along the low-relief basin margin opposite border fault, younger stratigraphic sequences commonly onlap older wedge-shaped sequences, although their internal geometry is often progradational. Intrabasinal faulting, flankuplift, and basaltic and felsic volcanism from centers at the northern end of the basin became more important at about 2.5 Ma when cross-rift transfer faults developed to link the Tukuyu-Karonga basin to the Rukwa basin. Local uplift and volcanic construction at the northern end of the basin led to a southeastward shift in the basin's depocenter. Sequence boundaries are commonly erosional along this low-relief (hanging wall) margin and conformable in the deep lake basin. The geometry of stratigraphic sequences and the distribution of the erosion indicate that horizontal and vertical crustal movements both across and along the length of the rift basin led to changes in levels of the lake, irrespective of paleoclimatic fluctuations.

Description: We report infrared heterodyne spectroscopy (lambda/delta lambda is approximately 10(exp 6)) of C2H6 emission at 11.9 microns from the northern Jovian auroral region, in observations conducted over December 2-7, 1989. Accurately measured line shapes provide information on C2H6 abundance as well as temperature and permit retrieval of the source pressure region. Enhanced emission was observed in the longitude range approximately 150-180 deg at approximately 60 deg north latitude, approximately corresponding to the CH4 7.8 micron hot spot and the region of brightest UV aurora. Significant brightness variations were observed in the hot spot emissions on a time scale of approximately 20 hours. Analysis of the brightest hot spot spectra indicates C2H6 mole fractions of approximately (6.3-6.8) x 10(exp -6) at temperatures of approximately 182-184 K at 1 mbar, compared to mole fractions of (3.8 +/- 1.4) x 10(exp -6) averaged over spectra outside the hot spot at a temperature of approximately 172 K at the same pressure. Fixing the mole fraction to the lower limit retrieved in the quiescent (non-hot spot) region allows the temperature at 1 mbar to be as high as approximately 200 K within the hot spot. These results provide upper limits to the temperature increase near the source of the C2H6 thermal infrared emission. Combined with results from similar measurements of ethylene emission probing the approximately 10-microbar region (Kostiuk et al., this issue), altitude information on the thermal structure of the Jovian auroral stratosphere can be obtained for the first time.

Description: We report on the observations of a number of quasi-dc electric field events associated with large-scale atmospheric weather formations. The observations were made by the electric field experiment onboard the San Marco D satellite, operational in an equatorial orbit from May to December 1988. Several theoretical studies suggest that electric fields generated by thunderstorms are present at high altitudes in the ionosphere. In spite of such favorable predictions, weather-related events are not often observed since they are relatively weak. We shall report here on a set of likely E field candidates for atmospheric-ionospheric causality, these being observed over the Indonesian Basin, northern South America, and the west coast of Africa; all known sites of atmospheric activity. As we shall demonstrate, individual events often be traced to specific active weather features. For example, a number of events were associated with spacecraft passages near Hurricane Joan in mid-October 1988. As a statistical set, the events appear to coincide with the most active regions of atmospheric weather.

Description: Doppler tracking data of three orbiting spacecraft have been reanalyzed to develop a new gravitational field model for the planet Mars, Goddard Mars Model 1 (GMM-1). This model employs nearly all available data, consisting of approximately 1100 days of S band tracking data collected by NASA's Deep Space Network from the Mariner 9 and Viking 1 and Viking 2 spacecraft, in seven different orbits, between 1971 and 1979. GMM-1 is complete to spherical harmonic degree and order 50, which corresponds to a half-wavelength spatial resolution of 200-300 km where the data permit. GMM-1 represents satellite orbits with considerably better accuracy than previous Mars gravity models and shows greater resolution of identifiable geological structures. The notable improvement in GMM-1 over previous models is a consequence of several factors: improved computational capabilities, the use of otpimum weighting and least squares collocation solution techniques which stabilized the behavior of the solution at high degree and order, and the use of longer satellite arcs than employed in previous solutions that were made possible by improved force and measurement models. The inclusion of X band tracking data from the 379-km altitude, nnear-polar orbiting Mars Observer spacecraft should provide a significant improvement over GMM-1, particularly at high latitudes where current data poorly resolve the gravitational signature of the planet.

Description: A series of surface-modified clays containing nanophase (np) iron/oxyhydroxides of extremely small particle sizes, with total iron contents as high as found in Mars soil, were prepared by iron deposition on the clay surface from ferrous chloride solution. Comprehensive studies of the iron mineralogy in these 'Mars-soil analogs' were conducted using chemical extractions, solubility analyses, pH and redox, x ray and electron diffractometry, electron microscopic imaging specific surface area and particle size determinations, differential thermal analyses, magnetic properties characterization, spectral reflectance, and Viking biology simulation experiments. The clay matrix and the procedure used for synthesis produced nanophase iron oxides containing a certain proportion of divalent iron, which slowly converts to more stable, fully oxidized iron minerals. The noncrystalline nature of the iron compounds precipitated on the surface of the clay was verified by their complete extractability in oxalate. Lepidocrocite (gamma-FeOOH) was detected by selected area electron diffraction. It is formed from a double iron Fe(II)/Fe(III) hydroxyl mineral such as 'green rust', or ferrosic hydroxide. Magnetic measurements suggested that lepidocrocite converted to the more stable meaghemite (gamma-Fe203) by mild heat treatment and then to nanophase hematite (aplha-Fe203) by extensive heat treatment. Their chemical reactivity offers a plausible mechanism for the somewhat puzzling observations of the Viking biology experiments. Their unique chemical reactivities are attributed to the combined catalytic effects of the iron oxide/oxyhydroxide and silicate phase surfaces. The mode of formation of these (nanophase) iron oxides on Mars is still unknown.

Description: A radiation model, together with National Meteorological Center temperature observations, was used to compute daily net heating rates in the northern hemisphere (NH) for the Arctic late fall and winter periods of both 1988-1989 and 1991-1992 and in the southern hemisphere (SH) for the Antarctic fall and winters of 1987 and 1992. The heating rates were interpolated to potential temperature (theta) surfaces between 400 K and 2000 K and averaged within the polar vortex, the boundary of which was determined by the maximum gradient in potential vorticity. The averaged heating rates were used in a one-dimensional vortex interior descent model to compute the change in potential temperature with time of air parcels initialized at various theta values, as well as to compute the descent in log pressure coodinates. In the NH vortex, air parcels which were initialized at 18 km on November 1, descended about 6 km by March 21, while air initially at 25 km descended 9 km in the same time period. this represents an average descent rate in the lower stratosphere of 1.3 to 2 km per month. Air initialized at 50 km descended 27 km between November 1 and March 21. In the SH vortex, parcels initialized at 18 km on March 1, descended 3 km, while air at 25 km descended 5-7 km by the end of October. This is equivalent to an average descent in the lower stratosphere of 0.4 to 0.9 km per month during this 8-month period. Air initialized at 52 km descended 26-29 km between March 1 and October 31. In both the NH and the SH, computed descent rates increased markedly with height. The descent for the NH winter of 1992-1993 and the SH winter of 1992 computed with a three-dimensional trajectory model using the same radiation code was within 1 to 2 km of that calculated by the one-dimensional model, thus validating the vortex averaging procedure. The computed descent rates generally agree well with observations of long-lived tracers, thus validating the radiative transfer model.

Description: In this paper we report radiometric measurements of tropospheric brightness temperatures obtained during the AASE 2 experiment. These measurements represent the first attempt to characterize effective radiative temperatures as seen from above the troposphere during the Arctic winter. The reported measurements include brightness temperatures at 6.7 and 10.5 microns as seen from the NASA DC-8 aircraft flying at about 11 km altitude. We also present radiative transfer calculations to estimate the effect of tropospheric brightness temperature on the lower stratospheric heating rates. Because of the recent massive eruption of the Pinatubo volcano, we also discuss the effects of a volcanic aerosol layer. It is concluded that small particles like the volcanic aerosol or polar stratospheric clouds (PSCs) type 1 do not affect stratospheric heating rates by much; on the other hand, larger particles, PSCs types 2 and 3, may have significant effects on heating rates and consequently on dynamics of the lower stratosphere. The dynamical effects of local stratospheric temperature variations are briefly discussed.

Description: We report tropospheric (altitudes greater than 5 km) observations of CO2, CO, CH4, and light hydrocarbons (C2-C4) over the latitude range from 90 deg N to 23 deg S recorded onboard the NASA DC-8 aircraft during the winter 1992 Second Airborne Arctic Stratospheric Expedition (AASE-2). Mixing ratios for these species exhibited significant north-south gradients with maximum values in subpolar and arctic regions and minima over the southern tropics. At latitudes greater than 40 deg N, the mixing ratios of most species increased significantly over the course of the 3-month measurement period. Also at high northern latitudes, the variations of all relatively long-lived reactive carbon species were linearly correlated with fluctuations of CO2 with CO, CH4, C2H6, C2H2, C3H8, and n-C4H10 exhibiting average enhancement ratios in terms of ppbv(X)/ppmb(CO2) of 13.8, 8.4, 0.21, 0.075, 0.085, and 0.037, respectively.

Description: Vertical column measurements of the gaseous composition of the tropical stratosphere were made from the NASA DC-8 aircraft early in 1992. As anticipated, the burdens of the stratospheric source gases (e.g., O3, HF, HCl, ClNO3, HNO3) were reduced from their mid-latitude values due to increased uplift and photolysis. The tracers revealed considerably more uplift near the equator than the sub-tropics. For example, the HF burdens at +/- 20 deg latitude were nearly double those at 5 deg N. This, together with results obtained from other long-lived gases (e.g. N2O, CH4, CF2Cl2) indicates that volume mixing ratios found at 22 km altitude at mid-latitudes occurred at 26 km in the sub-tropics and at 30 km in the equatorial zone. This zone of uplift was symmetrical about the equator even though the sun was overhead at 20 deg S.

Description: We welcome this opportunity which Drs. Ricard, Sabadini and Spada (RSS) have offered us to elaborate on some points concerning the connection between convective circulations and rotational dynamics. The relationship of the linear viscoelastic normal-mode theory to other approaches, such as viscous and the non-linear viscoelastic treatments, can be found in Moser et al. (1993). We will address the issues raised by RSS by first going over the reference system. This will be followed by a discussion concerning long-term polar wander and relative angular momentum. We will also go over the issues concerning the off-diagonal terms in the total moment of inertia and polar wander.

Description: In late July and early August of 1991, a major suborbital scientific campaign (NLC-91) involving scientists from eight countries was conducted as ESRANGE, Kiruna, Sweden and at Heiss Island, Russia. The purpose of the program was to investigate the chemical, dynamical, and electrodynamical properties of the polar summer mesosphere. Thirty one rocket flights were coordinated with two coherent radar facilities, EISCAT and CUPRI, and with other ground-based observatories and facilities. This permitted direct comparison between the in situ measurements and those obtained by remote sensing of the mesosphere via continuous ground-based monitoring. The primary objectives of the campaign were to study noctilucent clouds (NLCs) and polar mesospheric summer echoes (PMSEs), including their possible relationship to local aerosols and/or small scale turbulence. This overview describes the scientific program, discusses the geophysical conditions during launch activities, and reviews some of the preliminary results. More detailed results can be found in the papers which follow.

Description: The International Sun-Earth Explorer 3 (ISEE-3) magnetic field and plasma electron data from Jan - March 1983 have been searched to study thin current sheets in the deep tail region. 33 events were selected where the spacecraft crossed through the current sheet from lobe to lobe within 15 minutes. The average thickness of the observed current sheets was 2.45 R(sub E), and in 24 cases the current sheet was thinner than 3.0 R(sub E); 6 very thin current sheets (thickness lambda less than 0.5 R(sub E) were found. The electron data show that the very thin current sheets are associated with considerable temperature anisotropy. On average, the electron gradient current was about 17% of the total current, whereas the current arising from the electron temperature anisotropy varied between 8-45% of the total current determined from the lobe field magnitude.

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

Description: The 3-5 micrometer thermal emission of the nightside of Venus, recorded by the Near-Infrared Mapping Spectrometer (NIMS) instrument at the time of the Galileo flyby of Venus, is analysed to infer the properties of the upper cloud boundary. From the global maps of Venus at fixed wavelengths, the limb darkening of the flux is measured at several latitudes, within each infrared channel. By using the nominal Pioneer Venus thermal profile, these data give access to two parameters: the cloud deck temperature and the cloud scale height. It is verified independently, from the NIMS spectra, that this thermal profile is consistent with all the NIMS observations, and that the thermal structure does not vary significantly in the latitude range (25 deg S, 30 deg N). Within this range, the cloud scale height is found to be constant with latitude, and is H = 5.2 km, with an accuracy of about 15%, taking into account the various sources of theoretical and observational uncertainties. At higher latitudes, the temperature profile becomes more isothermal and the presented method to retrieve H is no longer valid.

Description: A large number of i.r. spectra of Venus was obtained using the Near-Infrared Mapping Spectrometer (NIMS) on the Galileo spacecraft, during the February 1990 encounter. Preliminary results show an apparent increase in the tropospheric CO volume mixing ratio (vmr) in the northern polar region. Other possible explanations of the observations are examined and rejected and an increase of the CO abundance north of 47 deg N of (35 +/- 15)% is inferred. Some possible causes of this enhancement are suggested.

Description: The spectroscopic data of the Near-Infrared Mapping Spectrometer (NIMS), recorded during the Galileo flyby of Venus, are analysed to retrieve the water vapor abundance variations in the lower atmosphere of Venus at night. The 1.18 micrometer spectral window, which probes altitude levels below 20 km, is used for this purpose. Constraints on the CO2 continuum and far-wing opacity from existing ground-based high-resolution observations are included in the modelling of the NIMS spectra. The NIMS measurements can be fitted with a water vapor mixing ratio of 30 +/- 15 ppm, in agreement with analyses of ground-based nightside observations. The water vapor abundance shows no horizontal variations exceeding 20% over a wide latitude range (40 deg S, 50 deg N) on the nightside of Venus. Within the same selection of NIMS spectra, a large enhancement in the O2 fluorescence emission at 1.27 micrometer is observed at a latitude of 40 deg S, over a spatial area about 100 km wide.

Description: Using Venus nightside data obtained by the Galileo Near-Infrared Mapping Spectrometer (NIMS), we have studied the correlation of 1.74 and 2.30 micrometer radiation which is transmitted through the clouds. Since the scattering and absorption properties of the cloud particles are different at these two wavelengths, one can distinguish between abundance variations and variations in the properties of the cloud particles themselves. The correlation of intensities shows a clustering of data into five distinct branches. Using radiative transfer calculations, we interpret these branches as regions of distinct but different mixes of Mode 2' and 3 particles. The data and calculations indicate large differences in these modal ratios, the active cloud regions varying in content from nearly pure Mode 2' particles to almost wholly Mode 3. The spatial distribution of these branches shows large scale sizes and both hemispheric symmetries and asymmetries. High-latitude concentrations of large particles are seen in both hemispheres and there is banded structure of small particles seen in both the North and South which may be related. The mean particle size in the Northern Hemisphere is greater than found in the South. If these different branch regions are due to mixing of vertically stratified source regions (e.g. photochemical and condensation source mechanisms), then the mixing must be coherent over very large spatial scales.

Description: The Radio and Plasma Wave Science (RPWS) experiment being built for the Cassini spacecraft will study a wide range of plasma and radio wave phenomena in the magnetosphere of Saturn and will also make valuable measurements during the cruise phase and at other encounters. A feature of data from wave receivers is the capability of producing vastly more data than the spacecraft telemetry link is capable of transmitting back to the Earth. Thus, techniques of on-board data compression and data reduction are important. The RPWS instrument has one processor dedicated to data compression tasks.

Description: The map of the coseismic displacement field generated by interferometric processing of synthetic aperture radar (SAR) images taken before and after the June 28, 1992, Landers earthquake sequence brings new insights into the nature of deformation caused by these earthquakes. We use the interferometric map generated by Massonnet et al. (1993) to analyze the surface displacement field in the vicinity of the fault trace. Complexities in the fringe pattern near the fault reflect short-wavelength variations of the surface rupture and slip distribution, and attest to large displacement gradients. Along two sections of the fault, characteristic fringe patterns can be recognized, contrasting in density and direction with patterns observed away from the rupture. In order to understand the observed fringe patterns, we compute synthetic interferograms in three simple cases: (1) rigid-body rotations about a vertical axis, (2) about a horizontal axis (tilt), and (3) distributed, simple shear. The orientation and spatial separation of interferometric fringes predicted by these models help constrain near-field deformation and rupture parameters. Where the Kickapoo fault connects with the Homestead Valley fault, the interferogram shows a clear pattern of parallel N20 deg W fringes separated by about 160 m. This pattern and vertical offsets measured along the Kickapoo fault suggest that the block between this fault and the Johnson Valley fault may have been tilted, down to the west. A 5-km block lifted by 1 m on one side would be tilted by an angle of 0.01 deg (190 microrad), producing fringes separated by about 160 m, parallel to the tilt axis. Such a tilt, parallel to a N20 deg W direction, would account for the gradual, northward increase of the vertical slip component observed along the Kickapoo fault. This tilt may also explain the 1 m of reverse slip observed along the 'slip gap' section of the Homestead Valley break. Between the southern end of the Johnson Valley fault and the Eureka Peak fault, where no surface rupture has been mapped, the dense pattern of fringes implies distributed shear, probably resulting from fault slip at depth. The density and direction of the fringes in the gap are consistent with a right-lateral slip of 1.2-3.8 m on a blind fault locked above the depth of 1.5-2 km. Such observations of small wavelength features in the SAR interferogram bring new insights into the near-field displacement gradient and thus on response of the uppermost crust to seismic rupture.

Description: Global distributions of surface and atmospheric cloud radiative forcing parameters have been derived using parameterized radiation models with satellite meteorological data from the International Satellite Cloud Climatology Project, and directly measured top-of-atmosphere radiative fluxes from the Earth Radiation Budget Experiment. Specifically, shortwave, longwave, and total cloud forcing at the surface, and column-averaged values of longwave cloud forcing of the atmosphere were derived for the midseasonal months of April, July, and October 1985 and January 1986, covering a complete annual cycle. Seasonal variability is illustrated by comparing the results for July 1985 and January 1986, which represent the seasonal extremes. Surface shortwave cloud forcing is always negative, representing a cooling of the surface, with strongest cooling (-120 to -180 W/sq m) occurring over midlatitude storm tracks of the summer hemisphere. Surface longwave cloud forcing is always positive, representing a warming of the surface, with strongest warming (60 to 75 W/sq m) occurring over storm tracks of the winter hemisphere. Zonal averages show the entire summer hemisphere dominated by shortwave cooling, the middle and high latitudes of the winter hemisphere dominated by longwave warming, and a broad zone of transition in between. The globally averaged total cloud forcing amounts to a cooling throughout the year, ranging from a low of about -12 W/sq m for July 1985 to a high of about -25 W/sq m for January 1986. The longwave cloud forcing of the atmosphere shows a strong warming over deep convective regions in the tropics and a moderate cooling outside the tropics, amounting to a weak cooling (-2 to -5 W/sq m) in the global average. Comparisons of the results with general circulation model simulations show broad qualitative agreement regarding the locations of prominent warming and cooling regions. Quantitative comparisons, on the other hand, show significant differences between the magnitudes of warming and cooling in these regions. Most of the larger differences can be attributed to known deficiences of the general circulation model simulations. Comparisons with satellite-derived results of other investigators show much better agreement.

Description: The acceleration of protons in a dynamically evolving magnetotail is investigated by tracing particles in the fields obtained from a three-dimensional resistive magnetohydrodynamic (MHD) simulation. The MHD simulation, representing plasmoid formation and ejection through a near-Earth reconnection process, leads to cross-tail electric fields of up to approximately 4 mV/m with integrated voltages across the tail of up to approximately 200 kV. Energization of particles takes place over a wide range along the tail, due to the large spatial extent of the increased electric field together with the finite cross-tail extent of the electric field region. Such accelerated particles appear earthward of the neutral line over a significant portion of the closed field line region inside of the separatrix, not just in the vicinity of the separatrix. Two different acceleration processes are identified: a 'quasi-potential' acceleration, due to particle motion in the direction of the cross-tail electric field, and a 'quasi-betatron' effect, which consists of multiple energy gains from repeated crossings of the acceleration region, mostly on Speiser-type orbits, in the spatially varying induced electric field. The major source region for accelerated particles in the hundreds of keV range is the central plasma sheet at the dawn flank outside the reconnection site. Since this source plasma is already hot and dense, its moderate energization by a factor of approximately 2 may be sufficient to explain the observed increases in the energetic particle fluxes. Particles from the tail are the source of beams at the plasma sheet/lobe boundary. The temporal increase in the energetic particle fluxes, estimated from the increase in energy gain, occurs on a fast timescale of a few minutes, coincident with a strong increase in B(sub z), despite the fact that the inner boundary ('injection boundary') of the distribution of energized particles is fairly smooth.

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

Description: A high-precision radiometric satellite tracking system, Doppler Orbitography and Radio-positioning Integrated by Satellite system (DORIS), has recently been developed by the French space agency, Centre National d'Etudes Spatiales (CNES). DORIS was designed to provide tracking support for missions such as the joint United States/French TOPEX/Poseidon. As part of the flight testing process, a DORIS package was flown on the French SPOT 2 satellite. A substantial quantity of geodetic quality tracking data was obtained on SPOT 2 from an extensive international DORIS tracking network. These data were analyzed to assess their accuracy and to evaluate the gravitational modeling enhancements provided by these data in combination with the Goddard Earth Model-T3 (GEM-T3) gravitational model. These observations have noise levels of 0.4 to 0.5 mm/s, with few residual systematic effects. Although the SPOT 2 satellite experiences high atmospheric drag forces, the precision and global coverage of the DORIS tracking data have enabled more extensive orbit parameterization to mitigate these effects. As a result, the SPOT 2 orbital errors have been reduced to an estimated radial accuracy in the 10-20 cm RMS range. The addition of these data, which encompass many regions heretofore lacking in precision satellite tracking, has significantly improved GEM-T3 and allowed greatly improved orbit accuracies for Sun-synchronous satellites like SPOT 2 (such as ERS 1 and EOS). Comparison of the ensuing gravity model with other contemporary fields (GRIM-4C2, TEG2B, and OSU91A) provides a means to assess the current state of knowledge of the Earth's gravity field. Thus, the DORIS experiment on SPOT 2 has provided a strong basis for evaluating this new orbit tracking technology and has demonstrated the important contribution of the DORIS network to the success of the TOPEX/Poseidon mission.

Description: The presence of anisotropic plasma distributions, trapped at the Earth's magnetic equator, has consequences for the electric field structure which must exist in equilibrium along the magnetic field line. Data from SCATHA and Dynamics Explorer 1 indicate that the core ion distributions at the magnetic equator can be well described as bi-Maxwellian distributions, with a perpendicular temperature an order of magnitude larger than the parallel temperature. A collisionless model is developed for the variation in plasma parameters, following the forms developed by Whipple (1977). If the core electron anisotropy is low, the resulting electric field of approximately 0.1 microV/m is pointed away from the equator. Under these conditions the self-consistent electric field will not overcome the effects of magnetic trapping. The resulting potential distribution results in a local maximum in total plasma density at the equator. Only when the electron distribution is primarily field-aligned can there be a density minimum at the equator. Comparisons are made between this model and the observed variations in DE 1 plasma parameters with latitude.

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

Description: A combined 3-dimensional circulation model and aerosol microphysical/transport model is used to simulate the dispersion of the Mt. Pinatubo volcanic cloud in the stratosphere for the first few months following the eruption. Radiative heating of the cloud due to upwelling infrared radiation from the troposphere is shown to be an important factor affecting the transport. Without cloud heating, cloud material stays mostly north of the equator, whereas with cloud heating, the cloud is transported southward across the equator within the first two weeks following the eruption. Generally the simulations agree with Total Ozone Mapping Spectrometer (TOMS), Advanced Very High Resolution Radiometer (AVHRR), and Stratospheric Aerosol and Gas Experiment (SAGE) satellite observations showing the latitude distribution of cloud material to be between about 20 deg S and 30 deg N within the first few months. Temperature perturbations in the stratosphere induced by the aerosol heating are generally 1-4 K, in the range of those observed.

Description: We use geodetic data from Very Long Baseline Interferometry (VLBI) to determine the pre- and postseismic velocities of two sites. We then place limits on variations in interseismic strain buildup. The 1987 and 1988 Gulf of Alaska earthquakes (each Ms = 7.6) broke the Pacific plate interior. During the earthquakes the Cape Yakataga site moved 78 mm toward southwest. During the 1989 Loma Prieta earthquake (Ms = 7.1) the Fort Ord site moved 48 mm toward north. Baselines (a) from Fairbanks to Cape Yakataga and (b) from Mojave to Fort Ord change at nearly the same rate before and after the earthquakes. Postseismic transients, which we determine from differences between post- and preseismic rates, are minor: at Cape Yakataga the transient is 3 +/- 4 mm in a postseismic interval of 23 months, and at Fort Ord the transient is 6 +/- 5 mm in 21 months. The slip beneath the Loma Prieta rupture needed to generate the Fort Ord transient is 0.22 +/- 0.19 m, one-tenth the coseismic slip (2 m). We analyze elastic lithosphere-viscous asthenosphere models to determine that the characteristic time describing exponential decay in deep fault slip is longer than 6 years. The VLBI measurements are consistent with uniform interseismic strain buildup. They disagree with fast postseismic rates caused by an asthenosphere with very low viscosity.

Description: The AEPI (Atmospheric Emissions Photometric Imager) experiment on the ATLAS-1 shuttle mission (launched March 24, 1992) imaged the earth night airglow emission in O2 Atmospheric (0,0) bands, at 762.0 nm. Earthward views of O2 A bands show structure from gravity waves which exhibit extended horizontal structure with horizontal wavelenghts on the order of 50-100 km. These observations of the O2 A (0,0) bands are particularly interesting since in this wavelength the lower atmosphere absorbs all the earth-reflected emissions and most of the spectrally diffuse backgrounds. Herein we present observations of gravity waves using a topside airglow imaging technique.

Description: Fujii et al. (1994) obtained characteristics of the electrodynamic parameters, that is, field-aligned currents, electric fields, and electron precipitation, which are associted with auroral substorm events in the nighttime sector, through a unique analysis that places the ionospheric measurements of these parameters into the context of a generic substorm determined from global auroral images. In this paper we investigate in considerably more detail the characteristics of the field-aligned currents using data from the same set of passes as the previuos study. We show for the first time that the net upward field-aligned currents throughout the surge and surge horn are sufficient to account for most if not all of the converging currents of the auroral electrojets. Current densities are largest in the surge and surge horn. Current region continuity does not appear to exist across the substorm bulge region. Much of the auroral substorm field-aligned current is composed of filamentary currents and finite current segments at large angles to each other. The westward electrojet may contain large gradients in intensity both in local time and latitude due to sets of localized field-aligned currents. The net downward current for several hours to the west of the surge is insufficient to account for the eastward electrojet, consistent with the concept that this electrojet originates primarily on the dayside. Our pattern of field-aligned currents associated with the surge has common features and also differs significantly from the patterns previously derived from data from radars and ground-based magnetometer arrays. Our pattern is considerably more complex, probably due to the much higher resolution in latitude of the satellite data. It is also larger in area, since our average substorm is much larger than those pertaining to the previous patterns, giving a substorm wedge considerably wider than that obtained from the radar and array data.

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

Description: An analytical model of the bidirectional reflectance for optically semi-infinite plant canopies has been extended to describe the reflectance of finite depth canopies contributions from the underlying soil. The model depends on 10 independent parameters describing vegetation and soil optical and structural properties. The model is inverted with a nonlinear minimization routine using directional reflectance data for lawn (leaf area index (LAI) is equal to 9.9), soybeans (LAI, 2.9) and simulated reflectance data (LAI, 1.0) from a numerical bidirectional reflectance distribution function (BRDF) model (Myneni et al., 1988). While the ten-parameter model results in relatively low rms differences for the BRDF, most of the retrieved parameters exhibit poor stability. The most stable parameter was the single-scattering albedo of the vegetation. Canopy albedo could be derived with an accuracy of less than 5% relative error in the visible and less than 1% in the near-infrared. Sensitivity were performed to determine which of the 10 parameters were most important and to assess the effects of Gaussian noise on the parameter retrievals. Out of the 10 parameters, three were identified which described most of the BRDF variability. At low LAI values the most influential parameters were the single-scattering albedos (both soil and vegetation) and LAI, while at higher LAI values (greater than 2.5) these shifted to the two scattering phase function parameters for vegetation and the single-scattering albedo of the vegetation. The three-parameter model, formed by fixing the seven least significant parameters, gave higher rms values but was less sensitive to noise in the BRDF than the full ten-parameter model. A full hemispherical reflectance data set for lawn was then interpolated to yield BRDF values corresponding to advanced very high resolution radiometer (AVHRR) scan geometries collected over a period of nine days. The resulting parameters and BRDFs are similar to those for the full sampling geometry, suggesting that the limited geometry of AVHRR measurements might be used to reliably retrieve BRDF and canopy albedo with this model.

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

Description: ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX computers to emulate Cray system time calls, which should be easy to modify for other machines as well. The standard distribution media for this version is a 9-track 1600 BPI ASCII Card Image format magnetic tape. The Cray version was developed in 1987. The IBM ES/3090 version is an IBM port of the Cray version. It is written in IBM VS FORTRAN and has the capability of executing in both vector and parallel modes on the MVS/XA operating system and in vector mode on the VM/XA operating system. Various options of the IBM VS FORTRAN compiler provide new features for the ES/3090 version, including 64-bit arithmetic and up to 2 GB of virtual addressability. The IBM ES/3090 version is available only as a 9-track, 1600 BPI IBM IEBCOPY format magnetic tape. The IBM ES/3090 version was developed in 1989. The DEC RISC ULTRIX version is a DEC port of the Cray version. It is written in FORTRAN 77 for RISC-based Digital Equipment platforms. The memory requirement is approximately 7Mb of main memory. It is available in UNIX tar format on TK50 tape cartridge. The port to DEC RISC ULTRIX was done in 1990. COS and UNICOS are trademarks and Cray is a registered trademark of Cray Research, Inc. IBM, ES/3090, VS FORTRAN, MVS/XA, and VM/XA are registered trademarks of International Business Machines. DEC and ULTRIX are registered trademarks of Digital Equipment Corporation.

Description: Panel method computer programs are software tools of moderate cost used for solving a wide range of engineering problems. The panel code PMARC_12 (Panel Method Ames Research Center, version 12) can compute the potential flow field around complex three-dimensional bodies such as complete aircraft models. PMARC_12 is a well-documented, highly structured code with an open architecture that facilitates modifications and the addition of new features. Adjustable arrays are used throughout the code, with dimensioning controlled by a set of parameter statements contained in an include file; thus, the size of the code (i.e. the number of panels that it can handle) can be changed very quickly. This allows the user to tailor PMARC_12 to specific problems and computer hardware constraints. In addition, PMARC_12 can be configured (through one of the parameter statements in the include file) so that the code's iterative matrix solver is run entirely in RAM, rather than reading a large matrix from disk at each iteration. This significantly increases the execution speed of the code, but it requires a large amount of RAM memory. PMARC_12 contains several advanced features, including internal flow modeling, a time-stepping wake model for simulating either steady or unsteady (including oscillatory) motions, a Trefftz plane induced drag computation, off-body and on-body streamline computations, and computation of boundary layer parameters using a two-dimensional integral boundary layer method along surface streamlines. In a panel method, the surface of the body over which the flow field is to be computed is represented by a set of panels. Singularities are distributed on the panels to perturb the flow field around the body surfaces. PMARC_12 uses constant strength source and doublet distributions over each panel, thus making it a low order panel method. Higher order panel methods allow the singularity strength to vary linearly or quadratically across each panel. Experience has shown that low order panel methods can provide nearly the same accuracy as higher order methods over a wide range of cases with significantly reduced computation times; hence, the low order formulation was adopted for PMARC_12. The flow problem is solved by modeling the body as a closed surface dividing space into two regions: the region external to the surface in which an unknown velocity potential exists representing the flow field of interest, and the region internal to the surface in which a known velocity potential (representing a fictitious flow) is prescribed as a boundary condition. Both velocity potentials are required to satisfy Laplace's equation. A surface integral equation for the unknown potential external to the surface can be written by applying Green's Theorem to the external region. Using the internal potential and zero flow through the surface as boundary conditions, the unknown potential external to the surface can be solved for. When the internal flow option, which allows the analysis of closed ducts, wind tunnels, and similar internal flow problems, is selected, the geometry is modeled such that the flow field of interest is inside the geometry and the fictitious flow is outside the geometry. Items such as wings, struts, or aircraft models can be included in the internal flow problem. The time-stepping wake model gives PMARC_12 the ability to model both steady and unsteady flow problems. The wake is convected downstream from the wake-separation line by the local velocity field. With each time step, a new row of wake panels is added to the wake at the wake-separation line. Time stepping can start from time t=0 (no initial wake) or from time t=t0 (an initial wake is specified). A wide range of motions can be prescribed, including constant rates of translation, constant rate of rotation about an arbitrary axis, oscillatory translation, and oscillatory rotation about any of the three coordinate axes. Investigators interested in a visual representation of the phenomenon they are studying with PMARC_12 may want to consider obtaining the program GVS (ARC-13361), the General Visualization System. GVS is a Silicon Graphics IRIS program which was created for the purpose of supporting the scientific visualization needs of PMARC_12. GVS is available separately from COSMIC. PMARC_12 is written in standard FORTRAN 77, with the exception of the NAMELIST extension used for input. This makes the code fairly machine independent. A compiler which supports the NAMELIST extension is required. The amount of free disk space and RAM memory required for PMARC_12 will vary depending on how the code is dimensioned using the parameter statements in the include file. The recommended minimum requirements are 20Mb of free disk space and 4Mb of RAM. PMARC_12 has been successfully implemented on a Macintosh II running System 6.0.7 or 7.0 (using MPW/Language Systems Fortran 3.0), a Sun SLC running SunOS 4.1.1, an HP 720 running HP-UX 8.07, an SGI IRIS running IRIX 4.0 (it will not run under IRIX 3.x.x without modifications), an IBM RS/6000 running AIX, a DECstation 3100 running ULTRIX, and a CRAY-YMP running UNICOS 6.0 or later. Due to its memory requirements, this program does not readily lend itself to implementation on MS-DOS based machines. The standard distribution medium for PMARC_12 is a set of three 3.5 inch 800K Macintosh format diskettes and one 3.5 inch 1.44Mb Macintosh format diskette which contains an electronic copy of the documentation in MS Word 5.0 format for the Macintosh. Alternate distribution media and formats are available upon request, but these will not include the electronic version of the document. No executables are included on the distribution media. This program is an update to PMARC version 11, which was released in 1989. PMARC_12 was released in 1993. It is available only for use by United States citizens.

Description: This program determines the supersonic flowfield surrounding three-dimensional wing-body configurations of a delta wing. It was designed to provide the numerical computation of three dimensional inviscid, flowfields of either perfect or real gases about supersonic or hypersonic airplanes. The governing equations in conservation law form are solved by a finite difference method using a second order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically. The flowfield between the body and outermost shock is treated in a shock capturing fashion and therefore allows for the correct formation of secondary internal shocks . The program operates in batch mode, is in CDC update format, has been implemented on the CDC 7600, and requires more than 140K (octal) word locations.

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

Description: MSIS is an empirical model of the thermosphere based on Mass Spectrometer and Incoherent Scatter data. It provides a description of atmospheric temperature, density, and composition for altitudes higher than 85 kilometers. There are coefficients in MSIS to account for yearly and daily variations, geodetic latitude and longitude, altitude, and solar activity. Variations due to magnetic storms are represented by three-hour magnetic ap indices. The MSIS model enables a more timely prediction of aeronomic densities for specific events, such as rocket flights. The database for this model is a comprehensive summary of rocket flight, satellite, incoherent scatter radar, grenade, and falling sphere measurements. Subsets of data were formed by random selection after sorting on altitude, latitude, time of day, etc. Curve fitting was done with four to five thousand data points at a time. The resulting coefficients are presented in subroutines which calculate thermospheric composition and temperature for a user-supplied position and time. MSIS is written in FORTRAN 77 for use with batch or interactive programs and has been implemented on a DEC VAX series computer operating under VMS 4.3 with a central memory requirement of approximately 25K of 8 bit bytes. MSIS is based on a 1977 thermosphere model and was last updated in 1987 to reflect the CIRA 1986 Neutral Thermosphere Model.

Description: This theoretical aerodynamics program, TAD, was developed to predict the aerodynamic characteristics of vehicles with sounding rocket configurations. These slender, axisymmetric finned vehicle configurations have a wide range of aeronautical applications from rockets to high speed armament. Over a given range of Mach numbers, TAD will compute the normal force coefficient derivative, the center-of-pressure, the roll forcing moment coefficient derivative, the roll damping moment coefficient derivative, and the pitch damping moment coefficient derivative of a sounding rocket configured vehicle. The vehicle may consist of a sharp pointed nose of cone or tangent ogive shape, up to nine other body divisions of conical shoulder, conical boattail, or circular cylinder shape, and fins of trapezoid planform shape with constant cross section and either three or four fins per fin set. The characteristics computed by TAD have been shown to be accurate to within ten percent of experimental data in the supersonic region. The TAD program calculates the characteristics of separate portions of the vehicle, calculates the interference between separate portions of the vehicle, and then combines the results to form a total vehicle solution. Also, TAD can be used to calculate the characteristics of the body or fins separately as an aid in the design process. Input to the TAD program consists of simple descriptions of the body and fin geometries and the Mach range of interest. Output includes the aerodynamic characteristics of the total vehicle, or user-selected portions, at specified points over the mach range. The TAD program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 123K of 8 bit bytes. The TAD program was originally developed in 1967 and last updated in 1972.

Description: This program, which is called 'AOFA', determines the complete viscous and inviscid flow around a body of revolution at a given angle of attack and traveling at supersonic speeds. The viscous calculations from this program agree with experimental values for surface and pitot pressures and with surface heating rates. At high speeds, lee-side flows are important because the local heating is difficult to correlate and because the shed vortices can interact with vehicle components such as a canopy or a vertical tail. This program should find application in the design analysis of any high speed vehicle. Lee-side flows are difficult to calculate because thin-boundary-layer theory is not applicable and the concept of matching inviscid and viscous flow is questionable. This program uses the parabolic approximation to the compressible Navier-Stokes equations and solves for the complete inviscid and viscous regions of flow, including the pressure. The parabolic approximation results from the assumption that the stress derivatives in the streamwise direction are small in comparison with derivatives in the normal and circumferential directions. This assumption permits the equation to be solved by an implicit finite difference marching technique which proceeds downstream from the initial data point, provided the inviscid portion of flow is supersonic. The viscous cross-flow separation is also determined as part of the solution. To use this method it is necessary to first determine an initial data point in a region where the inviscid portion of the flow is supersonic. Input to this program consists of two parts. Problem description is conveyed to the program by namelist input. Initial data is acquired by the program as formatted data. Because of the large amount of run time this program can consume the program includes a restart capability. Output is in printed format and magnetic tape for further processing. This program is written in FORTRAN IV and has been implemented on a CDC 7600 with a central memory requirement of approximately 35K (octal) of 60 bit words.