ALBERT

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: The Experimental Cloud Lidar Pilot Study (ECLIPS) was initiated to obtain statistics on cloud-base height, extinction, optical depth, cloud brokenness, and surface fluxes. Two observational phases have taken place, in October-December 1989 and April-July 1991, with intensive 30-day periods being selected within the two time intervals. Data are being archived at NASA Langley Research Center and, once there, are readily available to the international scientific community. This article describes the scale of the study in terms of its international involvement and in the range of data being recorded. Lidar observations of cloud height and backscatter coefficient have been taken from a number of ground-based stations spread around the globe. Solar shortwave and infrared longwave fluxes and infrared beam radiance have been measured at the surface wherever possible. The observations have been tailored to occur around the overpass times of the NOAA weather satellites. This article describes in some detail the various retrieval methods used to obtain results on cloud-base height, extinction coefficient, and infrared emittance, paying particular attention to the uncertainties involved.

Description: A diagnostic analysis of the VVP (volume velocity processing) retrieval method is presented, with emphasis on understanding the technique as a linear, multivariate regression. Similarities and differences to the velocity-azimuth display and extended velocity-azimuth display retrieval techniques are discussed, using this framework. Conventional regression diagnostics are then employed to quantitatively determine situations in which the VVP technique is likely to fail. An algorithm for preparation and analysis of a robust VVP retrieval is developed and applied to synthetic and actual datasets with high temporal and spatial resolution. A fundamental (but quantifiable) limitation to some forms of VVP analysis is inadequate sampling dispersion in the n space of the multivariate regression, manifest as a collinearity between the basis functions of some fitted parameters. Such collinearity may be present either in the definition of these basis functions or in their realization in a given sampling configuration. This nonorthogonality may cause numerical instability, variance inflation (decrease in robustness), and increased sensitivity to bias from neglected wind components. It is shown that these effects prevent the application of VVP to small azimuthal sectors of data. The behavior of the VVP regression is further diagnosed over a wide range of sampling constraints, and reasonable sector limits are established.

Description: Analysis of the version 16 Halogen Occultation Experiment (HALOE) CH4 data shows that this long-lived trace gas is well correlated with potential vorticity (PV) computed from National Meteorological Center balanced winds. Analyzing late September and October 1992 data, we show that very low CH4 values are confined to the interior of a vortex edge defined by the maximum gradient in PV. The CH4 and HF time tendency is used to estimate the descent rate in the Antarctic vortex. After removing a component of the trend correlated with the HALOE sampling pattern, we compute the lower stratosphere vertical descent rates and net heaing rates in the spring Antarctic vortex. Our computations of the spring Antarctic vortex heating rates give -0.5 to -0.1 K/day. Over the winter season, the overall lower stratospheric descent rate averages about 1.8-1.5 km/month. These computations are in line with radiative transfer estimates of the heating and descent rate. The HALOE data thus appear to be consistent with the picture of an isolated lower stratospheric Antarctic vortex.

Description: Using delta C-13 measurements in atmospheric CO2 from a cooperative global air sampling network, we determined the partitioning of the net uptake of CO2 between ocean and land as a function of latitude and time. The majority of delta C-13 measurements were made at the Institute of Arctic and Alpine Research (INSTAAR) of the University of Colorado. We perform an inverse deconvolution of both CO2 and delta C-13 observations, using a two-dimensional model of atmospheric transport. Also, the discrimination against C-13 by plant photosynthesis, as a function of latitude and time, is calculated from global runs of the simple biosphere (SiB) model. Uncertainty due to the longitudinal structure of the data, which is not represented by the model, is studied through a bootstrap analysis by adding and omitting measurement sites. The resulting error estimates for our inferred sources and sinks are of the order of 1 GTC (1 GTC = 10(exp 15) gC). Such error bars do not reflect potential systematic errors arising from our estimates of the isotopic disequilibria between the atmosphere and the oceans and biosphere, which are estimated in a separate sensitivity analysis. With respect to global totals for 1992 we found that 3.2 GTC of carbon dissolved into the ocean and that 1.5 GTC were sequestered by land ecosystems. Northern hemisphere ocean gyres north of 15 deg N absorbed 2.7 GTC. The equatorial oceans between 10 deg S and 10 deg N were a net source to the atmosphere of 0.9 GTC. We obtained a sink of 1.6 GTC in southern ocean gyres south of 20 deg S, although the deconvolution is poorly constrained by sparse data coverage at high southern latitudes. The seasonal uptake of CO2 in the northern gyres appears to be correlated with a bloom of phytoplankton in surface waters. On land, northern temperate and boreal ecosystems between 35 deg N and 65 deg N were found to be a major sink of CO2 in 1992, as large as 3.5 GTC. Northern tropical ecosystems (equator-30 deg N) appear to be a net source to the source to the atmosphere of 2 GTC which could reflect biomass burning. A small sink, 0.3 GTC, was inferred for southern tropical ecosystems (30 deg S-equator).

Description: Global sets of surface radiation budget (SRB) have been obtained from satellite programs. These satellite-based estimates need validation with ground-truth observations. This study validates the estimates of monthly mean surface insolation contained in two satellite-based SRB datasets with the surface measurements made at worldwide radiation stations from the Global Energy Balance Archive (GEBA). One dataset was developed from the Earth Radiation Budget Experiment (ERBE) using the algorithm of Li et al. (ERBE/SRB), and the other from the International Satellite Cloud Climatology Project (ISCCP) using the algorithm of Pinker and Laszlo and that of Staylor (GEWEX/SRB). Since the ERBE/SRB data contain the surface net solar radiation only, the values of surface insolation were derived by making use of the surface albedo data contained GEWEX/SRB product. The resulting surface insolation has a bias error near zero and a root-mean-square error (RMSE) between 8 and 28 W/sq m. The RMSE is mainly associated with poor representation of surface observations within a grid cell. When the number of surface observations are sufficient, the random error is estimated to be about 5 W/sq m with present satellite-based estimates. In addition to demonstrating the strength of the retrieving method, the small random error demonstrates how well the ERBE derives from the monthly mean fluxes at the top of the atmosphere (TOA). A larger scatter is found for the comparison of transmissivity than for that of insolation. Month to month comparison of insolation reveals a weak seasonal trend in bias error with an amplitude of about 3 W/sq m. As for the insolation data from the GEWEX/SRB, larger bias errors of 5-10 W/sq m are evident with stronger seasonal trends and almost identical RMSEs.

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

Description: 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: In this paper, we investigate the relative importance of local vs remote control on cloud radiative forcing using a cumulus ensemble model. It is found that cloud and surface radiation forcings are much more sensitive to the mean vertical motion assoicated with large scale tropical circulation than to the local SST (sea surface temperature). When the local SST is increased with the mean vertical motion held constant, increased surface latent and sensible heat flux associated with enhanced moisture recycling is found to be the primary mechanism for cooling the ocean surface. Large changes in surface shortwave fluxes are related to changes in cloudiness induced by changes in the large scale circulation. These results are consistent with a number of earlier empirical studies, which raised concerns regarding the validity of the cirrus-thermostat hypothesis (Ramanathan and Collins, 1991). It is argued that for a better understanding of cloud feedback, both local and remote controls need to be considered and that a cumulus ensemble model is a powerful tool that should be explored for such purpose.

Description: The probability of polar stratospheric cloud (PSC) occurrence in the Antarctic and Arctic has been estimated using Stratospheric Aerosol Measurement (SAM) II aerosol extinction data from 1978 to 1989. Antarctic PSCs are typically observed by SAM II from mid-May to early November, with a maximum zonal average probability of about 0.6 at 18-20 km in August. The typical Arctic PSC season extends only from late November to early March, with a peak zonal average probability of about 0.1 in early February at 20-22 km. There is considerable year-to-year variability in Arctic PSC sightings because of changes in the dynamics of the northern polar vortex. Year-to-year variability in Antarctic sightings is most prominent in the number of late season clouds. Maximum PSC sighting probabilities in both polar regions occur in the region from 90 deg W through the Greenwich meridian to 90 deg E, where temperatures are coldest on average. Arctic sighting probabilities approach zero outside this region, but clouds have been sighted in the Antarctic at all longitudes during most months. Inferred PSC formation temperatures remain constant throughout the Arctic winter and are similar to those in early Antarctic winter. PSC formation temperatures in the Antarctic drop markedly in the 15 to 20-km region by September, a pattern consistent with the irreversible loss of HNO3 and H2O vapor in sedimenting PSC particles.

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

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

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

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

Description: A one-speed Boltzmann transport theory, with diffusion approximations, is applied to study the radiative transfer properties of lightning in optically thick thunderclouds. Near-infrared (lambda = 0.7774 micrometers) photons associated with a prominent oxygen emission triplet in the lightning spectrum are considered. Transient and spatially complex lightning radiation sources are placed inside a rectangular parallelepiped thundercloud geometry and the effects of multiple scattering are studied. The cloud is assumed to be composed of a homogeneous collection of identical spherical water droplets, each droplet a nearly conservative, anisotropic scatterer. Conceptually, we treat the thundercloud like a nuclear reactor, with photons replaced by neutrons, and utilize standard one-speed neutron diffusion techniques common in nuclear reactor analyses. Valid analytic results for the intensity distribution (expanded in spherical harmonics) are obtained for regions sufficiently far from sources. Model estimates of the arrival-time delay and pulse width broadening of lightning signals radiated from within the cloud are determined and the results are in good agreement with both experimental data and previous Monte Carlo estimates. Additional model studies of this kind will be used to study the general information content of cloud top lightning radiation signatures.

Description: Participants in this workshop, which convened in Venice, Italy, 6-8 May 1993, met to consider the current state of climate monitoring programs and instrumentation for the purpose of climatological prediction on short-term (seasonal to interannual) timescales. Data quality and coverage requirements for definition of oceanographic heat and momentum fluxes, scales of inter- and intra-annual variability, and land-ocean-atmosphere exchange processes were examined. Advantages and disadvantages of earth-based and spaceborne monitoring systems were considered, as were the structures for future monitoring networks, research programs, and modeling studies.

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

Description: The work described here makes it possible to identify anomalous wind behavior such as the nighttime meridional wind abatements that occur at F-region heights. A new analysis technique uses a simple empirical wind model to simulate measurements of 'normal' winds (as measured by the Neutral Atmosphere and Temperature Experiment (NATE) that flew on the Atmosphere Explorer-E (AE-E)) to highlight anomalous wind measurements made by the satellite while in circular orbits at 270-290 km altitude. Our approach is based on the recognition that the 'in orbit' wind variation must show the combined effects of the diurnal wind variation as seen from the ground with the latitude variation of the satellite orbit. For the data period 77250-78035 examined thus far, the wind abatement always occurred with a corresponding pressure or temperature maximum, and was detected on 12 out of the 36 nights with data. This study has revealed that the wind abatement occur only during or shortly after increases in solar extreme ultraviolet (EUV) flux, as indicated by daily radio flux measurements. In the past, nighttime wind reversals at mid-latitudes have been associated with increased geomagnetic activity. This study indicates that intensified solar EUV heating may be responsible for anomalous thermospheric nighttime winds at mid-latitudes.

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

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

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

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

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

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

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

Description: Five years of meteorological and hydrological data from a typical New England watershed where winter snow cover is significant were used to drive and validate two off-line land surface schemes suitable for use in the Goddard Institute for Space Studies (GISS) general circulation model (GCM): a baseline scheme that does not model the physics of a snowpack and therefore, neglects the insulating properties of snow cover; and a modified scheme in which a three-layer snowpack is modeled. Comparing baseline model results with validation data reveals several model deficiencies. Surface radiation temperatures could not adequately be modeled and the ground froze to unreasonable depths. Furthermore, because of ground cooling resulting from large surface heat fluxes to the atmosphere from the uninsulated surface, deeper model layers did not unfreeze until midsummer. As such, the normal hydrologic processes of runoff, ground water infiltration, and movement, etc., are compromised for a good part of the year. With the inclusion of a simple three-layer snow model into the baseline model, not only are the ground and surface radiation temperatures adequately modeled but all the features of snowpack ripening that characterize pack growth/ablation are simulated.

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

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

Description: A three-year climatology of satellite-estimated rainfall for the warm season for the southwest United States and Mexico has been derived from data from the Special Sensor Microwave Imager (SSM/I). The microwave data have been stratified by month (June, July, August), year (1988, 1989, 1990), and time of day (morning and evening orbits). A rain algorithm was employed that relates 86-GHz brightness temperatures to rain rate using a coupled cloud-radiative transfer model. Results identify an early evening maximum in rainfall along the western slope of the Sierra Madre Occidental during all three months. A prominent morning rainfall maximum was found off the western Mexican coast near Mazatlan in July and August. Substantial differences between morning and evening estimates were noted. To the extent that three years constitute a climatology, results of interannual variability are presented. Results are compared and contrasted to high-resolution (8 km, hourly) infrared cloud climatologies, which consist of the frequency of occurrence of cloud colder than -38 C and -58 C. This comparison has broad implications for the estimation of rainfall by simple (cloud threshold) techniques. By sampling the infrared data to approximate the time and space resolution of the microwave, we produce ratios (or adjustment factors) by which we can adjust the infrared rain estimation schemes. This produces a combined microwave/infrared rain algorithm for monthly rainfall. Using a limited set of raingage data as ground truth, an improvement (lower bias and root-mean-square error) was demonstrated by this combined technique when compared to either method alone. The diurnal variability of convection during July 1990 was examined using hourly rain estimates from the Geostationary Operational Environmental Satellite (GOES) precipitation index and the convective stratiform technique, revealing a maximum in estimated rainfall from 1800 to 2100 local time. It is in this time period when the SSM/I evening orbit occurs. A high-resolution topographic database was available to aid in interpreting the influence of topography on the rainfall patterns.

Description: An advanced Microwave Precipitation Radiometer (AMPR) has been developed and flown in the NASA ER-2-high-altitude aircraft for imaging various atmospheric and surface processes, primarily the internal structure of rain clouds. The AMPR is a scanning four-frequency total power microwave radiometer that is externally calibrated with high-emissivity warm and cold loads. Separate antenna systems allow the sampling of the 10.7- and 19.35-GHz channels at the same spatial resolution, while the 37.1- and 85.5-GHz channels utilize the same multifrequency feedhorn as the 19.35-GHz channel. Spatial resolutions from an aircraft altitude of 20-km range from 0.6 km at 85.5 GHz to 2.8 km at 19.35 and 10.7 GHz. All channels are sampled every 0.6 km in both along-track and cross-track directions, leading to a contiguous sampling pattern of the 85.5-GHz 3-dB beamwidth footprints, 2.3X oversampling of the 37.1-GHz data, and 4.4X oversampling of the 19.35- and 10.7-GHz data. Radiometer temperature sensitivities range from 0.2 to 0.5 C. Details of the system are described, including two different calibration systems and their effect on the data collected. Examples of oceanic rain systems are presented from Florida and the tropical west Pacific that illustrate the wide variety of cloud water, rainwater, and precipitation-size ice combinations that are observable from aircraft altitudes.

Description: A general circulation model (GCM) is used to model global lightning distributions and frequencies. Both total and cloud-to-ground lightning frequencies are modeled using parameterizations that relate the depth of convective clouds to lightning frequencies. The model's simulations of lightning distributions in time and space show good agreement with available observations. The model's annual mean climatology shows a global lightning frequency of 77 flashes per second, with cloud-to-ground lightning making up 25% of the total. The maximum lightning activity in the GCM occurs during the Northern Hemisphere summer, with approximately 91% of all lightning occurring over continental and coastal regions.

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

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

Description: An automated scheme to characterize precipitation echoes within small windows in the radar field is presented and applied to previously subjectively classified tropical rain cloud systems near Darwin, Australia. The classification parameters are (a) E(sub e), effective efficiency, as determined by cloud-top and cloud-base water vapor saturation mixing ratios; (b) BBF, brightband fraction, as determined by the fraction of the radar echo area in which the maximal reflectivity occurs within +/- 1.5 km of the 0 C isotherm level; and (c) del(sub r) Z, radial reflectivity gradients (dB/km). These classification criteria were applied to tropical rain cloud systems near Darwin, Australia, and to winter convective rain cloud systems in Israel. Both sets of measurements were made with nearly identical networks of C-band radars and rain gauge networks. The results of the application of these objective classification criteria to several independently predetermined rain regimes in Darwin have shown that better organized rain systems have smaller del(sub r) Z and larger BBF. Similarly, smaller del(sub r)Z and larger BBF were also observed from maritime rain cloud systems, as compared to continental rain cloud systems with the same degree of organization. Continental rain cloud system, regardless of their degree of organization, have larger depths, as expressed by E(sub e). The rainfall analyses presented in this study are based exclusively on rain gauge measurements, while radar information was used only to classify the individual gauge measurements.

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

Description: This work is directed toward approximating the evolution of forecast error covariances for data assimilation. The performance of different algorithms based on simplification of the standard Kalman filter (KF) is studied. These are suboptimal schemes (SOSs) when compared to the KF, which is optimal for linear problems with known statistics. The SOSs considered here are several versions of optimal interpolation (OI), a scheme for height error variance advection, and a simplified KF in which the full height error covariance is advected. To employ a methodology for exact comparison among these schemes, a linear environment is maintained, in which a beta-plane shallow-water model linearized about a constant zonal flow is chosen for the test-bed dynamics. The results show that constructing dynamically balanced forecast error covariances rather than using conventional geostrophically balanced ones is essential for successful performance of any SOS. A posteriori initialization of SOSs to compensate for model - data imbalance sometimes results in poor performance. Instead, properly constructed dynamically balanced forecast error covariances eliminate the need for initialization. When the SOSs studied here make use of dynamically balanced forecast error covariances, the difference among their performances progresses naturally from conventional OI to the KF. In fact, the results suggest that even modest enhancements of OI, such as including an approximate dynamical equation for height error variances while leaving height error correlation structure homogeneous, go a long way toward achieving the performance of the KF, provided that dynamically balanced cross-covariances are constructed and that model errors are accounted for properly. The results indicate that such enhancements are necessary if unconventional data are to have a positive impact.

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

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

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

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

Description: A cumulus ensemble model is used to study the tropical water and energy cycles and their role in the climate system. The model includes cloud dynamics, radiative processes, and microphysics that incorporate all important production and conversion processes among water vapor and five species of hydrometeors. Radiative transfer in clouds is parameterized based on cloud contents and size distributions of each bulk hydrometeor. Several model integrations have been carried out under a variety of imposed boundary and large-scale conditions. In Part 1 of this paper, the primary focus is on the water and heat budgets of the control experiment, which is designed to simulate the convective - radiative equilibrium response of the model to an imposed vertical velocity and a fixed sea surface temperature at 28 C. The simulated atmosphere is conditionally unstable below the freezing level and close to neutral above the freezing level. The equilibrium water budget shows that the total moisture source, M(sub s), which is contributed by surface evaporation (0.24 M(sub s)) and the large-scale advection (0.76 M(sub s)), all converts to mean surface precipitation bar-P(sub s). Most of M(sub s) is transported verticaly in convective regions where much of the condensate is generated and falls to surface (0.68 bar-P(sub s)). The remaining condensate detrains at a rate of 0.48 bar-P(sub s) and constitutes 65% of the source for stratiform clouds above the melting level. The upper-level stratiform cloud dissipates into clear environment at a rate of 0.14 bar-P(sub s), which is a significant moisture source comparable to the detrained water vapor (0.15 bar-P(sub s)) to the upper troposphere from convective clouds. In the lower troposphere, stratiform clouds evaporate at a rate of 0.41 bar-P(sub s), which is a more dominant moisture source than surface evaporation (0.22 bar-P(sub s)). The precipitation falling to the surface in the stratiform region is about 0.32 bar-P(sub s). The associated latent heating in the water cycle is the dominant source in the heat budget that generates a net upward motion in convective regions, upper stratiform regions (above the freezing level), and a downward motion in the lower stratiform regions. The budgets reveal a cycle of water and energy resulted from radiation-dynamic-convection interactions that maintain equilibrium of the atmosphere.

Description: Radar or satellite observations of an area generate sequences of rain-rate maps. From a gridded map a histogram of rain rates can be obtained representing the relative areas occupied by rain rates of various strengths. The histograms vary with time as precipitating systems in the area evolve and decay and amounts of convective and stratiform rain in the area change. A method of decomposing the histograms into linear combinations of a few empirical distributions with time-dependent coefficients is developed, using principal component analysis as a starting point. When applied to a tropical Atlantic dataset (GATE), two distributions emerge naturally from the analysis, resembling stratiform and convective rain-rate distributions in that they peak at low and high rain rates, respectively. The two 'modes' have different timescales and only the high-rain-rate mode has a statistically significant diurnal cycle. The ability of just two modes to describe rain variabiltiy over an area can explain why methods of estimating area-averaged rain rate from the area covered by rain rates above a certain threshold are so successful.

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

Description: Carbon monoxide (CO) and methane (CH4) were measured in the 0.15- to 6-km portion of the troposphere over subarctic and boreal landscapes of midcontinent and eastern Canada during July - August 1990. In the mid-continent region, Arctic air entering the region was characterized by relatively uniform CO concentrations (86-108 parts per billion by volume (ppbv)) and CH4 concentrations (1729-1764 ppbv). Local biomass burning and long-range transport of CO into the area from industrial/urban sources and distant fires did frequently produce enhanced and variable concentrations. Emissions of CH4 from the Hudson Bay lowlands was the primary source for enhanced and variable concentrations, especially at altitudes of 0.15-1 km. In eastern Canada, most of the observed variability in CO and CH4 was similar in origin to the phenomena described for the midcontinent region. However, unexpectedly low concentrations of CO (51 ppbv) and CH4 (1688 ppbv) were measured in the midtroposphere on several flights. Combined meteorological and chemical data indicated that the low CO-CH4 events were the result of long-range transport of tropical Pacific marine air to subarctic latitudes.

Description: We report features of acidic gases in the troposphere from 9 to 5000 m altitude above ground over the Canadian taiga in the summer of 1990. The measurements were conducted at a 30-m meteorological tower and from the NASA Wallops Electra aircraft as part of the joint U.S.-Canadian Arctic Boundary Layer Expedition (ABLE) 3B Northern Wetland Studies (NOWES). We sampled air for acidic gases using the mist chamber collector coupled with subsequent analysis using ion chromatography. At the tower we collected samples at two heights during a 13-day period, including diurnal studies. Using eddy flux and profile data, we estimated the biosphere/troposphere fluxes of nitric, formic, and acetic acids and sulfur dioxide. For the organic acids, emissions from the taiga in the afternoon hours and deposition during the predawn morning hours were observed. The flux intensities alone were however not high enough to explain the observed changes in mixing ratios. The measured deposition fluxes of nitric acid were high enough to have a significant influence on its mixing ratio in the boundary layer. On three days we measured vertical profiles of nitric, formic, and acetic acids through the lower to midtroposphere. We found that the chemical composition of the troposphere was extremely heterogenous. Pronounced layers of polluted air were readily apparent from our measurements. Local photochemical production and episodic long-range transport of trace components, originating from biomass burning and possibly industrial emissions, appear to have a strong influence on the composition of the troposphere and biosphere/troposphere fluxes of acidic gases at this site.

Description: In this study we perform an error analysis for cloud-top pressure retrieval using the High-Resolution Infrared Radiometric Sounder (HIRS/2) 15-microns CO2 channels for the two-layer case of transmissive cirrus overlying an overcast, opaque stratiform cloud. This analysis includes standard deviation and bias error due to instrument noise and the presence of two cloud layers, the lower of which is opaque. Instantaneous cloud pressure retrieval errors are determined for a range of cloud amounts (0.1-1.0) and cloud-top pressures (850-250 mb). Large cloud-top pressure retrieval errors are found to occur when a lower opaque layer is present underneath an upper transmissive cloud layer in the satellite field of view (FOV). Errors tend to increase with decreasing upper-cloud effective cloud amount and with decreasing cloud height (increasing pressure). Errors in retrieved upper-cloud pressure result in corresponding errors in derived effective cloud amount. For the case in which a HIRS FOV has two distinct cloud layers, the difference between the retrieved and actual cloud-top pressure is positive in all cases, meaning that the retrieved upper-cloud height is lower than the actual upper-cloud height. In addition, errors in retrieved cloud pressure are found to depend upon the lapse rate between the low-level cloud top and the surface. We examined which sounder channel combinations would minimize the total errors in derived cirrus cloud height caused by instrument noise and by the presence of a lower-level cloud. We find that while the sounding channels that peak between 700 and 1000 mb minimize random errors, the sounding channels that peak at 300-500 mb minimize bias errors. For a cloud climatology, the bias errors are most critical.

Description: We examine the influence of aggregation errors on developing estimates of regional soil-CO2 flux from temperate forests. We find daily soil-CO2 fluxes to be more sensitive to changes in soil temperatures (Q(sub 10) = 3.08) than air temperatures (Q(sub 10) = 1.99). The direct use of mean monthly air temperatures with a daily flux model underestimates regional fluxes by approximately 4%. Temporal aggregation error varies with spatial resolution. Overall, our calibrated modeling approach reduces spatial aggregation error by 9.3% and temporal aggregation error by 15.5%. After minimizing spatial and temporal aggregation errors, mature temperate forest soils are estimated to contribute 12.9 Pg C/yr to the atmosphere as carbon dioxide. Georeferenced model estimates agree well with annual soil-CO2 fluxes measured during chamber studies in mature temperate forest stands around the globe.

Description: This study examines the role of synoptic-scale eddies during the development of persistent anticyclonic height anomalies over the central North Pacific in a general circulation model under perpetual January conditions. The General Circulation Model (GCM) replicates the basic characteristics of the evolution of the anomaly patterns found in observations. The life cycle is characterized by the rapid establishment of the major anomaly center and considerably longer maintenance and decay phases, which include the development of downstream anomaly centers. The simulation also shows a realistic evolution of synoptic-scale activity beginning with enhanced activity off the east coast of Asia prior to onset, followed by a northward shift of the Pacific storm track, which lasts throughout the maintenance phase. The initial enhancement of synoptic-scale eddy activity is associated with a large-scale cyclonic anomaly that developes over Siberia several days prior to the onset of the main anticyclonic anomaly over the central North Pacific. The observations, however, show considerable interdecadel variability in the details of the composite onset behavior; it is unclear whether this variability is real or whether it reflects differences in the data assimilation systems. The role of the time mean flow and synoptic-scale eddies in the development of the persistent Pacific anomalies is studied within the context of a kinetic energy budget in which the flow is decomposed into the time-mean, low-frequency (timescales longer than 10 days), and synoptic (timescales less than 6 days) components. The budget, which is carried out for the simulation at 500 mb, shows that the initial growth of the persistent anticyclonic anomalies is associated with barotropic conversions of energy, with approximately equal contributions coming from the mean flow and the synoptic-scale eddies. After onset the barotropic conversion from the mean flow dominates, whereas the decay phase is associated with baroclinic processes within the low-frequency flow.

Description: Winds derived from a stratospheric and tropospheric data assimilation system (STRATAN) are compared with balance winds derived from National Meteorological Center/Climate Analysis Center (NMC/CAC) heights. At middle latitudes in the lower stratosphere, the results show that STRATAN winds are comparable to the balance winds. In addition STRATAN winds provide useful horizontal divergence analyses, and hence, vertical velocity fields. More generally, the STRATAN winds are useful in a more extended domain than the balanced winds. In particular, they are useful in the Tropics and the upper stratosphere where the balanced winds fail. The assimilation also captures the quasi-biennial oscillation, but does not do a good job of representing tropical waves.

Description: The Institute for Naval Oceanography, in cooperation with Naval Research Laboratories and universities, executed the Data Assimilation and Model Evaluation Experiment (DAMEE) for the Gulf Stream region during fiscal years 1991-1993. Enormous effort has gone into the preparation of several high-quality and consistent datasets for model initialization and verification. This paper describes the preparation process, the temporal and spatial scopes, the contents, the structure, etc., of these datasets. The goal of DAMEE and the need of data for the four phases of experiment are briefly stated. The preparation of DAMEE datasets consisted of a series of processes: (1) collection of observational data; (2) analysis and interpretation; (3) interpolation using the Optimum Thermal Interpolation System package; (4) quality control and re-analysis; and (5) data archiving and software documentation. The data products from these processes included a time series of 3D fields of temperature and salinity, 2D fields of surface dynamic height and mixed-layer depth, analysis of the Gulf Stream and rings system, and bathythermograph profiles. To date, these are the most detailed and high-quality data for mesoscale ocean modeling, data assimilation, and forecasting research. Feedback from ocean modeling groups who tested this data was incorporated into its refinement. Suggestions for DAMEE data usages include (1) ocean modeling and data assimilation studies, (2) diagnosis and theoretical studies, and (3) comparisons with locally detailed observations.

Description: A new framework for interpreting the origin of the tropical intraseasonal oscillation (TISO), which avoids the speed and scale selection problems in the previous theories, is proposed in this study. In this interpretation TISO is viewed as an oscillation driven by an eastward moving convective region. This convective region consists of one or more super cloud clusters originating in the Indian Ocean and terminating in mid-Pacific, and is then followed by another convective region arising in the Indian Ocean in a period of 40-50 days. Additionally, a formal analogy is pointed out between super cloud clusters and the middle-latitude baroclinic wave packets. This study includes a simulation of TISO in a 2D model to support our interpretation. Experiments were conducted with four different convection schemes. The authors advocate that the successful simulation of TISO depends on the successful simulation of super clouds cluster, which in turn depends on the successful simulation of the life cycle of cloud clusters, which further in turn depends on the choice of cumulus convection scheme. What makes a cumulus convection scheme successful in simulating TISO is discussed.

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

Description: GEMPAK is a general meteorological software package developed at NASA/Goddard Space Flight Center. It includes programs to analyze and display surface, upper-air, and gridded data, including model output. There are very general programs to list, edit, and plot data on maps, to display profiles and time series, to draw and fill contours, to draw streamlines, to plot symbols for clouds, sky cover, and pressure tendency, and draw cross sections in the case of gridded data and sounding data. In addition, there are Barnes objective analysis programs to grid surface and upper-air data. The programs include the capabilities to derive meteorological parameters from those found in the dataset, to perform vertical interpolations of sounding data to different coordinate systems, and to compute an extensive set of gridded diagnostic quantities by specifying various nested combinations of scalars and vector arithmetic, algebraic, and differential operators. The GEMPAK 5.1 graphics/transformation subsystem, GEMPLT, provides device-independent graphics. GEMPLT also has the capability to display output in a variety of map projections or overlaid on satellite imagery. GEMPAK 5.1 is written in FORTRAN 77 and C-language and has been implemented on VAX computers under VMS and on computers running the UNIX operating system. During installation and normal use, this package occupies approximately 100Mb of hard disk space. The UNIX version of GEMPAK includes drivers for several graphic output systems including MIT's X Window System (X11,R4), Sun GKS, PostScript (color and monochrome), Silicon Graphics, and others. The VMS version of GEMPAK also includes drivers for several graphic output systems including PostScript (color and monochrome). The VMS version is delivered with the object code for the Transportable Applications Environment (TAE) program, version 4.1 which serves as a user interface. A color monitor is recommended for displaying maps on video display devices. Data for rendering regional maps is included with this package. The standard distribution medium for the UNIX version of GEMPAK 5.1 is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the VMS version of GEMPAK 5.1 is a 6250 BPI 9-track magnetic tape in DEC VAX BACKUP format. The VMS version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. This program was developed in 1985. The current version, GEMPAK 5.1, was released in 1992. The package is delivered with source code. An extensive collection of subroutine libraries allows users to format data for use by GEMPAK, to develop new programs, and to enhance existing ones.

Description: GEMPAK is a general meteorological software package developed at NASA/Goddard Space Flight Center. It includes programs to analyze and display surface, upper-air, and gridded data, including model output. There are very general programs to list, edit, and plot data on maps, to display profiles and time series, to draw and fill contours, to draw streamlines, to plot symbols for clouds, sky cover, and pressure tendency, and draw cross sections in the case of gridded data and sounding data. In addition, there are Barnes objective analysis programs to grid surface and upper-air data. The programs include the capabilities to derive meteorological parameters from those found in the dataset, to perform vertical interpolations of sounding data to different coordinate systems, and to compute an extensive set of gridded diagnostic quantities by specifying various nested combinations of scalars and vector arithmetic, algebraic, and differential operators. The GEMPAK 5.1 graphics/transformation subsystem, GEMPLT, provides device-independent graphics. GEMPLT also has the capability to display output in a variety of map projections or overlaid on satellite imagery. GEMPAK 5.1 is written in FORTRAN 77 and C-language and has been implemented on VAX computers under VMS and on computers running the UNIX operating system. During installation and normal use, this package occupies approximately 100Mb of hard disk space. The UNIX version of GEMPAK includes drivers for several graphic output systems including MIT's X Window System (X11,R4), Sun GKS, PostScript (color and monochrome), Silicon Graphics, and others. The VMS version of GEMPAK also includes drivers for several graphic output systems including PostScript (color and monochrome). The VMS version is delivered with the object code for the Transportable Applications Environment (TAE) program, version 4.1 which serves as a user interface. A color monitor is recommended for displaying maps on video display devices. Data for rendering regional maps is included with this package. The standard distribution medium for the UNIX version of GEMPAK 5.1 is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the VMS version of GEMPAK 5.1 is a 6250 BPI 9-track magnetic tape in DEC VAX BACKUP format. The VMS version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. This program was developed in 1985. The current version, GEMPAK 5.1, was released in 1992. The package is delivered with source code. An extensive collection of subroutine libraries allows users to format data for use by GEMPAK, to develop new programs, and to enhance existing ones.

Description: On an hourly time-scale the different roles of the solar wind and interplanetary magnetic field (IMF) parameters on ground micropulsation activity can be better investigated than at longer time-scales. A long-term comparison between ground measurements made at L'Aquila and IMP 8 observations confirms the solar wind speed as the key parameter for the onset of pulsations even at low latitudes, although additional control of the energy transfer from the interplanetary medium to the earth's magnetosphere is clearly exerted by the cone angle. Above about 20 mHz the frequency of pulsations is confirmed to be closely related to the IMF magnitude while, in agreement with model predictions, the IMF magnitude is related to the amplitude of the local fundamental resonant mode. We provide an interesting example in which high resolution measurements simultaneously obtained in the foreshock region and on the ground show that external transversal fluctuations do not penetrate deep into the low latitude magnetosphere.

Description: This paper discusses in situ measurements of O3, aerosol number density, and aerosol size obtained during NASA Arctic Boundary Layer Expedition aircraft flights over the Alaskan Arctic region. The major source of summer O3 for the troposphere in the intrusion of stratospheric air and subsequent transport to lower altitudes. Photochemistry of mixed layer emissions and O3 transported from high northern latitude urban/industrialized areas do not appear to play major roles as sources of O3 for the Alaska region. O3 gradients reflect the loss at the surface and supply from the stratosphere. Free tropospheric O3 averaged 74 ppbv compared to 32 ppbv for the mixed layer. O3 loss mechanisms are a combination of the destruction via photochemistry, chemical reaction with surface emissions, and direct loss through deposition to the surface. The boreal forest in the most efficient of the O3 sinks and has the largest increase in aerosol number density relative to the free troposphere.

Description: A meteorological overview of the Arctic Boundary Layer Expedition (ABLE 3A) flight series is presented. Synoptic analyses of mid-tropospheric circulation patterns are combined with isentropic back trajectory calculations to describe the long-range (400-3000 km) atmospheric transport mechanisms and pathways of air masses to the Arctic and sub-Arctic regions of North America during July and August 1988. Siberia and the northern Pacific Ocean were found to be the two most likely source areas for 3-day transport to the study areas in Alaska. Transport to the Barrow region was frequently influenced by polar vortices and associated short-wave troughs over the Arctic Ocean, while the Bethel area was most often affected by lows migrating across the Bering Sea and the Gulf of Alaska, as well as ridges of high pressure which built into interior Alaska. July 1988 was warmer and dryer than normal over much of Alaska. As a result, the 1988 Alaska fire season was one of the most active of the past decade. Airborne lidar measurements verified the presence of biomass burning plumes on many flights, often trapped in thin subsidence layer temperature inversions. Several cases of stratosphere/troposphere exchange were noted, based upon potential vorticity analyses and aircraft lidar data, especially in the Barrow region and during transit flights to and from Alaska.

Description: The aerosol optical depth over the Konza Prairie, near Manhattan, Kansas, was recorded at various locations by five separate teams. These measurements were made in support of the First ISLSCP Field Experiment (FIFE) and used to correct imagery from a variety of satellite and aircraft sensors for the effects of atmospheric scattering and absorption. The results from one instrument are reported here for 26 days in 1987 and for 7 in 1989. Daily averages span a range of 0.05 to 0.28 in the midvisible wavelengths. In addition, diurnal variations are noted in which the afternoon optical depths are greater than those of the morning by as much as 0.07. A comparison between instruments and processing techniques used to determine these aerosol optical depths is provided. The first comparisons are made using summer 1987 data. Differences of as much as 0.05 (midvisible) are observed. Although these data allow reasonable surface reflectance retrievals, they do not agree to within the performance limits typically associated with these types of instruments. With an accuracy goal of 0.02 a preseason calibration/comparison experiment was conducted at a mountain site prior to the final field campaign in 1989. Good calibration data were obtained, and good agreement (0.01, midvisible) was observed in the retrieved optical depth acquired over the Konza. By comparing data from the surface instruments at different locations, spatial inhomogeneities are determined. Then, data from the airborne tracking sunphotometer allow one to determine variations as a function of altitude. Finally, a technique is proposed for using the in situ data to establish an instrument calibration.

Description: The underlying mean and variance properties of surface net radiation, sensible-latent heat fluxes and soil heat flux are studied over the densely instrumented grassland region encompassing FIFE. Flux variability is discussed together with the problem of scaling up to area-averaged fluxes. Results are compared and contrasted for cloudy and clear situations and examined for the influence of surface-induced biophysical controls (burn and grazing treatments) and topographic controls (aspect ratios and slope factors).

Description: Airborne measurements of flux densities of carbon dioxide CO2, sensible heat, and latent heat (H2O) obtained over the First ISLSCP Field Experiment (FIFE) site during three intensive field campaigns in 1987 and one in 1989 are examined to characterize the spatial and temporal variations of CO2 and energy transfer processes. These data were collected by the National Research Council Twin Otter using low-level flight patterns, all flown at constant pressure altitude during relatively clear days. The spatial variations are larger in 1989 than in 1987 and a higher correlation is observed between the fluxes and the surface features. The temporal patterns are easier to characterize with the relatively homogeneous situation of 1987. Functional relationships obtained between fluxes of CO2 and latent heat, CO2 fluxes and greenness index, latent heat fluxes and greenness index, and between sensible heat fluxes and surface air temperature differences are presented for one day in 1987 and one in 1989 as an example of the kind of information that can be obtained from grid flights at constant pressure altitude.

Description: Radiometric measurements at 90 GHz and three sideband frequencies near the peak water vapor absorption line of 183.3 GHz were made with Advanced Microwave Moisture Sounder (AMMS) aboard the NASA DC-8 aircraft during the Global Aerosol Backscatter Experiment (GLOBE) mission over the Pacific Ocean in November 1989. Some of the measurements over the high-latitude regions (greater than 50 deg N or 50 deg S) were analyzed for the retrieval of total precipitable water less than 0.5 g/sq cm both over land and ocean surfaces. The results show that total precipitable water from a relatively dry atmosphere could be estimated with high sensitivity from these radiometric measurements. The retrieved values over ocean surface show a decrease toward the polar region as expected. The retrieved total precipitable water over land correlates positively with the aircraft radar altitude. This positive correlation is expected because the aircraft radar altitude provides a measure of atmospheric water vapor burden above the surface. Retrieved high reflectivities over land surface at 90 GHz and 183 GHz are presumably related to snow cover on the ground. This suggests that radiometric measurements at these frequencies could be used to map snow at high-latitude regions.

Description: The most important investigations leading to the International Reference Ionosphere 1990 (IRI-90) are overviewed, and the latest version of the model is described. The shortcomings and limitations of the IRI-90 are pointed out, together with the ways of overcoming them. The list of studies that the IRI group has yet to carry out includes the investigations of magnetic storm effects as the highest priority. This paper discusses determinations of and the available data on the electron density, plasma temperatures, ion composition, and ion drift in the ionosphere, together with future improvements needed on these parameters.

Description: In the aftermath of the eruption of Mt. Pinatubo, multiwavelength stratospheric aerosol extinction measurements by the satellite-borne Stratospheric Aerosol and Gas Experiment (SAGE II) revealed the presence of a previously unobserved mode of aerosol that exhibited high extinction but a small inferred particle size. This mode may represent a transitional phase between the very small aerosol created by gas-to-particle conversion and a quasi-steady state, post-volcanic aerosol that exhibits both large extinction and large particle size. The presence of a transitional small aerosol mode may have a significant impact on chemical and radiative processes in the stratosphere.

Description: CO2 exchange rates were measured at selected tundra sites near Bethel, Alaska using portable, climate-controlled, instrumented enclosures. The empirically modeled exchange rate for a representative area of vegetated tundra was 1.2 +/- 1.2 g/sq m/d, compared to a tower-measured exchange over the same time period of 1.1 +.0- 1.2 g/sq m/d. Net exchange in response to varying light levels was compared to wet meadow and dry upland tundra, and to the net exchange measured by the micrometeoroidal tower technique. The multispectral reflectance properties of the sites were measured and related to exchange rates in order to provide a quantitative foundation for the use of satellite remote sensing to monitor biosphere/atmosphere CO2 exchange in the tundra biome.

Description: Enhanced concentrations of CH4 in the unpolluted atmospheric mixed layer over both Arctic and subarctic tundra landscapes are documented here using data from the NASA Arctic Boundary Layer Expedition (ABLE 3A). The CH4 concentration gradients were determined mainly by interactions of biogenic emission from wet tundra and turbulent mixing proceses. The gradient were most frequently associated with intrusion of upper tropospheric or stratospheric air into the midtroposphere, emissions from forest and tundra fires, and long-range transport of enhanced concentration of these gases from unidentified sources. Summertime haze layers exhibited midtropospheric enhancements of CH4 similar to those measured in winter Arctic events. The observations confirm the importance of Arctic and Subarctic wetland environments as a regional source of global atmospheric CH4.

Description: The oxidation of the earth's crust and the increase in atmospheric oxygen early in earth history have been linked to the accumulation of reduced carbon in sedimentary rocks. Trends in the carbon isotope composition of sedimentary organic carbon and carbonate show that during the Proterozoic aeon (2.5-0.54 Gyr ago) the organic carbon reservoir grew in size, relative to the carbonate reservoir. This increase, and the concomitant release of oxidizing power in the environment, occurred mostly during episodes of global rifting and orogeny.

Description: Recently available spectral O3 line parameters are summarized and progress in remote sensing based on these parameters is discussed. Particular attention is given to covering line positions, line intensities and linewidths for both the main isotopic species (O-16)3 and the isotopic variants (O-16)(O-18)(O-16) and (O-16)(O-16)(O-18).

Description: The paper presents Alouette 2/ISIS 1 data which furnish evidence for the stimulation of high-order D(n) resonances (n greater than 4) by topside sounders, confirming the prediction of Osherovich (1990) based on an analogy with naturally occurring narrowband magnetospheric emissions. The results indicate that observations of stimulated ionospheric emissions and naturally occurring magnetospheric emissions can be used as complementary data sets to address such fundamental questions as the nature of the excitation mechanism for these emissions and the nature of the waves (i.e., whether the waves predominantly electrostatic or with a significant magnetic component present), and the question of the dominance among the Dn, Dn(+), and Dn(-) resonances and the conditions required for all three to be present at the same time.

Description: Ground magnetic field perturbations recorded by the CANOPUS magnetometer network in the 7 to 13 MLT sector are used to examine how reconfigurations of the dayside polar ionospheric flow take place in response to north-south changes of the IMF. During the 6-h interval in question, IMF Bz oscillates between +/- 7 nT with about a 1-h period. Corresponding variations in the ground magnetic disturbance are observed which we infer are due to changes in ionospheric flow. Cross correlation of the data obtained from two ground stations at 73.5 deg magnetic latitude, but separated by about 2 hours in MLT, shows that changes in the flow are initiated in the prenoon sector (about 10 MLT) and then spread outward toward dawn and dusk with a phase speed of about 5 km/s over the longitude range about 8 to 12 MLT, slowing to about 2 km/s outside this range. Cross correlating the data from these ground stations with IMP 8 IMF Bz records produces a MLT variation in the ground response delay relative to the IMF which is compatible with these deduced phase speeds.

Description: The particle scattering and current sheet stability features in the geomagnetic tail during the phase of substorm growth were investigated using Tsyganenko's (1989) magnetic field model. In a study of four substorm events which were observed both in the high-altitude nightside tail and in the auroral ionosphere, the model magnetic field was adjusted to each case so as to represent the global field development during the growth phase of the substorms. The model results suggest that the auroral brightenings are connected with processes taking place in the near-earth region inside about 15 earth radii. The results also suggest that there is a connection between the chaotization of the electrons and the auroral brightenings at substorm onset.

Description: Alternative geometries for the present-day configuration of plate boundaries in northeast Asia and Japan are tested using NUVEL-1 and 256 horizontal earthquake slip vectors from the Japan and northern Kuril trenches. Statistical analysis of the slip vectors is used to determine whether the North American, Eurasian, or Okhotsk plate overlies the trench. Along the northern Kuril trench, slip vectors are well-fit by the NUVEL-1 Pacific-North America Euler pole, but are poorly fit by the Pacific-Eurasia Euler pole. Results for the Japan trench are less conclusive, but suggest that much of Honshu and Hokkaido are also part of the North American plate. The simplest geometry consistent with the trench slip vectors is a geometry in which the North American plate extends south to 41 deg N, and possibly includes northern Honshu and southern Hokkaido. Although these results imply that the diffuse seismicity that connects the Lena River delta to Sakhalin Island and the eastern Sea of Japan records motion between Eurasia and North America, onshore geologic and seismic data define an additional belt of seismicity in Siberia that cannot be explained with this geometry. Assuming that these two seismic belts constitute evidence for an Okhotsk block, two published kinematic models for motion of the Okhotsk block are tested. The first model, which predicts motion of up to 15 mm/yr relative to North America, is rejected because Kuril and Japan trench slip vectors are fit more poorly than for the simpler geometry described above. The second model gives a good fit to the trench slip vectors, but only if Okhotsk-North America motion is slower than 5 mm/yr.

Description: This study is based on the NCAR King Air aircraft and radiosonde observations on 31 October 1986 during the FIRE in Wisconsin over Oshkosh. Aircraft step-up and spiral descent flights are used to obtain kinematic and thermodynamic data. In the step-up maneuver, six different penetrations were made between 1528 and 1616 UTC. Each penetration was about 30 km long separated in the vertical by about 300 m. The time difference between the two spiral soundings was about 43 min. The aircraft descended at a rate of 1.5 m/s during these spiral soundings. Kinematic, cloud physical, and radiometric observations from various instruments are used to estimate the different terms in the moisture- and heat-budget equations. The results show that the advection terms, estimated using the mean longitudinal wind and vertical velocities, and radiative fluxes are important in forming budgets for the cirrus layers. Ice-crystal growth is significant in the upper layers. The maintenance of cirrus can be attributed to relatively warm and moist air advection, radiative cooling at upper levels, and moisture advection in the vertical. Turbulent heat and moisture fluxes are found to be significant in the low levels of cirrus.

Description: Newell and Meng (1992) present maps of the occurrence probability of various classifications of particle precipitation as seen in the dayside topside ionosphere. It is argued that these are maps of the magnetospheric regions, a contention with which their critics disagree. The latter conclude that, because of convection, any one population of particles seen at low altitudes will have originated from a wide variety of locations, and particle characteristics cannot be mapped back to those in the magnetosphere without detailed knowledge of both the convection and magnetic field. Steplike boundaries between the regions will arise from nonsteady-state conditions and cannot be envisaged as steady-state magnetospheric boundaries between two plasma populations. In their reply Newell and Meng contend that convection does not move plasma from the LLBL into the cusp. Most of the LLBL plasma comes from the magnetosheath, so the direction of plasma transfer is in the other direction.

Description: The threshold method for estimation of area-average rain rate relies on determination of the fractional area where rain rate exceeds a preset level of intensity. Previous studies have shown that the optimal threshold level depends on the climatological rain-rate distribution (RRD). It has also been noted, however, that the climatological RRD may be composed of an aggregate of distributions, one for each of several distinctly different synoptic conditions, each having its own optimal threshold. In this study, the impact of RRD variations on the threshold method is shown in an analysis of 1-min rainrate data from a network of tipping-bucket gauges in Darwin, Australia. Data are analyzed for two distinct regimes: the premonsoon environment, having isolated intense thunderstorms, and the active monsoon rains, having organized convective cell clusters that generate large areas of stratiform rain. It is found that a threshold of 10 mm/h results in the same threshold coefficient for both regimes, suggesting an alternative definition of optimal threshold as that which is least sensitive to distribution variations. The observed behavior of the threshold coefficient is well simulated by assumption of lognormal distributions with different scale parameters and same shape parameters.

Description: To monitor future global temperature trends, it would be extremely useful if parameters nonlinearly related to surface temperature could be found, thereby amplifying any warming signal that may exist. Evidence that global thunderstorm activity is nonlinearly related to diurnal, seasonal and interannual temperature variations is presented. Since global thunderstorm activity is also well correlated with the earth's ionospheric potential, it appears that variations of ionospheric potential, that can be measured at a single location, may be able to supply valuable information regarding global surface temperature fluctuations. The observations presented enable a prediction that a 1 percent increase in global surface temperatures may result in a 20 percent increase in ionospheric potential.

Description: The largest known effusive eruptions during the Cenozoic and Mesozoic Eras, the voluminous flood basalts, have long been suspected as being associated with major extinctions of biotic species. Despite the possible errors attached to the dates in both time series of events, the significance level of the suspected correlation is found here to be 1 percent to 4 percent. Statistically, extinctions lag eruptions by a mean time interval that is indistinguishable from zero, being much less than the average residual derived from the correlation analysis. Oceanic flood basalts, however, must have had a different biological impact, which is still uncertain owing to the small number of known examples and differing physical factors. Although not all continental flood basalts can have produced major extinction events, the noncorrelating eruptions may have led to smaller marine extinction events that terminated at least some of the less catastrophically ending geologic stages. Consequently, the 26 Myr quasi-periodicity seen in major marine extinctions may be only a sampling effect, rather than a manifestation of underlying periodicity.

Description: The vicinity of KSC, where the primary ground truth site of the Tropical Rainfall Measuring Mission (TRMM) program is located, was the focal point of the Convection and Precipitation/Electrification (CaPE) experiment in July and Aug. 1991. In addition to several specialized radars, local coverage was provided by the C-band (5 cm) radar at Patrick AFB. Point measurements of rain rate were provided by tipping bucket rain gage networks. Besides these ground-based activities, airborne radar measurements with X- and Ka-band nadir-looking radars on board an aircraft were also recorded. A unique combination data set of airborne radar observations with ground-based observations was obtained in the summer convective rain regime of central Florida. We present a comparison of these data intending a preliminary validation. A convective rain event was observed simultaneously by all three instrument types on the evening of 27 July 1991. The high resolution aircraft radar was flown over convective cells with tops exceeding 10 km and observed reflectivities of 40 to 50 dBZ at 4 to 5 km altitude, while the low resolution surface radar observed 35 to 55 dBZ echoes and a rain gage indicated maximum surface rain rates exceeding 100 mm/hr. The height profile of reflectivity measured with the airborne radar show an attenuation of 6.5 dB/km (two way) for X-band, corresponding to a rainfall rate of 95 mm/hr.

Description: The Space Shuttle provides a unique perspective for cloud observation. From it, an observer can simulataneously view the grand scale and individual microscale components. An observer can differentiate between the exceptional and the routine and rapidly identify the phenomenon of interest. The capabilities include multiple aspect viewing, vision enhancement, instant data transmission, and the intelligence and decision-making ability required for meaningful selective observations. When these observations are used with other information, an improved understanding of the processes follows. Examples of these data and their application with other information sources are given along with suggestions for research.

Description: Flash rates measured from the Space Shuttle range from 27.8 flashes per minute to 77 flashes per minute. The cloud is an optically thick medium which effectively scatters the energy from a lightning discharge and thereby broadens the risetime and duration of each lightning pulse. Because of the small size, spacecraft sensors with resolutions of 1 km or more are unlikely to detect the individual lightning channels. Instead, the energy from the lightning channel is scattered within the cloud, thereby broadening the apparent area. All of these measurements of lightning flash area and flash rate have involved manual manipulation and analysis of the video or film data. Only a small percentage of the Space Shuttle lightning video has been analyzed. An attempt is made to combine the use of real-time digital disk system and an automated analysis routine in order to overcome this limitation and make processing of a sequence of video frames a much less labor-intensive task.

Description: A millimeter-wave imaging radiometer (MIR) developed by NASA Goddard Space Flight Center is described. The MIR is a nine-channel total power radiometer developed for atmospheric research. Three dual-pass band channels are centered about the strongly opaque 183-GHz water vapor absorption line; the frequencies are 183 +/- 1, +/- 3, and +/- 7 GHz. Another channel is located on the wing of this band at 150 GHz. These four channels have varying degrees of opacity from which the water vapor profile can be inferred. The design and salient characteristics of this instrument are discussed, together with its expected benefits.

Description: A new generation of automatic sunphotometers is used to systematically monitor clear sky total column aerosol concentration and optical properties, precipitable water and total column ozone diurnally and annually in West Africa and South America. The instruments are designed to measure direct beam sun, solar aureole and sky radiances in nine narrow spectral bands from the UV to the near infrared on an hourly basis. The instrumentation and the algorithms required to reduce the data for subsequent analysis are described.

Description: The monthly and daily samples of the Ap index for the interval from 1932 through 1982 were studied using the power spectrum technique. Results obtained for Bartel's period (about 27 days), the semiannual period, the dual-peak solar cycle distribution of geomagnetic storms, and certain other medium-scale periodicities are examined in detail. In addition, results on the cumulative occurrence number of storms per decade as a function of the Ap and Dst indices for the storm are presented.

Description: The circulation of polar wind ions from the high-latitude ionosphere to the plasma sheet is investigated. Single-particle trajectory simulations for the geomagnetic tail show, in addition to the breaking of adiabaticity, the possible creation of new high-altitude mirror points. This trajectory feature results from an intense parallel deceleration imparted by the magnetic field rotation during fast ExB transport. This centrifugal deceleration yields a critical parallel velocity which depends on the magnitude of the convection electric field and below which ions remain trapped inside the neutral sheet.

Description: High resolution in situ measurements of reactive nitrogen (NO(y)) and O3 were made in the upper troposphere and lower stratosphere at a variety of latitudes and seasons. In the lower stratosphere, NO(y) and O3 are very highly positively correlated at all times and spatial scales sampled. The ratio NO(y)/NO3 is much less variable than either species measured separately. The ratio has a much weaker gradient with altitude than the mixing ratios of O3 or NO(y). The ratio is smaller and decreases more rapidly with altitude in the tropics than at midlatitudes. In the upper troposphere NO(y) and O3 are only weakly correlated. Their ratio in the tropical upper troposphere is about 0.005-0.025 and the ratio in the midlatitude upper troposphere is about 0.004-0.010. The NO(y) in the upper troposhere is probably partly due to lightning.

Description: The chemistry and dynamics of a convective system observed on February 2, 1987 in the EMEX and STEP campaigns are analyzed. Chemical and thermodynamic profiles in undisturbed air near the EMEX 9 system indicate that the troposphere was well mixed by previous convection. As a consequence there was little direct transport from the boundary layer to the upper troposphere and air transported upward was detrained throughout the middle and upper troposphere. There was minimal effect on O3 production. Other convective complexes located 800-900 km upstream produced greater perturbations on trace gas profile immediately below the tropopause. Ozone production was reduced by about 0.25 ppbv/d at these altitudes, representing a reduction in P(O3) of 15-20 percent over the column from 14.5 to 17 km. P(O3) from 12 to 17 km in a region distant from active convection and subject to lightning was 2-3 times higher than it would have been without lightning.

Description: During the tropical experiment of NASA's Stratosphere-Troposphere Exchange Program (STEP), in situ radon and other trace constituent measurements were made aboard a NASA ER-2 high-altitude research aircraft to investigate the mechanisms of irreversible transfers from the troposphere into the tropical stratosphere. Observations made in and downwind of the cirrus shields of three large tropical cyclones and downwind of the cirrus anvil of a large cumulonimbus cloud cluster showed several clear instances of elevated radon activity occurring simultaneously with low total water mixing ratios. These observations are unambiguous evidence of an effective dehydration process, capable of reducing total water vapor mixing ratios to less than 2.5 ppmv, occurring in conjunction with troposphere-to-stratosphere transport and indicate that rapid localized convection, rather than slow regional mean motions, was responsible for the observed transports and associated with the accompanying dehydration. Radon activities measured in regions of active or recent troposphere-to-stratosphere transport were consistent with the 17 pCi/scm mean value needed to support the observed abundance of stratospheric 210 Pb.

Description: Measurements of stratospheric and upper tropospheric cloud water plus water vapor (total water) and water vapor were made with two Lyman alpha hygrometers as part of the STEP tropical experiment. The in situ measurements were made in the Darwin, Australia, area in January and February of 1987 on an ER-2 aircraft. Average stratospheric water vapor at a potential temperature of 375 K (the average value of Theta at the tropopause) was 2.4 parts per million by volume (ppmv). This water mixing ratio is below the 3.0 to 4.0 ppmv necessary to be consistent with the observed upper stratospheric dryness. Saturation with respect to ice and the potential for dehydration was observed up to Theta = 402 K.

Description: The STEP tropical objectives were successfully met during the flight experiments conducted from Darwin, Australia, January 16 to February 16, 1987. Necessary and sufficient measurements were made in, above, and downwind from very cold cirrus clouds, produced by three convective cloud types, to demonstrate irreversible mass transports into and dehydration in the lower tropical stratosphere. The three types are defined and described in terms of the physical processes that produce them and illustrated by examples derived from in situ and remote measurements. Intense solar heating is shown to produce, in addition to the usual vertical, sea breeze circulations normal to the coastline, an unusual pair of continental spanning, horizontal circulations. An upper tropospheric-lower stratospheric anticyclonic circulation, inclined upward toward the tropics, contributes to the dehydration of dissipating cirrus anvils and intensifies the upper level, tropical easterlies. The lower tropospheric cyclonic circulation with tropical westerlies and extratropical easterlies is in direct conflict with the normal tropical easterlies and extratropical westerlies. Impulsive switches between these two opposing lower-level wind systems create conditions favorable for each of these cloud types and explain the summer season's aperiodic variability.

Description: Overflights of a tropical cyclone during the Australian winter monsoon field experiment of the Stratosphere-Troposphere Exchange Project (STEP) show the presence of two mesoscale phenomena: a vertically propagating gravity wave with a horizontal wavelength of about 110 km and a feature with a horizontal scale comparable to that of the cyclone's entire cloud shield. The larger feature is fairly steady, though its physical interpretation is ambiguous. The 110-km gravity wave is transient, having maximum amplitude early in the flight and decreasing in amplitude thereafter. Its scale is comparable to that of 100-to 150-km-diameter cells of low satellite brightness temperatures within the overall cyclone cloud shield; these cells have lifetimes of 4.5 to 6 hrs. These cells correspond to regions of enhanced convection, higher cloud altitude, and upwardly displaced potential temperature surfaces. The temporal and spatial distribution of meteorological variables associated with the 110-km gravity wave can be simulated by a slowly moving transient forcing at the anvil top having an amplitude of 400-600 m, a lifetime of 4.5-6 hrs, and a size comparable to the cells of low brightness temperature.

Description: An overview is presented of the tropical component of STEP. The STEP cooperative experiments are described and summaries are presented of the STEP tropical ER-2 aircraft flights. STEP tropical results on dehydration and transfer and the mechanisms of upward transfer are summarized. Illustrations show flight paths for each sortie on satellite images and on 100 hPa synoptic flow charts, as well as the timing of flights with respect to overall cloudiness in the Australian region.

Description: The relationship between deep convection and tropopause cold trap conditions is examined for the tropical northern Australia region during the 1986-87 summer monsoon season, emphasizing the Australia Monsoon Experiment (AMEX) period when the NASA Stratosphere-Troposphere Exchange Project (STEP) was being conducted. The factors related to the spatial and temporal variability of the cold point potential temperature (CPPT) are investigated. A framework is developed for describing the relationships among surface average equivalent potential temperature in the surface layer (AEPTSL) the height of deep convection, and stratosphere-troposphere exchange. The time-mean pattern of convection, large-scale circulation, and surface AEPTSL in the Australian monsoon and the evolution of the convective environment during the monsoon period and the extended transition season which preceded it are described. The time-mean fields of cold point level variables are examined and the statistical relationships between mean CPPT, surface AEPTSL, and deep convection are described. Day-to-day variations of CPPT are examined in terms of these time mean relationships.

Description: Lateral variations of the temperature field in the lower mantle have been reconstructed using new results in mineral physics and seismic tomographic data. We show that, with the application of high-pressure experimental values of thermal expansivity and of sound velocities, the slow seismic anomalies in the lower mantle under the Pacific and Africa can be converted into realistic-looking plume structures with large dimensions of 0(1000 km). The outer fringes of the plumes have an excess temperature of around 400 K. In the core of the plumes are found tonguelike structures with extremely high thermal anomalies. These values can exceed 1200 K and are too high to be explained on the basis of thermal anomalies alone. We suggest that these major plumes in the deep mantle may be driven by both thermal and chemical buoyancies or that enhanced conductive heat-transfer may be important there.

Description: The largest known Cenozoic impact craters with the most accurately measured ages are found to correlate very closely with geologic stage boundaries. The level of confidence in this result is 98-99 percent even under the most pessimistic assumptions concerning dating errors. One or more large impacts may have led, in at least some cases, to the extinctions and first appearances of biotic species that mark many of the geologic stage boundaries.

Description: The study presents a methodology developed to apply a parameterization of radiative transfer calculations to satellite analyses of cirrus clouds. Cloud heights and optical depths are derived from visible and IR window measurements taken during FIRE when cirrus clouds were present. Geostationary satellite retrievals are compared to lidar-derived cloud heights and retrievals from a polar-orbiting satellite taken at different angles to determine which theoretical models of scattering phase function and single-scattering albedo best represent actual cirrus clouds. Models using small scattering phase function and single-scattering albedo best represent actual cirrus clouds. It is concluded that interpretation of cirrus reflectance with water-droplet models leads to biased results. The cloud-height and optical depth biases can be minimized with the aid of the C20 or cirrostratus models.

Description: Outgoing longwave radiation (OLR) data from Nimbus-7 ERB wide field-of-view instruments are compared with results from the ERBE instruments aboard the NOAA-9 and NOAA-10 satellites. Over most regions of the globe, the agreement between the two sets of OLR results is generally to within 8 W/sq m. There are larger differences at higher latitudes and regions concentrated over land and desert. Results of daytime and nighttime differences suggest that the shortwave channels may be at fault due to their different design for Nimbus-7 and NOAA-9. Some of the differences may also be related to different viewing geometry of the two satellites.

Description: Two simple techniques (Continuum Interpolated Band Ratio, CIBR, and Narrow/Wide ratio, N/W) to retrieve path precipitable water from AVIRIS high spectral resolution radiance data using the 940 nm water absorption band are compared. A sensitivity analysis was performed using the radiative transfer code LOWTRAN 7 to determine which one of these two approaches will provide a better estimate over land and water areas. The CIBR proved to be the technique less sensitive to perturbing effects, except for errors in visibility estimate. Both techniques were applied to AVIRIS radiance data acquired over Salton Sea, California. Resulting images confirmed that the use of a constant gray reflectance in the model led to a higher overestimation of the amount of water retrieved for N/W over vegetated areas. Validation was performed through comparison between an independent estimate of water vapor from concurrent Reagan sunphotometer measurements and AVIRIS estimates. Amounts retrieved using the N/W approach match more closely in situ measurements, even after adjusting model parameters for background reflectance, viewing geometry, and type of aerosol at the site.

Description: The viscous interfacial-sublayer model of Liu et al. (1979) is used to derive four bulk schemes (LKB, FG, D, and DB), with the flux-profile relationships of Lie et al., Francey and Garratt (1981), Dyer (1974), and Dyer and Bradley (1982). These schemes, with stability-dependent transfer coefficients, are tested against the eddy-correlation fluxes measured at the 50-m flight level above the western Atlantic Ocean during cold-air outbreaks. The bulk fluxes of momentum (tau), sensible heat (H), and latent heat (E) are found to increase with various von Karman constants. The dependence of transfer coefficients on wind speeds and roughness lengths is discussed. The transfer coefficients for tau and E agree excellently between LKB and FG. The ratio of the coefficent for H of LKB to that of FG, increasing with decreasing stability, is very sensitive to stability at low winds, but approaches the neutral value of 1.25 at high winds.

Description: Ozone profiles obtained with the Differential Absorption Lidar (DIAL) system at the National Institute for Environmental Studies (NIES) (Tsukuba, Japan) were compared with data provided by the satellite sensor SAGE II. The SAGE II data were selected based on criteria of spatial and temporal differences between the DIAL and the SAGE II measurements: five degrees in latitude and 15 degrees in longitude, within a latitudinal band from 31 deg to 41 deg N, and within one, three and five days after or before the DIAL measurements. Results show very good agreement for the individual and the zonal-mean profiles. The average mean difference between the DIAL and the SAGE II measurements over the altitudes 15-50 km was about 10 percent.