ALBERT

Description: Svante Arrhenius' research in atmospheric physics extended beyond the recent past and the near future states of the Earth, which today are at the center of sociopolitical attention. His plan encompassed all of the physical phenomena known at the time to relate to the formation and evolution of stars and planets. His two-volume textbook on cosmic physics is a comprehensive synopsis of the field. The inquiry into the possible cause of the ice ages and the theory of selective wavelength filter control led Arrhenius to consider the surface states of the other terrestrial planets, and of the ancient Earth before it had been modified by the emergence of life. The rapid escape of hydrogen and the equilibration with igneous rocks required that carbon in the early atmosphere prevailed mainly in oxidized form as carbon dioxide, together with other photoactive gases exerting a greenhouse effect orders of magnitude larger than in our present atmosphere. This effect, together with the ensuing chemical processes, would have set the conditions for life to evolve on our planet, seeded from spores spreading through an infinite Universe, and propelled, as Arrhenius thought, by stellar radiation pressure.

Description: The NASA/Goddard Space Flight Center Scanning Raman Lidar (SRL) had a highly successful deployment at the Department of Energy Cloud and Radiation Testbed (CART) Site in Billings, OK during April, 1994 for the first Intensive Operation Period (IOP) hosted there. During the IOP, the SRL operated from just after sundown to just before sunrise for all declared evenings of operation. The lidar acquired more than 123 hours of data over 15 nights with less than 1 hour of data lost due to minor system malfunction. The SRL acquired data both on the vertical and in scanning mode toward an instrumented 60 m tower during various meteorological conditions such as an intense cold frontal passage on April 15 which is the focus of this presentation.

Description: The first Atmospheric Radiation Measurement (ARM) Remote Cloud Study (RCS) Intensive Operations Period (IOP) was held during April 1994 at the Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) site near Lamont, Oklahoma. This experiment was conducted to evaluate and calibrate state-of-the-art, ground based remote sensing instruments and to use the data acquired by these instruments to validate retrieval algorithms developed under the ARM program. These activities are part of an overall plan to assess general circulation model (GCM) parameterization research. Since radiation processes are one of the key areas included in this parameterization research, measurements of water vapor and aerosols are required because of the important roles these atmospheric constituents play in radiative transfer. Two instruments were deployed during this IOP to measure water vapor and aerosols and study their relationship. The NASA/Goddard Space Flight Center (GSFC) Scanning Raman Lidar (SRL) acquired water vapor and aerosol profile data during 15 nights of operations. The lidar acquired vertical profiles as well as nearly horizontal profiles directed near an instrumented 60 meter tower. Aerosol optical thickness, phase function, size distribution, and integrated water vapor were derived from measurements with a multiband automatic sun and sky scanning radiometer deployed at this site.

Description: The Airborne Visible-Infrared Imaging Spectrometer (AVIRIS) measures the total upwelling spectral radiance from 400 to 2500 nm sampled at 10 nm intervals. The instrument acquires spectral data at an altitude of 20 km above sea level, as images of 11 by up to 100 km at 17x17 meter spatial sampling. We have developed a nonlinear spectral fitting algorithm coupled with a radiative transfer code to derive the total path water vapor from the spectrum, measured for each spatial element in an AVIRIS image. The algorithm compensates for variation in the surface spectral reflectance and atmospheric aerosols. It uses water vapor absorption bands centered at 940 nm, 1040 nm, and 1380 nm. We analyze data sets with water vapor abundances ranging from 1 to 40 perceptible millimeters. In one data set, the total path water vapor varies from 7 to 21 mm over a distance of less than 10 km. We have analyzed a time series of five images acquired at 12 minute intervals; these show spatially heterogeneous changes of advocated water vapor of 25 percent over 1 hour. The algorithm determines water vapor for images with a range of ground covers, including bare rock and soil, sparse to dense vegetation, snow and ice, open water, and clouds. The precision of the water vapor determination approaches one percent. However, the precision is sensitive to the absolute abundance and the absorption strength of the atmospheric water vapor band analyzed. We have evaluated the accuracy of the algorithm by comparing several surface-based determinations of water vapor at the time of the AVIRIS data acquisition. The agreement between the AVIRIS measured water vapor and the in situ surface radiometer and surface interferometer measured water vapor is 5 to 10 percent.

Description: In the DIAL technique, the water vapor concentration profile is determined by analyzing the lidar backscatter signals for laser wavelengths tuned 'on' and 'off' a water vapor absorption line. Desired characteristics of the on-line transmitted laser beam include: pulse energy greater than or equal to 100 mJ, high-resolution tuning capability (uncertainty less than 0.25 pm), good spectral stability (jitter less than 0.5 pm about the mean), and high spectral purity (greater than 99 percent). The off-line laser is generally detuned less than 100 pm away from the water vapor line. Its spectral requirements are much less stringent. In our past research, we developed and demonstrated the airborne DIAL technique for water vapor measurements in the 720-nm spectral region using a system based on an alexandrite laser as the transmitter for the on-line wavelength and a Nd:YAG laser-pumped dye laser for the off-line wavelength. This off-line laser has been replaced by a second alexandrite laser. Diode lasers are used to injection seed both lasers for frequency and linewidth control. This eliminates the need for the two intracavity etalons utilized in our previous alexandrite laser and thereby greatly reduces the risk of optical damage. Consequently, the transmitted pulse energy can be substantially increased, resulting in greater measurement range, higher data density, and increased measurement precision. In this paper, we describe the diode injection seed source, the two alexandrite lasers, and the device used to line lock the on-line seed source to the water vapor absorption feature.

Description: The differential absorption lidar (DIAL) technique was first applied to the remote measurement of atmospheric water vapor profiles from airborne platforms in 1981. The successful interpretation of the lidar profiles relies strongly on an accurate knowledge of specific water vapor absorption line parameters: line strength, pressure broadening coefficient, pressure-induced shift coefficient and the respective temperature-dependence factors. NASA Langley Research Center has developed and is currently testing an autonomous airborne water vapor lidar system: LASE (Lidar Atmospheric Sensing Experiment). This DIAL system uses a Nd:YAG-pumped Ti:Sapphire laser seeded by a diode laser as a lidar transmitter. The tunable diode has been selected to operate in the 813-818 nm wavelength region. This 5-nm spectral interval offers a large distribution of strengths for temperature-insensitive water vapor absorption lines. In support of the LASE project, a series of spectroscopic measurements were conducted for the 16 absorption lines that have been identified for use in the LASE measurements. Prior to this work, the experimental data for this water vapor absorption band were limited - to our knowledge - to the line strengths and to the line positions.

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 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: Antenna temperatures and the corresponding geolocation data from the five sources of the Special Sensor Microwave/Imager data from the Defense Meteorological Satellite Program F11 satellite have been characterized. Data from the Fleet Numerical Meteorology and Oceanography Center (FNMOC) have been compared with data from other sources to define and document the differences resulting from different processing systems. While all sources used similar methods to calculate antenna temperatures, different calibration averaging techniques and other processing methods yielded temperature differences. Analyses of the geolocation data identified perturbations in the FNMOC and National Environmental Satellite, Data and Information Service data. The effects of the temperature differences were examined by generating rain rates using the Goddard Scattering Algorithm. Differences in the geophysical precipitation products are directly attributable to antenna temperature differences.

Description: The edge technique is a new and powerful method for measuring small frequency shifts such as the Doppler shift of an atmospheric backscattered signal from a pulsed laser. The edge technique can be used for high spatial resolution, high accuracy ground and airborne wind measurements as well as high accuracy spaceborne wind measurements. We have recently made our first ground based wind measurements. These have a spatial resolution of 15 m and an accuracy of 25 cm/s and these measurements are presented in this paper. This is a unique capability and provides valuable information for studies of turbulent processes in the lower atmosphere. It could also be used for high sensitivity detection of wind shear and microbursts in the vicinity of airports. In addition, global wind measurements can be made with the edge technique from space with an accuracy of 1 m/s and a vertical resolution as high as 150 m in the boundary layer and 1 km through the troposphere. Such a system could make eyesafe wind measurements using well developed diode pumped solid state laser technology at 1.06 micron. Multi-pulse averaging would provide a spatially representative wind measurement.

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: The composition of the jovian atmosphere from 0.5 to 21 bars along the descent trajectory was determined by a quadrupole mass spectrometer on the Galileo probe. The mixing ratio of He (helium) to H2 (hydrogen), 0.156, is close to the solar ratio. The abundances of methane, water, argon, neon, and hydrogen sulfide were measured; krypton and xenon were detected. As measured in the jovian atmosphere, the amount of carbon is 2.9 times the solar abundance relative to H2, the amount of sulfur is greater than the solar abundance, and the amount of oxygen is much less than the solar abundance. The neon abundance compared with that of hydrogen is about an order of magnitude less than the solar abundance. Isotopic ratios of carbon and the noble gases are consistent with solar values. The measured ratio of deuterium to hydrogen (D/H) of (5 +/- 2) x 10(-5) indicates that this ratio is greater in solar-system hydrogen than in local interstellar hydrogen, and the 3He/4He ratio of (1.1 +/- 0.2) x 10(-4) provides a new value for protosolar (solar nebula) helium isotopes. Together, the D/H and 3He/4He ratios are consistent with conversion in the sun of protosolar deuterium to present-day 3He.

Description: We propose an index of climate change based on practical climate indicators such as heating degree days and the frequency of intense precipitation. We find that in most regions the index is positive, the sense predicted to accompany global warming. In a few regions, especially in Asia and western North America, the index indicates that climate change should be apparent already, but in most places climate trends are too small to stand out above year-to-year variability. The climate index is strongly correlated with global surface temperature, which has increased as rapidly as projected by climate models in the 1980s. We argue that the global area with obvious climate change will increase notably in the next few years. But we show that the growth rate of greenhouse gas climate forcing has declined in recent years, and thus there is an opportunity to keep climate change in the 21st century less than "business-as-usual" scenarios.

Description: The addition of the "dioxin-like" polychlorinated biphenyl (PCB) congeners to the assessment of risk associated with the 2,3,7,8-chlorine substituted dioxins and furans has dramatically increased the number of laboratories worldwide that are developing analytical procedures for their detection and quantitation. Most of these procedures are based on established sample preparation and analytical techniques employing high resolution gas chromatography/high resolution mass spectrometry (HRGC/HRMS), which are used for the analyses of dioxin/furans at low parts-per-trillion (ppt) levels. A significant and widespread problem that arises when using these sample preparation procedures for the analysis of coplanar PCBs is the presence of background levels of these congeners. Industrial processes, urban incineration, leaking electrical transformers, hazardous waste accidents, and improper waste disposal practices have released appreciable quantities of PCBs into the environment. This contamination has resulted in the global distribution of these compounds via the atmosphere and their ubiquitous presence in ambient air. The background presence of these compounds in method blanks must be addressed when determining the exact concentrations of these and other congeners in environmental samples. In this study reliable procedures were developed to accurately define these background levels and assess their variability over the course of the study. The background subtraction procedures developed and employed increase the probability that the values reported accurately represent the concentrations found in the samples and were not biased due to this background contamination.

Description: The composition of the Earth's early atmosphere is a subject of continuing debate. In particular, it has been suggested that elevated concentrations of atmospheric carbon dioxide would have been necessary to maintain normal surface temperatures in the face of lower solar luminosity in early Earth history. Fossil weathering profiles, known as palaeosols, have provided semi-quantitative constraints on atmospheric oxygen partial pressure (pO2) before 2.2 Gyr ago. Here we use the same well studied palaeosols to constrain atmospheric pCO2 between 2.75 and 2.2 Gyr ago. The observation that iron lost from the tops of these profiles was reprecipitated lower down as iron silicate minerals, rather than as iron carbonate, indicates that atmospheric pCO2 must have been less than 10(-1.4) atm--about 100 times today's level of 360 p.p.m., and at least five times lower than that required in one-dimensional climate models to compensate for lower solar luminosity at 2.75 Gyr. Our results suggest that either the Earth's early climate was much more sensitive to increases in pCO2 than has been thought, or that one or more greenhouse gases other than CO2 contributed significantly to the atmosphere's radiative balance during the late Archaean and early Proterozoic eons.

Description: The evolution of non-photosynthetic sulphide-oxidizing bacteria was contemporaneous with a large shift in the isotopic composition of biogenic sedimentary sulphides between 0.64 and 1.05 billion years ago. Both events were probably driven by a rise in atmospheric oxygen concentrations to greater than 5-18% of present levels--a change that may also have triggered the evolution of animals.

Description: Now that extrasolar planets have been found, it is timely to ask whether some of them might be suitable for life. Climatic constraints on planetary habitability indicate that a reasonably wide habitable zone exists around main sequence stars with spectral types in the early-F to mid-K range. However, it has not been demonstrated that planets orbiting such stars would be habitable when biologically-damaging energetic radiation is also considered. The large amounts of UV radiation emitted by early-type stars have been suggested to pose a problem for evolving life in their vicinity. But one might also argue that the real problem lies with late-type stars, which emit proportionally less radiation at the short wavelengths (lambda 〈 200 nm) required to split O2 and initiate ozone formation. We show here that neither of these concerns is necessarily fatal to the evolution of advanced life: Earth-like planets orbiting F and K stars may well receive less harmful UV radiation at their surfaces than does the Earth itself.

Description: Through analysis of spectral imaging data acquired with the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) from an ER-2 aircraft at 20 km altitude during several field programs, it was found that narrow channels near the center of the strong 1.38-micron water vapor band are very sensitive in detecting thin cirrus clouds. Based on this observation from AVIRIS data, a channel centered at 1.375 microns with a width of 30 nm was selected for the Moderate Resolution Imaging Spectrometer (MODIS) for remote sensing of cirrus clouds from space. The sensitivity of the 1.375-micron MODIS channel to detect thin cirrus clouds during the day time is expected to be one to two orders of magnitude better than the current infrared emission techniques. As a result, a larger fraction of the satellite data will likely be identified as containing cirrus clouds. In order to make better studies of surface reflectance properties, thin cirrus effects must be removed from satellite images. We have developed an empirical approach for removing/correcting thin cirrus effects in the 0.4 - 1.0 micron region using channels near 1.375 microns. This algorithm will be incorporated into the present MODIS atmospheric correction algorithms for ocean color and land applications and will yield improved MODIS atmospheric aerosol, land surface, and ocean color products.

Description: Thunderstorms separate charge. Most places they lift positive charge or lower negative, a few places they lift negative or lower positive. The electrical generator is stronger in some parts of the cloud than in others. Our long term goal is to map this generator. Cloud physicists tell us that uncharged ice and water particles become charged by collision, and that the charge transferred depends on size, temperature and humidity. There is still some disagreement about exactly how the charge transferred depends on size, temperature, and humidity. In principle, if we knew this ice physics, and also knew the distribution of particles everywhere in the storm, and the winds everywhere and the temperature and humidity everywhere, then we could compute everywhere the electrical power of the thunderstorm generator. In practice it is difficult to know all these things, particularly the distribution of particles, so it is difficult to use real thunderstorms to falsify cloud electrification theories. We here take one small step towards computing that map of electrical generator power, by relating radar reflectivity profiles of 2000 storms to lightning flash rates of those storms. This small step by itself doesn't falsify any existing electrification theories; it merely places weak constraints on the relation of electric generator power to cloud ice.

Description: Two primary detection techniques (optical and RF) have a proven capability for detecting lightning from low earth orbit. However, the lightning processes that generate the optical and RF signals are vastly different providing significantly different information content from each sensor type. Because of the intervening ionosphere, low frequency RF components do not reach satellite altitudes. As a consequence, many of the processes associated with the major energy release of a lightning event (i.e. return strokes, k-changes, recoil streamers, etc), in all likelihood contribute little to the RF signal arriving at the satellite. The optical output from lighting, on the other hand, has been shown to be highly correlated with the energetic, charge-transferring processes mentioned above. On the down side, the optical energy, while essentially unaffected by the atmosphere once it emerges from the cloud, is heavily scattered within the cloud. While there is little absorption by the cloud, the great optical depth makes the total light energy emerging from the cloud to be dependent on where in the cloud the lightning occurred. Analyses suggest that when lightning is confined to the lowest regions of the cloud, the light is strongly attenuated and detection becomes problematic. Fortunately, the vast majority of lightning flashes are comprised of channels that propagate through the middle of the cloud and higher. These flashes produce bright signals at the top of a cloud and are readily detectable. Presently, we have two optical instruments in orbit. The Optical Transient Detector (OTD) has been orbiting the earth since April, 1995, while the Lightning Imaging Sensor (LIS) was launched on the Tropical Rainfall Measuring Mission (TRMM) in November of 1997. Both instruments are relatively small, solid state optical imagers, designed specifically to detect and locate lightning activity from low earth orbit with high detection efficiency and location accuracy.

Description: For the past century, scientists have made quantitative measurements of lightning discharges. In the process, they refined the definition of a lightning unit, or basic quantum of lightning, in order to base it on observable parameters. In this paper, we will use cluster analysis to derive a basic spatial and temporal definition or scale length for the unit of lightning. We will use data from three different systems that detected pulses from the same storm complex over Central Oklahoma during June, 1998. Since the different instruments detect lightning in different ways with different resolutions, there may not be a single definition of the unit of lightning that can be applied to all three systems. However, common components can be found since all instrumentation are detecting aspects of the same phenomenon.

Description: The Lidar In-Space Technology Experiment (LITE) mission has demonstrated the utility of spaceborne lidar in observing multilayer clouds and has provided a dataset showing the distribution of tropospheric clouds and aerosols. These unambiguous observations of the vertical distribution of clouds will allow improved verification of current cloud climatologies and GCM cloud parameterizations. Although there is now great interest in cloud profiling radar, operating in the mm-wave region, for the spacebased observation of cloud heights the results of the LITE mission have shown that satellite lidars can also make significant contributions in this area.

Description: The use of unmanned aerial vehicles (UAVs), a platform for investigations in an environment hostile to manned spacecraft, is discussed. A program which includes the use of UAVs coupled with ground-based measurements to conduct scientific studies on the electrical state of the atmosphere during electrically active periods is proposed. The radiating power from alternate current and transient components of the storm electrification was investigated.

Description: Progress, future plans and publications regarding the following objectives are presented: (1) Determine the vertical and horizontal spatial distribution of hydrometeors in precipitating clouds; (2) Measure the spatial distribution of liquid water and ice in the clouds; and (3) Measure and determine the limits of measurement of the polarization characteristics related to the shapes and orientations of hydrometeors in precipitating clouds.

Description: Wind and water vapor are two major factors driving the Earth's atmospheric circulation, and direct measurement of these factors is needed for better understanding of basic atmospheric science, weather forecasting, and climate studies. Coherent lidar has proved to be a valuable tool for Doppler profiling of wind fields, and differential absorption lidar (DIAL) has shown its effectiveness in profiling water vapor. These two lidar techniques are generally considered distinctly different, but this paper explores an experimental combination of the Doppler and DIAL techniques for measuring both wind and water vapor with an eye-safe wavelength based on a solid-state laser material. Researchers have analyzed and demonstrated coherent DIAL water vapor measurements at 10 micrometers wavelength based on CO2 lasers. The hope of the research presented here is that the 2 gm wavelength in a holmium or thulium-based laser may offer smaller packaging and more rugged operation that the CO2-based approach. Researchers have extensively modeled 2 um coherent lasers for water vapor profiling, but no published demonstration is known. Studies have also been made, and results published on the Doppler portion, of a Nd:YAG-based coherent DIAL operating at 1.12 micrometers. Eye-safety of the 1.12 micrometer wavelength may be a concern, whereas the longer 2 micrometer and 10 micrometer systems allow a high level of eyesafety.

Description: Severe storms often have high flash rates (in excess of one flash per second) and are dominated by intracloud lightning activity. In addition to the extraordinary flash rates, there is a second distinguishing lightning characteristic of severe storms that seems to be important. When the total lightning history is examined, one finds sudden increases in the lightning rate, which we refer to as lightning "jumps," that precede the occurrence of severe weather by ten or more minutes. These jumps are typically 30-60 flashes/min, and are easily identified as anomalously large derivatives in the flash rate. This relationship is associated with updraft intensification and updraft strength is an important factor in storm severity (through the accumulation of condensate aloft and the stretching of vorticity). In several cases, evidence for diminishment of midlevel rotation and the descent of angular momentum from aloft is present prior to the appearance of the surface tornado. Based on our experience with severe and tornadic storms in Central Florida, we believe the total lightning may augment the more traditional use of NEXRAD radars and storm spotters. However, a more rigorous relation of these jumps to storm kinematics is needed if we are to apply total lightning in a decision tree that leads to improved warning lead times and decreased false alarm rates.

Description: The Lightning Imaging Sensor (LIS) is a NASA Earth Observing System (EOS) instrument on the Tropical Rainfall Measuring Mission (TRMM) platform designed to acquire and investigate the distribution and variability of total lightning (i.e., cloud-to-ground and intracloud) between q35' in latitude. Since lightning is one of the responses of the atmosphere to thermodynamic and dynamic forcing, the LIS data is being used to detect deep convection without land-ocean bias, estimate the precipitation mass in the mixed phased region of thunderclouds, and differentiate storms with strong updrafts from those with weak vertical motion.

Description: The Optical Transient Detector (OTD) is a space-based instrument specifically designed to detect and locate lightning discharges (intracloud and cloud-to-ground) as it orbits the Earth. A statistical examination of OTD lightning data reveals that nearly 1.2 billion flashes occurred over the entire earth during the one year period from September 1995 through August 1996. This translates to an average of 37 lightning flashes occurring around the globe every second, which is well below the traditional estimate of 100 flashes per second. An average of 75% of the global lightning activity during the year occurs between 30' S and 30' N. An analysis of the annual lightning distribution reveals that an average of 82% of the lightning flashes occur over the continents and 18% over the oceans, which translates to an average land-ocean flash density ratio of nearly 11.

Description: The lightning frequency model developed by Baker, Christian and Latham (1995) has been refined and extended, in an effort to provide a more realistic framework from which to examine computationally the relationships that might exist between lightning frequency f (which is now being routinely measured from satellites, using NASA/MSFC devices) and a variety of cloud physical parameters, including precipitation rate, updraught speed and non-precipitating ice content. Model results indicate the existence of a simple relationship between lightning frequency f and the upward flux of ice crystals into the thunderstorm anvil. It follows that, for a particular situation, one can assign a specific mass of non-precipitating ice to an individual lightning stroke. Therefore it may prove possible - using satellite measurements of global lightning - to estimate the atmospheric loading of ice crystals in thunderstorm anvils: a parameter of climatological importance. Early results from this work are presented, together with further studies of the relationships between f and other thundercloud parameters.

Description: A charge transfer process during the collision of a riming graupel pellet and an ice-crystal at low temperature is proposed. During riming, the surface structure of graupel deviates from perfect crystalline structure. A concept of quasi-solid layer (QSL) formation on the surface is introduced. This QSL contains defects formed during riming. In absence of impurities, positively charged X-defect abundance is considered in the outer layer. These defects are assumed to be the charge carriers during the charge transfer process. Some part of the QSL is stripped off by the colliding ice crystals, which thereby gain some positive charge, leaving the graupel pellet negatively charged. With the proposed model, fC to pC of charge transfer is observed per collision. A transition temperature between -10 C to -15 C is also noted beyond which the QSL concept does not hold. This transition temperature is dependent on the bulk liquid water content of the cloud.

Description: It is well known that most tropical cyclones (TCs) that make landfall along the Gulf coast of the United States spawn at least a few tornadoes. Although most landfalling TCs generate fewer than a dozen such tornadoes, a small proportion produce large swarm outbreaks, with as many as 25 or more tornadoes. Usually, these major outbreaks occur in large, intense hurricane-strength TCs, but on 15-17 August 1994 Tropical Storm Beryl spun off 37 tornadoes along its path from the Florida panhandle through the mid-Atlantic states. Some 32 of these tornadoes occurred on 16 August 1994 from eastern Georgia to southern Virginia, with most of these taking place in South Carolina. Beryl's 37 tornadoes moved it into what was at that time fifth place historically in terms of TC tornado productivity. The Beryl outbreak is especially noteworthy in that at least three of the tornadoes achieved peak intensity of F3 on the Fujita damage intensity scale. Although no fatalities resulted from the Beryl outbreak, at least 50 persons suffered injuries, and property damages totalled more than $50 million . The Beryl outbreak is a good example of a TC whose greatest danger to the public is its post-landfall severe weather. In this respect, and in the character of its swarm outbreak of tornadoes, it resembles another large tornado outbreak spawned by a relatively weak TC, Hurricane Danny of 1985). In the Danny outbreak, numerous shallow mini-supercell storms were found to have occurred, and it was noted that, because of the storms' relatively shallow depth, cloud-to-ground (CG) lightning was negligible. Better observations of future TC tornado outbreaks, especially with modern surveillance tools such as Doppler radars and the National Lightning Detection Network (NLDN), were recommended. Although the Beryl tornado outbreak is not the first set of TC-spawned tornado storms to be observed with the NLDN, it is one of the largest and likely the most intense such outbreak. The purpose of this paper is to document the NLDN-derived CG lightning characteristics of Beryl's tornadic storms, and to see how they compare with observations of CG lightning activity in other types of severe storms. In particular, we attempt to quantify the CG flash rates of TC tornadic cells, and to discover if there are any characteristics of their CG activity that may be useful to operational forecasters seeking to distinguish which cells are most likely to produce severe weather.

Description: The problem of retrieving ligntning, ground-strike location on a spherical Earth surface using a network of 4 or more time-of-arrival (TOA) sensors is considered, It is shown that this problem has an analytic solution and therefore does not require the use of nonlinear estimation theory (e.g., minimization). The mathematical robustness of the analytic solution is tested using computer-generated lightning sources and simulated TOA measurement errors. A summary of a quasi-analytic extension of the spherical Earth solution to an oblate spheroid Earth geometry is also provided.

Description: This study summarizes the results of an analysis of data from the LIS instrument on the TRMM platform. The data for the Indian summer monsoon season is examined to study the seasonal patterns of the geographic and diurnal distribution of lightning storms. The storms on the Tibetan plateau show a single large diurnal peak at about 1400 local solar time. A region of Northern Pakistan has two storm peaks at 0200 and 1400 local solar time. The morning peak is half the magnitude of the afternoon peak. The region south of the Himalayan Mountains has a combined diurnal cycle in location and time of storm occurrence.

Description: We have developed a new set of eight electric field mills that were flown on a NASA ER-2 high-altitude aircraft. During the Third Convection And Moisture EXperiment (CAMEX-3; Fall, 1998), measurements of electric field, storm dynamics, and ice microphysics were made over several hurricanes. Concurrently, the TExas-FLorida UNderflights (TEFLUN) program was being conducted to make the same measurements over Gulf Coast thunderstorms. Sample measurements are shown: typical flight altitude is 20km. Our new mills have an internal 16-bit A/D, with a resolution of 0.25V/m per bit at high gain, with a noise level less than the least significant bit. A second, lower gain channel gives us the ability to measure fields as high as 150 kV/m.

Description: The mapping of the lightning optical pulses detected by the Lightning Imaging Sensor (LIS) is compared with the radiation sources by Lightning Detection and Ranging (LDAR) and the National Lightning Detection Network (NLDN) for three thunderstorms observed during and overpasses on 15 August 1998. The comparison involves 122 flashes including 42 ground and 80 cloud flashes. For ground flash, the LIS recorded the subsequent strokes and changes inside the cloud. For cloud flashes, LIS recorded those with higher sources in altitude and larger number of sources. The discrepancies between the LIS and LDAR flash locations are about 4.3 km for cloud flashes and 12.2 km for ground flashes. The reason for these differences remain a mystery.

Description: Since April 1995, lightning activity around the globe has been monitored with the Optical Transient Detector (OTD). The OTD observations acquired during the one year period from September 1995 through August 1996 have been used to statistically determine the number of flashes that occur over the Earth during each hour of the diurnal cycle, expressed both as a function of local time and universal time. The globally averaged local [il,htnina activity displays a peak in late afternoon (1500-1800 local time) and a minimum in the morning hours (0600- 1000 local time) consistent with convection associated with diurnal heating. No diurnal variation is found for oceanic storms. The diurnal lightning distribution (universal time) for the globe displays a variation of about 35% about its mean as compared to the Carnegie curve which has a variation of only 15% above and below the mean.

Description: Total lightning observations made by the Optical Transient Detector (OTD) of a tornadic thunderstorm that occurred over Oklahoma on 17 April 1995 are presented. The average flash rate of the tornadic storm during the 3.2 min observation period was 45 flashes/min, with a flash rate density of 1.16 x 10(exp -4)/s sq km. The total flash rate was almost 18 times higher than the cloud-to-ground rate measured by the National Lightning Detection Network (NLDN). In addition, total lightning rates were observed to decrease prior to tornadic development.

Description: In recent years, atmospheric conductivity and electric field measurements over thunderstorms have been made at 20 km with a high altitude aircraft. After compensating for the effects of aircraft charging induced by external electric fields no significant variations in ambient conductivity above thunderstorms have been found. These Gerdien results contrast strongly with the large (and frequent) conductivity variations reported in studies using relaxation probe techniques.

Description: The El Nino Southern Oscillation (ENSO) is a climate anomaly responsible for world-wide weather impacts ranging from droughts to floods. In the United States, warm episode years are known to produce above normal rainfall along the Southeast US Gulf Coast and into the Gulf of Mexico, with the greatest response observed in the October-March period of the current warm-episode year. The 1997-98 warm episode, notable for being the strongest event since 1982-83, presents our first opportunity to examine the response to a major ENSO event and determine the variation of wintertime thunderstorm activity in this part of the world. Due to the recent launch of a lightning sensor on NASA's Tropical Rainfall Measuring Mission (TRMM) in November 1997 and the expanded coverage of the National Lightning Detection Network (NLDN), we are able to examine such year-to-year changes in lightning activity with far greater detail than ever before.

Description: The double-edge lidar technique for measuring the wind using molecular backscatter is described. Two high spectral resolution edge filters are located in the wings of the Rayleigh-Brillouin profile. This doubles the signal change per unit Doppler shift, the sensitivity, and gives nearly a factor of two improvement in measurement accuracy. The use of a crossover region is described where the sensitivity of a molecular and aerosol-based measurement are equal. This desensitizes the molecular measurement to the effects of aerosol scattering over a frequency range of +/- 100 m/s. We give methods for correcting for short-term frequency jitter and drift using a laser reference frequency measurement and methods for long-term frequency correction using a servo control system. The effects of Rayleigh-Brillouin scattering on the measurement are shown to be significant and are included in the analysis. Simulations for a conical scanning satellite-based lidar at 355 nm show an accuracy of 2-3 m/s for altitudes of 2 to 15 km for a 1 km vertical resolution, a satellite altitude of 400 km and a 200 km x 200 km spatial resolution. Results of ground based wind measurements are presented.

Description: The edge technique utilizes the edge of a high spectral resolution filter for high accuracy wind measurement using direct detection lidar. The signal is split between an edge filter channel and a broadband energy monitor channel. The energy monitor channel is used for signal normalization. The edge measurement is made as a differential frequency measurement between the outgoing laser signal and the atmospheric backscattered return for each pulse. As a result, the measurement is insensitive to laser and edge filter frequency jitter and drift at a level less than a few parts in 10(exp 10). We will discuss the methodology of the technique in detail, present a broad range of simulation results, and provide preprints of a journal article currently in press.

Description: Temperature measurements in the middle atmosphere using Rayleigh lidars have been performed for several decades now. The high accuracy and vertical resolution provided by lidars allow to study the temperature variability at various scales with high confidence levels. One of the numerous applications is the study of the middle atmospheric thermal tides. Although Rayleigh lidar measurements are basically possible only at nighttime, diurnal and semidiurnal components can often be extracted if the results are taken with care and correctly interpreted. Using results from more than 200 hours of nighttime measurements obtained by lidar in October 1996 and 1997 at Mauna Loa Observatory, Hawaii, a study of the middle atmospheric (25-90 km) thermal tides is presented in this paper. The amplitudes and phases of the diurnal and semidiurnal components were calculated for some altitudes where the fits converged significantly, and compared to that of the Global Scale Wave Model (GSWM).

Description: The primary constituents of the Earth's atmosphere are molecular nitrogen and molecular oxygen. Ozone is created when ultraviolet light from the sun photodissociates molecular oxygen into two oxygen atoms. The oxygen atoms undergo many collisions but eventually combine with a molecular oxygen to form ozone (O3). The ozone molecules absorb ultraviolet solar radiation, primarily in the wavelength region between 200 and 300 nanometers, resulting in the dissociation of ozone back into atomic oxygen and molecular oxygen. The oxygen atom reattaches to an O2 molecule, reforming ozone which can then absorb another ultraviolet photon. This sequence goes back and forth between atomic oxygen and ozone, each time absorbing a uv photon, until the oxygen atom collides with and ozone molecule to reform two oxygen molecules.

Description: We review the basic multiple scattering theory of off-beam lidar returns from optically thick clouds using the diffusion approximation. The shape of the temporal signal - the stretched pulse - depends primarily on the physical thickness of the cloud whereas its spatial counterpart - the diffuse spot - conveys specific information on the cloud's optical thickness, as do the absolute returns. This makes observation of the weak off-beam lidar returns an attractive prospect in remote sensing of cloud properties. By estimating the signal-to-noise ratio, we show that night-time measurements can be performed with existing technology. By the same criterion, day-time operation is a challenge that can only be met with a combination of cutting-edge techniques in filtering and in laser sources.

Description: In winter, active convective clouds frequently form along the coastline of the Hokuriku district, in association with strong advection of Siberian air masses over the Sea of Japan. On the other hand, in summer, many thunderclouds form in the Kanto region in the afternoon every day. Summer and winter thunderclouds were investigated by field works, operation of the C- and X-band weather radars and a car-borne fieldmill. The investigation found a very close relation between the temporal variation of 3-dimensional radar echo and surface electric field magnitude detected by a car-borne fieldmill in the case of summer thunderclouds and winter convective clouds or thunderclouds. The study probed the close relation among radar echoes, quantity of thunderclouds and surface electric field magnitude in the summer and winter seasons. We think that summer thundercloud activity can basically be equated with winter thundercloud lightning activity, except that the magnitude of surface electric field under summer thunderclouds in the case of the Kanto region cannot be equated with that under winter thunderclouds in the case of the Hokuriku district in winter.

Description: A six year record of optical observations of lightning-induced mesospheric transient luminous events (TLEs) is available from the Yucca Ridge Field Station (YRFS) near Ft. Collins, CO. Climatological analyses reveal sprites and elves occur in a variety of convective storm types, but principally mesoscale convective systems (MCSs) and squall lines. Severe supercell storms rarely produce TLEs, except during their dissipating stage. Few TLEs are observed during storms with radar echo areas 〈7,500 sq km. Above this size there is a modest correlation with radar areal coverage. A typical High Plains storm produces 45 TLEs over a 143 interval. Sprites and most elves are associated with +CGs. The probability of a TLE increases with peak current. In six storms, 5.1% of +CGs produced TLEs, the number increasing to 32% of +CGs with 〉75 kA and 52% of +CGs with 〉100 kA peak current.

Description: Significant differences are known to exist on a global scale between continental and oceanic total lightning regional flash rates, suggesting differences in the properties of convective storms in these regimes. Lightning properties observed by the Optical Transient Detector (OTD) and Lightning Imaging Sensor (LIS) over land and ocean are compared, limited to analysis over the tropics in order to simplify physical interpretation. We find that the mean flash rates of individual storms over tropical land only exceed those over ocean by a factor of 2 (far less than the observed differences in regional flash rates). However, the average nearest neighbor distance of continental thunderstorms is half that over oceans. Cloud-top lightning optical radiance in oceanic storms is also twice as large as over land, suggesting either more energetic flashes over the oceans or less intervening cloud particles.

Description: The danger of natural and triggered lightning significantly impacts space launch operations supported by the USAF. The lightning Launch Commit Criteria (LCC) are used by the USAF to avoid these lightning threats to space launches. This paper presents a brief overview of the LCC.

Description: Charged aerosol clouds allow to create strong electric fields which generate discharges. A character of appearance and development of the discharges is similar to thundercloud processes. That is why the charged aerosol clouds are used for experiments on studying of discharges in air. Experimental data on electromagnetic fields investigation produced by discharges between charged cloud and ground are considered.

Description: Scheduled for launch in 2001 as part of NASA's Earth Observing System (EOS), the Geoscience Laser Altimeter System (GLAS) will provide continuous laser sounding of the earth's atmosphere from space for the first time. From its polar orbit about 600 km above the surface, GLAS will employ a 40 Hz solid state laser operating at 1064 nm to measure topography to an accuracy of 10 cm. Simultaneously, the atmospheric channels (1064 and 532 nm) of GLAS will provide profiles of atmospheric backscatter from 40 km to the ground with 75 meter vertical resolution (Spinhirne and Palm, 1996). These measurements will give scientists an unprecedented global data set on the vertical structure of clouds and aerosols which will greatly aid research efforts aimed at understanding their effects on climate and their role in climate change (Hartman, 1994). To better understand and predict the performance of the GLAS atmospheric channels, a computer model was developed to simulate the type of signal that the instrument would likely produce. The model uses aircraft lidar data and provides realistic simulated GLAS data sets over large areas spanning a wide range of atmospheric conditions. These simulated GLAS datasets are invaluable for designing and testing algorithms for the retrieval of parameters such as cloud and aerosol layer height, optical depth and extinction cross section. This work is currently proceeding and in this paper we will present results of the cloud and aerosol detection algorithm with emphasis on the detection of Marine Atmospheric Boundary Layer (MABL) aerosol. In addition, we use a recently developed technique to ascertain the feasability of estimating MABL moisture and temperature structure from spaceborne systems such as GLAS.

Description: We have developed the theory for aerosol- and molecular-based lidar measurements of the wind using double edge versions of the edge technique. Aerosol-based wind measurements have been made at Goddard Space Flight Center and molecular-based wind measurements at the University of Geneva. We have demonstrated atmospheric measurements using these techniques for altitudes from 1 to more than 10 km. Measurement accuracies of better than 1.25 m/s have been obtained with integration times from 5 to 30 seconds. The measurements can be scaled to space and agree, within a factor of two, with satellite-based simulations of performance based on Poisson statistics.

Description: At mid and tropical latitudes, cirrus clouds are present more than 50% of the time in satellites observations. Due to their large spatial and temporal coverage, and associated low temperatures, cirrus clouds have a major influence on the Earth-Ocean-Atmosphere energy balance through their effects on the incoming solar radiation and outgoing infrared radiation. At present the impact of cirrus clouds on climate is well recognized but remains to be asserted more precisely, for their optical and radiative properties are not very well known. In order to understand the effects of cirrus clouds on climate, their optical and radiative characteristics of these clouds need to be determined accurately at different scales in different locations i.e. latitude. Lidars are well suited to observe cirrus clouds, they can detect very thin and semi-transparent layers, and retrieve the clouds geometrical properties i.e. altitude and multilayers, as well as radiative properties i.e. optical depth, backscattering phase functions of ice crystals. Moreover the linear depolarization ratio can give information on the ice crystal shape. In addition, the data collected with an airborne version of POLDER (POLarization and Directionality of Earth Reflectances) instrument have shown that bidirectional polarized measurements can provide information on cirrus cloud microphysical properties (crystal shapes, preferred orientation in space). The spaceborne version of POLDER-1 has been flown on ADEOS-1 platform during 8 months (October 96 - June 97), and the next POLDER-2 instrument will be launched in 2000 on ADEOS-2. The POLDER-1 cloud inversion algorithms are currently under validation. For cirrus clouds, a validation based on comparisons between cloud properties retrieved from POLDER-1 data and cloud properties inferred from a ground-based lidar network is currently under consideration. We present the first results of the validation.

Description: During the interval of 1944-1997, 120 intense hurricanes (category 3, 4, or 5 on the Saffir-Simpson hurricane scale) have been observed in the Atlantic basin. These intense hurricanes have had an observed annual frequency of 0-7 events per year (having a mean, mode, and median equal to about 2 events per year), being preferentially lower during El Ninio years and higher during non-El Ninio years. Also, it has recently been established that a long-term downward trend in the annual frequency of intense hurricanes, spanning about five decades, has taken place, although this trend can, alternatively, be explained as a shift from a more active state prior to the mid 1960's to a less active state thereafter (rather than as a simple linear decline). In this paper, on the basis of 10-yr moving averages, the long4erm trend of the frequency of intense hurricanes is compared against one for the annual mean temperature at Armagh Observatory, Northern Ireland (which serves as a proxy for climatic change). Interestingly, the two sets of 10-yr moving averages correlate extremely well, especially, when incorporating a slight 6-yr lag between them (with temperature leading; r = 0.90). This suggests that the current leading trend of temperature, which had been downward, but now is upward, may portend a return to the more active state for intense hurricanes. Thus, the 1998 season (presuming the abatement of El Ninio prior to the start of the hurricane season), and for several years thereafter (at least, into the early years of the next millennium), may have an annual frequency of intense hurricanes that is commensurate with the previously observed active state that was seen prior to the mid 1960's. If true, then, the shift to the more active state, probably, occurred in the mid-to-late 1980's, apparently, having gone undetected because of the masking, or modulating, effect of El Ninio, which has been rampant since the mid-to-late 1980's.

Description: Goals of the 3-month experiment GREENHOUSE using the equipment of greenhouse SVET (ECO-PSY-95) were to feature growth and development of wheat through the entire cycle of ontogeny under the maximally mimicked MIR environment, and to try out the procedures and timeline of space experiment GREENHOUSE-2 as a part of the fundamental biology investigations within the MIR/NASA space science program. Irradiation intensity (PAR) was 65 W/m2 and 38 W/m2 in the experiment and laboratory control, respectively. Values of other environmental parameters were MIR average (18-25 degrees C, relative air humidity in the interval between 40% and 75%, total gas pressure of about 660 to 860 mm Hg, partial oxygen pressure within the range from 140 to 200 mm Hg, partial carbon dioxide pressure up to 7 mm Hg). Experimental results showed that wheat cultivation in inhabited chamber under a modified lighting unit providing greater irradiation of the crop area produced more plant mass although seed production dropped. Low grain content in ears could be the aftermath of the gaseous trace contaminants in the chamber atmosphere.

Description: In 1996-1997 an experiment with super dwarf wheat (Greenhouse-2) was made aboard the orbital complex MIR as a part of the MIR-NASA space science program. The article deals with the main production and morphometric characteristics of plants that completed their vegetation cycle in the space flight. Lengths of the whole cycle of vegetation and its individual stages were essentially same as in ground control experiments. Dry mass of one plants equal, the number of headed shoots was in 2.7 times less in the flight harvest as compared with the control. The height of shoots was reduced by one half. No seeds were found in the heads formed in space. The architecture of heads was substantially different from what had been observed in the preceeding ground control experiments: mass of the heads was halved and lengths of inflorescence and palea awn shortened. The number of spikelets in a head reduced up to 8-10 vs. 13-14 in the controls, whereas the number of florets per a spikelet averaged 5 vs. 3 in the controls. The experiments showed that mainly the most profound changes in the productive and morphometric parameters of the super dwarf wheat plants were largely caused by the phytotoxic effects of ethylene rather than spaceflight specific factors as its concentrations in the MIR air amount to 0.3-1.8 mg/m3.

Description: The Galileo Probe entered the atmosphere of Jupiter on December 7, 1995. Measurements of the chemical and isotopic composition of the Jovian atmosphere were obtained by the mass spectrometer during the descent over the 0.5 to 21 bar pressure region over a time period of approximately 1 hour. The sampling was either of atmospheric gases directly introduced into the ion source of the mass spectrometer through capillary leaks or of gas, which had been chemically processed to enhance the sensitivity of the measurement to trace species or noble gases. The analysis of this data set continues to be refined based on supporting laboratory studies on an engineering unit. The mixing ratios of the major constituents of the atmosphere hydrogen and helium have been determined as well as mixing ratios or upper limits for several less abundant species including: methane, water, ammonia, ethane, ethylene, propane, hydrogen sulfide, neon, argon, krypton, and xenon. Analysis also suggests the presence of trace levels of other 3 and 4 carbon hydrocarbons, or carbon and nitrogen containing species, phosphine, hydrogen chloride, and of benzene. The data set also allows upper limits to be set for many species of interest which were not detected. Isotope ratios were measured for 3He/4He, D/H, 13C/12C, 20Ne/22Ne, 38Ar/36Ar and for isotopes of both Kr and Xe.

Description: We have performed high-resolution spectral observations at mid-infrared wavelengths of CH4 (8.14 micrometers), C2H6 (12.16 micrometers), and C2H2 (13.45 micrometers) on Jupiter. These emission features probe the stratosphere of the planet and provide information on the carbon-based photochemical processes taking place in that region of the atmosphere. The observations were performed using our cryogenic echelle spectrometer CELESTE, in conjunction with the McMath-Pierce 1.5-m solar telescope between November 1994 and February 1995. We used the methane observations to derive the temperature profile of the jovian atmosphere in the 1-10 mbar region of the stratosphere. This profile was then used in conjunction with height-dependent mixing ratios of each hydrocarbon to determine global abundances for ethane and acetylene. The resulting mixing ratios are 3.9(+1.9)(-1.3) x 10(-6) for C2H6 (5 mbar pressure level), and 2.3 +/- 0.5 x 10(-8) for C2H2 (8 mbar pressure level), where the quoted uncertainties are derived from model variations in the temperature profile which match the methane observation uncertainties. c1998 Academic Press.

Description: Water vapor plays an important role in many atmospheric processes related to radiation, climate change, atmospheric dynamics, meteorology, the global hydrologic cycle, and atmospheric chemistry, and yet our knowledge of the global distribution of water vapor is very limited. The differential absorption lidar (DIAL) technique has the potential of providing needed high resolution water vapor measurements from aircraft and from space, and the Lidar Atmospheric Sensing Experiment (LASE) is a key step in the development of this capability. The LASE instrument is the first fully engineered, autonomous DIAL system, and it is designed to operate from a high-altitude aircraft (ER-2) and to make water vapor and aerosol profile measurements across the troposphere. The LASE system was flown from the NASA Wallops Flight Facility in a series of engineering flights during September 1994. This paper discusses the characteristics of the LASE system and presents the first LASE measurements of water vapor and aerosol profiles.

Description: Low temperature Zn diffusion in GaSb, where the minimum temperature was 450 C, was studied. The pseudo-closed box (PCB) method was used for Zn diffusion into GaAs, AlGaAs, InP, InGaAs and InGaAsP. The PCB method avoids the inconvenience of sealed ampoules and proved to be simple and reproducible. The special design of the boat for Zn diffusion ensured the uniformality of Zn vapor pressure across the wafer surface, and thus the uniformity of the p-GaSb layer depth. The p-GaSb layers were studied using Raman scattering spectroscopy and the x-ray rocking curve method. As for the postdiffusion processing, an anodic oxidation was used for a precise thinning of the diffused GaSb layers. The results show the applicability of the PCB method for the large-scale production of the GaSb structures for solar cells.

Description: Indium phosphide is considered to be a strong contender for many photovoltaic space applications because of its radiation resistance and its potential for high efficiency. An overview of recent progress is presented, and possible future research directions for indium phosphide space solar cells are discussed. The topics considered include radiation damage studies and space flight experiments.

Description: Understanding solar cell response to pulsed laser outputs is important for the evaluation of power beaming applications. The time response of high efficiency GaAs and silicon solar cells to a 25 nS monochromatic pulse input is described. The PC-1D computer code is used to analyze the cell current during and after the pulse for various conditions.

Description: Four-terminal mechanically stacked solar cells were developed for advanced space arrays with line-focus reflective concentrators. The top cells are based on AlGaAs/GaAs multilayer heterostructures prepared by low temperature liquid phase epitaxy. The bottom cells are based on heteroepitaxial InP/InGaAs liquid phase epitaxy or on homo-junction GaSb, Zn-diffused structures. The sum of the highest reached efficiencies of the top and bottom cells is 29.4 percent. The best four-terminal tandems have an efficiency of 27 to 28 percent. Solar cells were irradiated with 1 MeV electrons and their performances were determined as a function of fluence up to 10(exp 16) cm(exp-2). It was shown that the radiation resistance of developed tandem cells is similar to the most radiative stable AlGaAs/GaAs cells with a thin p-GaAs photoactive layer.

Description: An investigation of solar cells based on AlGaAs/GaAs heterostructures with an internal Bragg reflector as the back-surface reflector is presented. The Bragg reflector is grown by low pressure metalorganic chemical vapor deposition on n-GaAs substrates in a horizontal resistively heated reactor. The Bragg reflector with its maximum reflectance centered at a wavelength of 860 nm consists of 12 pairs of AlAs/GaAs layers. The resulting Bragg reflector has a thickness of 0.072 micrometers for AlAs and 0.059 micrometers for GaAs. The multi-layered quasi-dielectric stack selectively reflects weakly absorbed photons with energies near to the GaAs band gap for a second pass through the photoactive region, thus increasing the photocurrent. The use of the Bragg reflector allows the external quantum efficiency to be increased in the long wavelength of the spectrum. The use of the Bragg reflector and an antireflective coating and prismatic cover allowed an efficiency of 23.4 percent to be obtained.

Description: The high efficiency solar cells based on multilayer AlGaAs/GaAs heterostructures, prepared by low temperature liquid phase epitaxy (LPE), were developed and tested. An investigation of the low temperature LPE process for the crystallization of AlGaAs heterostructures of as high as 24.0 to 24.7 percent under AMO conditions at concentration ratios of 20 to 100x, were reached. Developed solar cells show substantial radiation resistance to the damage induced by 3.75 MeV electrons.

Description: A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.

Description: Radiation resulting from interaction between the effluent cloud of a space shuttle thruster and the ambient atmosphere was observed with a spectograph aboard the shutttle. The spectral measurements were made between 400 and 800 nm with a resolutoion of 3 nm. The primary emissions are identified as NO2, HNO, O(1)D, and O(1)S. These are the first observations od O(1)S emission in the shuttle plume. These data are compared with the previous measurements, and possible excitation mechanisms are discussed. The results are also compared with a Monte Carlo simulation of thruster plume-atmosphere interaction radiation.

Description: We report multi-instrument observations during an isolated substorm on 17 October 1989. The European Incoherent Scatter (EISCAT) radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71 deg Lambda - 78 deg Lambda. Sub-Auroral Magnetometer Network (SAMNET) and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. Interplanetary Monitoring Platform-8 (IMP-8) magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux in the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71 deg Lambda by the time of the expansion phase onset. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the distant neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase.

Description: The impact of the extraterrestrial object that formed the Chicxulub crater in the northwestern Yucatan peninsula of Mexico is the leading suspect for the extinction of the dinosaurs. This article reports on a Planetary Society expedition to Albion Island in the Rio Hondo region of Belize to investigate evidence supporting the impact theory.

Description: Direct current (DC) electric field and ion density measurements near density depletion regions (that is, equatorial plasma bubbles) are used to estimate the vertical neutral wind speed. The measured zonal electric field in a series of density depletions crossed by the San Marco D satellite at 01.47-01.52 Universal Time (UT) on 25 October 1988, can be explained if a downward neutral wind of 15-30 m/s exists. Simultaneously, the F-region plasma was moving downward at a speed of 30-50 m/s. These events appear in the local time sector of 23.00-23.15 in which strong downward neutral winds may occur. Indeed, airglow measurements suggest that downward neutral velocities of 25-50 m/s are possible at time near midnight in the equatorial F-region.

Description: An overview of the observations of backstreaming electrons in the foreshock and the mechanisms that have been proposed to explain their properties will be presented. A primary characteristic of observed foreshock electrons is that their velocity distributions are spatially structured in a systematic way depending on distance from the magnetic field line which is tangent to the shock. There are two interrelated aspects to explaining the structure of velocity distributions in the foreshock, one involving the acceleration mechanism and the other, propagation from the source to the observing point. First, the source distribution of electrons energized by the shock must be determined along the shock surface. Proposed acceleration mechanisms include magnetic mirroring of incoming solar wind particles and mechanisms involving transmission of particles through the shock. Secondly, the kinematics of observable electrons streaming away from a curved shock with an initial parallel velocity and a downstream perpendicular velocity component due to the motional electric field must be determined. This is the context in which the observations and their explanations will be reviewed.

Description: This paper focuses on the comparison of cloud amounts derived from an atmospheric general circulation model (AGCM), Satellite-observed clouds, and Ground-based cloud observations. Unlike Earth Radiation Budget Experiment (ERBE)-type comparisons it does not mix potential errors in the cloud amount with those in the radiation code embedded in the model. Long term cloud climatologies were used to compare global cloud amounts and regional seasonal cycles. The AGCM successfully reproduced the signatures of the warm pool and North Pacific seasonal cycle cloudiness but failed in the low stratus region off the coast of South America, a known problem for AGCMs. The data sets also reproduced the anomaly signature associated with El Nino in the warm pool region, but the model amounts were lower. Global results had a similar success rate, with the model generally producing lower total cloud compared to the satellite and in situ measurements. To compare cloud vertical distributions the cloud height may need to be validated using the corresponding radiation fields. Unfortunately there were also some large discrepancies between the two observed cloud data sets. While tremendously improved over the last decade the character of the observed cloud data sets, must be substantially enhanced before they will be useful in validating AGCMs by any but the crudest levels of comparison.

Description: Data from both 27 sites in the Atlanta mesonet surface meteorological network and eight National Weather Service sites were analyzed for the period from 26 July to 3 August 1996. Analysis of the six precipitation events over the city during the period (each on a different day) showed that its urban heat island (UHI) induced a convergence zone that initiated three of the storms at different times of the day, i.e., 0630,0845, and 1445 EDT. Previous analysis has shown that New York City (NYC) effects summer daytime thunderstorm formation and/or movement. That study found that during nearly calm regional flow conditions the NYC UHI initiates convective activity. Moving thunderstorms, however, tended to bifurcate and to move around the city, due to its building barrier effect. The current Atlanta results thus agree with the NYC results with respect to thunderstorm initiation.