ALBERT

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

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

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

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

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

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

Description: Solar variability is examined in search of implications for global change. The topics covered include the following: solar variation modification of global surface temperature; the significance of solar variability with respect to future climate change; and methods of reducing the uncertainty of the potential amplitude of solar variability on longer time scales.

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

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

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

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

Description: We determined the infrared optical constants of nitric acid trihydrate, nitric acid dihydrate, nitric acid monohydrate, and solid amorphous nitric acid solutions which crystallize to form these hydrates. We have also found the infrared optical constants of H2O ice. We measured the transmission of infrared light throught thin films of varying thickness over the frequency range from about 7000 to 500/cm at temperatures below 200 K. We developed a theory for the transmission of light through a substrate that has thin films on both sides. We used an iterative Kramers-Kronig technique to determine the optical constants which gave the best match between measured transmission spectra and those calculated for a variety of films of different thickness. These optical constants should be useful for calculations of the infrared spectrum of polar stratospheric clouds.

Description: Meteorological conditions, extremely conducive to fire development and spread in the spring of 1987, resulted in forest fires burning over extremely large areas in the boreal forest zone in northeastern China and the southeastern region of Siberia. The great China fire, one of the largest and most destructive forest fires in recent history, occurred during this period in the Heilongjiang Province of China. Satellite imagery is used to examine the development and areal distribution of 1987 forest fires in this region. Overall trace gas emissions to the atmosphere from these fires are determined using a satellite-derived estimate of area burned in combination with fuel consumption figures and carbon emission ratios for boreal forest fires.

Description: In November to December 1991, a substantial number of remote sensors and in situ instruments were operated together in Coffeyville, Kansas, during the climate experiment First ISCCP Regional Experiment Phase 2 (FIRE 2). Includede in the suite of instruments were (1) the NOAA Environmental Technology Laboratory (ETL) three-channel microwave radiometer, (2) the NASA GSFC Raman lidar, (3) ETL radio acoustic sounding system (RASS), and (4) frequent, research-quality radiosondes. The Raman lidar operated only at night and the focus of this portion of the experiment concentrated on clear conditions. The lidar data, together with frequent radiosondes and measurements of temperature profiles (every 15 min) by RASS allowed profiles of temperature and absolute humidity to be estimated every minute. We compared 20 min measurements of brightness temperature (T(sub b) with calculations of T(sub b) that were based on the Liebe and Layton (1987) and Liebe et al. (1993) microwave propagation models, as well as the Waters (1976) model. The comparisons showed the best agreement at 20.6 GHz with the Waters model, with the Liebe et al. (1993) model being best at 31.65 GHz. The results at 90 GHz gave about equal success with the Liebe and Layton (1987) and Liebe et al. (1993) models. Comparisons of precipitable water vapor derived independently from the two instruments also showed excellent agreement, even for averages as short as 2 min. The rms difference between Raman and radiometric determinations of precipitable water vapor was 0.03 cm which is roughly 2%. The experiments clearly demonstrate the potential of simultaneous operation of radiometers and Raman lidars for fundamental physical studies of water vapor.

Description: This paper describes the validation of the International Satellite Cloud Climatology project (ISCCP)cloud detections by verifying the accuracy of the infrared clear-sky radiances. Comparison of retrieved surface temperatures to other measurements shows that bias errors are less than 2 K and random errors are about 2 K for sea surface (monthly means at 280-km scales) and that bias errors are less than 2 K and random errors are about 4 K for land surfaces (3 hourly at 280-km scales). Bias errors over a few persistently cloudy locations are sometimes -(2-4) K and over winter sea ice may be about +2 K. Surface reflectances are confirmed to be within 3% of other measurements and models for ocean, except for sun glint geometries, and to be within 3%-5% for land surfaces. Sufficiently accurate validation data are not avilable for visible reflectances of sea ice and snow-covered land, but some tests of specific cases suggest that errors are approximately 10%. These errors in clear-sky radiances suggest uncertainties in the ISCCP cloud detections of about 10% with a small (3%-6%) negative bias over land. Some specific regions exhibit both larger rms uncertainties and somewhat larger biases in cloud amount approaching 10%. ISCCP cloud detections are more in error over the polar regions than anywhere else. Based on comparisons with an anlysis of radiances measured at other wavelengths, the ISCCP analysis appears to miss 15%-25% of the clouds in summer but only 5%-10% of the winter clouds.

Description: Water vapor concentrations obtained by the Stratospheric Aerosol and Gas Experiment 2 (SAGE 2) and collocated temperatures provided by the National Meteorological Center (NMC) from 1986 to 1990 are used to deduce seasonally and zonally averaged acidity, density, and refractive index of stratospheric aerosols. It is found that the weight percentage of sulfuric acid in the aerosols increases from about 60 just above the tropopause to about 86 at 35 km. The density increases from about 1.55 to 1.85 g/cu cm between the same altitude limits. Some seasonal variations of composition and density are evident at high latitudes. The refractive indices at 1.02, 0.694, and 0.532 micrometers increase, respectively, from about 1.425, 1.430, and 1.435 just above the tropopause to about 1.445, 1.455, and 1.458 at altitudes above 27 km, depending on the season and latitude. The aerosol properties presented can be used in models to study the effectiveness of heterogeneous chemistry, the mass loading of stratospheric aerosols, and the extinction and backscatter of aerosols at different wavelengths. Computed aerosol surface areas, rate coefficients for the heterogeneous reaction ClONO2 + H2O yields HOCl + HNO3 and aerosol mass concentrations before and after the Pinatubo eruption in June 1991 are shown as sample applications.

Description: The global distribution of ozone during the October 31, 1978, to May 6, 1993, observing lifetime of the Nimbus 7 total ozone mapping spectrometer (TOMS) is described, with emphasis on the low ozone amounts observed during 1992 and 1993. Ozone amount time series are extended beyong May 6, 1993, to the end of July 1993 using preliminary Meteor 3 TOMS data. Time series for zonally averaged ozone amounts show that there has not been a significant shift in the seasonal patterns of ozone maxima and minima caused by the Mount Pinatubo eruption or by the onset of very low ozone values during 1992 and 1993. There has been a relatively slow, nearly linear decrease in the amount of ozone over the entire globe from 1979 to the end of 1991, with rates ranging from no change at the equator to a 4 - 6% decrease per decade at midlatitudes and a 10 - 12% decrease per decade at higher latitudes. After the eruption of Mount Pinatubo during June 1991, the ozone amount decreased in the equatorial latitudes (10 deg S to 10 deg N) for about 6 months (-10 Dobson units (DU) between 0 deg and 10 deg S and -3 DU between 0 deg and 10 deg N). During 1992 and continuing into 1993, the rate of ozone decrease deviated from the previously linear trend with the onset of changes that were large in comparison with the historical range of ozone values from 1979 to 1991. The first of the large decreases in ozone amount occurred earlier, in February 1990 to May 1990, at 50 deg - 70 deg N. At high northern latitudes, the 1993 decreased ozone amounts were about 12.5% below the envelope of historical values; at midlatitudes they were about 7% lower; and at low latitudes they were about 4% lower. Area-weighted averages in the northern and southern hemispheres show that most of the 1992 - 1993 ozone losses have occurred in the northern hemisphere. The 1993 global average (70 deg S to 70 deg N) ozone amount is 3% below the 1979 to 1991 minimum, 5% below the historical envelope in the northern hemisphere, and near the lower boundary of the historical envelope in the southern hemisphere. In the 70 deg - 60 deg S latitude band, the ozone losses between 1979 and 1993 have reduced the annual minimum amount to values below those seen in the equatorial regions.

Description: A new 8-year global cloud climatology has been produced by the International Satellite Cloud Climatology Project (ISCCP) that provides information every 3 h at 280-km spatial resolution covering the period from July 1983 through June 1991. If cloud detection errors and differences in area sampling are neglected, individual ISCCP cloud amounts agree with individual surface observations to within 15% rms with biases of only a few percent. When measurements of small-scale, broken clouds are isolated in the comparison, the rms differences between satellite and surface cloud amounts are about 25%, similar to the rms difference between ISCCP and Landsat determinations of cloud amount. For broken clouds, the average ISCCP cloud amounts are about 5% smaller than estimated by surface observers (difference between earth cover and sky cover), but about 5% larger than estimated from very high spatial resolution satellite observations (overestimate due to low spatial resolution offset by underestimate due to finite radiance thresholds). Detection errors caused by errors in the ISCCP clear-sky radiances or incorrect radiance threshold magnitudes, are the dominant source of error in monthly average cloud amounts. The ISCCP cloud amounts appear to be too low over land by about 10%, somewhat less in summer and somewhat more in winter, and about right (maybe slightly low) over oceans. In polar regions, ISCCP cloud amounts are probably too low by about 15%-25% in summer and 5%-10% in winter. Comparison of the ISCCP climatology to three other cloud climatologies shows excellent agreement in the geographic distribution and seasonal variation of cloud amounts; there is little agreement about day/night contrasts in cloud amount. Notable results from ISCCP are that the gobal annual mean cloud amount is about 63%, being about 23% higher over oceans than over land, that it varies by less than 1% rms from month to month, and that it has varied by about 4% on a time scale approximately equal to 2-4 years. The magnitude of interannual variations of local (280-km scale) monthly mean cloud amounts is about 7%-9%. Longitudinal contrasts in cloud amount are just as large as latitudinal contrasts. The largest seasonal variation of cloud amount occurs in the tropics, being larger in summer than in winter; the seasonal variation in the middle latitudes has the opposite phase. Polar regions may have little seasonable variability in cloud amount. The ISCCP results show slightly more nightime than daytime cloud amount over oceans and more daytime than nightime cloud amount over land.

Description: This paper, the first of three, describes the cloud detection part of the International Satellite Cloud Climatology Project (ISCCP) analysis. Key features of the cloud detection alogrithm are (1) use of space and time radiance variation tests over several different space and time domains to account for the global variety of cloudy and clear characteristics, (2) estimation of clear radiance values for every time and place, and, (3) use of radiance thresholds that vary with the type of surface and climate regime. Design of the detection algorithm was supported by global, multiyear surveys of the statistical behavior of satellite-measured infrared and visible radiances to determine those characteristics that differentiate cloudy and clear scenes and how these characteristics vary among climate regimes. A summary of these statistical results is presented to illustrate how the cloud detection method works in a variety of circumstances. The sensitivity of the results to changing test parameter values is determined to provide a first estimate of the uncertainty of ISCCP cloud amounts. These test results (which exclude polar regions) suggest detection uncertainties of about 10% with possible negative biases of 5% (especially at night).

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

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

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

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

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

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

Description: The Chicxulub impact crater in Mexico is the site of the impact purported to have caused mass extinctions at the Cretaceous/Tertiary (K/T) boundary. 2-D hydrocode modeling of the impact, coupled with studies of the impact site geology, indiate that between 0.4 and 7.0 x 10(exp 17) g of sulfur were vaporized by the impact into anhydrite target rocks. A small portion of the sulfur was released as SO3 or SO4, which converted rapidly into H2SO4 aerosol and fell as acid rain. A radiative transfer model, coupled with a model of coagulation indicates that the aerosol prolonged the initial blackout period caused by impact dust only if the aerosol contained impurities. A larger portion of sulfur was released as SO2, which converted to aerosol slowly, due to the rate-limiting oxidation of SO2. Our radiative transfer calculations, combined with rates of acid production, coagulation, and diffusion indicate that solar transmission was reduced to 10-20% of normal for a period of 8-13 yr. This reduction produced a climate forcing (cooling) of -300 W/sq.m, which far exceeded the +8 W/sq.m greenhouse warming, caused by the CO2 released through the vaporization of carbonates, and therefore produced a decade of freezing and near-freezing temperatures. Several decades of moderate warming followed the decade of severe cooling due to the long residence time of CO2. The prolonged impact winter may have been a major cause of the K/T extinctions.

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

Description: Clear-sky albedos and outgoing longwave radiation (OLR) determined from Earth Radiation Budget Experiment (ERBE) scanners on board the earth radiation budget satellite and NOAA-9 spacecraft were analyzed for three target sites for the months February 1985-January 1987. The targets were oceans, deserts, and a multiscene site covering half the earth's surface. Year-to-year ratios of the monthly albedos and OLR were within the 0.98-1.02 range with a standard error of about 1%. The data indicate that ERBE scanner measurements were stable to within a few tenths of a percent for the two-year periods.

Description: A set of system simulations has been performed to evaluate candidate scanner designs for an Earth Radiation Budget Instrument (ERBI) for the Earth Observing System (EOS) of the late 1990s. Five different instruments are considered: (1) the Active Cavity Array (ACA), (2) the Clouds and Earth's Radiant Energy System-Instrument (CERES-I), (3) the Conically Scanning Radiometer (CSR), (4) the Earth Radiation Budget Experiment Cross-Track Scanner (ERBE), and (5) the Nimbus-7 Biaxial Scanner (N7). Errors in instantaneous, top-of-the-atmosphere (TOA) satellite flux estimates are assumed to arise from two measurement problems: the sampling of space over a given geographic domain, and sampling in angle about a given spatial location. When angular sampling errors vanish due to the application of correct angular dependence models (ADMs) during inversion, the accuracy of each scanner design is determined by the instrument's ability to map the TOA radiance field in a uniform manner. In this regard, the instruments containing a cross-track scanning component (CERES-I and ERBE) do best. As errors in ADMs are encountered, cross-track instruments incur angular sampling errors more rapidly than biaxial instruments (N7, ACA, and CSR) and eventually overtake the biaxial designs in their total error amounts. A latitude bias (north-south error gradient) in the ADM error of cross-track instruments also exists. This would be objectionable when ADM errors are systematic over large areas of the globe. For instantaneous errors, however, cross-track scanners outperform biaxial or conical scanners for 2.5 deg latitude x 2.5 deg longitude target areas, providing that the ADM error is less than or equal to 30%. A key issue is the amount of systematic ADM error (departures from the mean models) that is present at the 2.5 deg resolution of the ERBE target areas. If this error is less than 30%, then the CERES-I, ERBE, and CSR, in order of increasing error, provide the most accurate instantaneous flux estimates, within 2-3 W/sq m of each other in reflected shortwave flux. The magnitude of this error is near the 10 W/sq m accuracy requirement of the user community. Longwave flux errors have been found to have the same space and time characteristics as errors in shortwave radiation, but only about 25% as large.

Description: This paper briefly reviews several single-frequency rain profiling methods for an airborne or spaceborne radar. The authors describe the different methods from a unified point of view starting from the basic differential equation. This facilitates the comparisons between the methods and also provides a better understanding of the physical and mathematical basis of the methods. The application of several methods to airborne radar data taken during the Convective and Precipitation/Electrification Experiment is shown. Finally, the authors consider a hybrid method that provides a smooth transition between the Hitschfeld-Bordan method, which performs well at low attenuations, and the surface reference method, for which the relative error decreases with increasing path attenuation.

Description: An analysis of a spontaneous sudden stratospheric warming that occurred during a 2-year integration of the Langley Research Center (LaRC) Atmospheric Simulation Model is presented. The simulated warming resembles observed 'wave 1' warmings in the Northern Hemisphere stratosphere and provides an opportunity to investigate the radiative and dynamical processes occurring during the warming event. Isentropic analysis of potential vorticity sources and sinks indicates that dynamically induced departures from radiative equilibrium play an important role in the warming event. Enhanced radiative cooling associated with a series of upper stratospheric warm pools leads to radiative dampening within the polar vortex. Within the 'surf zone' large-scale radiative cooling leads to diabatic advection of high potential vorticity air from aloft. Lagrangian area diagnostics of the simulated warming agree well with Limb Infrared Monitor of the Stratosphere (LIMS) analyses. Dynamical mixing is shown to account for the majority of the decrease in the size of the polar vortex during the simulated warming. An investigation of the nonlinear deformation of material lines that are initially coincident with diagnosed potential vorticity isopleths is conducted to clarify the relationship between the Lagrangian area diagnostics and potential vorticity advection during wave breaking events.

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

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

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

Description: Calculations of the quasi biennial oscillation (QBO) signal in Stratospheric Aerosol and Gas Experiment (SAGE) II O3 and NO2 data between 1984 and 1991 are presented and have been investigated by using a two-dimensional model. The isentropic 2D model is a fully interactive radiative-dynamical-chemical model in which the eddy fluxes of chemical species are calculated in a consistent manner. The QBO in the model has been forced by relaxing the equatorial zonal wind toward the observations at Singapore allowing the comparison of the model with observations from specific years. The model reproduces the observed vertical structure of the equatorial ozone anomaly with the well-known transition from dynamical to photochemical control at around 28km. The model also reproduces the observed vertical structure of the SAGE II observed NO2 anomaly. The model studies have shown that it is the QBO modulation of NO2 which the main cause of QBO signal in O3 above 30 km. The model also reproduces the observed latitudinal structure of the QBO signals in O3 and NO2. Due to the differing horizontal distribution of O3 and NO(y) the ozone signal shows a distinct phase change in the subtropics whereas the NO2 anomaly gives a broader signal.

Description: A generalized form of the second-order van Leer transport scheme is derived. Several constraints to the implied subgrid linear distribution are discussed. A very simple positive-definite scheme can be derived directly from the generalized form. A monotonic version of the scheme is applied to the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) for the moisture transport calculations, replacing the original fourth-order center-differencing scheme. Comparisons with the original scheme are made in idealized tests as well as in a summer climate simulation using the full GLA GCM. A distinct advantage of the monotonic transport scheme is its ability to transport sharp gradients without producing spurious oscillations and unphysical negative mixing ratio. Within the context of low-resolution climate simulations, the aforementioned characteristics are demonstrated to be very beneficial in regions where cumulus convection is active. The model-produced precipitation pattern using the new transport scheme is more coherently organized both in time and in space, and correlates better with observations. The side effect of the filling algorithm used in conjunction with the original scheme is also discussed, in the context of idealized tests. The major weakness of the proposed transport scheme with a local monotonic constraint is its substantial implicit diffusion at low resolution. Alternative constraints are discussed to counter this problem.

Description: A model is presented to compute the turbulent kinetic energy dissipation length scale l(sub epsilon) in a stably stratified shear flow. The expression for l(sub epsilon) is derived from solving the spectral balance equation for the turbulent kinetic energy. The buoyancy spectrum entering such equation is constructed using a Lagrangian timescale with modifications due to stratification. The final result for l(sub epsilon) is given in algebraic form as a function of the Froude number Fr and the flux Richardson number R(sub f), l(sub epsilon) = l(sub epsilon)(Fr, R(sub f). The model predicts that for R(sub f) less than R(sub fc), l(sub epsilon) decreases with stratification. An attractive feature of the present model is that it encompasses, as special cases, some seemingly different models for l(sub epsilon) that have been proposed in the past by Deardorff, Hunt et al., Weinstock, and Canuto and Minotti. An alternative form for the dissipation rate epsilon is also discussed that may be useful when one uses a prognostic equation for the heat flux. The present model is applicable to subgrid-scale models, which are needed in large eddy simulations (LES), as well as to ensemble average models. The model is applied to predict the variation of l(sub epsilon) with height z in the planetary boundary layer. The resulting l(sub epsilon) versus z profile reproduces very closely the nonmonotonic profile of l(sub epsilon) exhibited by many LES calculations, beginning with the one by Deardorff in 1974.

Description: Two stratospheric warmings during February and March 1993 are described using United Kingdom Meteorological Office (UKMO) analyses, calculated potential vorticity (PV) and diabetic heating, and N2O observed by the Cryogenic Limb Array Etalon Spectrometer (CLAES) instrument on the Upper Atmosphere Research Satellite (UARS). The first warming affected temperatures over a larger region, while the second produced a larger region of reversed zonal winds. Tilted baroclinic zones formed in the temperature field, and the polar vortex tilted westward with height. Narrow tongues of high PV and low N2O were drawn off the polar vortex, and irreversibly mixed. Tongues of material were drawn from low latitudes into the region between the polar vortex and the anticyclone; diabatic descent was also strongest in this region. Increased N2O over a broad region near the edge of the polar vortex indicates the importance of horizontal transport. N2O decreased in the vortex, consistent with enhanced diabatic descent during the warmings.

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

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

Description: In this paper, preliminary results in using orthogonal and continuous wavelet transform (WT) to identify period doubling and time-frequency localization in both synthetic and real data are presented. First, the Haar WT is applied to synthetic time series derived from a simple nonlinear dynamical system- a first-order quadratic difference equation. Second, the complex Morlet WT is used to study the time-frequency localization of tropical convection based on a high-resolution Japanese Geostationary Meteorological Satellite infrared (IR) radiance dataset. The Haar WT of the synthetic time series indicates the presence and distinct separation of multiple frequencies in a period-doubling sequence. The period-doubling process generates a multiplicity of intermediate frequencies, which are manifested in the nonuniformity in time with respect to the phase of oscillations in the lower frequencies. Wavelet transform also enables the detection of extremely weak signals in high-order subharmonics resulting from the period-doubling bifurcations. These signals are either undetected or considered statistically insignificant by traditional Fourier analysis. The Morlet WT of the IR radiance dataset indicates the presence of multiple timescales, which are localized in both frequency and time. There are two regimes in the variation of IR radiance, corresponding to the wet and dry periods. Multiple timescales, ranging from semidiurnal, diurnal, synoptic, to intraseasonal with embedding structures, are active in the wet regime. In particular, synoptic variability is more prominent during the wet phase of an intensive intraseasonal cycle. These are not only consistent with, but also show more details than, previous findings by using other techniques. The phase-locking relationships among the oscillations with different time-scales suggest that both synoptic and intraseasonal variations may be mixed oscillations due to the interaction of self-excited oscillations in the tropical atmosphere and external forcings such as annual and diurnal solar radiation variations. Both examples show that WT is a powerful tool for analysis of phenomena involving multiscale interactions that exhibit localization in both frequency and time. A discussion on the caveats in the use of WT in geophysical data analysis is also presented.

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

Description: Observations of upper tropospheric relative humidity obtained from Raman lidar and Cross-chain Loran Atmospheric Sounding System (CLASS) sonde instruments obtained during the First ISCCP Regional Experiment (FIRE) Cirrus-II field program are compared with satellite measurements from the GOES 6.7-micron channel. The 6.7-micron channel is sensitive to water vapor integrated over a broad layer in the upper troposphere (roughly 500-200 mbar). Instantaneous measurements of the upper tropospheric relative humidity from GOES are shown to agree to within roughly 6% of the nearest lidar observations and 9% of the nearest CLASS observations. The CLASS data exhibit a slight yet systematic dry bias in upper tropospheric humidity, a result which is consistent with previous radiosonde intercomparisons. Temporal stratification of the CLASS data indicates that the magnitude of the bias is dependent upon the time of day, suggesting a solar heating effect in the radiosonde sensor. Using CLASS profiles, the impact of vertical variability in relative humidity upon the GOES upper tropospheric humidity measurements is also examined. The upper tropospheric humidity inferred from the GOES 6.7-micron channel is demonstrated to agree to within roughly 5% of the relative humidity vertically averaged over the depth of atmosphere to which the 6.7-micron channel is sensitive. The results of this study encourage the use of satellite measurements in the 6.7-micron channel to quantitatively describe the distribution and temporal evolution of the upper tropospheric humidity field.

Description: Published estimates of cloud liquid water path (LWP) from satellite-measured microwave radiation show little agreement, even about the relative magnitudes of LWP in the tropics and midlatitudes. To understand these differences and to obtain more reliable estimate, optical and microwave LWP retrieval methods are compared using the International Satellite Cloud Climatology Project (ISCCP) and special sensor microwave/imager (SSM/I) data. Errors in microwave LWP retrieval associated with uncertainties in surface, atmosphere, and cloud properties are assessed. Sea surface temperature may not produce great LWP errors, if accurate contemporaneous measurements are used in the retrieval. An uncertainty of estimated near-surface wind speed as high as 2 m/s produces uncertainty in LWP of about 5 mg/sq cm. Cloud liquid water temperature has only a small effect on LWP retrievals (rms errors less than 2 mg/sq cm), if errors in the temperature are less than 5 C; however, such errors can produce spurious variations of LWP with latitude and season. Errors in atmospheric column water vapor (CWV) are strongly coupled with errors in LWP (for some retrieval methods) causing errors as large as 30 mg/sq cm. Because microwave radiation is much less sensitive to clouds with small LWP (less than 7 mg/sq cm) than visible wavelength radiation, the microwave results are very sensitive to the process used to separate clear and cloudy conditions. Different cloud detection sensitivities in different microwave retrieval methods bias estimated LWP values. Comparing ISCCP and SSM/I LWPs, we find that the two estimated values are consistent in global, zonal, and regional means for warm, nonprecipitating clouds, which have average LWP values of about 5 mg/sq cm and occur much more frequently than precipitating clouds. Ice water path (IWP) can be roughly estimated from the differences between ISCCP total water path and SSM/I LWP for cold, nonprecipitating clouds. IWP in the winter hemisphere is about 3 times the LWP but only half the LWP in the summer hemisphere. Precipitating clouds contribute significantly to monthly, zonal mean LWP values determined from microwave, especially in the intertropical convergence zone (ITCZ), because they have almost 10 times the liquid water (cloud plus precipitation) of nonprecipitating clouds on average. There are significant differences among microwave LWP estimates associated with the treatment of precipitating clouds.

Description: Cloud radiative forcing (CRF) is the radiative impact of clouds on the Earth's radiation budget. This study examines the diurnal variations of CRF using the Earth Radiation Budget Experiment (ERBE) monthly hourly flux data and the flux data derived from the International Satellite Cloud Climatology Project (ISCCP) using the Goddard Institute for Space Studies general circulation model radiation code. The results for the months of April, July, and October 1985 and January 1986 are analyzed. We found that, in general, two data sets agreed. For longwave (LW) CRF the diurnal range over land is generally greater than that observed over oceans. For the 4-month averages the ERBE values are 15.8 W/sq m and 6.8 W/sq m for land and ocean, respectively, compared with the ISCCP calculated values of 18.4 W/sq m and 8.0 W/sq m, respectively. The land/ocean contrast is largely associated with changes in cloud amount and the temperature difference between surface and cloud top. It would be more important to note that the clear-sky flux (i.e., surface temperature) variabilities are shown to be a major contributor to the large variabilities over land. The maximum diurnal range is found to be in the summer hemisphere, and the minimum values in the winter hemisphere. It is also shown that the daytime maximum and the nighttime minimum are seen over large portions of land, whereas they occur at any local hour over most oceans. For shortwave (SW) CRF the daytime maximum values are about twice as large as monthly averages, and their highest frequency occurs at local noon, indicating that solar insolation is a primary factor for the diurnal variation of SW CRF. However, the comparison of the ERBE data with the ISCCP results demonstrated that the largest differences in the diurnal range and monthly mean of LW CRF were associated with tropical convergence zones, where clear-sky fluxes could be easily biased by persistent cloudiness and the inadequate treatment of the atmospheric water vapor.

Description: Temperature variations in the stratosphere from 1979 to 1992 are investigated using 365-day running mean of the National Meteorological Center gridded analysis temperature data. Significant variations are seen at all levels between 70 and 1 mbar. The middle stratosphere shows temperature peaks during 1982 and 1983. The upper stratosphere has significant temperature declines between 1 and 10 mbar from 1981 to 1984. Temperatures at all levels recover to near their prior values after 1984, with the 5-mbar temperatures requiring the greatest time to fully recover. The temperature declines at 1 mbar occur in both hemispheres, over all longitudes, and in every month of the year. The decreases are largest in the middle latitudes and the polar regions and during the fall and the winter months. Such temperature variations, which appear to be of natural origin, must be taken into consideration when searching for temperature trends caused by the increase of CO2 or other greenhouse gases which affect the radiative balance of the Earth-atmosphere system or stratospheric ozone.

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

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

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

Description: This paper presents a comparison of the water vapor distribution obtained from two general circulation models, the European Centre for Medium-Range Weather Forecasts (ECMWF) model and the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM), with satellite observations of total precipitable water (TPW) from Special Sensor Microwave/Imager (SSM/I) and upper tropospheric relative humidity (UTH) from GOES. Overall, both models are successful in capturing the primary features of the observed water vapor distribution and its seasonal variation. For the ECMWF model, however, a systematic moist bias in TPW is noted over well-known stratocumulus regions in the eastern subtropical oceans. Comparison with radiosonde profiles suggests that this problem is attributable to difficulties in modeling the shallowness of the boundary layer and large vertical water vapor gradients which characterize these regions. In comparison, the CCM is more successful in capturing the low values of TPW in the stratocumulus regions, although it tends to exhibit a dry bias over the eastern half of the subtropical oceans and a corresponding moist bias in the western half. The CCM also significantly overestimates the daily variability of the moisture fields in convective regions, suggesting a problem in simulating the temporal nature of moisture transport by deep convection. Comparison of the monthly mean UTH distribution indicates generally larger discrepancies than were noted for TPW owing to the greater influence of large-scale dynamical processes in determining the distribution of UTH. In particular, the ECMWF model exhibits a distinct dry bias along the Intertropical Convergence Zone (ITCZ) and a moist bias over the subtropical descending branches of the Hadley cell, suggesting an underprediction in the strength of the Hadley circulation. The CCM, on the other hand, demonstrates greater discrepancies in UTH than are observed for the ECMWF model, but none that are as clearly correlated with well-known features of the large-scale circulation.

Description: In this paper we describe fundamental properties of an 'off-line' three-dimensional transport model, that is, a model which uses prescribed rather than predicted winds. The model is currently used primarily for studying problems of the middle atmosphere because we have not (yet) incorporated a formulation for the convective transport of trace species, a prerequisite for many tropospheric problems. The off-line model is simpler and less expensive than a model which predicts the wind and mass evolution (an 'on-line' model), but it is more complex than the two-dimensional (2-D) zonally averaged transport models often used in the study of chemistry and transport in the middle atmosphere. It thus serves as a model of intermediate complexity and can fill a useful niche for the study of transport and chemistry. We compare simulations of four tracers, released in the lower stratosphere, in both the on- and off-line models to document the difference resulting from differences in modeling the same problem with this intermediate model. These differences identify the price to be paid in going to a cheaper and simpler calculation. The off-line model transports a tracer in three dimensions. For this reason, it requires fewer approximations than 2-D transport model, which must parameterize the effects of mixing by transient and zonally asymmetric wind features. We compare simulations of the off-line model with simulations of a 2-D model for two problems. First, we compare 2-D and three-dimensional (3-D) models by simulating the emission of an NO(x)-like tracer by a fleet of high-speed aircraft. The off-line model is then used to simulate the transport of C-14 and to contrast its simulation properties to that of the host of 2-D models which participated in an identical simulation in a recent NASA model intercomparison. The off-line model is shown to be somewhat sensitive to the sampling strategy for off-line winds. Simulations with daily averaged winds are in very good qualitative agreement but are less diffusive than when driven with instantaneous winds sampled at half-hour intervals. Simulations with the off-line and 2-D models are quite similar in the middle and upper stratosphere but behave quite differently in the lower stratosphere, where the 3-D model has a substantially more vigorous circulation. The off-line model is quite realistic in its simulation of C-14. While there are still systematic differences between the 3-D calculation and the observations, the differences seem to be substantially reduced when compared with the body of 2-D simulations documented in the above mentioned NASA intercomparison, particularly at 31 deg N.

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

Description: SEASAT synthetic aperture radar (SAR) echoes from the sea have previously been shown to be the result of rain and winds produced by convective stroms; rain damps the surface waves and causes ech-free holes, while the diverging winds associated with downdraft generate waves and associated echoes surrounding the holes. Gust fronts are also evident. Such a snapshot from 8 July 1978 has been examined in conjunction with ground-based radar. This leads to the conclusion that the SAR storm footprints resulted from storm processes that occurred up to an hour or more prior to the snapshot. A sequence of events is discerned from the SAR imagery in which new cell growth is triggered in between the converging outflows of two preexisting cells. In turn, the new cell generates a mini-squall line along its expanding gust front. While such phenomena are well known over land, the spaceborne SAR now allows important inferences to be made about the nature and frequency of convective storms over the oceans. The storm effects on the sea have significant implications for spaceborne wind scatterometry and rainfall measurements. Some of the findings herein remain speculative because of the great distance to the Miami weather radar-the only source of corroborative data.

Description: Time average climatology and low-frequency variabilities of the global hydrologic cycle (GHC) in the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) were investigated in the present work. A 730-day experiment was conducted with the GLA GCM forced by insolation, sea surface temperature, and ice-snow undergoing climatological annual cycles. Ifluences of interactive soil moisture on time average climatology and natural variability of the GHC were also investigated by conducting 365-day experiments with and without interactive soil moisture. Insolation, sea surface temperature, and ice-snow were fixed at their July levels in the latter two experiments. Results show that the model's time average hydrologic cycle variables for July in all three experiments agree reasonably well with observations. Except in the case of precipitable water, the zonal average climates of the annual cycle experiment and the two perpetual July experiments are alike, i.e., their differences are within limits of the natural variability of the model's climate. Statistics of various components of the GHC, i.e., water vapor, evaporation, and precipitation, are significantly affected by the presence of interactive soil moisture. A long-term trend is found in the principal empirical modes of variability of ground wetness, evaporation, and sensible heat. Dominant modes of variability of these quantities over land are physically consistent with one another and with land surface energy balance requirements. The dominant mode of precipitation variability is found to be closely related to organized convection over the tropical western Pacific Ocean. The precipitation variability has timescales in the range of 2 to 3 months and can be identified with the stationary component of the Madden-Julian Oscillation. The precipitation mode is not sensitive to the presence of interactive soil moisture but is closely linked to both the rotational and divergent components of atmospheric moisture transport. The present results indicate that globally coherent natural variability of the GHC in the GLA GCM has two basic timescales in the absence of annual cycles of external forcings: a long-term trend associated with atmosphere-soil moisture interaction which affects the model atmosphere mostly over midlatitude continental regions and a large-scale 2- to 3-month variability associated with atmospheric moist processes over the western Pacific Ocean.

Description: Climatologies of convective precipitation were derived from passive microwave observations from the Special Sensor Microwave Imager using a scattering-based algorithm of Adler et al. Data were aggregated over periods of 3-5 months using data from 4 to 5 years. Data were also stratified by satellite overpass times (primarily 06 00 and 18 00 local time). Four regions (Mexico, Amazonia, western Africa, and the western equatorial Pacific Ocean (TOGA COARE area) were chosen for their meteorological interest and relative paucity of conventional observations. The strong diurnal variation over Mexico and the southern United States was the most striking aspect of the climatologies. Pronounced morning maxima occured offshore, often in concativities in the coastline, the result of the increased convergence caused by the coastline shape. The major feature of the evening rain field was a linear-shaped maximum along the western slope of the Sierra Madre Occidental. Topography exerted a strong control on the rainfall in other areas, particularly near the Nicaragua/Honduras border and in Guatemala, where maxima in excess of 700 mm/month were located adjacent to local maxima in terrain. The correlation between the estimates and monthly gage data over the southern United States was low (0.45), due mainly to poor temporal sampling in any month and an inadequate sampling of the diurnal cycle. Over the Amazon Basin the differences in morning versus evening rainfall were complex, with an alternating series of morning/evening maxima aligned southwest to northeast from the Andes to the northeast Brazilian coast. A real extent of rainfall in Amazonia was slightly higher in the evening, but a maximum in morning precipitation was found on the Amazon River just east of Manaus. Precipitation over the water in the intertropical convergence zone (ITCZ) north of Brazil was more pronounced in the morning, and a pronounced land-/sea-breeze circulation was found along the northeast coast of Brazil. Inter-comparison of four years revealed 1992 to be the driest over Amazonia, with about a 23% decrease in mean rate compared to the 4-year mean estimated rain rate.

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

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

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

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

Description: Ample evidence supports the significance of the high-latitude ionospheric contribution to magnetospheric plasma. Assuming flux conservation along a flux tube, the upward field-aligned ion flows observed in the magnetosphere require high-latitude ionospheric field-aligned ion upflows of the order of 10(exp 8) to 10(exp 9)/sq cm/s. Since radar and satellite observations of high-latitude F region flows at times exceed this flux requirement by an order of magnitude, the thermal ionospheric upflows are not simply the ionospheric response to a magnetospheric flux requirement. Several ionospheric ion upflow mechanisms have been proposed, but simulations based on fluid theory do not reproduce all the observed features of ionospheric ion upflows. Certain asymmetries in the statistical morphology of high-latitude F region ion upflows suggest that the ion upflows may be generated by ion-neutral frictional heating. We developed a single-component (O(+)), time-dependent gyro-kinetic model of the high-latitude F region response to frictional heating in which the neutral exobase is a discontinuous boundary between fully collisional and collisionless plasmas. The concept of a discontinuous neutreal exobase and the assumption of a constant and uniform polarization electric field reduce the ion velocity distribution function, from which we can compute the ion density, parallel velocity, parallel and perpendicular temperature, and parallel flux. Using our model, we simulated the response of a convecting flux tube between 500 km and 2500 km to various frictional heating inputs; the results were both qualitatively and quantitatively different from fluid model results, which may indicate an inadequacy of the fluid theory approach. The gyro-kinetic frictional heating model responses to the various simulations were qualitatively similar: (1) initial perturbations of all the modeled parameters propagated rapidly up the flux tube, (2) transient values of the ion parallel velocity, temperature, and flux exceeded 3 km/s, 2 x 10(exp 4) K, and 10(exp 9)/sq cm/s, respectively, (3) a second transient regime developed wherein the parallel temperature drops to very low values (a few hundred Kelvins), and (4) well after heating ceased, large parallel temperatures and large downward parallel velocities and fluxes developed as the flux tube slowly returned to diffusive equilibrium. The ion velocity distributions during the simulation are often non-Maxwellian and are sometimes composed of two distinct ion populations.

Description: The meteorological satellite program began in the United States as the result of the actions taken by a very small but dedicated group of people from the late 1940s to 1960. This paper provides firsthand accounts by two of these dedicated individuals. Their remarks provide an insight into the trials and tribulations they and the program encountered during these very early years. Those now active in the program, many of whom do not recall this time, might appreciate the effort of these pioneers and the legacy they left for us.

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: 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: 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: 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: 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: 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.