ALBERT

Description: The test capabilities of the Stability Wind Tunnel of the Virginia Polytechnic Institute and State University are described, and calibrations for curved and rolling flow techniques are given. Oscillatory snaking tests to determine pure yawing derivatives are considered. Representative aerodynamic data obtained for a current fighter configuration using the curved and rolling flow techniques are presented. The application of dynamic derivatives obtained in such tests to the analysis of airplane motions in general, and to high angle of attack flight conditions in particular, is discussed.

Description: The eruption of Mt. Pinatubo in June 1991 produced the largest enhancement of stratospheric aerosol loading ever observed by lidar over Hampton, Virginia. Low altitude layers (less than 20 km) were the first to arrive over Hampton in early August, the result of transport associated with a tropospheric anticyclonic cell over North America. The maximum peak scattering ratio, 34 at 22.4 km, and the maximum stratospheric integrated backscatter of 0.0053 sr(exp -1), both at 694 nm, observed since the eruption were measured on February 20, 1992. After decreasing during the spring and summer of 1992, the aerosol burden increased significantly during the winter of 1992-93, evidence of a poleward winter transport from the equatorial reservoir. Over the period from February 1992 to February 1994, the stratospheric aerosol loading decreased with an average 1/e decay time of 10.1 months. The vertical distribution, intensity, and transport of Pinatubo aerosols over this site are described and compared with similar measurements after El Chichon.

Description: T-matrix computations of light scattering by polydispersions of randomly oriented nonspherical aerosols and Mie computations for equivalent spheres are compared. Findings show that even moderate nonsphericity results in suubstantial errors in the retrieved aerosol optical thickness if satellite reflectance measurements are analyzed using Mie theory. On the other hand, the use of Mie theory for nonspherical aerosols produces negligible errors in the computation of albedo and flux related quantities, provided that the aerosol size distribution and optical thickness are known beforehand. The first result can be explained by large nonspherical-spherical differences in scattering phase function, while the second result follows from small nonspherical-spherical differences in single-scattering albedo and asymmetry parameter. No cancellation of errors occurs if one consistently uses Mie theory in the retrieval algorithm and then in computing the albedo for the retrieved aerosol optical thickness.

Description: Azimuthal asymmetries in the atmospheric refractive index can lead to errors in estimated vertical and horizontal station coordinates. Daily average gradient effects can be as large as 50 mm of delay at a 7 deg elevation. To model gradients, the constrained estimation of gradient paramters was added to the standard VLBI solution procedure. Here the analysis of two sets of data is summarized: the set of all geodetic VLBI experiments from 1990-1993 and a series of 12 state-of-the-art R&D experiments run on consecutive days in January 1994. In both cases, when the gradient parameters are estimated, the overall fit of the geodetic solution is improved at greater than the 99% confidence level. Repeatabilities of baseline lengths ranging up to 11,000 km are improved by 1 to 8 mm in a root-sum-square sense. This varies from about 20% to 40% of the total baseline length scatter without gradient modeling for the 1990-1993 series and 40% to 50% for the January series. Gradients estimated independently for each day as a piecewise linear function are mostly continuous from day to day within their formal uncertainties.

Description: Particles leving the neutral sheet in the distant magnetotail at times display adiabatic trajectory sequences characterized by an inflection toward the equator and subsequent mirroring in its vicinity. We demonstrate that this low-latitude mirroring results primarily from a centrifugal deceleration due to the fast direction-changing E x B drift. This effect which we refer to as 'centrifugal trapping' appears both in guiding centre and full particle treatments. It thus does not directly relate to nonadiabatic motion. However, pitch angle scattering due to nonadiabatic neutral sheet interaction does play a role in reducing the parallel speed of the particles. We show that centrifugal trapping is an important mechanism for the confinement of the slowest (typically below the equatorial E x B drift speed) plasma sheet populations to the midplane vicinity.

Description: Measurements of wing buffeting, using root strain gages, were made in the NASA Langley 0.3 m cryogenic wind tunnel to refine techniques which will be used in larger cryogenic facilities such as the United States National Transonic Facility (NTF) and the European Transonic Wind Tunnel (ETW). The questions addressed included the relative importance variations in frequency parameter and Reynolds number, the choice of model material (considering both stiffness and damping) and the effects of static aeroelastic distortion. The main series of tests was made on three half models of slender 65 deg delta wings with a sharp leading edge. The three delta wings had the same planform but widely differing bending stiffnesses and frequencies (obtained by varying both the material and the thickness of the wings). It was known that the steady flow on this configuration would be insensitive to variations in Reynolds number. On this wing at vortex breakdown the spectrum of the unsteady excitation is unusual, having a sharp peak at particular frequency parameter. Additional tests were made on one unswept half-wing of aspect ratio 1.5 with an NPL 9510 aerofoil section, known to be sensitive to variations in Reynolds number at transonic speeds. The test Mach numbers were M = 0.21 and 0.35 for the delta wings and to M = 0.30 for the unswept wing. On this wing the unsteady excitation spectrum is fairly flat (as on most wings). Hence correct representation of the frequency parameter is not particularly important.

Description: Measurements of superthermal electron fluxes in the solar wind indicate that field lines within coronal mass ejections, CMEs, near and beyond 1 AU are normally connected to the Sun at both ends. However, on occasion some field lines embedded deep within CMEs appear to be connected to the Sun at only one end. Here we propose an explanation for how such field lines arise in terms of 3-dimensional reconnection close to the Sun. Such reconnection also provides a natural explanation for the flux rope topology characteristic of many CMEs as well as the coronal loops formed during long-duration, solar X-ray events. Our consideration of the field topologies resulting from 3-dimensional reconnection indicates that field lines within and near CMEs may on occasion be connected to the outer heliosphere at both ends.

Description: We examined 11 cases when the interplanetary magnetic field (IMF) was intensely northward (greater than 10 nT) for long durations of time (greater than 3 hours), to quantitatively determine an uppler limit on the efficiency of solar wind energy injection into the magnetosphere. We have specifically selected these large B(sub N) events to minimize the effects of magnetic reconnection. Many of these cases occurred during intervals of high-speed streams associated with coronal mass ejections when viscous interaction effects might be at a maximum. It is found that the typical efficiency of solar wind energy injection into the magnetosphere is 1.0 x 10(exp -3) to 4.0 x 10(exp -3), 100 to 30 times less efficient than during periods of intense southward IMFs. Other energy sinks not included in these numbers are discussed. Estimates of their magnitudes are provided.

Description: The subject of this paper is a self-consistent, magnetohydrodynamic numerical realization for the Earth's magnetosphere which is in a quasi-steady dynamic equilibrium for a due northward interplanetary magnetic field (IMF). Although a few hours of steady northward IMF are required for this asymptotic state to be set up, it should still be of considerable theoretical interest because it constitutes a 'ground state' for the solar wind-magnetosphere interaction. Moreover, particular features of this ground state magnetosphere should be observable even under less extreme solar wind conditions. Certain characteristics of this magnetosphere, namely, NBZ Birkeland currents, four-cell ionospheric convection, a relatively weak cross-polar potential, and a prominent flow boundary layer, are widely expected. Other characteristics, such as no open tail lobes, no Earth-connected magnetic flux beyond 155 R(sub E) downstream, magnetic merging in a closed topology at the cusps, and a 'tadpole' shaped magnetospheric boundary, might not be expected. In this paper, we will present the evidence for this unusual but interesting magnetospheric equilibrium. We will also discuss our present understanding of this singular state.

Description: The stability of collisionless tearing modes is examined in the presence of curvature drift resonances and the trapped particle effects. A kinetic description for both electrons and ions is employed to investigate the stability of a two-dimensional equilibrium model. The main features of the study are to treat the ion dynamics properly by incorporating effects associated with particle trajectories in the tail fields and to include the linear coupling of trapped particle modes. Generalized dispersion relations are derived in several parameter regimes by considering two important sublayers of the reconnecting region. For a typical choice of parameters appropriate to the current sheet region, we demonstrate that localized tearing modes driven by ion curvature drift resonance effects are excited in the current sheet region with growth time of the order of a few seconds. Also, we examine nonlocal characteristics of tearing modes driven by curvature effects and show that modes growing in a fraction of a second arise when mode widths are larger than the current sheet width. Further, we show that trapped particle effects, in an interesting frequency regime, significantly enhance the growth rate of the tearing mode. The relevance of this theory for substorm onset phase and other features of the substorms is briefly discussed.

Description: The second pulsating aurora rocket (PULSAUR 2) rocket was launched into a pulsating aurora on 9 February 1994 and carried a broad range of instruments in order to perform a study of this type of aurora. The rocket measurements were complemented with a set of ground-based measurements. The particle measurements performed in the rocket are related to the ground-based optical measurements performed along the rocket trajectory. It is found that the high energy electrons are largely in phase with the measured luminosity. The EISCAT measurements carried out during the flight are reviewed. The PULSAUR 2 campaign is described. Results concerning the auroral conditions and the particle measurements are presented.

Description: The high-speed correction factor to the O(+)-O collision frequency, resulting from drift velocities between ions and neutrals, is calculated by solving the integral expression in this factor both numerically and analytically. Although the analytic solution is valid for either small or large drift velocities between ions and neutrals, for temperatures of interest and all drift velocities considered, agreement is found between analytic and detailed numerical integration results within less than 1% error. Let T(sub r) designate the average of the ion and neutral temperatures in K, and u = nu(sub d)/alpha, where nu(sub d) is the relative drift velocity in cm/s, and alpha = 4.56 x 10(exp 3) square root of T(sub r) cm/s is the thermal velocity of the O(+)-O system. Then, as u ranges from 0 to 2, the correction factor multiplying the collision frequency increases monotonically from 1 to about 1.5. An interesting result emerging from this calculation is that the correction factor for temperatures of aeronomical interest is to a good approximation independent of the temperature, depending only on the scaled velocity u.

Description: This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(exp -5).

Description: The relation between the solar wind input to the magetosphere, VB(sub South), and the auroral geomagnetic index AL is modeled with two linear moving-average filtering methods: linear prediction filters and a driven harmonic oscillator in the form of an electric circuit. Although the response of the three-parameter oscillator is simpler than the filter's, the methods yield similar linear timescales and values of the prediction-observation correlation and the prediction Chi(exp 2). Further the filter responses obtained by the two methods are similar in their long-term features. In these aspects the circuit model is equivalent to linear prediction filtering. This poses the question of uniqueness and proper interpretation of detailed features of the filters such as response peaks. Finally, the variation of timescales and filter responses with the AL activity level is discussed.

Description: Empirical data-based models of the magnetosphereic magnetic field have been widely used during recent years. However, the existing models (Tsyganenko, 1987, 1989a) have three serious deficiencies: (1) an unstable de facto magnetopause, (2) a crude parametrization by the K(sub p) index, and (3) inaccuracies in the equatorial magnetotail B(sub z) values. This paper describes a new approach to the problem; the essential new features are (1) a realistic shape and size of the magnetopause, based on fits to a large number of observed crossing (allowing a parametrization by the solar wind pressure), (2) fully controlled shielding of the magnetic field produced by all magnetospheric current systems, (3) new flexible representations for the tail and ring currents, and (4) a new directional criterion for fitting the model field to spacecraft data, providing improved accuracy for field line mapping. Results are presented from initial efforts to create models assembled from these modules and calibrated against spacecraft data sets.

Description: We have investigated the middle atmospheric response to the 27-day and 11-yr solar UV flux variations at low to middle latitudes using a two-dimensional photochemical model. The model reproduced most features of the observed 27-day sensitivity and phase lag of the profile ozone response in the upper stratosphere and lower mesosphere, with a maximum sensitivity of +0.51% per 1% change in 205 nm flux. The model also reproduced the observed transition to a negative phase lag above 2 mb, reflecting the increasing importance with height of the solar modulated HO(x) chemistry on the ozone response above 45 km. The model revealed the general anti-correlation of ozone and solar UV at 65-75 km, and simulated strong UV responses of water vapor and HO(x) species in the mesosphere. Consistent with previous 1D model studies, the observed upper mesospheric positive ozone response averaged over +/- 40 was simulated only when the model water vapor concentrations above 75 km were significantly reduced relative to current observations. In agreement with observations, the model computed a low to middle latitude total ozone phase lag of +3 days and a sensitivity of +0.077% per 1% change in 205 nm flux for the 27-day solar variation, and a total ozone sensitivity of +0.27% for the 11-yr solar cycle. This factor of 3 sensitivity difference is indicative of the photochemical time constant for ozone in the lower stratosphere which is comparable to the 27-day solar rotation period but is much shorter than the 11-yr solar cycle.

Description: The accuracy of Geosat satellite altimetry over the Greenland ice sheet is evaluated by comparing the measured heights to radar elevations from the airborne Greenland Aerogeophysics Project. At the center of the ice sheet where the ice surface is nearly level, surface comparisons show a fit at the 1 to 3 m level as expected, but even at moderately sloping ice regions (0.3 deg-0.6 deg), satellite altimetry mean errors in the range of 10 to 35 m are observed. These errors are found for slope-corrected and waveform-retracked data, so most previous accuracy estimates of current satellite altimetry ice sheet elevations in regions of slopping or undulating ice appear to be too optimistic.

Description: We report here on a number of examples of anomalous enhancements of eastward electric fields near sunrise in the equatorial ionospheric F-region. These examples were selected from the data base of the equatorial satellite, San Marco D (1988), which measured ionospheric electric fields during a period of solar minimum. The eastward electric fields reported correspond to vertical plasma drifts. The examples studied here are similar in signature and polarity to the pre-reversal electric field enhancements seen near sunset from ground-based radar systems. The morphology of these sunrise events, which are observed on about 14% of the morning-side satellite passes, are studied as a function of local zonal velocity, magnetic activity, geographic longitude and altitude. The nine events studied occur at locations where the zonal plasma flow is generally measured to be eastward, but reducing as a function of local time and at satellite longitudes where the magnetic declination has the opposite polarity as the declination of the sunrise terminator.

Description: A nonlinear filtering method is introduced for the study of the solar wind -- magnetosphere coupling and related to earlier linear techniques. The filters are derived from the magnetospheric state, a representation of the magnetospheric conditions in terms of a few global variables, here the auroral electrojet indices. The filters also couple to the input, a representation of the solar wind variables, here the rectified electric field. Filter-based iterative prediction of the indices has been obtained for up to 20 hours. The prediction is stable with respect to perturbations in the initial magnetospheric state; these decrease exponentially at the rate of 30/min. The performance of the method is examined for a wide range of parameters and is superior to that of other linear and nonlinear techniques. In the magnetospheric state representation the coupling is modeled as a small number of nonlinear equations under a time-dependent input.

Description: Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Description: Many problems in geophysical and astrophysical convection systems are characterized by fast rotation and spherical shell geometry. The combined effects of Coriolis forces and spherical shell geometry produce a unique spatial symmetry for the convection pattern in a rapidly rotating spherical shell. In this paper, we first discuss the general spatial symmetries for rotating spherical shell convection. A special model, a spherical shell heated from below, is then used to illustrate how and when the spatial symmetries are broken. Symmetry breaking occurs via a sequence of spatial transitions from the primary conducting state to the complex multiple-layered columnar structure. It is argued that, because of the dominant effects of rotation, the sequence of spatial transitions identified from this particular model is likely to be generally valid. Applications of the spatial symmetry breaking to planetary convection problems are also discussed.

Description: Initial observations of a newly documented type of optical emission above thunderstorms are reported. 'Blue jets', or narrowly collimated beams of blue light that appear to propagate upwards from the tops of thunderstorms, were recorded on B/W and color video cameras for the first time during the Sprites94 aircraft campaign, June-July, 1994. The jets appear to propagate upward at speeds of about 100 km/s and reach terminal altitudes of 40-50 m. Fifty six examples were recorded during a 22 minute interval during a storm over Arkansas. We examine some possible mechanisms, but have no satisfactory theory of this phenomenon.

Description: It is shown from flux transfer event (FTE) occurrence statistics, observed as a function of MLT by the ISEE satellites, that recent 2-dimensional analytic theories of the effects of pulsed Petschek reconnection predict FTEs to contribute between 50 and 200 kV to the total reconnection voltage when the magnetosheath field points southward. The upper limit (200 kV) allows the possibility that FTEs provide all the antisunward transport of open field lines into the tail lobe. This range is compared with the voltages associated with series of FTEs signatures, as inferred from ground-based observations, which are in the range 10-60 kV. We conclude that the contribution could sometimes be made by a series of single, large events; however, the voltage is often likely to be contributed by several FTEs at different MLT.

Description: New investigations of the core ion motion within high-latitude topside ionosphere and near-Earth magnetosphere, using data from the Dynamics Explorer (DE) retarding ion mass spectrometer (RIMS), reveal the existence of significant regions of downward moving O(+). The occurences of downgoing versus upgoing O(+) are not clearly separable in terms of either polar zenith angle or Kp but are well distinquished by the direction of the z component of the interplanetary magnetic field (IMF). On the average, down flow dominates when IMFB2 less than O, while upflow dominates for IMFB2 greater than 0. Combining cross-field convection velocities derived from an empirical convection electric field model with the observed parallel velocities yields a two dimensional ion velocity field. This velocity field is consistent with a senario which has O(+) of cusp/cleft and auroral zone origin concvecting into the polar cap and, because of the dominance of gravitational energy over the upward kinetic energy, falling back into the inonsphere. This provides additional confirmation of the results of studies of the cleft ion fountain. Estimates of the flux of O(+) in the upflowing and downflowing regions for Lambda greater than or equal to 60 deg give a total upflow of approximately 6 X 10(exp 25) ions/sec for IMFB2 greater than 0 and total upflow and downflow of approximately 4 X 10 (exp 25) ions/sec and 1 x 10(exp 25) ions/sec, respectively, for IMFB less than 0. In all cases the magnitude of the dayside outflow is consistent with previous work on upwelling ions. While the magnitudes vary for high and low Kp the ratios of upward to downward flow are roughly the same at approximately 1.7. The downflowing O(+) shows a correlation with the magnitude of the outflow of light ions in the same region but the cause and effect of this relationship is not distinquishable. hable.

Description: We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Description: The effects on double-probe electric field measurements induced by electron density and temperature gradients are investigated. We show that on some occasions such gradients may lead to marked spurious electric fields if the probes are assumed to lie at the same probe potential with repect to the plasma. The use of a proper bias current will decrease the magnitude of such an error. When the probes are near the plasma potential, the magnitude of these error signals, delta Epsilon, can vary as delta Epsilon approx. T(sub e)(Delta n(sub e)/n(sub e)) + 0.5 Delta T(sub e), where T(sub e) is the electron temperature, Delta n(sub e)/ n(sub e) the relative electron density variation between the two sensors, and Delta T(sub e) the electron temperature difference between the two sensors. This not only implies that the error signals will increase linearly with the density variations but also that such signatures grow with Delta T(sub e) i.e., such effects are 10 times larger in a 10-eV plasma than in a 1-eV plasma. This type of error is independent of the probe separation distance provided the gradient scale length is much larger than the distance. The largest errors occur when the probes are near to the plasma potential. During the crossing of a small structure (e.g, a double layer) the error signal appears as a bipolar signature. Our analysis shows that errors in double-probe measurements caused by plasma gradients are not significant at large scale (much greater than 1 km) plasma boundaries, and may only be important in cases where small-scale (less than 1 km), internal gradient structures exist. Bias currents tailored for each plasma parameter regime (i.e., variable bias current) would improve the double-probe response to gradient effects considerably.

Description: The evolution of the volcanic debris plume originating from the June 1991 eruption of Mt. Pinatubo has been monitored since its genesis using a ground-based backscatter lidar facility sited at the Jet Propulsion Laboratory (JPL). Both absolute and relative pre- and post-Pinatubo backscatter observations are in accord with Mie scattering projections based on measured aerosol particle size distributions reported in the literature. The post-Pinatubo column-integrated backscatter coefficient peaked approximately 400 days after the eruption, and the observed upper boundary of the aerosol column subsided at a rate of approximately 200 m/mon.

Description: Using Global Positioning System (GPS) receivers, we reoccupied several leveling benchmarks on the Kenai Peninsula of Alaska which had been surveyed by conventional leveling immediately following the March 27, 1964, Prince William Sound earthquake (M(sub w) = 9.3). By combining the two sets of measurements with a new, high-resolution model of the geoid in the region, we were able to determine the cumulative 1993-1964 postseismic vertical displacement. We find uplift at all of our benchmarks, relative to Seward, Alaska, a point that is stable according to tide gauge data. The maximum uplift of about 1 m occurs near the middle of the peninsula. The region of maximum uplift appears to be shifted northwest relative to the point of maximum coseismic subsidence. If we use tide gauge data at Nikishka and Seward to constrain the vertical motion, then the observed uplift has a trenchward tilt (down to the southeast) as well as an arching component. To explain the observations, we use creep-at-depth models. Most acceptable models require a fault slip of about 2.75 m, although this result is not unique. If the slip has been continuous since the 1964 earthquake, then the average slip rate is nearly 100 mm/yr, twice the plate convergence rate. Comparing the net uplift achieved in 29 years with that observed over 11 years in an adjacent region southeast of Anchorage, Alaska, we conclude that the rate of uplift is decreasing. A further decrease in the uplift rate is expected as the 29-year averaged displacement rate is about twice the plate convergence rate and therefore cannot be sustained over the entire earthquake cycle.

Description: Narrow enhancements of electron precipitation, with energy and flux well above typical values, have been observed with Dynamics Explorer 2 (DE 2) in the cusp/cleft region. The electron flux in the energy range 0.2-1 keV was 2 orders of magnitude higher in these structures than in the magnetosheath and were seen in approximaetly 80% of DE 2 cusp crossings at ionospheric altitudes. Typically, there was more than one electron structure in each cusp crossing. The position of these structures showed a systematic variation: for poleward ion dispersion (energy decreases with increasing latitude), electron structures were seen more often on the equatorial boundary of the cusp, while for equatorward ion dispersion (energy decreases with decreasing latitude), electron structures were more often seen on the poleward boundary. This suggests that the electron structures are associated with newly reconnected field lines. The electron spectra suggest that field-aligned acceleration processes could produce the electron structures, first near the boundary of the cusp/cleft during the reconnection of field lines and then in the cusp/cleft during the motion of reconnected flux tubes through the polar ionosphere.

Description: A permanent Global Positioning System (GPS) receiver at Casa Diablo Hot Springs, Long Valley Caldera, California was installed in January, 1993, and has operated almost continuously since then. The data have been transmitted daily to the Jet Propulsion Laboratory (JPL) for routine analysis with data from the Fiducial Laboratories for an International Natural sciences Network (FLINN) by the JPL FLINN analysis center. Results from these analyses have been used to interpret the on going deformation at Long Valley, with data excluded from periods when the antenna was covered under 2.5 meters of snow and from some periods when Anti Spoofing was enforced on the GPS signal. The remaining time series suggests that uplift of the resurgent dome of Long Valley Caldera during 1993 has been 2.5 +/- 1.1 cm/yr and horizontal motion has been 3.0 +/- 0.7 cm/yr at S53W in a no-net-rotation global reference frame, or 1.5 +/- 0.7 cm/yr at S14W relative to the Sierra Nevada block. These rates are consistent with uplift predicted from frequent horizontal strain measurements. Spectral analysis of the observations suggests that tidal forcing of the magma chamber is not a source of the variability in the 3 dimensional station location. These results suggest that remotely operated, continuously recording GPS receivers could prove to be a reliable tool for volcanic monitoring throughout the world.

Description: A comparison is made of the vertical distribution of high-level cloud tops derived from the Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements and from the International Satellite Cloud Climatology Project (ISCCP) for all Julys and Januarys in 1985 to 1990. The results suggest that ISCCP overestimates the pressure of high-level clouds by up to 50-150 mbar, particularly at low latitudes. This is caused by the frequent presence of clouds with diffuse tops (greater than 50% time when cloudy events are observed). The averaged vertical extent of the diffuse top is about 1.5 km. At midlatitudes where the SAGE II and ISCCP cloud top pressure agree best, clouds with distinct tops reach a maximum relative proportion of the total level cloud amount (about 30-40%), and diffuse-topped clouds are reduced to their minimum (30-40%). The ISCCP-defined cloud top pressure should be regarded not as the material physical height of the clouds but as the level which emits the same infrared radiance as observed. SAGE II and ISCCP cloud top pressures agree for clouds with distinct tops. There is also an indication that the cloud top pressures of optically thin clouds not overlying thicker clouds are poorly estimated by ISCCP at middle latitudes. The average vertical extent of these thin clouds is about 2.5 km.

Description: Global high-level clouds identified in Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements for January and July in the period 1985 to 1990 are compared with near-nadir-looking observations from the International Satellite Cloud Climatology Project (ISCCP). Global and zonal mean high-level cloud amounts from the two data sets agree very well, if clouds with layer extinction coefficients of less than 0.008/km at 1.02 micrometers wavelength are removed from the SAGE II results and all detected clouds are interpreted to have an average horizontal size of about 75 km along the 200 km transimission path length of the SAGE II observations. The SAGE II results are much more sensitive to variations of assumed cloud size than to variations of detection threshold. The geographical distribution of cloud fractions shows good agreement, but systematic regional differences also indicate that the average cloud size varies somewhat among different climate regimes. The more sensitive SAGE II results show that about one third of all high-level clouds are missed by ISCCP but that these clouds have very low optical thicknesses (less than 0.1 at 0.6 micrometers wavelength). SAGE II sampling error in monthly zonal cloud fraction is shown to produce no bias, to be less than the intraseasonal natural variability, but to be comparable with the natural variability at longer time scales.

Description: Enhancements in the fluxes of relativistic electrons trapped within the Earth's magnetosphere have been measured by the high-energy particle spectrometer, part of the particle environment monitor on the upper atmosphere research satellite (UARS). The largest increase in the electron fluxes with energies greater than 1 MeV observed on UARS from October 1991 through July 1994 was in early May 1992. The fluxes of trapped electrons in the drift loss cone and locally precipitating electrons showed differing buildup and decay rates as a function of invariant latitude. Increases of more than 2 orders of magnitude were observed in drift loss cone fluxes at magnetic latitudes of 40 deg-66 deg and in precipitating fluxes from 48 deg to 66 deg. The energy flux contained in the most intense local precipitation observed was approximately 0.1 erg/sq cm/s, entering the atmosphere and creating up to 1000 ion pairs/cu cm/s at 55-km altitude. The daily averaged energy flux from directly precipitating electrons with energies greater than 1 MeV deposited greater than 10(exp 20) erg/d worldwide into the atmosphere for the period May 12-21, 1992, producing greater than 10(exp 31) odd nitrogen molecules below 60-km altitude.

Description: Prandtl's theory is used to determine the airflow over bodies and wings adapted to supersonic flight. By making use of these results, and by incorporating in them an allowance for the probable skin friction, some estimates of expected lift-drag ratios are made for various flight speeds with the best configuration. At each speed a slender body and wings having the best angle of sweepback are considered. For the range of supersonic speeds shown an airplane of normal density and loading would be required to operate at an altitude of the order of 60,000 feet. The limiting value of 1-1/2 times the speed of sound corresponds to a flight speed of 1000 miles per hour. At this speed about 1.5 miles per gallon of fuel are expected. It is interesting to note that this value corresponds to a value of more than 15 miles per gallon when the weight is reduced to correspond to that of an ordinary automobile.

Description: In theory, the most efficient wing shape for transonic and low supersonic speeds is simply a long narrow straight subsonic wing turned at an oblique angle to the flight direction. This theory has been verified by tests at Mach numbers from .6 to 1.4 in supersonic wind tunnel and by comparative studies of transonic transport designs.

Description: Recent theoretical and experimental work in supersonic aerodynamics is reviewed with its practical application in mind. Several arrangements of supporting surfaces and bodies are discussed and in some cases comparisons of theory and experiment are made. Finally, certain phenomena connected with lift and drag in a rarefied medium are considered briefly.

Description: A method is reported for determining mathematically the combined disturbance field, and in certain cases the minimum drag, of wings at supersonic speeds. The simplest analytic example is provided by the wing of elliptic planform, which achieves its minimum drag when the lift is distributed uniformly over the surface. With a symmetrical distribution of thickness, the requirement of minimum drag for a given total volume is found to lead to profiles of constant curvature.

Description: The assumptions of the thin airfoil theory are found to provide certain necessary conditions for the minimum drag of airfoils having a given total lift, a given maximum thickness, or a given volume. The conditions are applicable to steady or unsteady motions and to subsonic or supersonic speeds without restriction on the planform. The computation of drag and the statement of the conditions for minimum drag depend on the consideration of a combined flow field, which is obtained by superimposing the disturbance velocities in forward and reversed motions. If the planform of the airfoil and its total lift are given, it is found that, for minimum drag, the lift must be distributed in such a way that the downwash in the combined field is constant over the entire planform. If the planform is given and the thickness of the airfoil is required to contain a specified volume, then the thickness must be distributed over the planform in such a way that the pressure gradient of the combined field in the direction of flight is constant at all points of the wing.

Description: The application of mathematical advances made in electricity and other branches to problems of airplane dynamics is demonstrated. The Heaviside-Bromwich methods of solution of linear differential equations are described and it is shown how these methods avoid the consideration of boundary conditions and of particular or complementary integrals. It is pointed out that if the solution of the differential equation is obtained for the case of a unit disturbance, the effect of varying disturbances may be found therefrom by Carson's theorem. A graphical solution of Carson's integral for irregular disturbances is given. The procedure of obtaining unit solutions of the equations is then taken up and the analogy between Heaviside's symbolic series solution and a physical procedure of approximation is shown. It is suggested that a fictitious impulsive disturbance be used in the treatment of initial motions.

Description: This paper presents the results of a series of total and spectral solar irradiance measurements made at ground surface (Table Mountain Facility, Calif., altitude 2.18 km). The spectral irradiance data are presented for the 0.3-3.0-micron spectral region for air mass 1.5.

Description: In linearized flow theory, certain very interesting extremal properties of wings can be derived under rather broad conditions without the use of a complicated mathematical apparatus. The present chapter reviews certain results of this theory and indicates some rather obvious extensions to incorporate various auxiliary conditions. Several examples illustrating the relation between the geometrical features of the wing and the lift distribution for minimum drag are given.

Description: The items discussed are: (1) a recently proposed correction formula for the effect of compressibility in two dimensional subsonic flow; (2) the equivalence rule and the area rule for transonic speeds; (3) reciprocal relations in linearized wing theory; and (4) some general results connected with the problem of minimum wave resistance. The paper concludes with an example showing indentation of the fuselage to obtain favorable interference with the wing at supersonic speeds.

Description: In the wing section theory the magnitude of the circulation, and hence of the lift, is determined by the velocity that would be induced near the trailing edge of the section in a non-lifting potential flow. In three dimensional flow the problem is complicated by the presence of the wake and no simple basic solution has been found. Treatment of the problem of a wing of finite span is reported on the basis of the two dimensional theory, corrected for the effect of the wake.

Description: In theory, antisymmetric arrangements of wings and bodies can have smaller wave drag than corresponding mirror-symmetric arrangements. Thus, a long narrow oblique wing which presents the same aspect for two opposite directions of flight is potentially more efficient than corresponding (i.e., structurally equivalent) swept wing. The single continuous wing panel also adapts itself more readily to varying angles of obliquity, and hence, to varying flight speeds. Previous work on the aerodynamics and flight stability of oblique wing combinations is reviewed and a possible mode of application to transport aircraft operating at moderate supersonic speeds is suggested.

Description: It is shown that the drag of any semi-infinite airfoil section in purely subsonic inviscid flow follows precisely the Prandtl-Glauert compressibility rule. The result for the parabola has application to leading edge corrections in thin airfoil theory.

Description: Comparisons of wing-body combinations may not disclose the full effect of a loss in aerodynamic efficiency. If the thrust needs to be increased at a given altitude then more or larger engines will have to be used and the possibility of concealing them becomes less. In this process the lift drag ratio of the complete airplane may become still more unfavorable than indicated by the comparison. Primarily aerodynamic and structural considerations point toward the development of turbojet engines specifically adapted to operation in an atmosphere of one tenth normal density. In addition to the numerous other technological problems associated with operation at these high altitudes, the problems of safe descent and effective limitation to low speeds at low altitudes seem important.

Description: If the thin airfoil theory is applied to an airfoil having a rounded leading edge, a certain error will arise in the determination of the pressure distribution around the nose. It is shown that the evaluation of the drag of such a blunt nosed airfoil by the thin airfoil theory requires the addition of a leading edge force, analogous to the leading edge thrust of the lifting airfoil. The method of calculation is illustrated by application to: (1) The Joukowski airfoil in subsonic flow; and (2) the thin elliptic cone in supersonic flow. A general formula for the edge force is provided which is applicable to a variety of wing forms.

Description: Some of the recent advances in the theory of thin airfoils are presented with particular reference to extensions of the theory to three dimensional flows and to supersonic speeds. The problem discussed herein is the calculation of the small disturbance velocities u, v, and w in the external field produced by the flight velocity V of the airfoil.

Description: Voyager 1 images show 14 volcanic centers wholly or partly within the Kane Patera quadrangle of Io, which are divided into four major classes: (1) shield with parallel flows; (2) shield with early radial fan shapd flows; (3) shield with radial fan shaped flows, surfaces of flows textured with longitudinal ridges; and (4) depression surrounded by plateau-forming scarp-bounded, untextured deposits. The interpretation attempted here hinges largely on the ability to distinguish lava flows from pyroclastic flows by remote sensing.

Description: The proposed general analytical model describes the anisotropic, elastoplastic, path-dependent, stress-strain properties of inviscid saturated clays under undrained conditions. Model parameters are determined by using results from strain-controlled simple shear tests on a saturated clay. The model's accuracy is evaluated by applying it to predict the results of other tests on the same clay, including monotonic and cyclic loading. The model explains the very anisotropic shear strength behavior observed for weak marine clays.

Description: A proposal for a multi-institutional investigation of the processes involved in the growth and maintenance of high level extended clouds is presented. Mapping of variability of the cloud and of its radiative characteristics in terms of the meteorological environment of the cloud; performance of case studies involving observation of the cloud microphysics and radiation characteristics; and investigation of the processes responsible for the generation, maintenance, and dissipation of the cloud system are recommended. Both modeling and monitoring activities are considered. The specific research projects which the author proposes to carry out are described. Suggestions for the administrative organization of the total effort are presented.

Description: The interplay of the various physical processes involved in the formation, maintenance, and decay of middle and upper tropospheric stratiform clouds is discussed. Ice phase fair weather cloud forms are considered. Simulations of cirriform clouds which attempt to incorporate the physical processes in an interactive manner were performed. A two dimensional time dependent Eulerian numerical model, which incorporates all of the important physical processes in a simplified way, is employed to investigate the role of these processes in the evolution of a cloud in an isolated atmospheric layer. Physical parameters considered are the eddy viscosity and the thermal, water vapor, and ice water eddy diffusivities.

Description: Structural and tectonic interpretations of planetary surfaces rely strongly on visual determination of regional structural grain. This grain can be very complex and confusing, and sorting out of discrete trends in time and space is of utmost importance. This study is a test of these techniques applied to a well known area having several discrete structural grains. In the Bighorn Basin region of Wyoming, a largely overlooked N10E structural grain has been verified with detailed structural analysis and indicates a significant change in stress orientation at the end of the Laramide orogeny.

Description: Tectoism in the Valles Marineris appears to have been accompanied by volcanism. The proposed volcanic features, though probably contemporaneous with the gigantic ones in the Tharsis area, are composed of small, mafic and, possibly, somewhat larger felsic flows. The size of these features is similar to that of volcanic flows on the Earth.

Description: Experiments were carried out in a steel pressure device using controlled amounts of water and thermite melt to examine the mechanical energy released on explosive mixing following the initial contact of the two materials. An experimental design was used to allow the direct calculation of the mechanical energy by the dynamic lift of the device as recorded both optically and physically. A large number of experiments were run to accurately determine the optimum mixture of water and melt for the conversion of thermal to mechanical energy. The maximum efficiency observed was about 12% at a water/thermite mass ratio of 0.50. These experiments are the basis for the development of models of hydroexplosions and melt fragmentation. Particles collected from the experimental products are similar in size and shape to pyroclasts produced by much larger hydrovolcanic explosions. Melt rupture at optimum ratios produces very fine particles whereas rupture at high or low water/melt ratios produces large melt fragments. Grain surface textures in the experimental products are also related to the water/melt ratio and the mechanism of explosive mixing. It is thus possible to have qualitative information about the nature of the explosion from the sizes and shapes of the fragments produced.

Description: An experimental study of slotted upper and lower walls in a two dimensional transonic wind tunnel with solid sidewalls is reported. Results are presented for several slot spacings and slot openness ratios. The experimental data were pressure measurements which were made on an airfoil model and on a sidewall near one of the slotted walls. The slotted-wall boundary condition coefficient, which related the pressure and streamline curvature near the wall, was determined from the wall pressure measurements. The measured wall-induced interference was correlated with the experimental values for the boundary condition coefficient. This correlation was compared with theory.

Description: A KC-135A aircraft equipped with wing tip winglets was flight tested to demonstrate and validate the potential performance gain of the winglet concept as predicted from analytical and wind tunnel data. Flight data were obtained at cruise conditions for Mach numbers of 0.70, 0.75, and 0.80 at a nominal altitude of 36,000 ft. and winglet configurations of 15 deg cant/-4 deg incidence, 0 deg cant/-4 deg incidence, and baseline. For the Mach numbers tested the data show that the addition of winglets did not affect the lifting characteristics of the wing. However, both winglet configurations showed a drag reduction over the baseline configuration, with the best winglet configuration being the 15 deg cant/-4 deg incidence configuration. This drag reduction due to winglets also increased with increasing lift coefficient. It was also shown that a small difference exists between the 15 deg cant/-4 deg incidence flight and wind tunnel predicted data. This difference was attributed to the pillowing of the winglet skins in flight which would decrease the winglet performance.

Description: A joint NASA/USAF program was conducted to accomplish the following objectives: (1) evaluate the benefits that could be achieved from the application of winglets to KC-135 aircraft; and (2) determine the ability of wind tunnel tests and analytical analysis to predict winglet characteristics. The program included wind-tunnel development of a test winglet configuration; analytical predictions of the changes to the aircraft resulting from the application of the test winglet; and finally, flight tests of the developed configuration. Pressure distribution, loads, stability and control, buffet, fuel mileage, and flutter data were obtained to fulfill the objectives of the program.

Description: A full-scale winglet flight test on a KC-135 airplane with an upper winglet was conducted. Data were taken at Mach numbers from 0.70 to 0.82 at altitudes from 34,000 feet to 39,000 feet at stabilized flight conditions for wing/winglet configurations of basic wing tip, 15/-4 deg, 15/-2 deg, and 0/-4 deg winglet cant/incidence. An analysis of selected pressure distribution and data showed that with the basic wing tip, the flight and wind tunnel wing pressure distribution data showed good agreement. With winglets installed, the effects on the wing pressure distribution were mainly near the tip. Also, the flight and wind tunnel winglet pressure distributions had some significant differences primarily due to the oilcanning in flight. However, in general, the agreement was good. For the winglet cant and incidence configuration presented, the incidence had the largest effect on the winglet pressure distributions. The incremental flight wing deflection data showed that the semispan wind tunnel model did a reasonable job of simulating the aeroelastic effects at the wing tip. The flight loads data showed good agreement with predictions at the design point and also substantiated the predicted structural penalty (load increase) of the 15 deg cant/-2 deg incidence winglet configuration.

Description: A joint NASA/U.S. industry program to test advanced technology airfoils in the Langley 0.3-meter Transonic Tunnel (TCT) was formulated under the Langley ACEE Project Office. The objectives include providing U.S. industry an opportunity to compare their most advanced airfoils to the latest NASA designs by means of high Reynolds number tests in the same facility. At the same time, industry would again experience in the design and construction of cryogenic test techniques. The status and details of the test program are presented. Typical aerodynamic results obtained, to date, are presented at chord Reynolds number up to 45 x 10(6) and are compared to results from other facilities and theory. Details of a joint agreement between NASA and the Deutsche Forschungs- und Versuchsantalt fur Luft- and Raumfahrt e.V. (DFVLR) for tests of two airfoils are also included. Results of these tests will be made available as soon as practical.

Description: Anorthosite massifs developed approximately 1.4 to 1.5 billion years ago along an arch which developed parallel to a zone of continental separation as a block which included North America, Europe, and probably Asia separated from a block which included parts of South America, Africa, India, and Australia. Anorthosite massifs also developed at the same time along a belt which runs through the continents which comprise Gondwanaland (South America), Africa, India, Australia, and Antarctica. This was a zone of continental separation which subsequently became a zone of continental collision about 1.2 billion years ago. The northern anorthosite belt also parallels an orogenic belt which was active between 1.8 and 1.7 billion years ago. Heat generated during this mountain building period helped in the formation of the anorthosites.

Description: The relative merits of pressure increment and partial derivative formulations of the ozone inversion problem are discussed briefly. The height range of validity of the retrieved ozone profile and the effects of adding wavelengths to or of dropping wave-length from the inversion system are indicated. Illustrative results are presented for profiles retrieved from BUV data using Backus-Gilbert, minimum information (Twomey), and quasi-optimum procedures.

Description: Stratospheric limb radiance profiles versus altitude of closest approach of the line of sight to the Earth's surface have been measured before and after the Mount St. Helens eruptions by means of photographs taken from a Sun-oriented balloon gondola floating above 35 km altitude over France. Preliminary data were reported for flights in October 1979 and in May and June 1980. The radiance integrated along the line of sight as in-situ radiance (R) can be derived taking into account absorption by ozone and air. The onion peeling inversion method was used to derive the vertical radiance (R) profiles respectively. The values of R were determined in the solar azimuth. The solar elevation angles are chosen larger for the backscattering observation than for the forward scattering observation to deal with as similar illumination conditions as possible despite the Earth's sphericity.

Description: Microscopical investigation of volcanic ash collected from ground stations during Mount St. Helens eruptions reveal a distinctive bimodel size distribution with high concentrations of particle ranges at (1) 200-100 microns and (2) 20-0.1 microns. Close examination of individual particles shows that most larger ones are solidified magma particles of porous pumice with numerous gas bubbles in the interior and the smaller ones are all glassy fragments without any detectable gas bubbles. Elemental analysis demonstrates that the fine fragments all have a composition similar to that of the larger pumice particles. Laboratory experiments suggest that the formation of the fine fragments is by bursting of glassy bubbles from a partially solidified surface of a crystallizing molten magma particle. The production of gas bubbles is due to the release of absorbed gases in molten magma particles when solubility decreases during phase transition. Diffusion cloud chamber experiments strongly indicate that sub-micron volcanic fragments are highly hygroscopic and extremely active as cloud condensation nuclei. Ice crystals also are evidently formed on those fragments in a supercooled (-20 C) cloud chamber. It has been reported that charge generation from ocean volcanic eruptions is due to contact of molten lava with sea water. This seems to be insufficient to explain the observed rapid and intense lightning activities over Mount St. Helens eruptions. Therefore, a hypothesis is presented here that highly electrically charged fine solid fragments are ejected by bursting of gas bubbles from the surface of a crystallizing molten magma particles.

Description: The direct solar radiation transmission record at Mauna Loa, dating from 1958 to the present, revealed with remarkable precision the presence of stratospheric aerosol from volcanic activity. This record can be used to quantify the intensity of the stratospheric volcanic aerosol perturbation following a significant eruption in reference to the Agung event in 1963. The Mount St. Helens' stratospheric cloud was first detected by lidar at 18 km over Mauna Loa on 17 July. The atmospheric transmission was seen to decrease slightly after that time, but only a few tenths of 1 percent. Although it is still fairly early to draw a definite conclusion on the ultimate magnitude of the Mount St. Helens stratospheric aerosol from the Mauna Loa results, it can be stated that the stratospheric aerosol optical depth presently observed is comparable with that observed from Fuego which erupted in 1974. At Boulder, Colorado, the atmospheric debris from Mount St. Helens was observed by lidar on a number of occasions. Also, observations of the diffuse, total and direct transmission of solar radiation were made on June 3 and 4. The latter set of observations is useful for deriving information on the scattering properties of the volcanic cloud. The lidar and solar radiation data are presented and some of their special features are discussed.

Description: Samples of stratospheric aerosols collected with U-2 aircraft for several months following the first three major eruptions of Mount St. Helens were analyzed for ash and liquid acid content. Ash grain sizes and compositions vary depending on collection altitude, location within the drifting cloud, and days following their injection. s computers Size distributions of ash particles vary with altitude. Generally small particles are depleted more rapidly at low altitudes (12 km) than at higher altitudes (17-18 km). Although samples collected 1 day after the first eruption of May 18, 1980, were dry, flow marks on the aircraft indicated parts of the cloud contained heavy acid concentrations. Indeed, all other samples obtained within 1 to 4 days after later eruptions (May 25 and June 12, 1980) were covered with copious amounts of liquid acid. Proportions of liquid to ash varied considerably depending on sampling location and cloud age. Because the acid-coated ash globules were large, they rapidly fell from the stratosphere until, by late June 1980, only a residue of acid droplets remained. Size distributions and concentrations of these droplets varied considerably.

Description: The inversions of multi-channel solar extinction measurements have been analyzed for the 0.35-1.0 micron wavelength region to retrieve stratospheric aerosol and ozone vertical profiles using both the constrained linear inversion scheme and the iterative scheme. The inversions of the multi-wavelength solar extinction data obtained from spacecraft have been analyzed based on the inversion of computer simulated data using various atmospheric models with differing amounts of aerosol and ozone in the stratosphere. The sensitivities of the inversion schemes to different experimental errors are discussed in terms of accuracy and resolution of the retrieved profiles.

Description: Polarization properties of the angularly scattered laser light from a volume of air are used to determine the size distribution of the aerosol particles within the volume by the use of appropriate inversion techniques. Similar techniques are employed to determine a mean size distribution of the particulates within a vertical column through the atmosphere from determinations of the aerosol optical depth as a function of wavelength. In both of these examples, a modification of an inversion technique originally described by Twomey has been employed. Details of this method are presented as well as results from actual measurements employing bistatic lidar and solar radiometer.

Description: Surface-based observations of downwelling microwave thermal emission are related to temperature and humidity profiles via a standard integral equation of radiative transfer. Both in clear and in cloudy atmospheres, statistical inversion techniques are used to retrieve profiles from a data vector of brightness observations and surface meteorological constraints. For the clear case, accuracy predictions and profile retrievals are illustrated for: (1) single frequency angular scanned data; (2) multi-frequency angular scanned data; and (3) multi-frequency zenith data. For the last case predicted and achieved accuracies were compared in a recently conducted radiometric experiment. Retrievals of cloud contaminated radiometric data are elaborated.

Description: Linear numerical inversion methods applied to atmospheric remote sounding generally can be categorized in two ways: (1) iterative, and (2) inverse matrix methods. However, these two categories are not unrelated; a duality exists between them. In other words, given an iterative scheme, a corresponding inverse matrix method exists, and conversely. This duality concept is developed for the more familiar linear methods. The iterative duals are compared with the classical linear iterative approaches and their differences analyzed. The importance of the initial profile in all methods is stressed. Calculations using simulated data are made to compare accuracies and to examine the dependence of the solution on the initial profile.

Description: An analytic model approach is applied to several simple atmospheric inversion problems. This method gives a sharp determination of aerosol size distribution parameters. It is shown that this analytic approach, together with ground level point sampling data measurements, can be used to infer information on the tropospheric ozone profile.

Description: Limb emission measurements are characterized by sharp weighting functions at high altitudes, and for temperature determinations, strongly nonlinear dependence of the weighting function on the temperature. Several methods for inverting this type of measurement have been described and used, including iterative, statistical, nonlinear and approximate direct approaches. These approaches are described and advantages and disadvantages of each are outlined.

Description: Nonlinear matrix inversion operators have been developed which, applied to observed radiances, infer maximal information regarding atmospheric scattering parameters and vertical distribution of radiant sources and sinks. The algorithm has the attractive feature of noise discrimination, attributing instrumental errors to extra-atmospheric sources.

Description: The inversion method provides a quantitative evaluation of the trade-off between vertical resolution of a retrieved profile and formal root-mean-square (rms) error due to measurement noise propagation. The problem of retrieving the top-side ozone profile from backscattered ultraviolet (BUV) measurements is considered. For measurements of the type currently being obtained with the Nimbus 4 and AE-E BUV experiments, it is found that a vertical resolution of approximately 0.75 scale height can be achieved for a formal volume mixing ratio profile error of 10%. Other examples include treatments of the retrieval of temperature profiles from measurements in the 15 micron CO2 absorption band for both the terrestrial and Martian atmospheres. Finally, the method is applied to the problem of retrieving temperature profiles of the Jovian planets from measurements in the far infrared pressure induced H2 lines to be obtained from the Mariner Jupiter/Saturn fly-by missions.

Description: Statistical methods are used to deal with the inverse problem of radiative transfer. All the available information about an unknown profile can be expressed in the form of values of functions of that profile and error estimates of these values. Estimation theory shows how these values are combined to give an estimate of the unknown profile and its error covariance. Many inversion methods are expressed in this form, although the error estimate is not usually carried out. Practical applications are described, both for inversion of individual profiles, and the global analysis of satellite data.

Description: Several commonly used methods for inversion--constrained linear inversion, synthesis (Backus-Gilbert) methods and nonlinear iterative techniques for the Chahine type--are discussed. It is demonstrated that a very close connection exists between Backus-Gilbert solutions and those given by constrained linear inversion. A number of examples of the application of such methods are presented, showing that resolution is not greatly different for quite different algorithms, a result quite in accord with general theoretical considerations: more resolution can be achieved at the expense of introducing greater a priori bias in the procedure.

Description: The radiative transfer in a scattering plane-parallel atmosphere is discussed, considering the exact analytical, the computational and the approximate methods. Some results of numerical comparisons are given. Finally, the difficulties of realistic atmospheric models are emphasized.

Description: Multiple scattering problems in a plane layer often permit the convenient use of different methods joined together. Sample numerical results to illustrate this point refer to X- and Y-functions, asymptotic fitting, the small-loss approximations, polarization in high orders, and photon path distribution.

Description: An instrument capable of observing the natural electron flux in the energy range from 0.1 to 12.0 kiloelectron volts is discussed for use in an experiment intended as a forerunner of a method that will utilize artificially accelerated electrons as tracer particles for electron fields parallel to the magnetic field. Effects that are of importance either as means of detecting the echo beam or as causes of beam perturbations (e.g., spacecraft charging effects and electron background) are to be studied. The use of electron accelerators as a tool to probe magnetospheric processes rather than to modify them is planned.

Description: A dual-channel video system mounted on a stabilized two-axis gimbal system (mounted on the pallet) with associated optics and data handling electronics described the low light flux observations are required for: (1) investigating ionospheric transport processes by observing Mg+ ions; (2) supporting magnetospheric electron bounce experiments; (3) measuring electron cross sections for selected atmospheric species; (4) detecting small particle contamination; and (5) studying natural auroras.

Description: A magnetometer experiment was designed to determine the local magnetic field by measuring the total of the Earth's magnetic field and that of an unknown spacecraft. The measured field vector components are available to all onboard experiments via the Spacelab command and data management system. The experiment consists of two parts, an electronic box and the magnetic field sensor. The sensor includes three independent measuring flux-gate magnetometers, each measuring one component. The physical background is the nonlinearity of the B-H curve of a ferrite material. Two coils wound around a ferrite rod are necessary. One of them, a tank coil, pumps the ferrite rod at approximately 20 kilohertz. As a consequence of the nonlinearity, many harmonics can be produced. The second coil (i.e., the detection coil) resonates to the first harmonic. If an unknown dc or low-frequency magnetic field exists, the amplitude of the first harmonic is a measure for the unknown magnetic field. The voltages detected by the sensors are to be digitized and transferred to the command and data management system.

Description: Study of sources of Lyman-alpha emission in the atmosphere, in the interplanetary medium, and perhaps in the galactic medium is planned. Sources of Lyman-alpha emission are described and a schematic of the instrument presented.

Description: The purpose of space experiments with particle accelerators (SEPAC) is to carry out active and interactive experiments on and in the Earth's ionosphere and magnetosphere. It is also intended to make an initial performance test for an overall program of Spacelab/SEPAC experiments. The instruments to be used are an electron beam accelerator, magnetoplasma dynamic arcjet, and associated diagnostic equipment. The accelerators are installed on the pallet, with monitoring and diagnostic observations being made by the gas plume release, beam-monitor TV, and particle-wave measuring instruments also mounted on the pallet. Command and display systems are installed in the module. Three major classes of investigations to be performed are vehicle charge neutralization, beam plasma physics, and beam atmosphere interactions. The first two are mainly onboard plasma physics experiments to measure the effect of phenomena in the vicinity of Spacelab. The last one is concerned with atmospheric modification and is supported by other Spacelab 1 investigations as well as by ground-based, remote sensing observations.

Description: Significant new geologic information has been revealed by comparing 1:5 million scale geologic maps of the equatorial zone quadrangles of Mercury (H-6, H-7 and H-8) to Earth-based elevation profiles and surface reflectivity maps of Mercury obtained in the early 1970's at the Arecibo (PR) and Goldstone (CA) radar facilities. These data consist of 23 Goldstone images and profiles of polarized return data at 12.5-cm wavelength and one Arecibo profile. Radar data with 150-m vertical accuracy and 10- to 20-km horizontal resolution are available for areas between latitudes 13 N. and 11 S. In general, these data sets show excellent correlation between: (1) relative elevation and roughness differences that are reflected by mapped geologic contacts; (2) mapped ridges and scarps that display distinctive radar signatures; and (3) position and morphology of crater-and-basin topographic elements. Inferences can also be drawn about topographic and geologic terrain beyond the area imaged by Mariner 10 cameras.

Description: On Mars, the association of gullied escarpments and chaotic terrain is evidence for failure and scarp retreat of poorly consolidated materials. Some martian gullies have no surface outlets and may have drained through subterranean channels. Similar features, though on a much smaller scale, can be seen in alluvium along terrestrial river banks in semiarid regions, such as the Rio Puerco Valley of central New Mexico. Many of the escarpments along the Rio Puerco are developing through formation of collapse gullies, which drain through soil pipes. Gully development can be monitored on aerial photographs taken in 1935, 1962, and 1980. A regression model was developed to quantify gully evolution over a known time span. Soil pipes and their associated collapse gullies make recognizable signatures on the air photos. The areal extent of this signature can be normalized to the scarp length of each pipe-gully system, which makes comparisons between systems possible.

Description: Geomorphic studies of impact structures in central Australia are being used to understand the complexities of fluvial dissection in the heavily cratered terrains of Mars. At Henbury, Northern Territory, approximately 12 small meteorite craters have interacted with a semiarid drainage system. The detailed mapping of the geologic and structural features at Henbury allowed this study to concentrate on degradational landforms. The breaching of crater rims by gullies was facilitated by the northward movement of sheetwash along an extensive pediment surface extending from the Bacon Range. South-facing crater rims have been preferentially breached because gullies on those sides were able to tap the largest amounts of runoff. At crater 6 a probable rim-gully system has captured the headward reaches of a pre-impact stream channel. The interactive history of impacts and drainage development is critical to understanding the relationships in the heavily cratered uplands of Mars. Whereas Henbury craters are younger than 4700 yrs. B.P., the Gosses Bluff structure formed about 130 million years ago. The bluff is essentially an etched central peak composed of resistant sandstone units. Fluvial erosion of this structure is also discussed.

Description: Eastern Acidalia Planitia contains a wide variety of terrain types on which are thousands of subkilometer volcanoes. Apparent morphometric variations were previously reported as a function of terrain type for the cones in the Cydonia area and extended to the rest of Acidalia for which high resolution Viking imagery exist. Crater counts are included for the six types of plains identified, density distributions of subkilometer cones found on each type of terrain, and orphometric data by morphological subclass as a function of terrain for more than 1400 cones.

Description: Preliminary mapping shows East Butte to be a single, large cumulo-dome composed dominantly of rhyolite which can be classified into three main groups based on color and structure. The rhyolite of East Butte is aphanitic with phenocrysts of sanidine and quartz which vary from 1 to 5 mm in length. Vesicular reddish black inclusions of basalt up to 10 cm in length, found in all varieties of the East Butte rhyolites are believed to have originated from fragmentation of the basalt walls of the conduit by rhyolitic magma as it was emplaced. Most of the inclusions contain plagioclase phenocrysts. These phenocrysts measure up to 1 to 2 cm in length and have a typical euhedral, tabular habit. A 250-m diameter depression which has the appearance of a crater is located at the top of East Butte. Evidence supporting the fact that the depression is a crater is displayed by three small (3 to 5 m in height) mounds of massive rhyolite which border the depression.

Description: Some of the geological relationships observed in the Mauna Loa sulfur flow may apply in considering volcanic processes on Io. Given the presence of sulfur/sulfur compounds in the eruption plumes and on the surface of Io, it is likely that extensive secondary deposits of sulfur exist, some of which may be of fumarolic origin and analogous to the Mauna Loa deposit. Given the likelihood of silicate volcanism of Io based on the inferred material properties of some flows, and the attendant high temperatures for silicate volcanism, it is likely that the secondary surface deposits of sulfur would have been mobilized without being heated to the high viscosity stage. Mobilized sulfur flows on Io may flow long distances as a result of: (1) low viscosities in the melting range; (2) sustained effusion resulting from continued heating source area; (3) continued remobilization within the flow as a consequence of surges from the source; and (4) extension via lava tubes, or similar conduits through which there is little heat loss. Sulfur flows may form a relatively thin veneer over silicate flows and other surface units, given their fluidity and low mobilization temperature. Active splashing and splattering may spread sulfur over a wider area contributing the bright blooms observed in association with some Ionian flows.

Description: Finite difference procedures were successfully used to solve the steady transonic flow about airfoils and appear to provide a practical means for calculating the corresponding unsteady flow. The purpose of the paper is to describe a finite difference procedure derived from the equations for the potential flow by assuming small perturbations and harmonic motion. The velocity potential is divided into steady and unsteady parts, and the resulting unsteady equation is linearized on the basis of small amplitudes of oscillation. The steady velocity potential, which must be calculated first, is described by the classical nonlinear transonic differential equation.

Description: An in-flight wing wake section drag investigation was conducted using traversing pitot and static probes. The primary objective was to develop measurement techniques and improve the accuracy of in-flight wing profile drag measurements for low values of dynamic pressure and Reynolds number. Data were obtained on a sailplane for speeds from about 40 knots to 125 knots at chord Reynolds numbers between 1,000,000 and 3,000,000. Tests were conducted with zero flap deflection, deflected flaps, and various degrees of surface roughness, and for smooth and rough atmospheric conditions. Several techniques were used to increase data reliability and to minimize certain bias errors. A discussion of the effects of a total pressure probe in a pressure gradient, and the effects of discrete turbulence levels, on the data presented and other experimental results is also included.

Description: A method for calculating the transonic flow over steady and oscillating airfoils was developed by Isogai. It solves the full potential equation with a semi-implicit, time-marching, finite difference technique. Steady flow solutions are obtained from time asymptotic solutions for a steady airfoil. Corresponding oscillatory solutions are obtained by initiating an oscillation and marching in time for several cycles until a converged periodic solution is achieved. In this paper the method is described in general terms, and results are compared with experimental data for both steady flow and for oscillations at several values of reduced frequency. Good agreement for static pressures is shown for subcritical speeds, with increasing deviation as Mach number is increased into the supercritical speed range. Fair agreement with experiment was obtained at high reduced frequencies with larger deviations at low reduced frequencies.

Description: The uses of laser Doppler velocimeter, hot wire, and surface hot film techniques in the study of turbulent flows are described, and data obtained in compressible flows are discussed. Applications are illustrated with measurements of wind tunnel freestream turbulence characteristics and with data obtained in transitional, turbulent, and separated shear flows. A new method which was developed for the study of time dependent and unsteady turbulent flows is also presented.

Description: Dynamic model verification is the process whereby an analytical model of a dynamic system is compared with experimental data, adjusted if necessary to bring it into agreement with the data, and then qualified for future use in predicting system response in a different dynamic environment. These are various ways to conduct model verification. The approach taken here employs Bayesian statistical parameter estimation. Unlike curve fitting, whose objective is to minimize the difference between some analytical function and a given quantity of test data (or curve), Bayesian estimation attempts also to minimize the difference between the parameter values of that funciton (the model) and their initial estimates, in a least squares sense. The objectives of dynamic model verification, therefore, are to produce a model which: (1) is in agreement with test data; (2) will assist in the interpretation of test data; (3) can be used to help verify a design; (4) will reliably predict performance; and (5) in the case of space structures, will facilitate dynamic control.

Description: Multidisciplinary analysis often requires optimization of nonlinear systems that are subject to constraints. Trajectory optimization is one example of this situation. The Program to Optimize Simulated Trajectories (POST) was used successfully for a number of problems. The purpose is to describe POST and a new optimization approach that has been incorporated into it. Typical uses of POST will also be illustrated. The projected-gradient approach to optimization is the preferred option in POST and is discussed. A new approach to optimization, the random-walk approach, is described, and results with the random-walk approach are presented.

Description: The purpose is not to provide a detailed discussion of several wall interference experiments, but rather to use these experiments (recently accomplished in the Boeing Transonic Wind Tunnel (BTWT) to illustrate the problems associated with many of the measurements required by current wall interference assessment/correction (WIAC) procedures. The wall correction to lift is emphasized. It is shown that, because conventional tunnels and relatively small models continue to be used, the flow field or flow boundary measurements to be made impose severe requirements on the experiment itself. In some cases, existing instrumentation and test techniques may not be adequate to obtain the data accuracies needed.

Description: Based upon limited, initial observations of wall interference corrections obtained for one airfoil test, there is a need for assessing the upstream flow direction. If there is no direct measurement then a two-pass correction procedure similar to the one described here is required. Questions have arisen pertaining to the correct interpretation of the pressure coefficients measured on the slats of a slotted tunnel wall, the interpretation of just what the calculated equivalent body encompasses or should include, and what can or should be considered as quantitative criteria for data correctability. Further studies using this modified procedure will address these questions. Hopefully, a meaningful WIAC procedure can be validated for the airfoil tests in the 0.3-m TCT.