ALBERT

Description: The development of an integrated approach to the modeling of forest dynamics encompassing submodels of forest growth and succession, soil processes and radiation interactions, is reported. Remote sensing technology is a key element of this study in that it provides data for developing, initializing, updating, and validating the models. The objectives are reviewed, the data collected and models in use are discussed, and a framework for studying interactions between the forest growth, soil process and energy interaction components, is described. Remote sensing technology used in the study includes optical and microwave field, aircraft and satellite borne instruments. The types of data collected during intensive field and aircraft campaigns included bidirectional reflectance, thermal emittance and multifrequency, multipolarization synthetic aperture radar backscatter. Synthetic imagery of derived products such as forest biomass and NDVI (Normalized Difference Vegetative Index), and collections of ground data are being assembled in a georeferenced data base. These data are used to drive or test multidiscipline simulations of forested ecosystems. Enhancements to the modeling environment permit considerable flexibility in configuring simulations and selecting results for reporting and graphical display.

Description: The quantitative interpretation of satellite observations requires the use of mathematical tools to extract the desired information on terrestrial environments from the radiation data collected in space. A whole range of approaches can be pursued, from the development of models capable of explaining the nature of the physical signal being measured and of characterizing the state of the system under observation, to the empirical correlations between the variables of interest and the space measurements. The premises and implications of these approaches are outlined, paying special attention to the mathematical and numerical requirements. The role and specific applications of empirical bidirectional reflectance models is also discussed, even though these models do not contribute to the understanding of the theory of radiation transfer or to the assessment of the variables of interest. The advantages and drawbacks of these various approaches and the research priorities for the next few years are discussed in the context of the planned availability of new sensors.

Description: Management of crop residues, the portion of a crop left in the field after harvest, is an important conservation practice for minimizing soil erosion and for improving water quality. Quantification of crop residue cover is required to evaluate the effectiveness of conservation tillage practices. Methods are needed to quantify residue cover that are rapid, accurate, and objective. The fluorescence of crop residue was found to be a broadband phenomenon with emission maxima at 420 to 495 nm for excitations of 350 to 420 nm. Soils had low intensity broadband emissions over the 400 to 690 nm region for excitations of 300 to 600 nm. The range of relative fluorescence intensities for the crop residues was much greater than the fluorescence observed of the soils. As the crop residues decompose their blue fluorescence values approach the fluorescence of the soil. Fluorescence techniques are concluded to be less ambiguous and better suited for discriminating crop residues and soils than reflectance methods. If properly implemented, fluorescence techniques can be used to quantify, not only crop residue cover, but also photosynthetic efficiency in the field.

Description: The importance of the measurement of wind fields is discussed. Wind regime data can be used to infer the amount and type of wind induced (aerolian) transport of sand and dust, or to establish global circulation models, for example on other planets. Since local measurements are costly and often impossible, it is desired to infer such data from remotely sensed information. A potential mechanism for remotely inferring the wind regime by using synthetic aperture radar data to describe the roughness of the surface is described. A project to estimate the practicality of using such a mechanism is described. An experiment that extends the mechanism to vegetated sites, where the goal is to measure potential for erosion, is reported.

Description: Surface reflectance is required to quantitatively investigate molecular absorption and particle scattering properties of materials on the Earth's surface. Atmospheric aerosol optical depth, surface pressure and water vapor are required to constrain a radiative transfer code for the inversion of measured spectral radiance to apparent surface reflectance. A suite of algorithms using nonlinear least squares fitting techniques are described that directly estimate these atmospheric parameters from spectral radiance measured by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). The derived atmospheric parameters are used to constrain a radiative transfer code for the inversion of the imaging spectrometer radiance to apparent reflectance. The derived apparent reflectance is validated with respect to in situ measurement on the same target.

Description: Investigations designed to study land surface hydrologic-atmospheric interactions, showing the potential of L band passive microwave radiometry for measuring surface soil moisture over large areas, are discussed. Satisfying the data needs of these investigations requires the ability to map large areas rapidly. With aircraft systems this means a need for more beam positions over a wider swath on each flightline. For satellite systems the essential problem is resolution. Both of these needs are currently being addressed through the development and verification of Electronically Scanned Thinned Array Radiometer (ESTAR) technology. The ESTAR L band radiometer was evaluated for soil moisture mapping applications in two studies. The first was conducted over the semiarid rangeland Walnut Gulch watershed located in south eastern Arizona (U.S.). The second was performed in the subhumid Little Washita watershed in south west Oklahoma (U.S.). Both tests showed that the ESTAR is capable of providing soil moisture with the same level of accuracy as existing systems.

Description: A weather resistant automatic scanning Sun photometer system is assessed and demonstrated as practical for measurements of aerosol concentrations and properties at remote sites. Interfaced with a transmitter using the Geostationary Data Collection System (GDCS), the data are processed in near real time. The processing allows a time dependence of the aerosols and water vapor and an ongoing assessment of the health and calibration of the instruments. The system's automatic data acquisition, transmission, and processing offer immediate application to atmospheric monitoring and modeling on a regional to global scale and validation of satellite retrievals. It is estimated that under normal circumstances the retrieved aerosol optical thickness has a network wide accuracy of +/- 0.02 from 340 nm to 1020 nm, water vapor +/- 0.2 cm and size distribution from 0.1 to 3 micrometers.

Description: As part of a global program to validate the ocean surface sensors on board ERS-1, a joint experiment on the Grand Banks of Newfoundland was carried out in Nov. 1991. The principal objective was to provide a field validation of ERS-1 Synthetic Aperture Radar (SAR) measurement of ocean surface structure. The NASA-P3 aircraft measurements made during this experiment provide independent measurements of the ocean surface along the validation swath. The Radar Ocean Wave Spectrometer (ROWS) is a radar sensor designed to measure direction of the long wave components using spectral analysis of the tilt induced radar backscatter modulation. This technique greatly differs from SAR and thus, provides a unique set of measurements for use in evaluating SAR performance. Also, an altimeter channel in the ROWS gives simultaneous information on the surface wave height and radar mean square slope parameter. The sets of geophysical parameters (wind speed, significant wave height, directional spectrum) are used to study the SAR's ability to accurately measure ocean gravity waves. The known distortion imposed on the true directional spectrum by the SAR imaging mechanism is discussed in light of the direct comparisons between ERS-1 SAR, airborne Canadian Center for Remote Sensing (CCRS) SAR, and ROWS spectra and the use of the nonlinear ocean SAR transform.

Description: Terrain slopes, which can be measured with Synthetic Aperture Radar (SAR) interferometry either from a height map or from the interferometric phase gradient, were used to calculate the local incidence angle and the correct pixel area. Both are required for correct thematic interpretation of SAR data. The interferometric correlation depends on the pixel area projected on a plane perpendicular to the look vector and requires correction for slope effects. Methods for normalization of the backscatter and interferometric correlation for ERS-1 SAR are presented.

Description: The relationship between the gradient of the interferometric phase and the terrain slope, which, it is thought, would allow a derivation of the terrain slopes without phase unwrapping, is presented. A linear relationship between the interferometric phase gradient and the terrain slopes was found. A quantitative error analysis showed that only very small errors are introduced by these approximations for orbital Synthetic Aperture Radar (SAR) geometries. An example of a slope map for repeat pass interferometry from ERS-1 SAR data is given. A number of direct and indirect applications of the terrain slope are indicated: erosion and avalanche hazard studies, radiometric calibration of SAR data, and normalization of the interferometric correlation coefficient.

Description: Research on the use of active microwaves in remote sensing, presented during plenary and poster sessions, is summarized. The main highlights are: calibration techniques are well understood; innovative modeling approaches have been developed which increase active microwave applications (segmentation prior to model inversion, use of ERS-1 scatterometer, simulations); polarization angle and frequency diversity improves characterization of ice sheets, vegetation, and determination of soil moisture (X band sensor study); SAR (Synthetic Aperture Radar) interferometry potential is emerging; use of multiple sensors/extended spectral signatures is important (increase emphasis).

Description: Crop residues, the portion of the crop left in the field after harvest, can be an important management factor in controlling soil erosion. Methods to quantify residue cover are needed that are rapid, accurate, and objective. Scenes with known amounts of crop residue were illuminated with long wave ultraviolet (UV) radiation and fluorescence images were recorded with an intensified video camera fitted with a 453 to 488 nm band pass filter. A light colored soil and a dark colored soil were used as background for the weathered soybean stems. Residue cover was determined by counting the proportion of the pixels in the image with fluorescence values greater than a threshold. Soil pixels had the lowest gray levels in the images. The values of the soybean residue pixels spanned nearly the full range of the 8-bit video data. Classification accuracies typically were within 3(absolute units) of measured cover values. Video imaging can provide an intuitive understanding of the fraction of the soil covered by residue.

Description: Global study of land surface properties uses AVHRR channels 1 and 2, but channel 3 may be of interest, although its use requires preprocessing. It consists of both a reflective part and an emissive part, the former can be derived from T3, T4 and T5. Since the water vapor affects channel 3, its content is retrieved from the channel 4 and 5 using the split window technique. A formula of reflective part retrieval at 3.75 micrometers is tested in the case of sunglint observations where the emissivities of channels 4 and 5 can be set to the unity. The formula is adapted and validated to land surface using the FIFE-87 data set. Preliminary applications of the reflectance at 3.75 micrometers to the studies of surface properties retrieval, aerosol retrieval over land, and desertic aerosol retrieval, are addressed.

Description: The polarization of the sunlight scattered by atmospheric aerosols or cloud droplets and reflected from ground surfaces or plant canopies may convey much information when used for remote sensing purposes. The typical polarization features of aerosols, cloud droplets, and plant canopies, as observed by ground based and airborne sensors, are investigated, looking especially for those invariant properties amenable to description by simple models when possible. The question of polarization measurements from space is addressed. The interest of such measurements for remote sensing purposes is investigated, and their feasibility is tested by using results obtained during field campaigns of the airborne POLDER instrument, a radiometer designed to measure the directionality and polarization of the sunlight scattered by the ground atmosphere system.

Description: Papers focused on land surface, atmospheric, and ocean properties are reported. Specific comments pertaining to polarization, models and inversion, and measurements, are given. Recommendations are: continued research into the application potential of the BRDF (Bidirectional Reflectance Distribution Function) and polarization properties of ground surface and atmospheric targets; three dimensional models, which account for the statistical behavior of remotely sensed data, should be extended and inverted in order to support analysis of data potentially covering rolling terrain such that pixels represent heterogeneous mixtures of surface cover types and project ground footprints with sizes between 10 to 6 km, the ground pixel sizes of planned future sensors; available reflectance models should be further validated by means of multi dimensional (directional, spectral, temporal) field data and existing models should be intercompared in more depth to evaluate their performance and limitations; existing methods for model inversion should be validated in more depth in order to quantify the practical limitations and the expected accuracy of the parameters retrieved and new approaches should be developed based upon apriori knowledge of plant canopy development and spectral BRDF properties; there is a need to establish a protocol of validation and intercomparison of the indices and compositing techniques which have been proposed during these last years.

Description: Traditionally, the remote sensing community has relied totally on spectral knowledge to extract vegetation characteristics. However, there are other knowledge bases (KB's) that can be used to significantly improve the accuracy and robustness of inference techniques. Using AI (artificial intelligence) techniques a KB system (VEG) was developed that integrates input spectral measurements with diverse KB's. These KB's consist of data sets of directional reflectance measurements, knowledge from literature, and knowledge from experts which are combined into an intelligent and efficient system for making vegetation inferences. VEG accepts spectral data of an unknown target as input, determines the best techniques for inferring the desired vegetation characteristic(s), applies the techniques to the target data, and provides a rigorous estimate of the accuracy of the inference. VEG was developed to: infer spectral hemispherical reflectance from any combination of nadir and/or off-nadir view angles; infer percent ground cover from any combination of nadir and/or off-nadir view angles; infer unknown view angle(s) from known view angle(s) (known as view angle extension); and discriminate between user defined vegetation classes using spectral and directional reflectance relationships developed from an automated learning algorithm. The errors for these techniques were generally very good ranging between 2 to 15% (proportional root mean square). The system is designed to aid scientists in developing, testing, and applying new inference techniques using directional reflectance data.

Description: Most earth surfaces, particularly those supporting natural vegetation ecosystems, constitute structurally and spectrally complex surfaces that are distinctly non-Lambertian reflectors. Obtaining meaningful measurements of the directional radiances of landscapes and obtaining estimates of the complete bidirectional reflectance distribution functions of ground targets with complex and variable landscape and radiometric features are challenging tasks. Reasons for the increased interest in directional radiance measurements are presented, and the issues that must be addressed when trying to acquire directional radiances for vegetated land surfaces from different types of remote sensing platforms are discussed. Priority research emphases are suggested, concerning field measurements of directional surface radiances and reflectances for future research. Primarily, emphasis must be given to the acquisition of more complete and directly associated radiometric and biometric parameter data sets that will empower the exploitation of the 'angular dimension' in remote sensing of vegetation through enabling the further development and rigorous validation of state of the art plant canopy models.

Description: Synthetic aperture radar (SAR) images of the Greenland ice sheet collected by an airborne system clearly reveal the four melting facies of this sheet defined 30 years ago from snow stratigraphy studies by glaciologists. In particular, the radar echoes from the percolation facies have radiometric and polarimetric characteristics that are unique among terrestrial surfaces, but that resemble the exotic radar echoes recorded from the icy Galilean satellites. There, the radar signals interact with subsurface, massive ice features created in the cold, dry snow by seasonal melting and refreezing events. The subsurface features act as efficient reflectors of the incident radiation most likely via internal reflections. In the soaked-snow facies, the radar reflectivity is much lower because radar signals are attenuated by the wetter snow before they can interact with subsurface structures. Inversion algorithms to derive geophysical information from the SAR data are developed in both cases to estimate snow wetness in the soaked-snow facies and the mass of ice water retained in the percolation facies.

Description: An operational stratospheric correction scheme used after the Mount Pinatubo (Phillipines) eruption (Jun. 1991) is presented. The stratospheric aerosol distribution is assumed to be only variable with latitude. Each 9 days the latitudinal distribution of the optical thickness is computed by inverting radiances observed in the NOAA AVHRR channel 1 (0.63 micrometers) and channel 2 (0.83 micrometers) over the Pacific Ocean. This radiance data set is used to check the validity of model used for inversion by checking consistency of the optical thickness deduced from each channel as well as optical thickness deduced from different scattering angles. Using the optical thickness profile previously computed and radiative transfer code assuming Lambertian boundary condition, each pixel of channel 1 and 2 are corrected prior to computation of NDVI (Normalized Difference Vegetation Index). Comparison between corrected, non corrected, and years prior to Pinatubo eruption (1989 to 1990) NDVI composite, shows the necessity and the accuracy of the operational correction scheme.

Description: Aspects of aerosol studies and remote sensing are reviewed. Aerosol scatters solar radiation before it reaches the surface and scatters and absorbs it again after it is reflected from the surface and before it reaches the satellite sensor. The effect is spectrally and spatially dependent. Therefore atmospheric aerosol (dust, smoke and air pollution particles) has a significant effect on remote sensing. Correction for the aerosol effect was never achieved on an operational basis though several case studies were demonstrated. Correction can be done in a direct way by deriving the aerosol loading from the image itself and correcting for it using the appropriate radiative transfer model or by an indirect way, by defining remote sensing functions that are less dependent on the aerosol loading. To some degree this was already achieved in global remote sensing of vegetation where a composite of several days of NDVI (Normalized Difference Vegetation Index) measurements, choosing the maximal value, was used instead of a single cloud screened value. The Atmospheric Resistant Vegetation Index (ARVI) introduced recently for the NASA Earth Observing System EOS-MODIS is the most appropriate example of indirect correction, where the index is defined in such a way that the atmospheric effect in the blue spectral channel cancels to a large degree the atmospheric in the red channel in computations of a vegetation index. Atmospheric corrections can also use aerosol climatology and ground based instrumentation.

Description: A session dedicated to high spectral resolution in the solar spectrum, covering topics of calibration, atmospheric correction, geology/pedology, inland water, and vegetation, is reported. The session showed a high degree of diversity in the topics and the approaches used. It was highlighted that high spectral resolution data could provide atmospherically corrected ground level calibrated reflectance values. Important advances were shown in the use of radiative transfer models applied either on water bodies or vegetation. Several studies highlighted the high degree of redundancy contained in high spectral resolution data.

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: 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: 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: Spectral absorption-coefficients (cross-sections) kappa(sub nu) (/cm/atm) have been measured in the 7.62, 8.97, and 12.3 micrometer bands of HCFC-22 (CHClF2) and the 10.6 micrometer bands of SF6 employing a high-resolution Fourier-transform spectrometer. Temperature and total pressure have been varied to simulate conditions corresponding to tropospheric and stratospheric layers in the atmosphere. The kappa(sub nu) are compared with values measured by us previously using a tunable diode laser spectrometer and with the appropriate entries in HITRAN and GEISA, two of the databases known to the atmospheric scientist. The measured absolute intensities of the bands are compared with previously published values.

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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: Microwave radar and radiometer measurements of grasslands indicate a substantial reduction in sensor sensitivity to soil moisture in the presence of a thatch layer. When this layer is wet it masks changes in the underlying soil, making the canopy appear warm in the case of passive sensors (radiometer) and decreasing backscatter in the active case (scatterometer). A model for a grass canopy with thatch will be presented in this paper to explain this behavior and to compare with observations. The canopy model consists of three layers: grass, thatch, and the underlying soil. The grass blades are modeled by elongated elliptical discs and the thatch is modeled as a collection of disk shaped water droplets (i.e., the dry matter is neglected). The ground is homogeneous and flat. The distorted Born approximation is used to compute the radar cross section of this three layer canopy and the emissivity is computed from the radar cross section using the Peake formulation for the passive problem. Results are computed at L-band (1.4 GHz) and C-band (4.75 GHz) using canopy parameters (i.e., plant geometry, soil moisture, plant moisture, etc.) representative of Konza Prairie grasslands. The results are compared to C-band scatterometer measurements and L-band radiometer measurements at these grasslands.

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

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

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

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

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

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

Description: There have been many significant improvements in the public access to the Space Shuttle Earth Observations Photography Database. New information is provided for the user community on the recently released videodisc of this database. Topics covered included the following: earlier attempts; our first laser videodisc in 1992; the new laser videodisc in 1994; and electronic database access.

Description: The Challenge Awards are designed to provide a unique perspective to students gifted in the arts and humanities from which to understand scientific endeavor by giving students an opportunity to participate in an ongoing research project. In the graduate program, seven students who had participated in previous Challenge Awards programs were selected to help develop the tools for Earth observations for the astronauts on the Space Radar Laboratory (SRL) missions. The goal of the Challenge Awards program was to prepare a training manual for the astronauts on the SRL missions. This paper describes the observations to be made by the astronauts on the SRL missions. The emphasis is on the dynamic seasonal features of the Earth's surface and atmosphere which justify the need for more than one flight of the SRL. Complete notebooks of the sites, global seasonal patterns, examples of radar and the Measurement of Air Pollution from Satellites data, and shuttle photographs have been given to each of the SRL crews.

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

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

Description: MacSigma0 is an interactive tool for the Macintosh which allows you to display and make computations from radar data collected by the following sensors: the JPL AIRSAR, ERS-1, JERS-1, and Magellan. The JPL AIRSAR system is a multi-polarimetric airborne synthetic aperture radar developed and operated by the Jet Propulsion Laboratory. It includes the single-frequency L-band sensor mounted on the NASA CV990 aircraft and its replacement, the multi-frequency P-, L-, and C-band sensors mounted on the NASA DC-8. MacSigma0 works with data in the standard JPL AIRSAR output product format, the compressed Stokes matrix format. ERS-1 and JERS-1 are single-frequency, single-polarization spaceborne synthetic aperture radars launched by the European Space Agency and NASDA respectively. To be usable by MacSigma0, The data must have been processed at the Alaska SAR Facility and must be in the "low-resolution" format. Magellan is a spacecraft mission to map the surface of Venus with imaging radar. The project is managed by the Jet Propulsion Laboratory. The spacecraft carries a single-frequency, single-polarization synthetic aperture radar. MacSigma0 works with framelets of the standard MIDR CD-ROM data products. MacSigma0 provides four basic functions: synthesis of images (if necessary), statistical analysis of selected areas, analysis of corner reflectors as a calibration measure (if appropriate and possible), and informative mouse tracking. For instance, the JPL AIRSAR data can be used to synthesize a variety of images such as a total power image. The total power image displays the sum of the polarized and unpolarized components of the backscatter for each pixel. Other images which can be synthesized are HH, HV, VV, RL, RR, HHVV*, HHHV*, HVVV*, HHVV* phase and correlation coefficient images. For the complex and phase images, phase is displayed using color and magnitude is displayed using intensity. MacSigma0 can also be used to compute statistics from within a selected area. The statistics computed depend on the image type. For JPL AIRSAR data, the HH, HV, VV, HHVV* phase, and correlation coefficient means and standard deviation measures are calculated. The mean, relative standard deviation, minimum, and maximum values are calculated for all other data types. A histogram of the selected area is also calculated and displayed. The selected area can be rectangular, linear, or polygonal in shape. The user is allowed to select multiple rectangular areas, but not multiple linear or polygonal areas. The statistics and histogram are displayed to the user and can either be printed or saved as a text file. MacSigma0 can also be used to analyze corner reflectors as a measure of the calibration for JPL AIRSAR, ERS-1, and JERS-1 data types. It computes a theoretical radar cross section and the actual radar cross section for a selected trihedral corner reflector. The theoretical cross section, measured cross section, their ratio in dBs, and other information are displayed to the user and can be saved into a text file. For ERS-1, JERS-1, and Magellan data, MacSigma0 simultaneously displays pixel location in data coordinates and in latitude, longitude coordinates. It also displays sigma0, the incidence angle (for Magellan data), the original pixel value (for Magellan data), and the noise power value (for ERS-1 and JERS-1 data). Grey scale computed images can be saved in a byte format (a headerless format which saves the image as a string of byte values) or a PICT format (a standard format readable by other image processing programs for the Macintosh). Images can also be printed. MacSigma0 is written in C-language for use on Macintosh series computers. The minimum configuration requirements for MacSigma0 are System 6.0, Finder 6.1, 1Mb of RAM, and at least a 4-bit color or grey-scale graphics display. MacSigma0 is also System 7 compatible. To compile the source code, Apple's Macintosh Programmers Workbench (MPW) 3.2 and the MPW C language compiler version 3.2 are required. The source code will not compile with a later version of the compiler; however, the compiled application which will run under the minimum hardware configuration is provided on the distribution medium. In addition, the distribution media includes an executable which runs significantly faster but requires a 68881 compatible math coprocessor and a 68020 compatible CPU. Since JPL AIRSAR data files can be very large, it is often desirable to reduce the size of a data file before transferring it to the Macintosh for use in MacSigma0. A small FORTRAN program which can be used for this purpose is included on the distribution media. MacSigma0 will print statistics on any output device which supports QuickDraw, and it will print images on any device which supports QuickDraw or PostScript. The standard distribution medium for MacSigma0 is a set of five 1.4Mb Macintosh format diskettes. This program was developed in 1992 and is a copyrighted work with all copyright vested in NASA. Version 4.2 of MacSigma0 was released in 1993.

Description: VICAR (Video Image Communication and Retrieval) is a general purpose image processing software system that has been under continuous development since the late 1960's. Originally intended for data from the NASA Jet Propulsion Laboratory's unmanned planetary spacecraft, VICAR is now used for a variety of other applications including biomedical image processing, cartography, earth resources, and geological exploration. The development of this newest version of VICAR emphasized a standardized, easily-understood user interface, a shield between the user and the host operating system, and a comprehensive array of image processing capabilities. Structurally, VICAR can be divided into roughly two parts; a suite of applications programs and an executive which serves as the interfaces between the applications, the operating system, and the user. There are several hundred applications programs ranging in function from interactive image editing, data compression/decompression, and map projection, to blemish, noise, and artifact removal, mosaic generation, and pattern recognition and location. An information management system designed specifically for handling image related data can merge image data with other types of data files. The user accesses these programs through the VICAR executive, which consists of a supervisor and a run-time library. From the viewpoint of the user and the applications programs, the executive is an environment that is independent of the operating system. VICAR does not replace the host computer's operating system; instead, it overlays the host resources. The core of the executive is the VICAR Supervisor, which is based on NASA Goddard Space Flight Center's Transportable Applications Executive (TAE). Various modifications and extensions have been made to optimize TAE for image processing applications, resulting in a user friendly environment. The rest of the executive consists of the VICAR Run-Time Library, which provides a set of subroutines (image I/O, label I/O, parameter I/O, etc.) to facilitate image processing and provide the fastest I/O possible while maintaining a wide variety of capabilities. The run-time library also includes the Virtual Raster Display Interface (VRDI) which allows display oriented applications programs to be written for a variety of display devices using a set of common routines. (A display device can be any frame-buffer type device which is attached to the host computer and has memory planes for the display and manipulation of images. A display device may have any number of separate 8-bit image memory planes (IMPs), a graphics overlay plane, pseudo-color capabilities, hardware zoom and pan, and other features). The VRDI supports the following display devices: VICOM (Gould/Deanza) IP8500, RAMTEK RM-9465, ADAGE (Ikonas) IK3000 and the International Imaging Systems IVAS. VRDI's purpose is to provide a uniform operating environment not only for an application programmer, but for the user as well. The programmer is able to write programs without being concerned with the specifics of the device for which the application is intended. The VICAR Interactive Display Subsystem (VIDS) is a collection of utilities for easy interactive display and manipulation of images on a display device. VIDS has characteristics of both the executive and an application program, and offers a wide menu of image manipulation options. VIDS uses the VRDI to communicate with display devices. The first step in using VIDS to analyze and enhance an image (one simple example of VICAR's numerous capabilities) is to examine the histogram of the image. The histogram is a plot of frequency of occurrence for each pixel value (0 - 255) loaded in the image plane. If, for example, the histogram shows that there are no pixel values below 64 or above 192, the histogram can be "stretched" so that the value of 64 is mapped to zero and 192 is mapped to 255. Now the user can use the full dynamic range of the display device to display the data and better see its contents. Another example of a VIDS procedure is the JMOVIE command, which allows the user to run animations interactively on the display device. JMOVIE uses the concept of "frames", which are the individual frames which comprise the animation to be viewed. The user loads images into the frames after the size and number of frames has been selected. VICAR's source languages are primarily FORTRAN and C, with some VAX Assembler and array processor code. The VICAR run-time library is designed to work equally easily from either FORTRAN or C. The program was implemented on a DEC VAX series computer operating under VMS 4.7. The virtual memory required is 1.5MB. Approximately 180,000 blocks of storage are needed for the saveset. VICAR (version 2.3A/3G/13H) is a copyrighted work with all copyright vested in NASA and is available by license for a period of ten (10) years to approved licensees. This program was developed in 1989.

Description: The Interactive Image Display Program (IMDISP) is an interactive image display utility for the IBM Personal Computer (PC, XT and AT) and compatibles. Until recently, efforts to utilize small computer systems for display and analysis of scientific data have been hampered by the lack of sufficient data storage capacity to accomodate large image arrays. Most planetary images, for example, require nearly a megabyte of storage. The recent development of the "CDROM" (Compact Disk Read-Only Memory) storage technology makes possible the storage of up to 680 megabytes of data on a single 4.72-inch disk. IMDISP was developed for use with the CDROM storage system which is currently being evaluated by the Planetary Data System. The latest disks to be produced by the Planetary Data System are a set of three disks containing all of the images of Uranus acquired by the Voyager spacecraft. The images are in both compressed and uncompressed format. IMDISP can read the uncompressed images directly, but special software is provided to decompress the compressed images, which can not be processed directly. IMDISP can also display images stored on floppy or hard disks. A digital image is a picture converted to numerical form so that it can be stored and used in a computer. The image is divided into a matrix of small regions called picture elements, or pixels. The rows and columns of pixels are called "lines" and "samples", respectively. Each pixel has a numerical value, or DN (data number) value, quantifying the darkness or brightness of the image at that spot. In total, each pixel has an address (line number, sample number) and a DN value, which is all that the computer needs for processing. DISPLAY commands allow the IMDISP user to display all or part of an image at various positions on the display screen. The user may also zoom in and out from a point on the image defined by the cursor, and may pan around the image. To enable more or all of the original image to be displayed on the screen at once, the image can be "subsampled." For example, if the image were subsampled by a factor of 2, every other pixel from every other line would be displayed, starting from the upper left corner of the image. Any positive integer may be used for subsampling. The user may produce a histogram of an image file, which is a graph showing the number of pixels per DN value, or per range of DN values, for the entire image. IMDISP can also plot the DN value versus pixels along a line between two points on the image. The user can "stretch" or increase the contrast of an image by specifying low and high DN values; all pixels with values lower than the specified "low" will then become black, and all pixels higher than the specified "high" value will become white. Pixels between the low and high values will be evenly shaded between black and white. IMDISP is written in a modular form to make it easy to change it to work with different display devices or on other computers. The code can also be adapted for use in other application programs. There are device dependent image display modules, general image display subroutines, image I/O routines, and image label and command line parsing routines. The IMDISP system is written in C-language (94%) and Assembler (6%). It was implemented on an IBM PC with the MS DOS 3.21 operating system. IMDISP has a memory requirement of about 142k bytes. IMDISP was developed in 1989 and is a copyrighted work with all copyright vested in NASA. Additional planetary images can be obtained from the National Space Science Data Center at (301) 286-6695.

Description: PC-SEAPAK is a user-interactive satellite data analysis software package specifically developed for oceanographic research. The program is used to process and interpret data obtained from the Nimbus-7/Coastal Zone Color Scanner (CZCS), and the NOAA Advanced Very High Resolution Radiometer (AVHRR). PC-SEAPAK is a set of independent microcomputer-based image analysis programs that provide the user with a flexible, user-friendly, standardized interface, and facilitates relatively low-cost analysis of oceanographic satellite data. Version 4.0 includes 114 programs. PC-SEAPAK programs are organized into categories which include CZCS and AVHRR level-1 ingest, level-2 analyses, statistical analyses, data extraction, remapping to standard projections, graphics manipulation, image board memory manipulation, hardcopy output support and general utilities. Most programs allow user interaction through menu and command modes and also by the use of a mouse. Most programs also provide for ASCII file generation for further analysis in spreadsheets, graphics packages, etc. The CZCS scanning radiometer aboard the NIMBUS-7 satellite was designed to measure the concentration of photosynthetic pigments and their degradation products in the ocean. AVHRR data is used to compute sea surface temperatures and is supported for the NOAA 6, 7, 8, 9, 10, 11, and 12 satellites. The CZCS operated from November 1978 to June 1986. CZCS data may be obtained free of charge from the CZCS archive at NASA/Goddard Space Flight Center. AVHRR data may be purchased through NOAA's Satellite Data Service Division. Ordering information is included in the PC-SEAPAK documentation. Although PC-SEAPAK was developed on a COMPAQ Deskpro 386/20, it can be run on most 386-compatible computers with an AT bus, EGA controller, Intel 80387 coprocessor, and MS-DOS 3.3 or higher. A Matrox MVP-AT image board with appropriate monitor and cables is also required. Note that the authors have received some reports of incompatibilities between the MVP-AT image board and ZENITH computers. Also, the MVP-AT image board is not necessarily compatible with 486-based systems; users of 486-based systems should consult with Matrox about compatibility concerns. Other PC-SEAPAK requirements include a Microsoft mouse (serial version), 2Mb RAM, and 100Mb hard disk space. For data ingest and backup, 9-track tape, 8mm tape and optical disks are supported and recommended. PC-SEAPAK has been under development since 1988. Version 4.0 was updated in 1992, and is distributed without source code. It is available only as a set of 36 1.2Mb 5.25 inch IBM MS-DOS format diskettes. PC-SEAPAK is a copyrighted product with all copyright vested in the National Aeronautics and Space Administration. Phar Lap's DOS_Extender run-time version is integrated into several of the programs; therefore, the PC-SEAPAK programs may not be duplicated. Three of the distribution diskettes contain DOS_Extender files. One of the distribution diskettes contains Media Cybernetics' HALO88 font files, also licensed by NASA for dissemination but not duplication. IBM is a registered trademark of International Business Machines. MS-DOS is a registered trademark of Microsoft Corporation. HALO88 is a registered trademark of Media Cybernetics, but the product was discontinued in 1991.

Description: The Land Analysis System (LAS) is an image analysis system designed to manipulate and analyze digital data in raster format and provide the user with a wide spectrum of functions and statistical tools for analysis. LAS offers these features under VMS with optional image display capabilities for IVAS and other display devices as well as the X-Windows environment. LAS provides a flexible framework for algorithm development as well as for the processing and analysis of image data. Users may choose between mouse-driven commands or the traditional command line input mode. LAS functions include supervised and unsupervised image classification, film product generation, geometric registration, image repair, radiometric correction and image statistical analysis. Data files accepted by LAS include formats such as Multi-Spectral Scanner (MSS), Thematic Mapper (TM) and Advanced Very High Resolution Radiometer (AVHRR). The enhanced geometric registration package now includes both image to image and map to map transformations. The over 200 LAS functions fall into image processing scenario categories which include: arithmetic and logical functions, data transformations, fourier transforms, geometric registration, hard copy output, image restoration, intensity transformation, multispectral and statistical analysis, file transfer, tape profiling and file management among others. Internal improvements to the LAS code have eliminated the VAX VMS dependencies and improved overall system performance. The maximum LAS image size has been increased to 20,000 lines by 20,000 samples with a maximum of 256 bands per image. The catalog management system used in earlier versions of LAS has been replaced by a more streamlined and maintenance-free method of file management. This system is not dependent on VAX/VMS and relies on file naming conventions alone to allow the use of identical LAS file names on different operating systems. While the LAS code has been improved, the original capabilities of the system have been preserved. These include maintaining associated image history, session logging, and batch, asynchronous and interactive mode of operation. The LAS application programs are integrated under version 4.1 of an interface called the Transportable Applications Executive (TAE). TAE 4.1 has four modes of user interaction: menu, direct command, tutor (or help), and dynamic tutor. In addition TAE 4.1 allows the operation of LAS functions using mouse-driven commands under the TAE-Facelift environment provided with TAE 4.1. These modes of operation allow users, from the beginner to the expert, to exercise specific application options. LAS is written in C-language and FORTRAN 77 for use with DEC VAX computers running VMS with approximately 16Mb of physical memory. This program runs under TAE 4.1. Since TAE 4.1 is not a current version of TAE, TAE 4.1 is included within the LAS distribution. Approximately 130,000 blocks (65Mb) of disk storage space are necessary to store the source code and files generated by the installation procedure for LAS and 44,000 blocks (22Mb) of disk storage space are necessary for TAE 4.1 installation. The only other dependencies for LAS are the subroutine libraries for the specific display device(s) that will be used with LAS/DMS (e.g. X-Windows and/or IVAS). The standard distribution medium for LAS is a set of two 9~track 6250 BPI magnetic tapes in DEC VAX BACKUP format. It is also available on a set of two TK50 tape cartridges in DEC VAX BACKUP format. This program was developed in 1986 and last updated in 1992.

Description: Spaceborne Synthetic Aperture Radar (SAR) images are useful for planetary mapping and Earth sciences investigations. However, swath widths rarely exceed 100 Kilometers, and images must be patched together to create a mosaic in order to analyze larger areas. The primary function of this program is to generate large digital mosaics of SAR imagery without manually marked tiepoints. MOSK can produce multiframe mosaics by combining images in the along-track, adjacent cross-track swaths, or ascending and descending passes. Geocoded map registered images, such as the ones produced by MAPJTC (NPO-17718), are required as input. The output is a geocoded mosaic on a standard map grid which permits easy registration with other geocoded data sets. Mosaicking of geocoded SAR imagery involves three steps. First, a match point is selected at the center of the overlapping area, then an image patch around the match point is extracted from both images and cross-correlation is done on this area. Then, images with their refined match points are merged together to form a mosaic. To handle the large data volume of overlapping intermediate stages, large mosaics are divided into equal size quadrants with each quadrant cut from an intermediate mosaic. The full mosaic can then be assembled from the individual quadrants. Finally, radiometric disparities at the image seams are smoothed by a "feathering" technique. The automatic mosaic system generates output with minimal operator interaction. However, manual tiepointing is required in cases of a large registration error or two images with smooth surfaces such as ocean images. MOSK is implemented on a DEC VAX 11/785 running VMS 4.5. Most subroutines are in FORTRAN, but three are in MAXL and one is in APAL. The program requires 1 Mb of memory and a Floating Point Systems AP-5210 array processor. The system memory usage is approximately 1000 pages and the requirement of page file size is 2000 blocks. MOSK was developed in 1988.

Description: MAPJTC was designed to rectify and transform the standard image output of the digital Synthetic Aperture Radar (SAR) correlator into a geocoded map registered image without operator interaction or manual tiepointing. This is accomplished by modeling the distortions and predicting the pixel displacements based on platform and radar parameters. The map projection implemented in MAPJTC is the Universal Transverse Mercator. Since the re-sampling operation is independent of the transformation data generation, other cartographic projections can be implemented with few software modifications. MAPJTC makes a precise determination of the geodetic location of an arbitrary pixel within the image frame based on the simultaneous solution of a set of earth model equations, SAR Doppler equations, and SAR range equations that identify the slant range from the sensor to the target at a specific image pixel. Based on a table of geodetic coordinates of the image pixels, the image is then mapped onto the desired cartographic projection by applying the appropriate transformation equations. Typically, mapping involves a two-dimensional re-sampling and is very computationally intensive. MAPJTC reduces the procedure to two one-dimensional passes, which saves computer time. Geocoding transforms the rectified image into a grid defined by a specific map projection. (The image is rotated and rectified to match the map projection.) Again, the two dimensional re-sampling process can be separated into two one-dimensional re-sampling processes. Optionally, MAPJTC can correct terrain-induced distortions in SAR imagery when a digital elevation map is available. MAPJTC was developed on a DEC VAX 11/785 under VMS 4.5. The program is written in FORTRAN (84%), APAL (2%), and MAXL (14%). It requires 6Mb of memory and a Floating Point Systems AP-5210 Array Processor equipped with 1Mb of memory. MAPJTC can run interactively or as a batch job. MAPJTC was developed in 1987.

Description: The Spectral Analysis Manager (SPAM) was developed to allow easy qualitative analysis of multi-dimensional imaging spectrometer data. Imaging spectrometers provide sufficient spectral sampling to define unique spectral signatures on a per pixel basis. Thus direct material identification becomes possible for geologic studies. SPAM provides a variety of capabilities for carrying out interactive analysis of the massive and complex datasets associated with multispectral remote sensing observations. In addition to normal image processing functions, SPAM provides multiple levels of on-line help, a flexible command interpretation, graceful error recovery, and a program structure which can be implemented in a variety of environments. SPAM was designed to be visually oriented and user friendly with the liberal employment of graphics for rapid and efficient exploratory analysis of imaging spectrometry data. SPAM provides functions to enable arithmetic manipulations of the data, such as normalization, linear mixing, band ratio discrimination, and low-pass filtering. SPAM can be used to examine the spectra of an individual pixel or the average spectra over a number of pixels. SPAM also supports image segmentation, fast spectral signature matching, spectral library usage, mixture analysis, and feature extraction. High speed spectral signature matching is performed by using a binary spectral encoding algorithm to separate and identify mineral components present in the scene. The same binary encoding allows automatic spectral clustering. Spectral data may be entered from a digitizing tablet, stored in a user library, compared to the master library containing mineral standards, and then displayed as a timesequence spectral movie. The output plots, histograms, and stretched histograms produced by SPAM can be sent to a lineprinter, stored as separate RGB disk files, or sent to a Quick Color Recorder. SPAM is written in C for interactive execution and is available for two different machine environments. There is a DEC VAX/VMS version with a central memory requirement of approximately 242K of 8 bit bytes and a machine independent UNIX 4.2 version. The display device currently supported is the Raster Technologies display processor. Other 512 x 512 resolution color display devices, such as De Anza, may be added with minor code modifications. This program was developed in 1986.

Description: The Spectral Analysis Manager (SPAM) was developed to allow easy qualitative analysis of multi-dimensional imaging spectrometer data. Imaging spectrometers provide sufficient spectral sampling to define unique spectral signatures on a per pixel basis. Thus direct material identification becomes possible for geologic studies. SPAM provides a variety of capabilities for carrying out interactive analysis of the massive and complex datasets associated with multispectral remote sensing observations. In addition to normal image processing functions, SPAM provides multiple levels of on-line help, a flexible command interpretation, graceful error recovery, and a program structure which can be implemented in a variety of environments. SPAM was designed to be visually oriented and user friendly with the liberal employment of graphics for rapid and efficient exploratory analysis of imaging spectrometry data. SPAM provides functions to enable arithmetic manipulations of the data, such as normalization, linear mixing, band ratio discrimination, and low-pass filtering. SPAM can be used to examine the spectra of an individual pixel or the average spectra over a number of pixels. SPAM also supports image segmentation, fast spectral signature matching, spectral library usage, mixture analysis, and feature extraction. High speed spectral signature matching is performed by using a binary spectral encoding algorithm to separate and identify mineral components present in the scene. The same binary encoding allows automatic spectral clustering. Spectral data may be entered from a digitizing tablet, stored in a user library, compared to the master library containing mineral standards, and then displayed as a timesequence spectral movie. The output plots, histograms, and stretched histograms produced by SPAM can be sent to a lineprinter, stored as separate RGB disk files, or sent to a Quick Color Recorder. SPAM is written in C for interactive execution and is available for two different machine environments. There is a DEC VAX/VMS version with a central memory requirement of approximately 242K of 8 bit bytes and a machine independent UNIX 4.2 version. The display device currently supported is the Raster Technologies display processor. Other 512 x 512 resolution color display devices, such as De Anza, may be added with minor code modifications. This program was developed in 1986.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: IMAGEP is a FORTRAN computer algorithm containing various image processing, analysis, and enhancement functions. It is a keyboard-driven program organized into nine subroutines. Within the subroutines are other routines, also, selected via keyboard. Some of the functions performed by IMAGEP include digitization, storage and retrieval of images; image enhancement by contrast expansion, addition and subtraction, magnification, inversion, and bit shifting; display and movement of cursor; display of grey level histogram of image; and display of the variation of grey level intensity as a function of image position. This algorithm has possible scientific, industrial, and biomedical applications in material flaw studies, steel and ore analysis, and pathology, respectively. IMAGEP is written in VAX FORTRAN for DEC VAX series computers running VMS. The program requires the use of a Grinnell 274 image processor which can be obtained from Mark McCloud Associates, Campbell, CA. An object library of the required GMR series software is included on the distribution media. IMAGEP requires 1Mb of RAM for execution. The standard distribution medium for this program is a 1600 BPI 9~track magnetic tape in VAX FILES-11 format. It is also available on a TK50 tape cartridge in VAX FILES-11 format. This program was developed in 1991. DEC, VAX, VMS, and TK50 are trademarks of Digital Equipment Corporation.

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

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

Description: Analyses incorporating future data, as well as current and past data, are termed retrospective analyses. In this paper, the fixed-lag Kalman smoother (FLKS) is proposed as a means of providing retrospective analysis capability in data assimilation. The FLKS is a direct generalization of the Kalman filter. It incorporates all data observed up to and including some fixed amount of time past each analysis time. A computationally efficient form of the FLKS is derived. A simple scalar examination of the FLKS demonstrates that incorporating future data improves analyses the most in the presence of dynamical instabilities, for accurate models and for accurate observations. An implementation of the FLKS for two-dimensional linear shallow-water model corroborates the scalar analysis. The numerical experiments also demonstrate the ability of the FLKS to propagate information upstream as well as downstream, thus improving analysis quality substantially in data voids.

Description: An assimilation technique is described in which satellite-observed surface skin temperature tendencies are used in a model surface energy budget so that the predicted rate of temperature change in the model more closely agrees with the satellite observations. Both visible and infrared GOES satellite data are used in the assimilation. The technique is based on analytically recovering surface moisture from similarity expressions derived from an evapotranspiration residual obtained as a difference between the unadjusted model evapotranspiration and the satellite-inferred evapotranspiration. The technique has application in regional-scale models where surface parameters such as root zone moisture, soil moisture, etc., are unknown. It is assumed that the largest error in the surface energy budget is in the evapotranspiration term. Two tests are given for the technique, first, a one-dimensional test against FIFE data and, second, a three-dimensional test over Oklahoma. In these cases the technique appears to correctly adjust the model response to agree better with observations.

Description: Trajectory calculations using horizontal winds from the U.K. Meteorological Office data assimilation system and vertical velocities from a radiation calculation are used to simulate the three-dimensional motion of air through the stratospheric polar vortex for Northern Hemisphere (NH) and Southern Hemisphere (SH) winters since the launch of the Upper Atmosphere Research Satellite (UARS). Throughout the winter, air from the upper stratosphere moves poleward and descends into the middle stratosphere. In the SH lower to middle stratosphere, strongest descent occurs near the edge of the polar vortex, with that edge defined by mixing characteristics. The NH shows a similar pattern in late winter, but in early winter strongest descent is near the center of the vortex, except when wave activity is particularly strong. Strong barriers to latitudinal mixing exist above about 420 K throughout the winter. Below this, the polar night jet is weak in early winter, so air descending below that level mixes between polar and middle latitudes. In late winter, parcels descend less and the polar night jet moves downward, so there is less latitudinal mixing. The degree of mixing in the lower stratosphere thus depends strongly on the position and evolution of the polar night jet and on the amount of descent experienced by the air parcels; these characteristics show considerable interannual variability in both hemispheres. The computed trajectories provide a three-dimensional picture of air motion during the final warming. Large tongues of air are drawn off the vortex and stretched into increasingly long and narrow tongues extending into low latitudes. This vortex erosion process proceeds more rapidly in the NH than in he SH. In the lower stratosphere, the majority of air parcels remain confined within a lingering region of strong potential vorticity gradients into December in the SH and April in the NH, well after the vortex breaks up in the midstratosphere.

Description: Meteorological analyses, produced at the U.K. Meteorological Office by data assimilation, are used to study the circulation of the stratosphere in the Northern Hemisphere during winter 1991/92. The analyses are supplemented by Lagrangian visualizations of the circulation. The main features discussed are (1) the changes in vertical structure of the circulation, (2) the merger of anticyclones that precipitated a strong stratospheric warming, (3) vortex roll up in the upper stratosphere, (4) the entrainment of air into the polar vortex in the middle and upper stratosphere, and (5) the influence of tropospheric blocking on the lower stratosphere. The study provides a meteorological basis for the interpretation of data from the Upper Atmosphere Research Satellite (UARS).

Description: The evolution of zonal wind and zonal wavenumber one (wave 1) in the Southern Hemisphere subpolar middle atmosphere is described for the period December 1978 - May 1979 using temperature and ozone measurements from the Limb Infrared Monitor of the Stratosphere (LIMS) experiment. In late December maximum zonal easterlies of approx. -70 m/s are observed at 0.1 mb, 60 deg S. A zonal flow reversal occurs during late February and westerlies subsequently increase to 60-70 m/s in the upper stratosphere by April - May. LIMS zonal winds are compared with rocketsonde measurements and nadir sounder (derived) winds for summer and autumn. Although quantitative agreement is found at stratospheric levels, substantial discrepancies are evident in the mesosphere, most likely a reflection of sampling and resolution differences in the respective datasets. Stationary and traveling wave 1 temperature disturbances (amplitudes approx. 1 - 2 K at 60 deg S) are observed by LIMS during summer. The stationary wave is confined to the lower stratosphere near the level of zero zonal- mean wind flow, whereas the traveling wave is prominent in the middle stratosphere moves west at a rate similar to the zonal-mean wind, and exhibits a vertical - meridional structure similar to a P(sub 4)(sup 1) normal mode Rossby wave. A substantial intensification of wave 1 activity occurs during autumn (amplitudes approx. 5 - 10 K), which is found to be associated with an upward-directed Eliasse - Palm flux near the subpolar tropopause level. Evidence relating wave 1 activity in the lower - middle stratosphere to the occurrence of zonal ozone perturbations of 10% - 20% amplitude is presented for summer and autumn.

Description: The evolution of the vertical electric field profile deduced from simultaneous field measurements at several levels below a thundercloud shows the development of a space charge layer at least up to 600 m. The average charge density in the whole layer from 0 m to 600 m can reach about 1 nC m(exp -3). The ions are generated at the ground by corona effect and the production rate is evaluated with a new method from the comparison of field evolutions at the ground and at altitude after a lightning flash. The modeling of the relevant processes shows tht ground corona accounts for the observed field evolutions and that the aerosol particles concentration has a very large effect on the evolution of corona ions. However, with a realistic value for this concentration a large amount of ground corona ions reach the level of 600 m.

Description: The one-dimensional cirrus model described in part 1 of this issue has been used to study the sensitivity of simulated cirrus microphysical and radiative properties to poorly known model parameters, poorly understood physical processes, and environmental conditions. Model parameters and physical processes investigated include nucleation rate, mode of nucleation (e.g., homogeneous freezing of aerosols and liquid droplets or heterogeneous deposition), ice crystal shape, and coagulation. These studies suggest that the leading sources of uncertainty in the model are the phase change (liquid-solid) energy barrier and the ice-water surface energy which dominate the homogeneous freezing nucleation rate and the coagulation sticking efficiency at low temperatures which controls the production of large ice crystals (radii greater than 100 mcirons). Environmental conditions considered in sensitivity tests were CN size distribution, vertical wind speed, and cloud height. We found that (unlike stratus clouds) variations in the total number of condensation nuclei (NC) have little effect on cirrus microphysical and radiative properties, since nucleation occurs only on the largest CN at the tail of the size distribution. The total number of ice crystals which nucleate has little or no relationship to the number of CN present and depends primarily on the temperature and the cooling rate. Stronger updrafts (more rapid cooling) generate higher ice number densities, ice water content, cloud optical depth, and net radiative forcing. Increasing the height of the clouds in the model leads to an increase in ice number density, a decrease in effective radius, and a decrease in ice water content. The most prominent effect of increasing cloud height was a rapid increase in the net cloud radiative forcing which can be attributed to the change in cloud temperature as well as change in cloud ice size distributions. It has long been recognized that changes in cloud height or cloud area have the greatest potential for causing feedbacks on climate change. Our results suggest that variations in vertical velocity or cloud microphysical changes associatd with cloud height changes may also be important.

Description: Special Sensor Microwave/Imager (SSM/I) observations were used to examine the spatial and temporal changes of the precipitation characteristics of tropical cyclones. SSM/I observations were also combined with the results of a tropical cyclone numerical model to examine the role of inner-core diabatic heating in subsequent intensity changes of tropical cyclones. Included in the SSM/I observations were rainfall characteristics of 18 named western North Atlantic tropical cyclones between 1987 and 1989. The SSM/I rain-rate algorithm that employed the 85-GHz channel provided an analysis of the rain-rate distribution in greater detail. However, the SSM/I algorithm underestimated the rain rates when compared to in situ techniques but appeared to be comparable to the rain rates obtained from other satellite-borne passive microwave radiometers. The analysis of SSM/I observations found that more intense systems had higher rain rates, more latent heat release, and a greater contribution from heavier rain to the total tropical cyclone rainfall. In addition, regions with the heaviest rain rates were found near the center of the most intense tropical cyclones. Observational analysis from SSM/I also revealed that the greatest rain rates in the inner-core regions were found in the right half of fast-moving cyclones, while the heaviest rain rates in slow-moving tropical cyclones were found in the forward half. The combination of SSM/I observations and an interpretation of numerical model simulations revealed that the correlation between changes in the inner core diabetic heating and the subsequent intensity became greater as the tropical cyclones became more intense.

Description: The intensity, spatial, and temporal changes in precipitation were examined in three North Atlantic hurricanes during 1989 (Dean, Gabrielle, and Hugo) using precipitation estimates made from Special Sensor Microwave/Imager (SSM/I) measurements. In addition, analyses from a barotropic hurricane forecast model and the European Centre for Medium-Range Weather Forecast model were used to examine the relationship between the evolution of the precipitation in these tropical cyclones and external forcing. The external forcing parameters examined were (1) mean climatological sea surface temperatures, (2) vertical wind shear, (3) environmental tropospheric water vapor flux, and (4) upper-tropospheric eddy relative angular momentum flux convergence. The analyses revealed that (1) the SSM/I precipitation estimates were able to delineate and monitor convective ring cycles similar to those observed with land-based and aircraft radar and in situ measurements; (2) tropical cyclone intensification was observed to occur when these convective rings propagated into the inner core of these systems (within 111 km of the center) and when the precipitation rates increased; (3) tropical cyclone weakening was observed to occur when these inner-core convective rings dissipated; (4) the inward propagation of the outer convective rings coincided with the dissipation of the inner convective rings when they came within 55 km of each other; (5) in regions with the combined warm sea surface temperatures (above 26 C) and low vertical wind shear (less than 5 m/s), convective rings outside the region of strong lower-tropospheric inertial stability could be initiated by strong surges of tropospheric moisture, while convective rings inside the region of strong lower-tropospheric inertial stability could be enhanced by upper-tropospheric eddy relative angular momentum flux convergence.

Description: Measurements of ozone (O3) and aerosol distributions were made with an airborne lidar system in the lowland and boreal forest regions of eastern Canada during July - August 1990 as part of the NASA Global Tropospheric Experiment/Arctic Boundary Layer Expedition (ABLE) 3B. Aerosol and O3 profiles were measured simultaneously above and below the Electra aircraft from near the surface to above the tropopause on long-range flights over these important ecosystems. A broad range of atmospheric conditions were encountered during repeated flights over intensive study sites in the Hudson Bay lowlands near Moosonee, Ontario, and over the boreal forest near Schefferville, Quebec. The tropospheric composition in this high-latitude region was found to be strongly influenced by stratospheric intrusions. Regions of low aerosol scattering and enhanced O3 mixing ratios were correlated with descending air from the lower stratosphere. Over 33% of the troposphere (0-12 km) along our flight track at latitudes from about 45 deg to 55 deg N had significantly enhanced O3 due to stratospheric intrusions, and in the middle to upper troposphere the extent of the enhanced O3 gnerally exceeded 40%. Ozone mixing ratios of 80 parts per billion by volume (ppbv) near 6 km were common in strong intrusions. In the boundary layer over the lowlands, O3 was in the 20-30 ppbv range with a vertical O3 gradient of 6.7 ppbv/km to about 45 ppbv at 3 km. Above 6 km the background tropospheric O3 profile was nearly constant with an average value of 53 ppbv. Due to forest fires in Canada and Alaska, plumes from biomass-burning sources were observed on many flights. Biomass-burning plumes influenced about 25% of the free troposphere below 4 km, and in some of the plumes, O3 was enhanced by 10-20 ppbv over ambient levels of 30-45 ppbv. Several air masses transported from the tropical Pacific were observed over Canada in the middle to upper troposphere with O3 levels 10-20 ppbv below background values of 50-55 ppbv.

Description: Airborne heat, moisture, O3, CO, and CH4 flux measurements were obtained over the Hudson Bay lowlands (HBL) and northern boreal forest regions of Canada during July - August 1990. The airborne flux measurements were an integral part of the NASA/Arctic Boundary Layer Expedition (ABLE) 3B field experiment executed in collaboration with the Canadian Northern Wetlands Study (NOWES). Airborne CH4 flux measurements were taken over a large portion of the HBL. The surface level flux of CH4 was obtained from downward extrapolations of multiple-level CH4 flux measurements. Methane source strengths ranged from -1 to 31 mg m(exp -2)/d, with the higher values occurring in relatively small, isolated areas. Similar measurements of the CH4 source strength in the boreal forest region of Schefferville, Quebec, ranged from 6 to 27 mg m(exp -2)/d and exhibited a diurnal dependence. The CH4 source strengths found during the ABLE 3B expedition were much lower than the seasonally averaged source strength of 51 mg m(exp -2)/d found for the Yukon-Kuskokwim delta region of Alaska during the previous ABLE 3A study. Large positive CO fluxes (0.31 to 0.53 parts per billion by volume (ppbv) m/s) were observed over the inland, forested regions of the HBL study area, although the mechanism for the generation of these fluxes was not identified. Repetitive measurements along the same ground track at various times of day near the Schefferville site also suggested a diurnal dependence for CO emissions. Measurements of surface resistance to the uptake of O3 (1.91 to 0.80 s/cm) for the HBL areas investigated were comparable to those observed near the Schefferville site (3.40 to 1.10 s/cm). Surface resistance values for the ABLE 3B study area were somewhat less than those observed over the Yukon-Kuskokwim delta during the previous ABLE 3A study. The budgets for heat, moisture, O3, CO, and CH4 were evaluated. The residuals from these budget studies indicated, for the cases selected, a moderate net photochemical production of O3 present in the boundary layer over the HBL that coincided with an in situ destruction of CO, although the mechanism responsible for the destruction of CO was not identified. Results from the O3 budget analysis indicate the importance of in situ photochemical production and its possible dominance over surface deposition to the local O3 budget at the Schefferville site. Measurements of the in situ production of O3 indicated a direct relationship between the presence of biomass burning or large-scale pollution effects. Residuals from budget calculations for conserved quantities (heat, moisture, and CH4) were compared with their respective surface fluxes to provide a measure of the internal self-consistency of the flux measurements.

Description: The Arctic Boundary Layer Expedition (ABLE) 3B was conducted to determine the summertime tropospheric distribution, sources, and sinks of important trace gas and aerosol species over the wetlands and boreal forests of central and eastern Canada. Isentropic trajectories and analyzed midtropospheric circulation patterns were used to group flights according to the transport histories of polar, midlatitude, or tropical air masses which were sampled. These data were then divided into bands of potential temperature levels representing the low, middle, and maximum aircraft altitudes to assess the effects of both local and long distance transport and natural and man-made pollutants to the measured chemical species. Detailed case studies are provided to depict the complex three-dimensional airflow regimes that transported air with differing chemical signatures to the study area. Mission 6 details the large-scale movement of smoke in the generally prevailing west to northwesterly airflow that was observed on the majority of flights. Mission 1 analyzes the horizontal and vertical motions of maritime Pacific air in the upper troposphere that was routinely mixed downward to the aircraft altitude. Finally, mission 14 tracks the far northward excursion of tropical air that had been associated with a Pacific typhoon. The following three factors all had important influences on the collected chemical data sets: (1) local and distant stratospheric in puts into the upper and middle troposphere; (2) biomass-burning plumes from active fires in Alaska and Canada; (3) a band of 'low ozone' upper tropospheric air that was observed by airborne differential absorption lidar (DIAL) above the aircraft maximum altitude. Other modification factors observed on some flights included urban pollution from U.S. and Canadian cities, tropical air that had been associated with a Pacific typhoon, and precipitation scavenging by clouds and rain. Many flights were affected by several of the above factors which led to complex chemical signatures that will be discussed in other companion papers.

Description: The use of an airborne CO2 lidar to obtain cloud backscatter and extinction data at a thermal infrared wavelength is described. The extinction in this spectral region is proportional to the cloud liquid water content. The use of coherent detection results in high sensitivity and narrow field of view, the latter property greatly reducing multiple-scattering effects. Backscatter measurements in absolute units are obtained through a hard target calibration methodology. For clouds of low to moderate optical thickness at the lidar wavelength, both geometric thickness and optical thickness can be measured. The sea surface reflectance signal is used to obtain estimates of the cloud optical thickness. The utility of this technique results from studies that indicate that the spatial scale of variability of the sea surface reflectance is generally large compared with that of cloud optical thickness. Selected results are presented from data taken during flights over the Pacific Ocean.

Description: Airborne dual-wavelength and dual-polarization radar data are analyzed for measurements taken in stratiform rain in the western Pacific during September 1990. The focus of the paper is on the vertical profiles of the linear depolarization ratio, LDR (10 GHz); the reflectivity factor, dBZ (10 GHz); and the dual-frequency ratio, DFR (10, 34.45 GHz). Statistical characterizations of the maxima of these quantities and the relative locations at which they occur suggest that the eccentricity of the melting particles is fairly large and that the shape and size of the particles are correlated. To try to explain these features, two types of simulation are presented. In the first, a set of measured drop size distributions is used in the context of a standard model of the melting layer. Variations in snow density, as well as shape, size, and orientation distributions are used to study the relationship between these parameters and the radar measurements. To reduce the amount of ambiguity in the estimation, a second type of simulation is described in which the size distribution of the snow is estimated. Comparisons between the simulated and measured profiles indicate that radar measurements can be used to derive certain characteristics of the particle size and shape distributions in the melting layer.

Description: A brief mission overview of STS-49 is given, and some of the pictorially outstanding and scientifically interesting photographs obtained during the mission are presented. The Earth observations are described and include the following: the Southwestern Pacific Ocean -- wind and water; the Southwestern Pacific Ocean -- coasts and volcanoes; the US; Cuba and the Bahamas; South America; Africa; the Red Sea and Western Indian Ocean; and the Indian Subcontinent.

Description: This paper provides a description of the variability of global atmospheric angular momentum (GAM) and its relationship with principal modes of three-dimensional atmospheric circulation anomalies. The data used are 5-day mean global wind fields from the European Centre for Medium-Range Weather Forecasts initialized dataset for 1980-1989. Significant seasonal variation of GAM is observed with maxima in April and November and a minimum during late July. The amplitude of the annual cycle is largest in the upper troposphere and decreases toward the surface. Although the lower tropospheric contribution to the total angular momentum is relatively small, its annual cycle is out of phase with those of the upper atmosphere and GAM. Also identified is a distinct semiannual component, with double peaks appearing in April and November. This signal is most noticeable in the upper troposphere above the 300-mb level. The principal modes of zonal-mean angular momentum and meridional circulation anomalies and their coupled modes are obtained by using empirical orthogonal function analysis and singular value decomposition. It is shown that the leading modes of the angular momentum and meridional circulation are coupled with each other and are responsible for much of the variability in GAM. The coupled modes represent fluctuations of upper-level subtropical zonal flow, which are linked to the modulation of Hadley circulation intensity in both hemispheres. It is found that GAM is highly correlated with the first eigenvector of upper-level streamfunction anomalies, which consists of a superrotational flow in the tropics and subtropics, except over the central Pacific where a 'blocked' flow with two subtropical anticyclonic circulation cells straddling the equator is found. Much of the blocked flow is due to the establishment of dipole anomalies in the velocity potential with centers over the central Pacific and the Maritime Continent on the interannual time scale. On the intraseasonal time scale, GAM fluctuation is dominated by superrotational flow in the tropics, with the blocked flow present to a much lesser extent. The associated velocity potential anomaly has a weak dipole structure with centers over the Indian Ocean and the eastern Pacific. The implications of the above results on the total angular momentum balance of the earth-atmosphere system are also discussed.

Description: Simulations with the UCLA atmospheric general circulation model (AGCM) using two different global sea surface temperature (SST) datasets for January 1979 are compared. One of these datasets is based on Comprehensive Ocean-Atmosphere Data Set (COADS) (SSTs) at locations where there are ship reports, and climatology elsewhere; the other is derived from measurements by instruments onboard NOAA satellites. In the former dataset (COADS SST), data are concentrated along shipping routes in the Northern Hemisphere; in the latter dataset High Resolution Infrared Sounder (HIRS SST), data cover the global domain. Ensembles of five 30-day mean fields are obtained from integrations performed in the perpetual-January mode. The results are presented as anomalies, that is, departures of each ensemble mean from that produced in a control simulation with climatological SSTs. Large differences are found between the anomalies obtained using COADS and HIRS SSTs, even in the Northern Hemisphere where the datasets are most similar to each other. The internal variability of the circulation in the control simulation and the simulated atmospheric response to anomalous forcings appear to be linked in that the pattern of geopotential height anomalies obtained using COADS SSTs resembles the first empirical orthogonal function (EOF 1) in the control simulation. The corresponding pattern obtained using HIRS SSTs is substantially different and somewhat resembles EOF 2 in the sector from central North America to central Asia. To gain insight into the reasons for these results, three additional simulations are carried out with SST anomalies confined to regions where COADS SSTs are substantially warmer than HIRS SSTs. The regions correspond to warm pools in the northwest and northeast Pacific, and the northwest Atlantic. These warm pools tend to produce positive geopotential height anomalies in the northeastern part of the corresponding oceans. Both warm pools in the Pacific produce large-scale circulation anomalies with a pattern that resembles that obtained using COADS SSTs as well as EOF 1 of the control simulation; the warm pool in the Atlantic does not. These results suggest that the differences obtained with COADS SSTs and HIRS SSTs are mostly due to the differences in the datasets over the northern Pacific. There was a blocking episode near Greenland in late January 1979. Both simulations with warm SST anomalies over the northwest and northeast Pacific show a tendency toward increased incidence of North Atlantic blocking; the simulation with warm SST anomalies over the northwest Atlantic shows a tendency toward decreased incidence. These results suggest that features in both SST datasets that do not have a counterpart in the other dataset contribute signficantly to the differences between the simulated and observed fields. The results of this study imply that uncertainties in current SST distributions for the world oceans can be as important as the SST anomalies themselves in terms of their impact on the atmospheric circulation. Caution should be exercised, therefore, when linking anomalous circulation and SST patterns, especially in long-range prediction.

Description: The hydrologic cycle is a major part of the global climate system. There is an atmospheric flux of water from the ocean surface to the continents. The cycle is closed by return flow in rivers. In this paper a river routing model is developed to use with grid box climate models for the whole earth. The routing model needs an algorithm for the river mass flow and a river direction file, which has been compiled for 4 deg x 5 deg and 2 deg x 2.5 deg resolutions. River basins are defined by the direction files. The river flow leaving each grid box depends on river and lake mass, downstream distance, and an effective flow speed that depends on topography. As input the routing model uses monthly land source runoff from a 5-yr simulation of the NASA/GISS atmospheric climate model (Hansen et al.). The land source runoff from the 4 deg x 5 deg resolution model is quartered onto a 2 deg x 2.5 deg grid, and the effect of grid resolution is examined. Monthly flow at the mouth of the world's major rivers is compared with observations, and a global error function for river flow is used to evaluate the routing model and its sensitivity to physical parameters. Three basinwide parameters are introduced: the river length weighted by source runoff, the turnover rate, and the basinwide speed. Although the values of these parameters depend on the resolution at which the rivers are defined, the values should converge as the grid resolution becomes finer. When the routing scheme described here is coupled with a climate model's source runoff, it provides the basis for closing the hydrologic cycle in coupled atmosphere-ocean models by realistically allowing water to return to the ocean at the correct location and with the proper magnitude and timing.

Description: The authors analyze the influence of Sea Surface Temperature (SST) and surface wind divergence on atmospheric thermodynamic structure and the resulting effects on the occurrence of deep convection using National Meteorological Center radiosonde data and International Satellite Cloud Climatology Program data for July 1983-July 1985. The onset of deep convection requires not only the existence of positive convective available potential energy (CAPE), but also an unstable planetary boundary layer (PBL). A stable PBL is observed to suppress deep convection even when CAPE is positive. Variations of SST have a major effect on CAPE, but surface wind divergence can also affect deep convection by changing the lapse rate in the lower troposphere and humidity in the PBL. Specifically, when SST is greater than or equal to 28 C, CAPE is always positive, and surface wind divergence does not qualitatively change the buoyancy profile above the PBL. Strong surface wind divergence, however, stabilizes the PBL so as to suppress the initiation of deep convection. In warm SST regions, CAPE is greater than 0 regardless of assumptions about condensate loading, although the pseudoadiabatic limit is more consistent with the observed deep convection than the reversible moist-adiabatic limit under these circumstances. When SST is less than 27 C, CAPE is usually negative and inhibits convection, but strong surface wind convergence can destabilize the inversion layer and moisten the PBL enough to make the atmosphere neutrally stable in the mean. As a result, deep convection is generally enhanced either when SST is greater than or equal to 28 C in the absence of strong surface wind divergence or when strong surface wind convergence occurs even if SST is less than 27 C. The anomalous suppression of deep convection in the warm area of the equatorial west Pacific lying between the intertropical convergence zone (ITCZ) and south Pacific convergence zone (SPCZ) is probably caused by dryness in the PBL and an inversion in that area. The seasonal cycles of deep convection and surface wind divergence are in phase with the maximum solar radiation and lead SST for one to three months in the central Pacific. The change of PBL relative humidity plays a critical role in the changeover to convective instability in this case. The seasonal change of deep convection and associated clouds seems not to have important effects on the seasonal change of local SST in the central Pacific.