ALBERT

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: 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: How much of the interannual variation in the satellite derived radiation balance can be purely attributed to changes taking place at the land surface, was examined. The role of surface latent heating was examined in relation to its control of the precipitation pattern from one year to the next.

Description: As a check on structure safety aspects, two approaches in seismic analysis for the large 70-m antennas are presented. The first approach, commonly used by civil engineers, utilizes known recommended design response spectra. The second approach, which is the full transient analysis, is versatile and applicable not only to earthquake loading but also to other dynamic forcing functions. The results obtained at the fundamental structural frequency show that the two approaches are in good agreement with each other and both approaches show a safe design. The results also confirm past 64-m antenna seismic studies done by the Caltech Seismology Staff.

Description: The major criterion for the Atmospheric General Circulation Experiment (AGCE) design is that it be possible to realize strong baroclinic instability in the spherical configuration chosen. A configuration was selected in which a hemispherical shell of fluid is subjected to latitudinal temperature gradients on its spherical boundaries and the latitudinal boundaries are insulators. Work in the laboratory with a cylindrical version of this configuration revealed more instabilities than baroclinic instability. Since researchers fully expect these additional instabilities to appear in the spherical configuration also, they decided to continue the laboratory cylindrical annulus studies. Four flow regimes were identified: an axisymmetric Hadley circulation, boundary layer convection, baroclinic waves and deep thermal convection. Regime diagrams were prepared.

Description: A two-layer truncated baroclinic spectral model was developed to study the long-term evolution of disturbances to a baroclinically unstable mean flow. Topography and crudely-parameterized radiative processes were accounted for. As a result of Robert Schlaak's discovery of the underlying barotropic nature of the index oscillation as well as reviewers suggestions about the original manuscript, the model has been revised to allow for barotropic as well as baroclinic wave-mean flow interactions. The form-drag exerted by the topography on the barotropic part of the mean flow is larger than on the baroclinic part and thus researchers anticipate significant changes from the original calculations on the index oscillation when it is strongly modulated by topography. Researchers believe that since the index oscillation accounts for a significant portion of atmospheric temporal variance, the long term predictability could be improved if reliable forecasts of the index oscillation were available. Two spectral models of the index oscillation, one barotropic and the other baroclinic, have been developed. The latter allows for moisture, radiation, land-sea temperature countrasts, and energy exchanges with the underlying surface.

Description: Recent atmospheric calculation suggest that the prebiological atmosphere was most probably composed of nitrogen, carbon dioxide, and water vapor, resulting from volatile outgassing, as opposed to the older view of a strongly reducing early atmosphere composed of methane, ammonia, and hydrogen. Photochemical calculations indicate that methane would have been readily destroyed via reaction with the hydroxyl radical produced from water vapor and that ammonia would have been readily lost via photolysis and rainout. The rapid loss of methane and ammonia, coupled with the absence of a significant source of these gases, suggest that atmospheric methane and ammonia were very short lived, if they were present at all. An early atmosphere of N2, CO2, and H2O is stable and leads to the chemical production of a number of atmospheric species of biological significance, including oxygen, ozone, carbon monoxide, formaldehyde, and hydrogen cyanide. Using a photochemical model of the early atmosphere, the chemical productionof these species over a wide range of atmospheric parameters were investigated. These calculations indicate that early atmospheric levels of O3 were significantly below the levels needed to provide UV shielding. The fate of volcanically emitted sulfur species, e.g., sulfur dioxide and hydrogen sulfide, was investigated in the early atmosphere to assess their UV shielding properties. The photochemical calculations show that these species were of insufficient levels, due in part to their short photochemical lifetimes, to provide UV shielding.

Description: A false-color multipolarization version of one of the images of Owens Valley area acquired by the JPL Synthetic Aperture Radar (SAR) is given. A geologic map of the alluvial fans there (Gillespie, 1982) is also given for comparison. In general, brightness in the multipolarization images can be seen to be inversely proportional to the age of the surfaces. A more detailed investigation of the relationship between backscatter and age of the surfaces was undertaken with calibrated aircraft SAR data. The quantitative relationship between backscatter coefficient and age for the three polarizations is shown. The straight lines connecting the measured data points imply a steady-state process, although the process or processes leading to this relationship may have operated at rates that varied with climate fluctuations, such as the glacial ages. It is expected that the relationship between radar brightness and age is a consistent one, and that with the wider availability of calibrated radar backscatter data, these relationships can be less well-known areas. The effect of variable such as past climate fluctuations, tectonic disturbance, and rock type must be understood before extension beyond the Mojave Desert region can be attempted.

Description: Mobile Very Long Base Interferometry (VLBI) and Global Positioning System (GPS) geodetic measurements have many error sources in common. Calibration of the effects of water vapor on signal transmission through the atmosphere, however, remains the primary limitation to the accuracy of vertical crustal motion measurements made by either technique. The two primary methods of water vapor calibration currently in use for mobile VLBI baseline measurements were evaluated: radiometric measurements of the sky brightness near the 22 GHz emission line of free water molecules and surface meteorological measurements used as input to an atmospheric model. Based upon a limited set of 9 baselines, it is shown that calibrating VLBI data with water vapor radiometer measurements provides a significantly better fit to the theoretical decay model than calibrating the same data with surface meteorological measurements. The effect of estimating a systematic error in the surface meteorological calibration is shown to improve the consistency of the vertical baseline components obtained by the two calibration methods. A detailed error model for the vertical baseline components obtained indicates current mobile VLBI technology should allow accuracies of order 3 cm with WVR calibration and 10 cm when surface meteorological calibration is used.

Description: Meteorologists and astrophysicists interested in large scale planetary and solar circulations have come to recognize the importance of rotation and stratification in determining the character of these flows. In particular, the effect of latitude-dependent Coriolis force on nonlinear convection is thought to play a crucial role in such phenomena as differential rotation on the Sun, cloud band orientation on Jupiter, and the generation of magnetic fields in thermally driven dynamos. The continuous low-gravity environment of the orbiting space shuttle offers a unique opportunity to make laboratory studies of such large-scale thermally driven flows under the constraint imposed by rotation and sphericity. This is possible because polarization forces in a dielectric liquid, which are linearly dependent on fluid temperature, give rise to an effectively radial buoyancy force when a radial electrostatic field is imposed. The Geophysical Fluid Flow Cell (GFFC) is an implementation of this ideal in which fluid is contained between two rotating hemispheres that are differentially heated and stressed with a large a-c voltage. The experiment, to be flown on Spacelab III (currently set for launch April 29, 1985), will explore non-linear mode selection and high Rayleigh number turbulence in a rotating convecting spherical shell of liquid. Experiments will be carried out in a low driving parameter range where some limited numerical experimentation is currently feasible, as well as in a parameter range significantly beyond numerical computation for many years.

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: Comparisons are presented of satellite, rocket, and balloon ozone profiles near Natal, Brazil (5.9 deg S, 35.2 deg W). The low variability of stratospheric ozone at Natal during March and April of 1985 has allowed intercomparisons of reasonably large data sets, rather than a small number of paired satellite/in situ comparisons. There are sharp differences between the profile from the SBUV instrument on Nimbus 7 and the in situ measurements. These results support the conclusions of the NASA Ozone Trends Panel that there is an instrumental cause for the very large changes in upper stratospheric ozone seen by SBUV. Along with other comparisons, these results are being used in a reassessment of the SBUV instrument and its data reduction procedures. The agreement between the ozone profiles from the SAGE II instrument on the ERBS satellite and the rocket values is excellent over the full range of comparisons. Both SAGE II and ROCOZ-A must convert from altitude to pressure for intercomparisons with SME and with SBUV-type instruments. The conversion between pressure and altitude is as important as the ozone measurements, especially in the upper stratosphere where the scale height for ozone is approximately half that for pressure.

Description: A recent reanalysis of the International Latitude Serivce (ILS) polar motion data-day has been processed using Kalman filtering techniques to generate the polar motion excitation function over the time-span from 1960 to 1965. The resulting excitation function has been examined for the effects of 1960 Chile in an attempt to determine experimentally how large earthquake affect polar motion. The resulting upper bound of about 75 x 10 to the 22nd N-m for a 10-deg dip (about 36 x 10 to the 22nd N-m for 20-deg dip) is consistent with results obtained from previous seismic studies, including a recent normal mode excitation result. Following future great earthquakes, monitoring of polar motion by space-based techniques such as VLBI should continue at high temporal resolution for several weeks in order to directly measure the rheological parameters of the upper mantle.

Description: The Clark Mountains in eastern California form a rugged, highly dissected area nearly 5000 ft above sea level, with Clark Mountain rising to 8000 ft. The rocks of the Clark Mountains and the Mescal Range just to the south are Paleozoic carbonate and clastic rocks, and Mesozoic clastic and volcanic rocks standing in pronounced relief above the fractured Precambrian gneisses to the east. The Permian Kaibab Limestone and the Triassic Moenkopi and Chinle Formations are exposed in the Mescal Range, which is the only place in California where these rocks, which are typical of the Colorado Plateau, are found. To the west, the mountains are bordered by the broad alluvial plains of Shadow Valley. Cima Dome, which is an erosional remnant carved on a batholithic intrusion of quartz monzonite, is found at the south end of the valley. To the east of the Clark and Mescal Mountains is found the Ivanpah Valley, in the center of which is located the Ivanpah Play. Studies of the Clark Mountains with the airborne visible/infrared imaging spectrometer are briefly described.

Description: Microwave limb-sounding can be used to improve understanding of earth's upper atmosphere. This paper summarizes general features of the technique. An experiment being developed for the Upper Atmosphere Research Satellite is described. Plans for the future Earth Observing System are discussed.

Description: Seasonal and diurnal emissions of NO and N2O from agricultural sites in Jamestown, Virginia and Boulder, Colorado are estimated in terms of soil temperature; percent moisture; and exchangeable nitrate, nitrite, and ammonium concentrations. The techniques and procedures used to analyze the soil parameters are described. The spatial and temporal variability of the NO and N2O emissions is studied. A correlation between NO fluxes in the Virginia sample and nitrate concentration, temperature, and percent moisture is detected, and NO fluxes for the Colorado site correspond with temperature and moisture. It is observed that the N2O emissions are only present when percent moisture approaches or exceeds the field capacity of the soil. The data suggest that NO is produced primarily by nitrification in aerobic soils, and N2O is formed by denitrification in anaerobic soils.

Description: Mass spectrometer measurements of ozone made during two balloon flights included its heavy isotopes at mass 49 and 50. Both flights were flown during the day and during summer from Palestine, TX. At float altitudes above 42 km the enrichments in heavy ozone were 41 percent and 23 percent, respectively. The enrichment appears to be mass independent since, at high altitudes, both 49 and 50 show the same enhancement. During the descent the enrichment in heavy ozone decreased, faster during the first flight than during the second, reaching values between 15 and 20 percent above 30 km. Near and below this altitude another increase is observed. During a night flight, previously reported, an enhancement in heavy ozone of over 40 percent at 32 km was found, decreasing both toward higher and lower altitudes.

Description: The observation of a geomagnetic signature in the zonal eastward plasma flow, which is a striking feature of the equatorial ionosphere in the evening quadrant is reported. These observations were derived fronm (E x B)/B-squared measurements made with the cylindrical double-floating-probe experiment carried on the Dynamics Explorer 2 satellite. The signature consists of a crest-trough-crest effect in the latitude dependence of the eastward plasma flow with the crests at + or - 8 dip latitude and the trough nearly centered at the dip equator at all geographic longitudes. This phenomenon can be readily interpreted in terms of the altitude dependence of the F region dynamo electric field, and it is related to dip equator signatures in the plasma density and the magnetic declination which have been reported earlier.

Description: The ROCOZ-A radiometer measures ozone by long pathlength photometry in the stratosphere and lower mesosphere. After a rocket launch to an apogee of 70 km, the instrument measures the solar ultraviolet irradiance over its four filter wavelengths as it descends on a parachute. The fundamental values from ROCOZ-A are ozone overburdens versus radar altitude from 53 to 20 km. The slope of these values gives ozone number density. At one standard deviation the repeatability of the ozone overburden measurements averages 2.4 percent. For ozone number density the repeatability averages 3.2 percent with a significant increase at altitudes below the ozone number density maximum. The accuracy limits for overburden and number density are estimated at 5-7 percent. With auxiliary measurements of pressure and temperature, ozone results are also produced in terms of ozone mixing ratio, albeit with a slight broadening of the estimated accuracy limits. The vertical response of ROCOZ-A ozone measurements (full width at half maximum) is 4 km. The assembly of ROCOZ-A profiles can be used to compare with measurements from each of the current NASA and NOAA satellite ozone instruments. In addition, the repeatability of ROCOZ-A allows the use of this instrument as a transfer standard between satellite instruments with different fundamental ozone measurements.

Description: A possible cause of the large plasma flow velocities parallel to the magnetic field (which were observed in spacecraft experiments) near the boundary of the plasma sheet in the earth's magnetotail is considered in the framework of a magnetohydrodynamic model. It is shown for steady-state configurations that high parallel flow velocities can be expected to exist on field lines connecting to a region of weak magnetic field. The physical mechanism causing large values of the parallel velocity component can be visualized as a strong imbalance of perpendicular mass flux into and out of magnetic flux tubes passing through regions where the magnetic field is weak and inhomogeneous. The value of the parallel velocity component is evaluated, and it is found that it can substantially exceed the perpendicular velocity (by as much as a factor of 40). The results are applied to the earth's magnetotail; it is found that this mechanism is able to explain the parallel flow velocities near the boundary of the plasma sheet in the range of several hundreds of km/s.

Description: This paper presents a representative example of an enhancement in energetic ion flux associated with the International Sun-Earth Explorer 3 (ISEE 3) spacecraft's encounter with a traveling compression region (TCR). Data from the energetic particle anisotropy spectrometer (EPAS) instrument on ISEE 3 are studied, along with magnetic field data from the vector helium magnetometer. It is concluded that the ion enhancements seen are spatial in nature, thus supporting the idea that TCRs are the lobe signatures of plasmoids moving along the magnetotail, away from earth.

Description: The spatial and temporal characteristics of ozone density measured from the SBUV (Solar Backscatter Ultraviolet) spectrometer on Nimbus-7 and the UV and the UV and the IR spectrometers on SME (Solar Mesosphere Explorer) are compared in the altitude region near 50 km where the three data sets overlap. Their temporal characteristics, when averaged over the same longitude range, are remarkably similar with respect to seasonal variations and short term fluctuations induced by transient planetary waves. The long term trends in the three data sets, however, differ significantly with each other. Over the three year period after 1982 ozone mixing ratio at 1 mb decreased by about 10 percent based on SEUV measurements but increased by 12 and 30 percent respectively based on SME-IR and SME-UV measurements. None of these estimates are consistent with the predicted decrease of about 2 percent based on solar UV flux and temperature changes during this period.

Description: The phase diagram of FeO has been experimentally determined to pressures of 155 GPa and temperatures of 4000 K using shock-wave and diamond-cell techniques. A metallic phase of FeO is observed at pressures greater than 70 GPa and temperatures exceeding 1000 K. The metallization of FeO at high pressures implies that oxygen can be present as the light alloying element of the earth's outer core, in accord with the geochemical predictions of Ringwood (1977 and 1979). The high pressures necessary for this metallization suggest that the core has acquired its composition well after the initial stages of the earth's accretion. Direct experimental observations at elevated pressures and temperatures indicate that core-forming alloy can react chemically with oxides such as those forming the mantle. The core and mantle may never have reached complete chemical equilibrium, however. If this is the case, the core-mantle boundary is likely to be a zone of active chemical reactions.

Description: A method for calculating the photodissociation rates needed for photochemical modeling of the stratosphere, which includes the effects of molecular scattering, is described. The procedure is based on Sokolov's method of averaging functional correction. The radiation model and approximations used to calculate the radiation field are examined. The approximated diffuse fields and photolysis rates are compared with exact data. It is observed that the approximate solutions differ from the exact result by 10 percent or less at altitudes above 15 km; the photolysis rates differ from the exact rates by less than 5 percent for altitudes above 10 km and all zenith angles, and by less than 1 percent for altitudes above 15 km.

Description: Data from the BOIC which consisted of three balloon missions conducted in Palestine, Texas from June 1983 to March 1984 are presented. The BOIC was to assess the ability to perform ozone measurements from balloon platforms. The accuracy and precision of the various ozone measurement systems, which were composed of a photometer, a mass spectrometer, and solar UV absorption sensors, are evaluated. The ozone observations obtained with the instruments on the three flight missions are analyzed and intercompared. The flight in situ data are also compared to the National Bureau of Standards reference photometer, satellite measurements, and under simulated stratospheric pressure and ozone concentrations.

Description: An analysis procedure has been developed for derivation of information about the photochemical behavior of ozone near 1 mbar by explicitly accounting for the dynamical transport terms in the continuity equation for perturbations from the zonal mean. The procedure is shown to be valid using data from a numerical transport model and is then applied to LIMS ozone and temperature data, using geostrophic winds to estimate the transport terms. The data study is restricted to March at 2, 1, and 0.7 mbar. Because the temperature deviations are dynamically produced, large temperature deviations are associated with significant ozone transport terms. The anticorrelation between the deviations of ozone and temperature disappears when the transport terms are small. The derived photochemical information is compared to photochemical theory. Although there is overall agreement in the magnitude as well as the latitude, altitude, and time dependencies, discrepancies are suggested which may be related to the long-standing failure of photochemical models to calculate ozone accurately near 1 mbar. The theory suggests that the addition of Cl(x) to the stratosphere will affect the relationship of ozone and temperature at 2 and 1 mbar. Comparison of the photochemical information derived from future ozone and temperature measurements with the results of the present analysis should provide a critical test of the photochemical scheme thought to describe the behavior of odd chlorine in the stratosphere.

Description: A detailed statistical analysis of monthly averages of ozonesond readings is performed to assess trends in ozone in the troposphere and the lower to midstratosphere. Regression time series models, which include seasonal and trend factors, are estimated for 13 stations located mainly in the midlatitudes of the Northern Hemisphere. At each station, trend estimates are calculated for 14 'fractional' Umkehr layers covering the altitude range from 0 to 33 km. For the 1970-1982 period, the main findings indicate an overall negative trend in ozonesonde data in the lower stratosphere (15-21 km) of about -0.5 percent per year, and some evidence of a positive trend in the troposphere (0-5 km) of about 0.8 percent per year. An in-depth sensitivity study of the trend estimates is performed with respect to various correction procedures used to normalize ozonesonde readings to Dobson total ozone measurements. The main results indicate that the negative trend findings in the 15- to 21-km altitude region are robust to the normalization procedures considered.

Description: ISEE particle and wave data are noted to furnish substantial support for the basic features of the velocity dispersed model at the foreshock boundary that was proposed by Filbert and Kellogg (1979). Among many remaining discrepancies between this model and observation, it is noted that unstable reduced velocity distributions have been discovered behind the thin boundary proposed by the model, and that these are at suprathermal energies lying far below those explainable in terms of an oscillating, two-stream instability. Although the long-theorized unstable beam of electrons has been found in the foreshock, there is still no ready explanation of the means by which it could have gotten there.

Description: Theoretical models of planetary-atmosphere tidal fields are examined analytically, comparing models based on geometric-height coordinates with those employing log-pressure coordinates. The relationship between the linearized meteorological variables in the two systems is explored for classical tidal theory and its reduction to the Laplace tidal equation, and it is shown that identical horizontal and vertical equations are obtained. Also considered are the tidal zonal-mean bilinear flux convergences and their Eliassen-Palm formulations. Numerical results for problems involving the earth and Mars atmospheres are presented in graphs and discussed.

Description: The analytical results of Zurek (1986) are applied to compare the vertical-structure equations obtained by Dickinson and Geller (1968) and Lindzen and McKenzie (1967) in solving the classical tidal equations for planetary atmospheres when Newtonian cooling (NC) is included as part of diabatic tidal forcing. The two diabatic forcings are found to be different when the basic state temperature varies with height, and it is recommended that a log-pressure NC formulation be used whenever temperature is made a function of pressure in the radiative-damping calculation. Scaling arguments are presented to show that the choice of NC formulation has significant effects only in cases where the radiative time constants are short (such as the thin CO2 atmosphere of Mars). Numerical computations for that case indicate a difference of 20 percent in amplitude and 15 deg in phase for the (probably meteorologically significant) 4.1-d wavenumber-one Rossby tidal mode.

Description: Recent ground-based infrared solar spectra at 0.02 per cm resolution in the 3000 per cm region have been analyzed for the atmospheric content of HCl. Nonlinear spectral least-squares fitting applied to spectra obtained at several zenith angles shows little sensitivity of the results to tropospheric HCl but provides an accurate measurement of the total column amount.

Description: This review presents a summary of past work on the ISEE-3 distant tail magnetic field observations. An attempt has been made to bring the many results together as a coherent whole, in the hope that the reader can envision the direction of future research necessary to achieve an understanding of the dynamics of the magnetotail from 60 to 240 earth radii and perhaps beyond.

Description: Ozone measurements from 1970 to 1984 from the Nimbus 4 backscattered ultraviolet and the Nimbus 7 solar backscattered ultraviolet spectrometers show significant decrease in total ozone only after 1979. The downward trend is most apparent in October south of 70 deg S in the longitude zone 0 to 30 deg W where planetary wave activity is weak. Outside this longitude region, the trend in total ozone is much smaller due to strong interannual variability of wave activity. This paper gives a phenomenological description of ozone depletion in the Antarctic region based on vertical advection and transient planetary waves.

Description: A more detailed examination of the TOMS observations reveals a number of important aspects of the Antarctic ozone depletion. First, it is noted that the presence of large scale disturbances around the edge of the ozone hole can lead to highly variable station observations. Second, an examination of the zonal mean total ozone for the seven Octobers shows that the large systematic decline is not simply confined to the polar region but exists at midlatitudes as well. Integrations of the total ozone from the South Pole northwards show that a portion of the systematic trend of decreasing Antarctic total ozone (prior to 1985) seems to be due to a redistribution of total ozone to subpolar and midlatitude regions. That is, decreases at high latitudes are compensated by increases at lower latitudes. The correlation between zonal mean total ozone and 70 mb zonal mean temperatures from polar to midlatitudes shows that the systematic decreases in total ozone is well correlated with a systematic decrease in stratospheric temperatures.

Description: The ozone profiles measured by SBUV in the region of the Antarctic ozone minimum are in error as archived for altitudes below the 7 mb level in the atmosphere because of the lack of a reasonable climatology to use for the initial guess profile. The ozone profile error in this region is examined, and it is shown that use of a reasonable climatology in this unusual region results in ozone profiles that agree substantially better with the balloon observations available from the Syowa station. The corrected profiles show that by 1984 there was 35 percent less ozone in the 15 to 30 km region than in 1979, and 14 percent less even in the 40 to 50 km region.

Description: Daytime and nighttime vertical profiles of the tropospheric trace gas N2O were determined from grab sample collections off the Atlantic and Gulf coasts of Florida. The grab samples were collected during the week of October 7-13, 1984, from a Lear jet during descent spirals over an altitude range of 12.5-0.3 km in approximately 1.2-km intervals. During this period there were two distinct airflow regimes sampled: (1) the surface boundary layer (less than 2 km), in which the wind direction was typically easterly; and (2) the regime above the boundary layer, which was predominantly characterized by westerly flow. N2O mixing ratios, normalized to dry air, were determined from 148 daytime and nighttime samplings. N2O was found to be uniformly mixed at all altitudes at 301.9 + or - 2.4 parts per billion by volume.

Description: Mobile Very Long Base Interferometry (VLBI) and Global Positioning System (GPS) geodetic measurements have many error sources in common. Calibration of the effects of water vapor on signal transmission through the atmosphere, however, remains the primary limitation to the accuracy of vertical crustal motion measurements made by either technique. The two primary methods of water vapor calibration currently in use for mobile VLBI baseline measurements were evaluated: radiometric measurements of the sky brightness near the 22 GHz emission line of free water molecules and surface meteorological measurements used as input to an atmospheric model. Based upon a limited set of 9 baselines, it is shown that calibrating VLBI data with water vapor radiometer measurements provides a significantly better fit to the theoretical decay model than calibrating the same data with surface meteorological measurements. The effect of estimating a systematic error in the surface meteorological calibration is shown to improve the consistency of the vertical baseline components obtained by the two calibration methods. A detailed error model for the vertical baseline components obtained indicates current mobile VLBI technology should allow accuracies of order 3 cm with WVR calibration and 10 cm when surface meteorological calibration is used.

Description: Techniques employed to map the global distribution of stratospheric aerosols using SAGE II data at four wavelengths (385, 453, 525 and 1020 microns) are described. The methods were devised to integrate the 900 data points scanned by the instrument with account taken of the different times the vertical profiles were obtained. The data have been used to calculate the O3, H2O, NO2 and aerosol optical thicknesses. Sample details of data gathered during April 1985 are discussed, with emphasis on difficulties being encountered in analyzing the differences being observed in the dawn and dusk data.

Description: Data on the characteristics of the stratospheric aerosol as measured with sensors on the SAM II and SAGE I satellites and with ground-based and airborne lidar are discussed. Emphasis is placed on the impact of the El Chichon eruptions. The volcanic cloud was tracked to an altitude of 30 km, and was observed to travel around the earth in 3 weeks. The maximum stratospheric loading is estimated at 12 Mtons, which increased the stratospheric optical depth to 0.15-2.0 at the peak period. The particulate loading was predicted to lower the Northern Hemisphere average temperatures by 0.4-0.5 C in 1984-85.

Description: Alongtrack data collected with the Earth Radiation Budget Experiment (ERBE) scanner instrument were used to study an observed limb-darkening phenomenon. A numerical model developed for the longwave exitance as a function of the SZA has agreed well with ERBE data, indicating that the atmosphere is in adiabatic, rather than radiative, equilibrium. A corrected form of the model has been defined for SZA over 60 deg. The model was used to parameterize diurnal data for the Sahara-Saudi and Australian deserts.

Description: Design features and the performance envelope of the SAGE II stratospheric aerosol monitoring instrument on the Earth Radiation Budget Satellite are described. SAGE II was designed to obtain vertical profiles of stratospheric aerosols, monitor global seasonal changes in aerosols, provide data on stratospheric circulation and the behavior of transient events such as volcanic particulate injections, and to investigate atmospheric chemistry. The mmeasurements are centered on extinctions due to aerosols, NO2, O3 and water vapor.

Description: Instruments on board the (presently two) satellites of the Earth Radiation Budget Experiment (ERBE) are collecting data for determining monthly averaged radiation exitances at the top of the atmosphere (TOA). To achieve the accuracy desired of the mission, radiances at the satellite are first calculated, with allowance made of sensor optical properties and the directionality of the TOA radiation field. The subsatellite surface type is classified to adjust for albedo changes and correction values are added for the types of cloud cover detected.

Description: Early results are reported from measurements of the diurnal variability of total and clear-sky regional radiative parameters by the ERBE instruments on one dedicated satellite and the polar-orbiting NOAA-9 satellite. Attention is focused on November 1984, the first complete data set. The scene is identified in terms of longwave and shortwave radiances (daytime) or longwave radiation (night) and maximum likelihood estimates carried out with the addition of Earth Radiation Budget data from Nimbus-7. Analysis of the first data set revealed significant differences between total and clear-sky albedo. The clear-sky and LRE both reach maximum around noon and minimum values at midnight.

Description: The Earth Radiation Budget Experiment (ERBE) is to be a three-satellite complement, of which two satellites were in orbit as of October 1984. Each satellite carries scanners for earth radiance sensing and one nonscanner for solar constant sensing. Of the four channels for earth scanning, two are dedicated to limb-to-limb observations and two view a swath 1000 km wide. Comparisons are made among the channel readings, thus far producing agreement to within 1 percent. Inversion techniques which are being applied to data gathered by the sensors are checked against equivalent procedures and results with the Nimbus-7 and GOES satellites.

Description: The effects changes in stratospheric ozone concentrations have on photochemical models of the atmosphere are evaluated. The study was spurred by the appearance of excessive stratosphere heating above 40 km in radiative transfer computations. A 20 percent reduction in ozone in the 40-55 m altitude interval would offset the anomalous heating values. Comparisons were made between calculated heating rates and in-situ data that included ozone concentrations. Heating rates decreased 0.4 C/day when clouds were present at the 300 mb level, and increased by 0.3 C/day with clouds at the 80 mb level.

Description: The SAGE II limb-scanning radiometer carried on the Earth Radiation Budget Satellite functions at wavelengths of 0.385, 0.45, 0.525, and 1.02 microns to identify vertical profiles of aerosol density by atmospheric extinction measurements from cloud tops upward. The data are being validated by correlating the satellite data with data gathered with, e.g., lidar, sunphotometer, and dustsonde instruments. Work thus far has shown that the 1 micron measurements from the ground and satellite are highly correlated and are therefore accurate to within measurement uncertainty.

Description: Nimbus-7 wide-FOV irradiance data collected from 1981-1985 are used to evaluate the effects of the El Chichon eruptions of 1982 on the earth radiation budget. The north polar region displayed a maximum response of 20 percent in the winter of 1982-1983, with the variation being most apparent in the near-IR 2.8 micron and 0.2-3.8 microns shortwave bands. The data indicate that the particle size distribution was constant for a year after eruptions.

Description: Laboratory kinetic studies suggest that an appreciable fraction of the ozone produced by recombination of atomic oxygen is vibrationally excited in the nu3 mode. Further, laboratory data are available for the physical quenching of nu3, its radiative relaxation, and its radiative excitation by resonant absorption. It is shown that these chemical and physical processes are likely to result in substantial departures from local thermodynamic equilibrium for the nu3 mode of ozone in the mesosphere and therefore have important effects on infrared mesospheric ozone measurments by emission in the 9.6-micron band. Implications of these effects on data from the Limb Infrared Monitor of the Stratosphere (Remsberg et al, 1984), particularly their night/day ratio, are discussed.

Description: Sm-Nd isotopic studies of anorthosites can be used to provide information on their ages of crystallization and metamorphism, contamination history, and mantle sources. Proterozoic anorthosites in the Grenville and Nain Provinces of eastern North America crystallized between about 1100 and 1600 Ma, and some were metamorphosed at about 1000 Ma. Grenville Province anorthosite massifs were derived from depleted mantle. It is not clear whether massifs and related mafic intrusions throughout the Nain Province of Labrador were derived from enriched mantle, or were contaminated by early Archean (greater than 3500 Ma) silicic crustal materials, heretofore thought to be restricted to coastal Labrador.

Description: The relation between rotational modulation of the ultraviolet solar irradiance and variations in atmospheric ozone has been investigated using Fourier transform harmonic analysis and cross-correlations. Ozone variations with the same period and phase as 13.5 day or 27-day solar flux variations occur at tropical and subtropical latitudes over a range of pressure levels centered about 3 mbar. The solar-forced oscillation is stronger in the summer hemisphere; as temperature-related variations would be stronger in winter. Changes in solar irradiance over the 11-year cycle can be estimated by scaling rotational modulation. Using this estimate and the ozone-sun relation obtained for rotational modulation yields solar cycle changes of 3.5 percent in 3 percent mixing ratio comparable to that predicted from halocarbons and 0.7 mbar in total ozone.

Description: Ozone profile data from the Solar Backscattered Ultraviolet Instrument on Nimbus 7 from 1979 to the present and clear cases of ozone destruction associated with five sudden proton events (SPEs) on June 7, 1979, August 21, 1979, October 13-14, 1981, July 13, 1982, and December 8, 1982 are found. During the SPE on July 13, 1982, the largest of this solar cycle, no depletion at all at 45 km is observed, but there is a 15 percent ozone depletion at 50 km increasing to 27 percent at 55 km, all at a solar zenith angle of 85 deg. A strong variation of the observed depletion with solar zenith angle is found, with maximum depletion occurring at the largest zenith angles (near 85 deg) decreasing to near zero for angles below about 70 deg. The observed depletion is short lived, disappearing within hours of the end of the SPE.

Description: Absorption cross-sections of ozone have been measured over the range 230 nm to 350 nm, and for temperatures 200 K to 300 K, with improved photometric accuracy and spectral resolution. These measurements are referred to the cross-section at the 253.65 nm mercury line by the Hearn value (1961), 1147 x 10 to the -20th/sq cm, and show an internal consistency of + or - 1 percent. Tables of ozone absorption cross-section in the ultraviolet have been prepared for intervals of 0.05 nm over the range 245 to 340 nm. at each wavelength entry in the table a set of coefficients has been derived that permits the cross-section to be computed as a function of temperature, between 200 K and 300 K, with an accuracy of 1 percent.

Description: Intercomparisons of remote and in-situ techniques used to measure stratospheric ozone are made using results obtained on the Balloon Intercomparison Campaign of 1982 and 1983. Two in-situ and four remote instruments participated. These included ECC ozonesondes, a UV absorption photometer, and microwave emission, IR emission, and absorption spectrometers. Differences are generally less than 15 percent, and are within the quoted error bars. Flights which involved different sets of instruments were made on four separate days, and results are intercompared in plots of ozone density versus altitude. A careful assessment of errors was made for each instrument, and a plot of absolute errors versus altitude is given.

Description: In the Balloon Ozone Intercomparison Campaign (BOIC), several in situ UV absorption photometers, two solar UV absorption photometers, electrochemical sondes, and a mass spectrometer were intercompared in three flight missions. Concurrent data from Umkehr and satellite observations are also intercompared. The National Bureau of Standards provided a 'standard' ozone source for intercomparing the in situ instruments and ground pressure. Preliminary findings indicate that the standard deviation of the sensitivities among 17 instruments against the NBS reference was about 11 percent. These differences appear in flight at the lower levels and change at higher altitudes, indicating height-dependent errors. The difference among five in-situ UV photometers flown together ranged by plus or minus 8 percent during ascent to about 41 km. During float at 42 km, the difference nearly doubled. During descent, the difference decreased to about 4 percent, which is much closer to the expected accuracy of these instruments. Results from UV solar radiometers have been systematically higher than those from UV photometers by 15 to 20 percent - a very important disagreement that needs to be resolved.

Description: Utilizing independent estimates of ozone and temperature fields from the SBUV (Nimbus 7) and NOAA operational satellites, respectively, for the period 1978-1981, the coefficient of variation between the two parameters is determined. This coefficient is defined as A = Delta-O3 x (T)/Delta T x (O3) wehre Delta is an incremental change in either temperature or ozone and the bracket is a mean state. In practice, A is determined on a daily basis by regression of ozone mixing ratio versus temperature around a latitude circle during the winter season and the bracket value is the daily zonal average. This has the advantage of keeping the solar zenith angle fixed for a daily value while allowing it to change during the season. This is done at 30, 10, 5, 2, and 1 mb from 20 deg to 60 deg latitude in both hemispheres. The results are summarized and compared with those determined from a one-dimensional photochemical model applied to different latitudes.

Description: Until recently, Umkehr data taken by 20 Dobson stations around the world have been the principal source of information about the behavior of upper stratospheric ozone. Umkehr results are used also for detecting drifts in satellite instruments and for determining intersatellite biases. However, a systematic evaluation of the quality of Umkehr data taken by the various stations has been lacking. Five years of ozone profile data from the Solar Backscattered Ultraviolet (SBUV) have been used to examine and intercompare the quality of Umkehr stations, and to assess the degradation of their performance after the El Chichon volcano eruption in southern Mexico. In contrast to Umkehr, the SBUV ozone measurments in layers 7 through 9 (1-8 mb) were unaffected by the massive amounts of dust and gases ejected by El Chichon.

Description: Ozone profile data are intercompared with those from the LIMS and SBUV flown on Nimbus 7, SAGE flown on Atmospheric Explorer Mission 2, and the Ultraviolet and Infrared Spectrometers flown on the SME. Ozone data were derived from the measurements with independently derived processing algorithms. The data cover different time periods and have different spatial and temporal resolutions. The similarities and differences between the individual data sets are determined with effort focused on directly comparing the spacecraft data sets in the form of individual profiles, zonal mean profiles and global analyses. Comparisons with ground data from balloons and Umkehr stations are made.

Description: A study is reported based on four years of SBUV data and five years of published ozonesonde data. The Solar Backscatter Ultraviolet Experiment (SBUV), has been operating continuously since November 1978, providing some 1000 profiles per day. The data are of excellent quality from about 50 km down to the tropopause. By combining the satellite-derived data with ozonesonde data, a data base of standard ozone profiles and variance/covariance matrices has been developed, defining the observed variation of the atmospheric ozone around the mean profiles. Highlights of the data base are presented along with parameters of an analytical fit that describes the essential features of the data set. Although the data base was created specifically for use as a priori profiles for SBUV and Umkehr retrievals, it should be useful in other remote sensing applications and in modeling the UV radiation fields in the stratosphere. Another suggested application is in estimating ozone above the peak altitude of the ozone balloons.

Description: An overview is given of the Nimbus 7 Total Ozone Mapping Spectrometer (TOMS) which was designed to observe the spatial characteristics of total ozone that were not resolved by the nadir-viewing Nimbus BUV and SBUV instruments. At the wavelengths suitable for total ozone measurements, the radiance is large enough that the entire daytime atmosphere could be surveyed with about 50-km resolution from a polar orbiting satellite. The resulting high spatial resolution TOMS ozone images are found to reflect the internal dynamic structure of the lower atmosphere. Features which can be identified and tracked include: planetary wave scale troughs and ridges, mesoscale cutoff lows and rapidly moving troughs, jet stream confluence and difluence areas, hurricanes, and polar night lows. These features control the ozone above any given location and account for nearly all the variance in the total ozone. The instrument has been used to track the volcanic eruption clouds from El Chichon, Mount St. Helens, Alaid, and smaller eruptions such as Galunggung. It would be feasible to use a similar instrument on a geostationary platform to obtain half-hourly maps. Determination of the vertical ozone distribution in the lower stratosphere using Radon transform principles would be of importance in measuring jet stream folds and the related troposphere-stratosphere exchange.

Description: Results are given of SBUV (solar backscattered ultraviolet instrument) measurements from Nimbus 7, when - one day per month - it is operated in a spectral scan mode, scanning from 160 nm to 400 nm in 0.2-nm steps. By measuring the intensity of a series of nitric oxide (NO) gamma band fluorescence features in this wavelength range, it has been possible to estimate the amount of NO in the upper stratosphere and mesosphere. The background of atmospherically scattered sunlight normally masks the much weaker NO gamma band emission, but these emission features are discriminated by subtracting a synthetic spectrum calculated for a model atmosphere that includes only Rayleigh scattering and absorption by ozone and oxygen. The resulting difference plot clearly reveals features resulting from processes not included in the simple model, such as NO gamma band emission. Nitric oxide is inferred by measuring the absolute intensity of various bands relative to the adjacent background and relating this intensity to total NO above an altitude determined by the backscattering contribution function for that band. NO observations near the solstice and at various latitudes are reported.

Description: Analyses of a three-year time series of rocket ozone measurements at Wallops Island, VA, a set of rocket ozone soundings across the Southern Hemisphere, and rocket soundings at Fort Churchill, Manitoba are reported. Evidence is obtained that the NOx budget is not simply explained by oxidation of biospheric nitrous oxide. I 1-D time-dependent photochemical model is used to compute the amount of NO2 required to maintain odd oxygen in a steady state after accounting for Chapman, odd hydrogen, and odd chlorine reactions. At Wallops Island, a mid-latitude station, the inferred seasonal variation of NOx is small with the fall and winter mixing ratios about 20 percent greater than the spring and summer values. The soundings at Fort Churchill require about the same NOx amount as at Wallops Island in the spring and summer months but more than twice this amount in late fall and winter. Results indicate that the nitrous oxide source of NOx is supplemented by a polar source during the fall and winter months. This is consistent with the descent of thermospheric air with its high nitric oxide content during the period of strong cooling in the polar night.

Description: An analysis is given of some results of the LIMS experiment, launched on the Nimbus 7 satellite on October 24, 1978 for the purpose of sounding the middle atmosphere composition and structure. One of the LIMS channels was centered spectrally in the 6 micron region to measure vertical profiles and the global distribution of NO2 mixing ratio. The data reveal that NO2 is highly variable with altitude, latitude, longitude and time. Steep latitudinal gradients are observed in high latitude winter ('the NO2 cliff') when a wave number one pattern exists in the geopotential height field. The daytime gradient under these conditions is less than that for the night. The summer hemisphere column amount is considerably greater (by a factor of about 2) than in the winter hemisphere. The largest mixing ratios occur at 4 mb at night (about 20 ppbv) and 9 mb in the day (about 7 ppbv) and the latitude region of the peak mixing ratio is skewed toward the Southern Hemisphere. The mixing ratio variability is greater at night (by a factor of about 2 to 3) than in the day and there are significant long-term changes revealed by mixing ratio time series analyses on a pressure vs time grid.