ALBERT

Description: The Interactive Controls Analysis (INCA) program was developed to provide a user friendly environment for the design and analysis of linear control systems, primarily feedback control systems. INCA is designed for use with both small and large order systems. Using the interactive graphics capability, the INCA user can quickly plot a root locus, frequency response, or time response of either a continuous time system or a sampled data system. The system configuration and parameters can be easily changed, allowing the INCA user to design compensation networks and perform sensitivity analysis in a very convenient manner. A journal file capability is included. This stores an entire sequence of commands, generated during an INCA session into a file which can be accessed later. Also included in INCA are a context-sensitive help library, a screen editor, and plot windows. INCA is robust to VAX-specific overflow problems. The transfer function is the basic unit of INCA. Transfer functions are automatically saved and are available to the INCA user at any time. A powerful, user friendly transfer function manipulation and editing capability is built into the INCA program. The user can do all transfer function manipulations and plotting without leaving INCA, although provisions are made to input transfer functions from data files. By using a small set of commands, the user may compute and edit transfer functions, and then examine these functions by using the ROOT_LOCUS, FREQUENCY_RESPONSE, and TIME_RESPONSE capabilities. Basic input data, including gains, are handled as single-input single-output transfer functions. These functions can be developed using the function editor or by using FORTRAN- like arithmetic expressions. In addition to the arithmetic functions, special functions are available to 1) compute step, ramp, and sinusoid functions, 2) compute closed loop transfer functions, 3) convert from S plane to Z plane with optional advanced Z transform, and 4) convert from Z plane to W plane and back. These capabilities allow the INCA user to perform block diagram algebraic manipulations quickly for functions in the S, Z, and W domains. Additionally, a versatile digital control capability has been included in INCA. Special plane transformations allow the user to easily convert functions from one domain to another. Other digital control capabilities include: 1) totally independent open loop frequency response analyses on a continuous plant, discrete control system with a delay, 2) advanced Z-transform capability for systems with delays, and 3) multirate sampling analyses. The current version of INCA includes Dynamic Functions (which change when a parameter changes), standard filter generation, PD and PID controller generation, incorporation of the QZ-algorithm (function addition, inverse Laplace), and describing functions that allow the user to calculate the gain and phase characteristics of a nonlinear device. The INCA graphic modes provide the user with a convenient means to document and study frequency response, time response, and root locus analyses. General graphics features include: 1) zooming and dezooming, 2) plot documentation, 3) a table of analytic computation results, 4) multiple curves on the same plot, and 5) displaying frequency and gain information for a specific point on a curve. Additional capabilities in the frequency response mode include: 1) a full complement of graphical methods Bode magnitude, Bode phase, Bode combined magnitude and phase, Bode strip plots, root contour plots, Nyquist, Nichols, and Popov plots; 2) user selected plot scaling; and 3) gain and phase margin calculation and display. In the time response mode, additional capabilities include: 1) support for inverse Laplace and inverse Z transforms, 2) support for various input functions, 3) closed loop response evaluation, 4) loop gain sensitivity analyses, 5) intersample time response for discrete systems using the advanced Z transform, and 6) closed loop time response using mixed plane (S, Z, W) operations with delay. A Graphics mode command was added to the current version of INCA, version 3.13, to produce Metafiles (graphic files) of the currently displayed plot. The metafile can be displayed and edited using the QPLOT Graphics Editor and Replotter for Metafiles (GERM) program included with the INCA package. The INCA program is written in Pascal and FORTRAN for interactive or batch execution and has been implemented on a DEC VAX series computer under VMS. Both source code and executable code are supplied for INCA. Full INCA graphics capabilities are supported for various Tektronix 40xx and 41xx terminals; DEC VT graphics terminals; many PC and Macintosh terminal emulators; TEK014 hardcopy devices such as the LN03 Laserprinter; and bit map graphics external hardcopy devices. Also included for the TEK4510 rasterizer users are a multiple copy feature, a wide line feature, and additional graphics fonts. The INCA program was developed in 1985, Version 2.04 was released in 1986, Version 3.00 was released in 1988, and Version 3.13 was released in 1989. An INCA version 2.0X conversion program is included.

Description: Expensive analysis programs are often combined with optimization procedures to solve engineering problems. An optimal solution requires numerous iterations between the analysis program and an optimizer. This often becomes prohibitive due to cost and amount of computer time needed to converge to an optimal solution. NETS/PROSSS was developed to provide a system for combining NETS (MSC-21588), a neural network program developed at NASA's Johnson Space Center, and the optimization program CONMIN (Constrained Function Minimization, ARC-10836) developed at Ames Research Center. After training, NETS approximates the results from the analysis program, possibly allowing the user to reach a near-optimal solution in much less time than before. These results can then be used as a starting point in a normal optimization process, possibly allowing the user to converge to an optimal solution in significantly fewer iterations. NETS/PROSSS is written in C-language and FORTRAN 77 for Sun series computers running SunOS. The required CONMIN and NETS v3.0 files are included in this package. The documentation for CONMIN and NETS are included with the documentation of NETS/PROSSS. The program requires 342K of RAM for execution. The standard distribution medium for this program is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. NETS/PROSSS was developed in 1991.

Description: AESOP was developed to solve a number of problems associated with the design of controls and state estimators for linear time-invariant systems. The systems considered are modeled in state-variable form by a set of linear differential and algebraic equations with constant coefficients. Two key problems solved by AESOP are the linear quadratic regulator (LQR) design problem and the steady-state Kalman filter design problem. AESOP is designed to be used in an interactive manner. The user can solve design problems and analyze the solutions in a single interactive session. Both numerical and graphical information are available to the user during the session. The AESOP program is structured around a list of predefined functions. Each function performs a single computation associated with control, estimation, or system response determination. AESOP contains over sixty functions and permits the easy inclusion of user defined functions. The user accesses these functions either by inputting a list of desired functions in the order they are to be performed, or by specifying a single function to be performed. The latter case is used when the choice of function and function order depends on the results of previous functions. The available AESOP functions are divided into several general areas including: 1) program control, 2) matrix input and revision, 3) matrix formation, 4) open-loop system analysis, 5) frequency response, 6) transient response, 7) transient function zeros, 8) LQR and Kalman filter design, 9) eigenvalues and eigenvectors, 10) covariances, and 11) user-defined functions. The most important functions are those that design linear quadratic regulators and Kalman filters. The user interacts with AESOP when using these functions by inputting design weighting parameters and by viewing displays of designed system response. Support functions obtain system transient and frequency responses, transfer functions, and covariance matrices. AESOP can also provide the user with open-loop system information including stability, controllability, and observability. The AESOP program is written in FORTRAN IV for interactive execution and has been implemented on an IBM 3033 computer using TSS 370. As currently configured, AESOP has a central memory requirement of approximately 2 Megs of 8 bit bytes. Memory requirements can be reduced by redimensioning arrays in the AESOP program. Graphical output requires adaptation of the AESOP plot routines to whatever device is available. The AESOP program was developed in 1984.

Description: The Automation Technology Branch of NASA's Langley Research Center is employing increasingly complex degrees of operator/robot cooperation (telerobotics). A good relationship between the operator and computer is essential for smooth performance by a telerobotic system. ESG (Expert Script Generator) is a software package that automatically generates high-level task objective commands from the NASA Intelligent Systems Research Lab's (ISRL's) complex menu-driven language. ESG reduces errors and makes the telerobotics lab accessible to researchers who are not familiar with the comprehensive language developed by ISRL for interacting with the various systems of the ISRL testbed. ESG incorporates expert system technology to capture the typical rules of operation that a skilled operator would use. The result is an operator interface which optimizes the system's capability to perform a task remotely in a hazardous environment, in a timely manner, and without undue stress to the operator, while minimizing the chance for operator errors that may damage equipment. The intricate menu-driven command interface which provides for various control modes of both manipulators and their associated sensors in the TeleRobotic System Simulation (TRSS) has a syntax which is both irregular and verbose. ESG eliminates the following two problems with this command "language": 1) knowing the correct command sequence to accomplish a task, and 2) inputting a known command sequence without typos and other errors. ESG serves as an additional layer of interface, working in conjunction with the menu command processor, not supplanting it. By specifying task-level commands, such as GRASP, CONNECT, etc., ESG will generate the appropriate menu elements to accomplish the task. These elements will be collected in a script file which can then be executed by the ISRL menu command processor. In addition, the operator can extend the list of task-level commands to include customized tasks composed of sub-task commands. This mechanism gives the operator the ability to build a task-hierarchy tree of increasingly powerful commands. ESG also provides automatic regeneration of scripts based on system knowledge of telerobotic environment updates. The commands generated by ESG may be displayed at the terminal screen and/or stored. ESG is implemented as a rule-based expert system written in CLIPS (C Language Integrated Production System). The system consists of a knowledge-base of task heuristics, a static (unchanged during execution) database which describes the physical features of objects, and a dynamic (may change as a result of task achievement) database which maintains changes in the environment. Capabilites are provided for adding new environmental objects and for modifying existing objects and configuration data. Options are available for interactively viewing both the static and dynamic attribute values of database items. Execution of the ESG may be suspended to allow access to system-level functions. ESG was implemented on a VAX 11/780 with the VMS 4.7 operating system using a VT100 compatible terminal. Its source code is 47% CLIPS and 53% C-language, with a memory requirement of approximately 205 KB. The program was developed in 1988.

Description: The expert system called EXADS was developed to aid users of the Automated Design Synthesis (ADS) general purpose optimization program. Because of the general purpose nature of ADS, it is difficult for a nonexpert to select the best choice of strategy, optimizer, and one-dimensional search options from the one hundred or so combinations that are available. EXADS aids engineers in determining the best combination based on their knowledge of the problem and the expert knowledge previously stored by experts who developed ADS. EXADS is a customized application of the AESOP artificial intelligence program (the general version of AESOP is available separately from COSMIC. The ADS program is also available from COSMIC.) The expert system consists of two main components. The knowledge base contains about 200 rules and is divided into three categories: constrained, unconstrained, and constrained treated as unconstrained. The EXADS inference engine is rule-based and makes decisions about a particular situation using hypotheses (potential solutions), rules, and answers to questions drawn from the rule base. EXADS is backward-chaining, that is, it works from hypothesis to facts. The rule base was compiled from sources such as literature searches, ADS documentation, and engineer surveys. EXADS will accept answers such as yes, no, maybe, likely, and don't know, or a certainty factor ranging from 0 to 10. When any hypothesis reaches a confidence level of 90% or more, it is deemed as the best choice and displayed to the user. If no hypothesis is confirmed, the user can examine explanations of why the hypotheses failed to reach the 90% level. The IBM PC version of EXADS is written in IQ-LISP for execution under DOS 2.0 or higher with a central memory requirement of approximately 512K of 8 bit bytes. This program was developed in 1986.

Description: CFORM was developed by the Kennedy Space Center Robotics Lab to assist in linear control system design and analysis using closed form and transient response mechanisms. The program computes the closed form solution and transient response of a linear (constant coefficient) differential equation. CFORM allows a choice of three input functions: the Unit Step (a unit change in displacement); the Ramp function (step velocity); and the Parabolic function (step acceleration). It is only accurate in cases where the differential equation has distinct roots, and does not handle the case for roots at the origin (s=0). Initial conditions must be zero. Differential equations may be input to CFORM in two forms - polynomial and product of factors. In some linear control analyses, it may be more appropriate to use a related program, Linear Control System Design and Analysis (KSC-11376), which uses root locus and frequency response methods. CFORM was written in VAX FORTRAN for a VAX 11/780 under VAX VMS 4.7. It has a central memory requirement of 30K. CFORM was developed in 1987.

Description: This control theory design package, called Optimal Regulator Algorithms for the Control of Linear Systems (ORACLS), was developed to aid in the design of controllers and optimal filters for systems which can be modeled by linear, time-invariant differential and difference equations. Optimal linear quadratic regulator theory, currently referred to as the Linear-Quadratic-Gaussian (LQG) problem, has become the most widely accepted method of determining optimal control policy. Within this theory, the infinite duration time-invariant problems, which lead to constant gain feedback control laws and constant Kalman-Bucy filter gains for reconstruction of the system state, exhibit high tractability and potential ease of implementation. A variety of new and efficient methods in the field of numerical linear algebra have been combined into the ORACLS program, which provides for the solution to time-invariant continuous or discrete LQG problems. The ORACLS package is particularly attractive to the control system designer because it provides a rigorous tool for dealing with multi-input and multi-output dynamic systems in both continuous and discrete form. The ORACLS programming system is a collection of subroutines which can be used to formulate, manipulate, and solve various LQG design problems. The ORACLS program is constructed in a manner which permits the user to maintain considerable flexibility at each operational state. This flexibility is accomplished by providing primary operations, analysis of linear time-invariant systems, and control synthesis based on LQG methodology. The input-output routines handle the reading and writing of numerical matrices, printing heading information, and accumulating output information. The basic vector-matrix operations include addition, subtraction, multiplication, equation, norm construction, tracing, transposition, scaling, juxtaposition, and construction of null and identity matrices. The analysis routines provide for the following computations: the eigenvalues and eigenvectors of real matrices; the relative stability of a given matrix; matrix factorization; the solution of linear constant coefficient vector-matrix algebraic equations; the controllability properties of a linear time-invariant system; the steady-state covariance matrix of an open-loop stable system forced by white noise; and the transient response of continuous linear time-invariant systems. The control law design routines of ORACLS implement some of the more common techniques of time-invariant LQG methodology. For the finite-duration optimal linear regulator problem with noise-free measurements, continuous dynamics, and integral performance index, a routine is provided which implements the negative exponential method for finding both the transient and steady-state solutions to the matrix Riccati equation. For the discrete version of this problem, the method of backwards differencing is applied to find the solutions to the discrete Riccati equation. A routine is also included to solve the steady-state Riccati equation by the Newton algorithms described by Klein, for continuous problems, and by Hewer, for discrete problems. Another routine calculates the prefilter gain to eliminate control state cross-product terms in the quadratic performance index and the weighting matrices for the sampled data optimal linear regulator problem. For cases with measurement noise, duality theory and optimal regulator algorithms are used to calculate solutions to the continuous and discrete Kalman-Bucy filter problems. Finally, routines are included to implement the continuous and discrete forms of the explicit (model-in-the-system) and implicit (model-in-the-performance-index) model following theory. These routines generate linear control laws which cause the output of a dynamic time-invariant system to track the output of a prescribed model. In order to apply ORACLS, the user must write an executive (driver) program which inputs the problem coefficients, formulates and selects the routines to be used to solve the problem, and specifies the desired output. There are three versions of ORACLS source code available for implementation: CDC, IBM, and DEC. The CDC version has been implemented on a CDC 6000 series computer with a central memory of approximately 13K (octal) of 60 bit words. The CDC version is written in FORTRAN IV, was developed in 1978, and last updated in 1989. The IBM version has been implemented on an IBM 370 series computer with a central memory requirement of approximately 300K of 8 bit bytes. The IBM version is written in FORTRAN IV and was generated in 1981. The DEC version has been implemented on a VAX series computer operating under VMS. The VAX version is written in FORTRAN 77 and was generated in 1986.

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 paper presents the synthesis of neural network based feedback laws for dynamic systems using the computed optimal and time histories of the state and control variables. The efficacy of the proposed approach has been successfully demonstrated on a minimum time orbit injection problem. If the method is found to be effective to real life problems with many state and control variables, it can used for a variety of guidance and control problems.

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.

Description: Sea surface temperature (SST) in roughly 50% of the tropical Pacific Ocean is warm enough (SST greater than 300 K) to permit deep convection. This paper examines the effects of deep convection on the climatological mean vertical distributions of water vapor and its greenhouse effect over such warm oceans. The study, which uses a combination of satellite radiation budget observations, atmospheric soundings deployed from ships, and radiation model calculations, also examines the link between SST, vertical distribution of water vapor, and its greenhouse effect in the tropical oceans. Since the focus of the study is on the radiative effects of water vapor, the radiation model calculations do not include the effects of clouds. The data are grouped into nonconvective and convective categories using SST as an index for convective activity. On average, convective regions are more humid, trap significantly more longwave radiation, and emit more radiation to the sea surface. The greenhouse effect in regions of convection operates as per classical ideas, that is, as the SST increases, the atmosphere traps the excess longwave energy emitted by the surface and reradiates it locally back to the ocean surface. The important departure from the classical picture is that the net (up minus down) fluxes at the surface and at the top of the atmosphere decrease with an increase in SST; that is, the surface and the surface-troposphere column lose the ability to radiate the excess energy to space. The cause of this super greenhouse effect at the surface is the rapid increase in the lower-troposphere humidity with SST; that of the column is due to a combination of increase in humidity in the entire column and increase in the lapse rate within the lower troposphere. The increase in the vertical distribution of humidity far exceeds that which can be attributed to the temperature dependence of saturation vapor pressure; that is, the tropospheric relative humidity is larger in convective regions. The positive coupling between SST and the radiative warming of the surface by the water vapor greenhouse effect is also shown to exist on interannual time scales.

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: MacMultiview is an interactive tool for the Macintosh II family which allows one to display and make computations utilizing polarimetric radar data collected by the Jet Propulsion Laboratory's imaging SAR (synthetic aperture radar) polarimeter system. The system 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. MacMultiview provides two basic functions: computation of synthesized polarimetric images and computation of polarization signatures. The radar data can be used to compute a variety of images. The total power image displays the sum of the polarized and unpolarized components of the backscatter for each pixel. The magnitude/phase difference image displays the HH (horizontal transmit and horizontal receive polarization) to VV (vertical transmit and vertical receive polarization) phase difference using color. Magnitude is displayed using intensity. The user may also select any transmit and receive polarization combination from which an image is synthesized. This image displays the backscatter which would have been observed had the sensor been configured using the selected transmit and receive polarizations. MacMultiview can also be used to compute polarization signatures, three dimensional plots of backscatter versus transmit and receive polarizations. The standard co-polarization signatures (transmit and receive polarizations are the same) and cross-polarization signatures (transmit and receive polarizations are orthogonal) can be plotted for any rectangular subset of pixels within a radar data set. In addition, the ratio of co- and cross-polarization signatures computed from different subsets within the same data set can also be computed. Computed images can be saved in a variety of formats: byte format (headerless format which saves the image as a string of byte values), MacMultiview (a byte image preceded by an ASCII header), and PICT2 format (standard format readable by MacMultiview and other image processing programs for the Macintosh). Images can also be printed on PostScript output devices. Polarization signatures can be saved in either a PICT format or as a text file containing PostScript commands and printed on any QuickDraw output device. The associated Stokes matrices can be stored in a text file. MacMultiview is written in C-language for Macintosh II series computers. MacMultiview will only run on Macintosh II series computers with 8-bit video displays (gray shades or color). The program also requires a minimum configuration of System 6.0, Finder 6.1, and 1Mb of RAM. MacMultiview is NOT compatible with System 7.0. It requires 32-Bit QuickDraw. Note: MacMultiview may not be fully compatible with preliminary versions of 32-Bit QuickDraw. Macintosh Programmer's Workshop and Macintosh Programmer's Workshop C (version 3.0) are required for recompiling and relinking. The standard distribution medium for this package is a set of three 800K 3.5 inch diskettes in Macintosh format. This program was developed in 1989 and updated in 1991. MacMultiview is a copyrighted work with all copyright vested in NASA. QuickDraw, Finder, Macintosh, and System 7 are trademarks of Apple Computer, Inc.

Description: The accurate partitioning of available energy into sensible and latent heat flux is crucial to the understanding of surface atmosphere interactions. This issue is more complicated in arid and semi arid regions where the relative contribution to surface fluxes from the soil and vegetation may vary significantly throughout the day and throughout the season. A three component model to estimate sensible heat flux over heterogeneous surfaces is presented. The surface was represented with two adjacent compartments. The first compartment is made up of two components, shrubs and shaded soil, the second of open 'illuminated' soil. Data collected at two different sites in Nevada (U.S.) during the Summers of 1991 and 1992 were used to evaluate model performance. The results show that the present model is sufficiently general to yield satisfactory results for both sites.

Description: Since September 1991, the ERS-1 SAR (Synthetic Aperture Radar) has collected approximately 170 frames of ocean imagery on 28 passes over the western Gulf Stream in support of ESA experiment US8-2c. SAR signatures of the north wall are seen on nearly all passes, with modulation depth varying from 3% to 35% for 100 m samples. Many small and mesoscale circulation features associated with the Stream are evident. The detailed form of the signature varies considerably, however. Narrow bright or dark linear features appear to follow streamlines, and on one occasion a strong dark line was associated with an in situ measurement of a sharp current shear. Similarly, larger spatial scale changes in backscatter over the Stream were associated with in situ measurements of atmospheric stability transitions. Physical explanations for the narrow features are not so obvious. However, the accumulation of surfactants along converging current boundaries or local short wave straining and breaking appear plausible. These preliminary results strongly suggest that a wide swath (approximately 500 km) SAR with at least 100 m resolution would be a useful adjunct to existing satellite AVHRR (Advanced Very High Resolution Radiometer) imagery.

Description: ERS-1 SAR (Synthetic Aperture Radar) data were used to specify Winter physical processes on the surface of the Greenland Sea, and SSM/I (Special Sensor Microwave/Imager) data were used to characterize the regional behavior of the ice cover. Examination of the SSM/I data indicated that the convective water was likely to be confined to small (less than 100 km) domains near the ice edge. Using ERS-1 SAR data crossing the ice edge, it was possible to identify ice edge features that are very similar to modeled plumes; they have a diameter of about 100 m, a spacing of about 300 m, and cover an area about 20 by 90 km. This plume observation is the first such identification of convection. The plumes are topped by ice, and the return-water areas are open, indicating a freshened top layer of the sea and emphasizing the importance of this layer to convection.

Description: Global water isotope geochemisty during the last glacial maximum (LGM) is simulated with an 8 deg x 10 deg atmospheric general circulation model (GCM). The simulation results suggest that the spatial delta O-18/temperature relationships observed for the present day and LGM climates are very similar. Furthermore, the temporal delta O-18/temperature relationship is similar to the present-day spatial relationship in regions for which the LGM/present-day temperature change is significant. This helps justify the standard practice of applying the latter to the interpretation of paleodata, despite the possible influence of other factors, such as changes in the evaportive sources of precipitation or in the seasonality of precipitation. The model suggests, for example, that temperature shifts inferred from ice core data may differ from the true shifts by only about 30%.

Description: There is renewed interest in using optical remote sensing techniques for measuring and monitoring airborne toxins. Many instruments have been used to measure traces of atmospheric molecular species. Some optical remote sensing technology is in use for measuring or monitoring toxic or pollutant gases. The examples in this article provide a number of lessons that may be of use in the new field of optical remote sensing of toxic gases. Topics covered include the following: atmospheric properties, Fourier transform spectrometry, gas-filter correlation, differential absorption laser systems, aircraft missions, laser long-path instruments, differential optical absorption spectrometry, and ozone measurement.

Description: Satellite observations indicate that the humidity of the upper troposphere is higher in summer than in winter. We use general circulation model (GCM) simulations to explore the processes that maintain upper troposphere water vapor and determine its seasonal cycle. In the subtropics, drying by Hadley cell subsidence and stratiform condensation is offset primarily by moistening by eddies, with moist convection playing a minor role. Elsewhere, both mean meridional circulation and eddies moisten the upper troposphere and are balanced primarily by stratiform condensation drying. The effect of the seasonal shift of the Hadley cell is limited to latitudes equatorward of 30 deg. At higher latitudes where the largest observed summer moistening occurs, eddy moisture fluxes are primarily responsible despite the eddies being weaker in summer than winter. The same mechanism causes upper level humidity to increase in GCM climate warming simulations. The observed seasonal variation may thus be a good proxy for decadal climate change. This suggests that upper troposphere water vapor feedback is positive at all latitudes, consistent with GCM predictions.

Description: The response of radiation budgets to changes in water vapor and clouds in an El Nino episode is investigated using the analyzed sea surface temperature (SST) and satellite-derived clouds and the earth radiation budgets for the tropical Pacific (30 deg N-30 deg S, 100 deg E-100 deg W). Analyses are performed for April 1985 and April 1987. The former is a non-El Nino year and the latter is an El Nino year. Compared to April 1985, when the SST over the central and eastern equatorial Pacific is approximately 2 C lower, the high-level cloudiness in April 1987 increases in the central and eastern equatorial Pacific. Corresponding to the increase in cloudiness, the outgoing longwave radiation and the net downward solar radiation at the top of the atmosphere decrease. The patterns of these changes are reversed in the western tropical Pacific and the Northern Hemispheric (NH) subsidence region centered at approximately 20 deg N, indicating an eastward shift of the convection center from the maritime continents to the central equatorial Pacific and a strengthened NH Hadley circulation. The earth-atmosphere system in the region receives less radiative energy by 4 W/sq m in the warmer month of April 1987 than in the month of April 1985, which is primarily caused by a reduced atmospheric clear sky greenhouse effect in the NH tropical Pacific in April 1987. Clouds have strong effects on both the IR and solar radiation, but the net effect on the radiation budget at the top of the atmopshere changes only slightly between April 1985 and April 1987. The results are consistent with Lindzen's hypothesis that reduced upper-tropospheric water vapor in the vicinity of the enhanced convection region produces cooling that counteracts warming in the Tropics.

Description: The impact of time compositing on infrared profiling from geosynchronous orbit was evaluated for two convective outbreak cases. Time compositing is the accumulation of the data from several successive images taken at short intervals to provide a single field of measurements with the temporal resolution equal to the time to take all of the images. This is especially effective when the variability of the measurements is slow compared to the image interval. Time compositing should be able to reduce the interference of clouds for infrared measurments since clouds move and change. The convective outbreak cases were on 4 and 21 May 1990 over the eastern Midwest and southeastern United States, respectively. Geostationary Operational Environmental (GOES) Satellite imagery was used to outline clear areas at hourly intervals by two independent analysts. Time compositing was done every 3 h (1330-1530 UTC; 1630-1830 UTC) and over the full 5-h period. For both cases, a significant increase in coverage was measured with each 3-h compositing (about a factor of 2) and a further increase over the full period (approximately a factor of 3). The increase was especially useful in areas of broken cloud cover where large gaps between potential profiling areas on each image were reduced. To provide information on measurement variability over local areas, the regions where the clear-area analyses were done were subdivided into 0.5 deg latitude-longitude boxes, and if some portion of each box was clear, it was assumed that at least one profile could be obtained within the box. In the largest clear areas, at least some portion was clear every hour. Even in the cloudier regions, multiple clear looks possible during the entire period.

Description: National Data Buoy Center (NDBC) observations from the western Gulf of Mexico Provided several early indications of more rapid cyclogenesis on 12 March 1993 than was forecast by numerical weather prediction (NWP) models. Observations demostrated a closed circulation with gale-force winds shortly after the storm entered the Gulf of Mexico. Pressure measurements at two buoys off the Texas coast were 4-6 hPa lower than the 12-h NWP forecasts, a significant forecast error. Observations from NDBC's moored buoys and Coastal-Marine Automated Network stations revealed that the developing storm was significantly deeper than was indicated on the National Meteorological Center's automated surface analyses. Ocean wave observations reveal some of the steepest waves NDBC has ever measured, indicating phenomenonal wave growth and a high potential for damage to vessels and structures. A warm eddy caused sea surface temperatures (SST) to be several degrees above normal under the track of the storm, creating a strong SST gradient to the north. This provided ample energy and strengthened the baroclinity. NDBC observations showed the eddy to be somewhat larger and warmer than indicated by the most recent National Hurricane Center analysis. This event demonstrates the tremendous value of NDBC marine observations for accurately detecting the occurrence and strength of coastal cyclogenesis events.

Description: A global mountain wave parameterization for prediction of wave-related displacements and turbulence is described. The parameterization is used with input from National Meteorological Center (NMC) analyses of wind and temperature to examine small-scale disturbances encountered by the National Aeronautics and Space Administration high-altitude ER-2 during the Second Airborne Arctic Stratosphere Experiment (AASE-II). The magnitude and location of observed large wave events are well reproduced. A strong correlation is suggested between patches of moderate turbulence encountered by the ER-2 and locations where breaking mountain waves are predicted by the parameterization. These facts suggest that useful forecasts of global mountain wave activity, including wave-related clear-air turbulence, can be made quickly and inexpensively using our mountain wave parameterization with input from current numerical forecast models.

Description: The first Special Sensor Microwave Imager (SSMI) was launched on the Defense Meteorological Satellite Program (DMSP) F8 spacecraft in July 1987, and wind speed was no longer retrieved after December 1991. A second SSMI was launched on DMSP F10 in December 1990. Interpretation of the 1987-1993 (or longer) SSMI wind speed time series is dependent upon the space and time characteristics of the differences between F8 and F10 SSMI measurements. The 10 deg-zonal averaged monthly mean F8-F10 wind speed difference was negative (positive) for wind speeds less (greater) than 7.9 m/s, reaching -0.43 (0.32) m/s at 5(10) m/s. Between 60 deg S and 60 deg N the 10 deg-zonal averaged monthly mean F8-F10 wind speed bias was greater than +/- 0.5 m/s on several occasions. From 60 deg S - 60 deg N the 1991 average value of the monthly mean root-mean-square difference between daily F8 and F10 wind speeds in 10 deg-longitudinal bands was 2.0 m/s.In the 60 deg S - 60 deg N region, about 50% of the daily F8 and F10 wind speed differences was caused by measurement non-simultaneity and about 50% of the difference was attributed to other factors, such as instrument noise and the different azimuthal orientations of each SSMI.

Description: The evolution of Southern Hemisphere air masses observed by the Halogen Occultation Experiment (HALOE) during September 21 through October 15, 1992, is investigated using isentropic trajectories computed from United Kingdom Meteorological Office (UKMO) assimilated winds and temperatures. Maps of constituent concentrations are obtained by accumulation of air masses from previous HALOE occultations. Lagged correlations between initial and subsequent HALOE observations of the same air mass are used to validate the air mass trajectories. High correlations are found for lag times as large as 10 days. Frequency distributions of the air mass constituent concentrations are used to examine constituent distributions in and around the Southern Hemisphere polar vortex.

Description: A new model is presented for feedback between rotational and climatic variations, operative on time scales of 10(exp 4) - 10(exp 7) years. Due to the combined effect of planetary perturbations to the Earth's orbit plane and luni-solar torques on the oblate figure of the Earth, the obliquity varies by approximately 1 deg on a 4 x 10(exp 4) year time scale. Associated changes in the seasonal and latitudinal pattern of incident solar radiation cause major glaciations. Mass transport from the oceans to the polar ice sheets during these glaciations can change the gravitational oblateness of the Earth by amounts approaching 1%. As the rate of spin axis precession is directly proportional to the oblateness, the climatically forced mass transport can by dynamically significant. A simple parameterization of this coupled orbital-rotational-climatic system suggests that there is a strong tendency for the system to evolve away from climatically sensitive values of the obliquity. This may explain the mid-Pleistocene transition from an obliquity dominated regime to the present regime in which most climatic variability is concentrated at longer (10(exp 5) year) periods.

Description: Potential vorticity (PV) gradients defining the lower stratospheric vortex during the 1992-1993 winter were anomalously strong and persistent compared to those during the last 16 Arctic winters. For approximately equal to 3 months PV gradients were closer to typical Antarctic values than to most Arctic values. Air motion diagnostics computed for 3-dimensional air parcels confirm that the 1992-1993 Arctic lower stratospheric vortex was substantially more isolated than is typical. Such isolation will delay and reduce the export of the higher ozone typical of the winter lower stratospheric vortex to mid-latitudes. This may have contributed to the record-low total ozone amounts observed in northern mid-latitudes in 1993.

Description: Lagrangian material line simulations are performed using U.K. Meteorological Office simulated winds and temperatures to examine mixing processes in the middle- and lower-stratospheric polar night jet during the 1992 Southern Hemisphere spring and Northern Hemisphere winter. The Lagrangian simulations are undertaken to provide insight into the effects of mixing within the polar night jet on observations of the polar vortex made by instruments onboard the Upper Atmosphere Research Satellite (UARS) during these periods. A moderate to strong kinematic barrier to large-scale isentropic exchange, similar to the barrier identified in General Circulation Model (GCM) simulations, is identified during both of these periods. Characteristic timescales for mixing by large-scale isentropic motions within the polar night jet range from 20 days in the Southern Hemisphere lower stratosphere to years in the Northern Hemisphere middle stratosphere. The long mixing timescales found in the Northern Hemisphere polar night jet do not persist. Instead, the Northern Hemisphere kinematic barriers are broken down as part of the large-scale stratospheric response to a strong tropospheric blocking event. A series of Lagrangian experiments are conducted to investigate the sensitivity of the kinematic barrier to diabatic effects and to small-scale inertial gravity wave motions. Differential diabatic descent is found to have a significant impact on mixing processes within the Southern Hemisphere middle-stratospheric jet core. The interaction between small-scale displacements by idealized, inertial gravity waves and the large-scale flow is found to have a significant impact on mixing within the polar night jet in both hemispheres. These sensitivity experiments suggest that scales of motion that are unresolved in global assimilated datasets may contribute to mass exchange across the kinematic barrier to large-scale isentropic motion.

Description: The degree to which the Southern Hemisphere polar vortex is isolated against horizontal (isentropic) mixing is investigated using data from the Halogen Occultation Experiment (HALOE), U.K. Meteorological Office (UKMO) potential vorticity (PV), and contour advection diagnostics. Measurements of methane and water vapor taken by HALOE during a disturbed period in the Southern Hemisphere springtime (21 September - 15 October 1992) are interpreted in light of the prevailing synoptic meteorology. Daily fields of winds and PV are shown to be essential in the interpretation of the data. A climatological high pressure region is responsible for a distorted vortex, and a substantial 'vortex stripping' event is present, associated with the early stages of vortex breakdown. This leads to significant temporal, zonal, and altitudinal variations in the distribution of tracers. The authors point out the difficulties this presents for the interpretation of solar occultation data, especially with regard to the use of zonal average time series. Longitude-height methane distributions from two days during the period are examined. Both days show substantial variations in abundance around a latitude circle. In particular, the authors investigate HALOE measurements at 77 deg S on 15 October 1992, which indicate an abundance of methane in the height region 600-2000 K (approximately 30-1 mb) that is more typical of midlatitude air. Similar distributions, observed in the 1991 HALOE data, have previously been interpreted as evidence for the penetration of midlatitude air into the vortex. Gradients of potential vorticity and contour advection diagnositcs are employed to examine whether the UKMO winds are consistent with this hypothesis in 1992. Although midlatitude air is able to penetrate poleward of the main jet core by advection processes alone, an essentially intact inner core of vortex air remains, which does not mix to any great extent with air from lower latitudes. The authors show that the high-latitude HALOE abundances that are typical of midlatitude air were observed in a region of extensive filamentation and mixing, rather than within the inner, more isolated, core.

Description: The distribution of dehydrated air in the middle and lower stratosphere during the 1992 Southern Hemisphere spring is investigated using Halogen Occultation Experiment (HALOE) observations and trajectory techniques. Comparisons between previously published Version 9 and the improved Version 16 retrievals on the 700-K isentropic surface show very slight (0.05 ppmv) increases in Version 16 CH4 relative to Version 9 within the polar vortex. Version 16 H2O mixing ratios show a reduction of 0.5 ppmv relative to Version 9 within the polar night jet and a reduction of nearly 1.0 ppmv in middle latitudes when compared to Version 9. The version 16 HALOE retrievals show low mixing ratios of total hydrogen (2CH4 + H2O) within the polar vortex on both 700 and 425 K isentropic surfaces relative to typical middle-stratospheric 2CH4 + H2O mixing ratios. The low 2CH4 + H2O mixing ratios are associated with dehydration. Slight reductions in total hydrogen, relative to typical middle-stratospheric values, are found at these levels throughout the Southern Hemisphere during this period. Trajectory calculations show that middle-latitude air masses are composed of a mixture of air from within the polar night jet and air from middle latitudes. A strong kinematic barrier to large-scale exchange is found on the poleward flank of the polar night jet at 700 K. A much weaker kinematic barrier is found at 425 K. The impact of the finite tangent pathlength of the HALOE measurements is investigated using an idealized tracer distribution. This experiment suggests that HALOE should be able to resolve the kinematic barrier, if it exists.

Description: An aircraft locally distorts the ambient thundercloud electric field. In order to determine the field in the absence of the aircraft, an aircraft calibration is required. In this work a matrix inversion method is introduced for calibrating an aircraft equipped with four or more electric field sensors and a high-voltage corona point that is capable of charging the aircraft. An analytic, closed form solution for the estimate of a (3 x 3) aircraft calibration matrix is derived, and an absolute calibration experiment is used to improve the relative magnitudes of the elements of this matrix. To demonstrate the calibration procedure, we analyze actual calibration date derived from a Lear jet 28/29 that was equipped with five shutter-type field mill sensors (each with sensitivities of better than 1 V/m) located on the top, bottom, port, starboard, and aft positions. As a test of the calibration method, we analyze computer-simulated calibration data (derived from known aircraft and ambient fields) and explicitly determine the errors involved in deriving the variety of calibration matrices. We extend our formalism to arrive at an analytic solution for the ambient field, and again carry all errors explicitly.

Description: This paper is on the control of nonlinear-nonstationary vibration of a frame-stringer structure resulting from high levels of excitation from a nearby supersonic jet exhaust. The structure exhibits periodic, chaotic, or random behaviors when forced by high-intensity sound from a supersonic jet exhaust with 'shock' loading superimposed on a broadband response. The time history of the pressure, showing the rotation and flapping of the shock structure in the jet column due to large-scale instabilities, indicates that the response is not only nonlinear but also nonstationary. The acoustic pressure radiated by the structure also contains shocks and the formation of harmonics with distance. Control of the structural response is achieved by actively forcing the structure with an actuator at the shock oscillation frequency whose amplitude is locked into a self-control cycle. Results show that the peak power level is reduced by a factor of 63, or 18 dB. As a result, new broadband components emerge with at least four harmonics. At accelerating and decelerating supersonic speeds, the exhaust from the jet induces higher transient loading on the nearby flexible structure due to the occurrence of multiple shock from the jet.

Description: This paper presents a novel method to design decentralized controllers for large complex flexible structures by using the idea of joint decoupling. Decoupling of joint degrees of freedom from the interior degrees of freedom is achieved by setting the joint actuator commands to cancel the internal forces exerting on the joint degrees of freedom. By doing so, the interactions between substructures are eliminated. The global structure control design problem is then decomposed into several substructure control design problems. Control commands for interior actuators are set to be localized state feedback using decentralized observers for state estimation. The proposed decentralized controllers can operate successfully at the individual substructure level as well as at the global structure level. Not only control design but also control implementation is decentralized. A two-component mass-spring-damper system is used as an example to demonstrate the proposed method.

Description: Simulations of natural variability by two general circulation models (GCMs) are examined. One GCM is a sector model, allowing relatively rapid integration without simplification of the model physics, which would potentially exclude mechanisms of variability. Two mechanisms are found in which tropical surface temperature and sea surface temperature (SST) vary on interannual and longer timescales. Both are related to changes in cloud cover that modulate SST through the surface radiative flux. Over the equatorial ocean, SST and surface temperature vary on an interannual timescale, which is determined by the magnitude of the associated cloud cover anomalies. Over the subtropical ocean, variations in low cloud cover drive SST variations. In the sector model, the variability has no preferred timescale, but instead is characterized by a 'red' spectrum with increasing power at longer periods. In the terrestrial GCM, SST variability associated with low cloud anomalies has a decadal timescale and is the dominant form of global temperature variability. Both GCMs are coupled to a mixed layer ocean model, where dynamical heat transports are prescribed, thus filtering out El Nino-Southern Oscillation (ENSO) and thermohaline circulation variability. The occurrence of variability in the absence of dynamical ocean feedbacks suggests that climatic variability on long timescales can arise from atmospheric processes alone.

Description: The Goddard Institute for Space Studies (GISS) general circulation model (GCM) is used to study the possible implications of past and future climate change on global lightning frequencies. Two climate change experiments were conducted: one for a 2 x CO2 climate (representing a 4.2 degs C global warming) and one for a 2% decrease in the solar constant (representing a 5.9 degs C global cooling). The results suggest at 30% increase in global lightning activity for the warmer climate and a 24% decrease in global lightning activity for the colder climate. This implies an approximate 5-6% change in global lightning frequencies for every 1 degs C global warming/cooling. Both intracloud and cloud-to-ground frequencies are modeled, with cloud-to-ground lightning frequencies showing larger sensitivity to climate change than intracloud frequencies. The magnitude of the modeled lightning changes depends on season, location, and even time of day.

Description: We have used a one-dimensional model of cirrus formation to study the development of cirrus clouds during the 1986 First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment (FIRE) intensive field observations (IFO). The cirrus model includes microphysical, dynamical, and radiative processes. Sulfate aerosols, solution drops, ice crystals, and water vapor are all treated as interactive elements in the model. Ice crystal size distributions are fully resolved based on calculations of homogeneous freezing nucleation, growth by water vapor deposition, evaporation, coagulation, and vertical transport. We have focused on the cirrus observed on November 1, 1986. Vertical wind speed for the one-dimensional simulation is taken from a mesoscale model simulation for the appropriate time period. The mesoscale model simulation suggested that strong upward motions over Wyoming and subsequent horizontal transport of upper level moisture were responsible for the cirrus observed over Wisconsin on this date. We assumed that our one-dimensional model could be used to represent a vertical column moving from Wyoming to Wisconsin over a period of several hours. Ice crystal nucleation occurs in our model in the 8 to 10-km region as a result of the strong updrafts (and cooling) early in the simulation. Growth, coagulation, and sedimentation of these ice crystals result in a broad cloud region (5-10 km thick) with an optical depth of 1-2 after a few hours, in agreement with the FIRE measurements. Comparison with aircraft microphysical measurements made over Wisconsin indicates that the simulation generated reasonable ice water content, but the predicted ice number densities are too low, especially for radii less than about 50 microns. Sensitivity tests suggest that better agreement between simulated and observed microphysical properties is achieved if the nucleation rate is higher or stronger vertical mixing (perhaps associated with multidimensional motions) is present.

Description: In this paper the influences of recent technology developments in the areas of lasers, detectors, andoptical filters of a differential absorption lidar (DIAL) system on the measurenent of tropospheric water vapor (H2O) profiles are discussed. The lidar parameters selected are based upon a diode-seeded Ti:sapphire laser that is locked to an H2O line in the 820- or 930-nm band of H2O. To assess the influence of the mode of deployment on the measurement of tropospheric H2O, DIAL performance is evaluated for operation from a medium-altitude (12 km) aircraft, the ground, and space-based systems. It is found that incorporation of these developments could greatly enhance DIAL measurement capability.

Description: Ozone measurements were obtained between the surface and the 6-km altitude on aircraft flights over central and eastern Canada during the summer 1990 NASA Global Tropospheric Experiment Arctic Boundary Layer Expedition (GTE/ABLE 3B). Tropospheric O3 budgets for these regions were observed to be highly variable and significantly impacted by long-range transport and regional scale air mass modification processes. For example, integrated O3 abundance below 5-km altitude averaged 40% and 30% greater in air masses influenced by anthropogenic sources and biomass burning, respectively, than in background (polar) air. Conversely, aged air transported from subtropical areas of the Pacific at times reduced O3 abundance in this height interval by up to 20%. Though intrusion of anthropogenic air was infrequent during the experiment period, the influence of biomass-burning emissions was particularly notable as two thirds of the flights sampled air influenced by plumes from fires burning in Alaska and western Canada. The impinging pollution, both natural and anthropogenic, not only elevated O3 levels directly but also was a source of reactive nitrogen (and nonmethane hydrocarbons) which generally increases the tropospheric lifetime of O3 via moderation of photochemical destruction rates.

Description: NASA's Atmospheric Boundary Layer Experiment conducted during the summer of 1990 focused on the distribution of trace species in central and northeastern Canada (altitudes less than 6 km) and the importance of surface sources/sinks, local emissions, distant transport, tropospheric/stratospheric exchange. Aircraft flights were based from North Bay, Ontario, and Goose Bay, Labrador, Canada. As part of the aircraft measurements, aerosol number density (0.2- to 3-micrometers diameter) was measured using an optical laser technique. Results show that summertime aerosol budgets of central and northeastern Canada can be significantly impacted by the transport of pollutants from distant source regions. Biomass burning in Alaska and western and central Canada exerts major influences on regional aerosol budgets. Urban emissions transported from the U.S./Canadian border regions are also important. Aerosol enhancements (mixed layer and free troposphere) were most prevalent in air with carbon monoxide mixing ratios greater than 110 parts per billion by volume (ppbv). When data were grouped as to the source of the air (5-day back trajectories) either north or south of the polar jet, aerosol number density in the mixed layer showed a tendency to be enhanced for air south of the jet relative to north of the jet. However, this difference was not observed for measurements at the higher altitudes (4 to 6 km). For some flights, mixed layer aerosol number densities were greater than 100 higher than free-tropospheric values (3- to 6-km altitude). The majority of the observed mixed layer enhancement was associated with transport of effluent-rich air into the Canadian regions. Aerosol emissions from natural Canadian ecosystems were relatively small when compared to transport.

Description: Based on point measurements of methane flux from wetlands in the boreal and subarctic regions, northern wetlands are a major source of atmospheric methane. However, measurements have not been carried out in large continuous peatlands such as the Hudson Bay Lowland (HBL) (320,000 sq km) and the Western Siberian lowland (540,000 sq km), which together account for over 30% of the wetlands north of 40 deg N. To determine the role the Hudson Bay Lowland as a source of atmospheric methane, fluxes were measured by enclosures throughout the 1990 snow-free period in all the major wetland types and also by an aircraft in July. Two detailed survey areas were investigated: one (approximately 900 sq km) was in the high subarctic region of the northern lowland and the second area (approximately 4,800 sq km) straddled the Low Subarctic and High Boreal regions of the southern lowland. The fluxes were integrated over the study period to produce annual methane emissions for each wetland type. The fluxes were then weighted by the area of 16 different habitats for the southern area and 5 habitats for the northern area, as determined from Landsat thematic mapper to yield an annual habitat-weighted emission. On a per unit area basis, 1.31 +/- 0.11 and 2.79 +/- 0.39 g CH4 m(exp -2)/yr were emitted from the southern and northern survey areas, respectively. The extrapolated enclosure estimates for a 3-week period in July were compared to within 10% of the flux derived by airborne eddy correlation measurements made during the same period. The aircraft mean flux of 10 +/- 9 mg CH4 m(exp -2)/d was not statistically different from the extrapolated mean flux of 20 +/- 16 mg CH4 m(exp -2)/d. The annual habitat-weighted emission for the entire HBL using six wetland classes is estimated as 0.538 +/- 0.187 Tg CH4/yr (range of extreme cases is 0.057 to 2.112 Tg CH4/yr). This value is much lower than expected, based on previous emission estimates from northern wetlands.

Description: A method for detecting cirrus clouds in terms of brightness temperature differences between narrowbands at 8, 11, and 12 microns has been proposed by Ackerman et al. In this method, the variation of emissivity with wavelength for different surface targets was not taken into consideration. Based on state-of-the-art laboratory measurements of reflectance spectra of terrestrial materials by Salisbury and D'Aria, it is found that the brightness temperature differences between the 8- and 11-microns bands for soils, rocks, and minerals, and dry vegetation can vary between approximately -8 and +8 K due solely to surface emissivity variations. The large brightness temperature differences are sufficient to cause false detection of cirrus clouds from remote sensing data acquired over certain surface targets using the 8-11-12-microns method directly. It is suggested that the 8-11-12-microns method should be improved to include the surface emissivity effects. In addition, it is recommended that in the future the variation of surface emissivity with wavelength should be taken into account in algorithms for retrieving surface temperatures and low-level atmospheric temperature and water vapor profiles.

Description: On the basis of 30 samples from near-simultaneous overwater measurements by pairs of anemometers located at different heights in the Gulf of Mexico and off the Chesapeake Bay, Virginia, the mean and standard deviation for the exponent of the power-law wind profile over the ocean under near-neutral atmospheric stability conditions were determined to be 0.11 +/- 0.03. Because this mean value is obtained from both deep and shallow water environments, it is recommended for use at sea to adjust the wind speed measurements at different heights to the standard height of 10 m above the mean sea surface. An example to apply this P value to estimate the momentum flux or wind stress is provided.

Description: A new airborne rain-mapping radar (ARMAR) has been developed by NASA and the Jet Propulsion Laboratory for operation on the NASA Ames DC-8 aircraft. The radar operates at 13.8 GHz, the frequency to be used by the radar on the Tropical Rainfall Measuring Mission (TRMM). ARMAR simulates the TRMM radar geometry by looking downward and scanning its antenna in the cross-track direction. This basic compatibility between ARMAR and TRMM allows ARMAR to provide information useful for the TRMM radar design, for rain retrieval algorithm development, and for postlaunch calibration. ARMAR has additional capabilities, including multiple polarization, Doppler velocity measurement, and a radiometer channel for brightness temperature measurement. The system has been tested in both ground-based and airborne configurations. This paper describes the design of the system and shows results of field tests.

Description: The Earth Radiation Budget Experiment (ERBE) consisted of a scanning radiometer and non-scanning radiometers on each of three spacecraft. These instruments began flying in October 1984. The nonscanning radiometers continue to operate, providing broadband radiation measurements of the Earth's outgoing longwave radiation and reflected solar radiation, in addition to measurements of the solar output. The Clouds and Earth Radiant Energy System (CERES) features a scanning radiometer, which is an improved version of the ERBE scanning radiometer, and will fly on the Tropical Rainfall Measurement Mission and Earth Observation System platforms in the late nineties. The CERES project will provide not only radiant fluxes at the 'top of the atmosphere' (TOA), but also at the surface and will compute radiant flux divergence through the atmosphere.

Description: Global and regional temperature variations in the lower troposphere and lower stratosphere are examined for the period 1979-92 from Microwave Sounder Unit (MSU) data obtained by the Television Infrared Observation Satellite (TIROS)-N series of National Oceanic and Atmospheric Administration (NOAA) operational satellites. In the lower troposphere, globally-averaged temperature variations appear to be dominated by tropical El Nino (warm) and La Nina (cool) events and volcanic eruptions. The Pinatubo volcanic eruption in June 1991 appears to have initiated a cooling trend which persisted through the most recent data analyzed (July, 1992), and largely overwhelmed the warming from the 1991-92 El Nino. The cooling has been stronger in the Northern Hemisphere than in the Southern Hemisphere. The temperature trend over the 13.5 year satellite record is small (+0.03 C) compared to the year-to-year variability (0.2-0.4 C), making detection of any global warming signal fruitless to date. However, the future global warming trend, currently predicted to be around 0.3 C/decade, will be much easier to discern should it develop. The lower stratospheric temperature record is dominated by warm episodes from the Pinatubo eruption and the March 1982 eruption of El Chichon volcano.

Description: Several types of algorithms are generally used to process digital imagery such as Landsat data. The most commonly used algorithms perform the task of registration, compression, and classification. Because there are different techniques available for performing registration, compression, and classification, imagery data users need a rationale for selecting a particular approach to meet their particular needs. This collection of registration, compression, and classification algorithms was developed so that different approaches could be evaluated and the best approach for a particular application determined. Routines are included for six registration algorithms, six compression algorithms, and two classification algorithms. The package also includes routines for evaluating the effects of processing on the image data. This collection of routines should be useful to anyone using or developing image processing software. Registration of image data involves the geometrical alteration of the imagery. Registration routines available in the evaluation package include image magnification, mapping functions, partitioning, map overlay, and data interpolation. The compression of image data involves reducing the volume of data needed for a given image. Compression routines available in the package include adaptive differential pulse code modulation, two-dimensional transforms, clustering, vector reduction, and picture segmentation. Classification of image data involves analyzing the uncompressed or compressed image data to produce inventories and maps of areas of similar spectral properties within a scene. The classification routines available include a sequential linear technique and a maximum likelihood technique. The choice of the appropriate evaluation criteria is quite important in evaluating the image processing functions. The user is therefore given a choice of evaluation criteria with which to investigate the available image processing functions. All of the available evaluation criteria basically compare the observed results with the expected results. For the image reconstruction processes of registration and compression, the expected results are usually the original data or some selected characteristics of the original data. For classification processes the expected result is the ground truth of the scene. Thus, the comparison process consists of determining what changes occur in processing, where the changes occur, how much change occurs, and the amplitude of the change. The package includes evaluation routines for performing such comparisons as average uncertainty, average information transfer, chi-square statistics, multidimensional histograms, and computation of contingency matrices. This collection of routines is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 662K of 8 bit bytes. This collection of image processing and evaluation routines was developed in 1979.

Description: A common approach to supervised classification and prediction in artificial intelligence and statistical pattern recognition is the use of decision trees. A tree is "grown" from data using a recursive partitioning algorithm to create a tree which has good prediction of classes on new data. Standard algorithms are CART (by Breiman Friedman, Olshen and Stone) and ID3 and its successor C4 (by Quinlan). As well as reimplementing parts of these algorithms and offering experimental control suites, IND also introduces Bayesian and MML methods and more sophisticated search in growing trees. These produce more accurate class probability estimates that are important in applications like diagnosis. IND is applicable to most data sets consisting of independent instances, each described by a fixed length vector of attribute values. An attribute value may be a number, one of a set of attribute specific symbols, or it may be omitted. One of the attributes is delegated the "target" and IND grows trees to predict the target. Prediction can then be done on new data or the decision tree printed out for inspection. IND provides a range of features and styles with convenience for the casual user as well as fine-tuning for the advanced user or those interested in research. IND can be operated in a CART-like mode (but without regression trees, surrogate splits or multivariate splits), and in a mode like the early version of C4. Advanced features allow more extensive search, interactive control and display of tree growing, and Bayesian and MML algorithms for tree pruning and smoothing. These often produce more accurate class probability estimates at the leaves. IND also comes with a comprehensive experimental control suite. IND consists of four basic kinds of routines: data manipulation routines, tree generation routines, tree testing routines, and tree display routines. The data manipulation routines are used to partition a single large data set into smaller training and test sets. The generation routines are used to build classifiers. The test routines are used to evaluate classifiers and to classify data using a classifier. And the display routines are used to display classifiers in various formats. IND is written in C-language for Sun4 series computers. It consists of several programs with controlling shell scripts. Extensive UNIX man entries are included. IND is designed to be used on any UNIX system, although it has only been thoroughly tested on SUN platforms. The standard distribution medium for IND is a .25 inch streaming magnetic tape cartridge in UNIX tar format. An electronic copy of the documentation in PostScript format is included on the distribution medium. IND was developed in 1992.

Description: The program AUTOCLASS III, Automatic Class Discovery from Data, uses Bayesian probability theory to provide a simple and extensible approach to problems such as classification and general mixture separation. Its theoretical basis is free from ad hoc quantities, and in particular free of any measures which alter the data to suit the needs of the program. As a result, the elementary classification model used lends itself easily to extensions. The standard approach to classification in much of artificial intelligence and statistical pattern recognition research involves partitioning of the data into separate subsets, known as classes. AUTOCLASS III uses the Bayesian approach in which classes are described by probability distributions over the attributes of the objects, specified by a model function and its parameters. The calculation of the probability of each object's membership in each class provides a more intuitive classification than absolute partitioning techniques. AUTOCLASS III is applicable to most data sets consisting of independent instances, each described by a fixed length vector of attribute values. An attribute value may be a number, one of a set of attribute specific symbols, or omitted. The user specifies a class probability distribution function by associating attribute sets with supplied likelihood function terms. AUTOCLASS then searches in the space of class numbers and parameters for the maximally probable combination. It returns the set of class probability function parameters, and the class membership probabilities for each data instance. AUTOCLASS III is written in Common Lisp, and is designed to be platform independent. This program has been successfully run on Symbolics and Explorer Lisp machines. It has been successfully used with the following implementations of Common LISP on the Sun: Franz Allegro CL, Lucid Common Lisp, and Austin Kyoto Common Lisp and similar UNIX platforms; under the Lucid Common Lisp implementations on VAX/VMS v5.4, VAX/Ultrix v4.1, and MIPS/Ultrix v4, rev. 179; and on the Macintosh personal computer. The minimum Macintosh required is the IIci. This program will not run under CMU Common Lisp or VAX/VMS DEC Common Lisp. A minimum of 8Mb of RAM is required for Macintosh platforms and 16Mb for workstations. The standard distribution medium for this program is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format and a 3.5 inch diskette in Macintosh format. An electronic copy of the documentation is included on the distribution medium. AUTOCLASS was developed between March 1988 and March 1992. It was initially released in May 1991. Sun is a trademark of Sun Microsystems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories. DEC, VAX, VMS, and ULTRIX are trademarks of Digital Equipment Corporation. Macintosh is a trademark of Apple Computer, Inc. Allegro CL is a registered trademark of Franz, Inc.

Description: A multispectral, multiresolution (MSMR) method is developed for analyzing scenes of overlapping cloud layers. The MSMR method is applied to data from the NOAA 11 advanced very high resolution radiometer (AVHRR) and the high-resolution infrared radiometer sounder (HIRS-2). The data are from a nighttime oceanic scene in which a semitransparent cirrus veil overlays a large-scale stratus cloud. Low-cloud and clear-sky radiances are determined using a spatial coherence technique. Middle to upper level cloud pressures and radiances are estimated from HIRS-2 15 micrometer CO2 band radiometric data. The MSMR method improves the interpretation of a nighttime, oceanic scene containing thin cirrus over a large-scale stratiform cloud. If, for example, the same scene is analyzed using only the AVHRR 10.8 micrometer channel, the accompanying retrieved cloud heights are found to be between the cirrus and stratus cloud heights and are incorrectly identified as midlevel altostratus clouds. Theoretical radiative transfer model results for both water droplet spheres and randomly oriented hexagonal ice crystals are compared to observed AVHRR brightness temperature differences (BTD) between the 3.7- and 10.8 micrometer channels (BTD(sup 34)) and between the 10.8- and 12- micrometer channels (BTD(sup 45)) to distinguish among the effects of cloud optical depth, particle size, and phase for both single-layer clouds and overlapping two-layer clouds. Theoretical BTD calculations are used to estimate the range of effective particle sizes for eac h cloud layer. The data for the cirrus in the case study region near Bermuda are consistent with theoretical results for relatively small randomly oriented hexagonal ice crystals. The observed BTD(sup 34) and BTD(sup 45) values are lower for the cirrus above a lower-level cloud than for single-level cirrus with no underlying cloud. In certain cases the BTD analysis provides a way to distinguish between clouds composed of supercooled water droplets rather than ice particles. Analysis of nighttime data permits determination of stratus infrared optical depths smaller than 4.

Description: The lower boundary of a spectral mechanistic model is prescribed with 100 hPa geopotentials, and its performance during a November 1989 through March 1990 integration is compared with National Meteorological Center observations. Although the stratopause temperatures quickly become biased near the pole in both hemispheres, the model develops a residual mean circulation which shows significant descent over the winter pole and ascent in the tropics and over the summer pole at pressures less than 10 hPa. The daily correspondence of observed to modeled features in the upper stratosphere and mesosphere degrades after one month. However, the long-term variability qualitatively follows the observations. The results of off-line transport experiments are also described. A passive tracer is instantaneously injected into the flow over the poles and evolves in a manner which is consistent with the residual mean circulation. It demonstrates a significant cross-equatorial flux in the mesosphere near solstice, and air which originates in the southern hemisphere polar mesosphere can be found descending deep into the nothern polar stratosphere at the end of the integration. Nitrous oxide is also transported, and its ability to act as a dynamical tracer is evaluated by comparison to the evolution of the passive tracer.

Description: Aircraft-based observations of turbulence fields of velocity, moisture, and temperature are used to study coherent turbulent structures that dominate turbulent transfer of moisture and heat above three different eco-systems. Flux traces are defragmented, to reconstruct the presumed full size (along the sampled transect) of these structures, and flux traces are simplified by elimination of those that contribute negligibly to the flux estimate. Structures are analyzed in terms of size, spatial distribution, and contribution to the flux, in the four 'quadrant' modes of eddy-covariance transfer (excess up/down and deficit up/down). The effect of nonlinear detrending of moisture and temperature data on this 'structural analysis,' over surfaces with heterogeneous surface wetness, is also examined. Results over grassland, wetland, and moist and dry agricultural land, show that nonlinear detrending may provide a more physically realistic description of structures. Significant differences are observed between structure size and associated relative flux contribution, between moist and dry areas, with smaller structures playing a more important role over the moist areas. Structure size generally increases with height, as spatial reorganization from smaller structures into larger ones takes place. This coincides with a gradual loss of surface 'signature' (position and clustering of plumes above localized source areas). The data are expected to provide a basis for an eventual statistical description of boundary-layer transfer events , and help to interpret the link between boundary-layer transfer and hydrological surface conditions.

Description: An ocean general circulation model (OGCM) is used to study the response of ocean heat and mass transport to positive and negative heat flux anomalies at the ocean surface. As expected, tropical and low-latitude mixed layers respond rapidly (e-folding time about 50-70 years) to external forcing, while the response of the high-latitude mixed layer, especially the Southern Ocean and northern North Atlantic, is very slow (e-folding time greater than 300 yr). The overall response is faster for negative than positive heat flux anomaly at the surface. The meridional heat transport changes by 15% in the first 50 yr in the southern high latitudes. Surprisingly, for the next 400-500 yr the change is very small. The analysis shows that the meridional mass transport intensifies in response to a negative surface heat flux anomaly but weakens in response to a positive heat flux anomaly. For example, at model year 100 the North Atlantic Deep Water (NADW) is reduced from about 18 Sv to about 10 Sv for the positive heat flux experiment but increased to about 26 Sv for the negative heat flux experiment.

Description: Data from the Earth Radiation Budget Experiment scanning radiometer aboard the NOAA-9 operational meteorological satellite are used to investigate the spatial variability of Outgoing Longwave Radiation (OLR). Daily and monthly radiation maps at 2.5 deg latitude-longitude scale are used as a basis for the study. The regions of greatest variability are in the tropics and subtropics. Storm tracks such as the South Pacific convergence zone appear as regions of high OLR variability. Spatial spectra in longitude show two regimes of OLR. At large scales (wavenumbers less than 6), the spatial spectrum is flat. For wavenumbers greater than 10, the spectra decrease as wavenumber to the -3 power. The spatial spectrum of daily anomalies from the mean is a strong function of latitude and season, with interesting features. Correlations of daily anomalies from the monthly mean decrease exponentially in latitude but have a damped-wave structure in longitude. The spatial variability of the daily maps, as measured by degree variance, have 10 times the power at degree 24 than the monthly maps, but at scales between 1 and 10, the degree variance is practically the same for daily as for monthly.

Description: This paper contains two parts. In the first part, a new set of diagnostic equations is derived for the third-order moments for a buoyancy-driven flow, by exact inversion of the prognostic equations for the third-order moment equations in the stationary case. The third-order moments exhibit a universal structure: they all are a linear combination of the derivatives of all the second-order moments, bar-w(exp 2), bar-w theta, bar-theta(exp 2), and bar-q(exp 2). Each term of the sum contains a turbulent diffusivity D(sub t), which also exhibits a universal structure of the form D(sub t) = a nu(sub t) + b bar-w theta. Since the sign of the convective flux changes depending on stable or unstable stratification, D(sub t) varies according to the type of stratification. Here nu(sub t) approximately equal to wl (l is a mixing length and w is an rms velocity) represents the 'mechanical' part, while the 'buoyancy' part is represented by the convective flux bar-w theta. The quantities a and b are functions of the variable N(sub tau)(exp 2), where N(exp 2) = g alpha derivative of Theta with respect to z and tau is the turbulence time scale. The new expressions for the third-order moments generalize those of Zeman and Lumley, which were subsequently adopted by Sun and Ogura, Chen and Cotton, and Finger and Schmidt in their treatments of the convective boundary layer. In the second part, the new expressions for the third-order moments are used to solve the ensemble average equations describing a purely convective boundary laye r heated from below at a constant rate. The computed second- and third-order moments are then compared with the corresponding Large Eddy Simulation (LES) results, most of which are obtained by running a new LES code, and part of which are taken from published results. The ensemble average results compare favorably with the LES data.

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: Calibration of polarimetric radar systems is a field of research in which great progress has been made over the last few years. POLCAL (Polarimetric Radar Calibration) is a software tool intended to assist in the calibration of Synthetic Aperture Radar (SAR) systems. In particular, POLCAL calibrates Stokes matrix format data produced as the standard product by the NASA/Jet Propulsion Laboratory (JPL) airborne imaging synthetic aperture radar (AIRSAR). POLCAL was designed to be used in conjunction with data collected by the NASA/JPL AIRSAR system. AIRSAR is a multifrequency (6 cm, 24 cm, and 68 cm wavelength), fully polarimetric SAR system which produces 12 x 12 km imagery at 10 m resolution. AIRSTAR was designed as a testbed for NASA's Spaceborne Imaging Radar program. While the images produced after 1991 are thought to be calibrated (phase calibrated, cross-talk removed, channel imbalance removed, and absolutely calibrated), POLCAL can and should still be used to check the accuracy of the calibration and to correct it if necessary. Version 4.0 of POLCAL is an upgrade of POLCAL version 2.0 released to AIRSAR investigators in June, 1990. New options in version 4.0 include automatic absolute calibration of 89/90 data, distributed target analysis, calibration of nearby scenes with calibration parameters from a scene with corner reflectors, altitude or roll angle corrections, and calibration of errors introduced by known topography. Many sources of error can lead to false conclusions about the nature of scatterers on the surface. Errors in the phase relationship between polarization channels result in incorrect synthesis of polarization states. Cross-talk, caused by imperfections in the radar antenna itself, can also lead to error. POLCAL reduces cross-talk and corrects phase calibration without the use of ground calibration equipment. Removing the antenna patterns during SAR processing also forms a very important part of the calibration of SAR data. Errors in the processing altitude or in the aircraft roll angle are possible causes of error in computing the antenna patterns inside the processor. POLCAL uses an altitude error correction algorithm to correctly remove the antenna pattern from the SAR images. POLCAL also uses a topographic calibration algorithm to reduce calibration errors resulting from ground topography. By utilizing the backscatter measurements from either the corner reflectors or a well-known distributed target, POLCAL can correct the residual amplitude offsets in the various polarization channels and correct for the absolute gain of the radar system. POLCAL also gives the user the option of calibrating a scene using the calibration data from a nearby site. This allows precise calibration of all the scenes acquired on a flight line where corner reflectors were present. Construction and positioning of corner reflectors is covered extensively in the program documentation. In an effort to keep the POLCAL code as transportable as possible, the authors eliminated all interactions with a graphics display system. For this reason, it is assumed that users will have their own software for doing the following: (1) synthesize an image using HH or VV polarization, (2) display the synthesized image on any display device, and (3) read the pixel locations of the corner reflectors from the image. The only inputs used by the software (in addition to the input Stokes matrix data file) is a small data file with the corner reflector information. POLCAL is written in FORTRAN 77 for use on Sun series computers running SunOS and DEC VAX computers running VMS. It requires 4Mb of RAM under SunOS and 3.7Mb of RAM under VMS for execution. The standard distribution medium for POLCAL is a .25 inch streaming magnetic tape cartridge in UNIX tar format. It is also available on a 9-track 1600 BPI magnetic tape in DEC VAX FILES-11 format or on a TK50 tape cartridge in DEC VAX FILES-11 format. Other distribution media may be available upon request. Documentation is included in the price of the program. POLCAL 4.0 was released in 1992 and is a copyrighted work with all copyright vested in NASA.

Description: A significant decline in Antarctic total column ozone and upper air temperatures has been observed in recent years. Furthermore, high correlations between monthly mean values of ozone and stratospheric temperature have been measured above Syowa, Antarctica. For the observations reported here, data from TOMS (Total Ozone Mapping Spectrometer) aboard the Nimbus 7 satellite have been used to examine the 1980 to 1990 decrease in total column ozone above the South African Antarctic base of SANAE (70 deg 18 min S, 2 deg 21 min W). The cooling of the Antarctic stratosphere above SANAE during this period has been investigated by examining upper air temperatures at the 150, 100, 70, 50, and 30 hPa levels obtained from daily radiosonde balloon launches. Furthermore, these two data sets have been used to examine long-term, medium-term, and short-term correlations between total column ozone and the temperatures at each of the five levels. The trend in SANAE total column ozone has been found to be -4.9 DU/year, while upper air temperatures have been found to decrease at around 0.3 C/year. An analysis of monthly average SANAE total column ozone has shown the decrease to be most severe during the month of September with a trend of -7.7 DU/year. A strong correlation (r(exp 2) = 0.92) has been found between yearly average total column ozone and temperature at the 100 hPa level. Daily ozone and temperature correlations show high values from September to November, at a time when the polar vortex is breaking down.

Description: Vertical profiles of ozone and temperature have been measured at McMurdo Station, Antarctica, during the springs of 1986 to 1991, roughly every two days from 25 August to 31 October. Comparisons of temporal histories and average vertical structure for these years reveals some striking consistency in the ozone depletion process. Ozone depletion generally begins in early September, and with a half-life of 20-30 days, reaches its maximum in mid-October. The depletion occurs almost exclusively between 12 and 20 km. At the time of maximum depletion total ozone has been decreased roughly 40 percent while ozone between 12 and 20 km has been reduced 80 percent. Recovery generally begins in late October with the influx, above 20 km, of ozone rich air from the lower latitudes. From this record the worst years for ozone depletion were 1987, 1989, and 1990. A new region of ozone depletion, below 12 km, was observed in 1991, coinciding with the entrainment of a volcanic cloud into the polar vortex.

Description: In this paper we study the robustness with respect to stability of the closed-loop system with collocated rate sensor using LQG (mean square rate) optimized compensators. Our main result is that the transmission zeros of the compensator are precisely the structure modes when the actuator/sensor locations are 'pinned' and/or 'clamped': i.e., motion in the direction sensed is not allowed. We have stability even under parameter mismatch, except in the unlikely situation where such a mode frequency of the assumed system coincides with an undamped mode frequency of the real system and the corresponding mode shape is an eigenvector of the compensator transfer function matrix at that frequency. For a truncated modal model - such as that of the NASA LaRC Phase Zero Evolutionary model - the transmission zeros of the corresponding compensator transfer function can be interpreted as the structure modes when motion in the directions sensed is prohibited.

Description: In optimal placement of actuators for stochastic systems, it is commonly assumed that the actuator noise variances are not related to the feedback matrix and the actuator locations. In this paper, we will discuss the limitation of that assumption and develop a more practical noise variance model. Various properties associated with optimal actuator placement under the assumption of this noise variance model are discovered through the analytical study of a second order system.

Description: Optimal regulation of hyperbolic systems in the presence of unknown disturbances is considered. Necessary conditions for determining the optimal control that tracks a desired trajectory in the presence of the worst possible perturbations are developed. The results also characterize the worst possible disturbance that the system will be able to tolerate before any degradation of the system performance. Numerical results on the control of a vibrating beam are presented.

Description: This paper presents the development of a general-purpose fuzzy logic (FL) control methodology for isolating the external vibratory disturbances of space-based devices. According to the desired performance specifications, a full investigation regarding the development of an FL controller was done using different scenarios, such as variances of passive reaction-compensating components and external disturbance load. It was shown that the proposed FL controller is robust in that the FL-controlled system closely follows the prespecified ideal reference model. The comparative study also reveals that the FL-controlled system achieves significant improvement in reducing vibrations over passive systems.

Description: Remote viewing is critical for teleoperations, but the inherent limitations of standard video reduce the operator's effectiveness. These limitations have been compensated for in many ways, from using the operator's adaptability, to augmenting his capability with feedback from a variety of sensors and simulations. Omniview can overcome some of these limitations and improve the operator's efficiency without adding additional sensors or computational burden. It can minimize the potential collisions with facility equipment, provide peripheral vision, and display multiple images simultaneously from a single input device. The Omniview technology provides electronic pan, tilt, magnify, and rotational orientation within a hemispherical field-of-view without any moving parts. Image sizes, viewing directions, scale, offset, etc., may be adjusted to fit operator needs. This paper discusses the derivation of the image transformation, the design of the electronics, and two applications to telepresence that are under development. These are Video Emulated Tweening (VET), and Manipulator Guidance and Positioning (ManGAP). The VET effort uses Omniview to compensate for time-delayed video in teleoperation of remote vehicles. In ManGAP two Omniview systems are used to provide two sets of orientation vectors to points in the field-of-view (FOV). These vectors then provide absolute position information to both control the position of the telerobot, and to avoid collisions with the work sight equipment.

Description: A fuzzy-neural control system simulation was developed for the control of a camera platform used to observe aircraft on final approach to an aircraft carrier. The fuzzy-neural approach to control combines the structure of a fuzzy knowledge base with a supervised neural network's ability to adapt and improve. The performance characteristics of this hybrid system were compared to those of a fuzzy system and a neural network system developed independently to determine if the fusion of these two technologies offers any advantage over the use of one or the other. The results of this study indicate that the fuzzy-neural approach to control offers some advantages over either fuzzy or neural control alone.

Description: A fuzzy classifier system that discovers rules for controlling a mathematical model of a pH titration system was developed by researchers at the U.S. Bureau of Mines (USBM). Fuzzy classifier systems successfully combine the strengths of learning classifier systems and fuzzy logic controllers. Learning classifier systems resemble familiar production rule-based systems, but they represent their IF-THEN rules by strings of characters rather than in the traditional linguistic terms. Fuzzy logic is a tool that allows for the incorporation of abstract concepts into rule based-systems, thereby allowing the rules to resemble the familiar 'rules-of-thumb' commonly used by humans when solving difficult process control and reasoning problems. Like learning classifier systems, fuzzy classifier systems employ a genetic algorithm to explore and sample new rules for manipulating the problem environment. Like fuzzy logic controllers, fuzzy classifier systems encapsulate knowledge in the form of production rules. The results presented in this paper demonstrate the ability of fuzzy classifier systems to generate a fuzzy logic-based process control system.

Description: Neural networks are an outgrowth of interdisciplinary studies concerning the brain. These studies are guiding the field of Artificial Intelligence towards the, so-called, 6th Generation Computer. Enormous amounts of resources have been poured into R/D. Wavelet Transforms (WT) have replaced Fourier Transforms (FT) in Wideband Transient (WT) cases since the discovery of WT in 1985. The list of successful applications includes the following: earthquake prediction; radar identification; speech recognition; stock market forecasting; FBI finger print image compression; and telecommunication ISDN-data compression.

Description: Artificial neural systems (ANS), also known as neural networks, are an attempt to develop computer systems that emulate the neural reasoning behavior of biological neural systems (e.g. the human brain). As such, they are loosely based on biological neural networks. The ANS consists of a series of nodes (neurons) and weighted connections (axons) that, when presented with a specific input pattern, can associate specific output patterns. It is essentially a highly complex, nonlinear, mathematical relationship or transform. These constructs have two significant properties that have proven useful to the authors in signal processing and process modeling: noise tolerance and complex pattern recognition. Specifically, the authors have developed a new network learning algorithm that has resulted in the successful application of ANS's to high speed signal processing and to developing models of highly complex processes. Two of the applications, the Weld Bead Geometry Control System and the Welding Penetration Monitoring System, are discussed in the body of this paper.

Description: Neural network capabilities include automatic and organized handling of complex information, quick adaptation to continuously changing environments, nonlinear modeling, and parallel implementation. This viewgraph presentation presents Bellcore work on applications, learning chip computational function, learning system block diagram, neural network equalization, broadband access control, calling-card fraud detection, software reliability prediction, and conclusions.

Description: Over the past three years, our group has concentrated on the application of neural network methods to the training of controllers for real-world systems. This presentation describes our approach, surveys what we have found to be important, mentions some contributions to the field, and shows some representative results. Topics discussed include: (1) executing model studies as rehearsal for experimental studies; (2) the importance of correct derivatives; (3) effective training with second-order (DEKF) methods; (4) the efficacy of time-lagged recurrent networks; (5) liberation from the tyranny of the control cycle using asynchronous truncated backpropagation through time; and (6) multistream training for robustness. Results from model studies of automotive idle speed control serve as examples for several of these topics.

Description: The research mission is the development of computer assisted diagnostic (CAD) methods for improved diagnosis of medical images including digital x-ray sensors and tomographic imaging modalities. The CAD algorithms include advanced methods for adaptive nonlinear filters for image noise suppression, hybrid wavelet methods for feature segmentation and enhancement, and high convergence neural networks for feature detection and VLSI implementation of neural networks for real time analysis. Other missions include (1) implementation of CAD methods on hospital based picture archiving computer systems (PACS) and information networks for central and remote diagnosis and (2) collaboration with defense and medical industry, NASA, and federal laboratories in the area of dual use technology conversion from defense or aerospace to medicine.

Description: This viewgraph presentation presents four working analog VLSI vision chips: (1) time-derivative retina, (2) zero-crossing chip, (3) resistive fuse, and (4) figure-ground chip; work in progress on computing motion and neuromorphic systems; and conceptual and practical lessons learned.

Description: The problem under consideration in this viewgraph presentation is to understand, predict, and control the fluid mechanics of dynamic maneuvers, unsteady boundary layers, and vortex dominated flows. One solution is the application of neural networks demonstrating closed-loop control. Neural networks offer unique opportunities: simplify modeling of three dimensional, vortex dominated, unsteady separated flow fields; are effective means for controlling unsteady aerodynamics; and address integration of sensors, controllers, and time lags into adaptive control systems.

Description: Irvine Sensors Corporation (ISC), working closely with JPL under BMDO/ONR sponsorship, is developing a radically new neural computing technology. Primarily aimed at discrimination and target recognition for BMDO missile interceptor applications, it appears to have near term commercial applicability to such problems as handwriting and face recognition, just to name two. In its earliest form it will be able to perform inner product computation using 262 thousand 64x64 templates (weighted synapse arrays) where the 64(exp 5) weights can all be changed every millisecond. Internal switching provides an inherent capability to zoom, translate, or rotate the templates. The 3D silicon architecture is manufactured on a commercial, high volume DRAM production line at very low cost, enabling its commercialization. Two technology thrusts are beginning: in the first, the 64 layer capability of 3DANN-I will be extended to 1024 layers and beyond. In the second layer size will be shrunk to 2-3 millimeters to reduce layer costs. Our workshop goal is to expose this technology to the neural network community as an emerging tool for their use and to obtain their desire for its future development.

Description: Electronic and optoelectronic hardware implementations of highly parallel computing architectures address several ill-defined and/or computation-intensive problems not easily solved by conventional computing techniques. The concurrent processing architectures developed are derived from a variety of advanced computing paradigms including neural network models, fuzzy logic, and cellular automata. Hardware implementation technologies range from state-of-the-art digital/analog custom-VLSI to advanced optoelectronic devices such as computer-generated holograms and e-beam fabricated Dammann gratings. JPL's concurrent processing devices group has developed a broad technology base in hardware implementable parallel algorithms, low-power and high-speed VLSI designs and building block VLSI chips, leading to application-specific high-performance embeddable processors. Application areas include high throughput map-data classification using feedforward neural networks, terrain based tactical movement planner using cellular automata, resource optimization (weapon-target assignment) using a multidimensional feedback network with lateral inhibition, and classification of rocks using an inner-product scheme on thematic mapper data. In addition to addressing specific functional needs of DOD and NASA, the JPL-developed concurrent processing device technology is also being customized for a variety of commercial applications (in collaboration with industrial partners), and is being transferred to U.S. industries. This viewgraph p resentation focuses on two application-specific processors which solve the computation intensive tasks of resource allocation (weapon-target assignment) and terrain based tactical movement planning using two extremely different topologies. Resource allocation is implemented as an asynchronous analog competitive assignment architecture inspired by the Hopfield network. Hardware realization leads to a two to four order of magnitude speed-up over conventional techniques and enables multiple assignments, (many to many), not achievable with standard statistical approaches. Tactical movement planning (finding the best path from A to B) is accomplished with a digital two-dimensional concurrent processor array. By exploiting the natural parallel decomposition of the problem in silicon, a four order of magnitude speed-up over optimized software approaches has been demonstrated.

Description: We have previously reported on the use of neural networks for detection and identification of faults in complex microprocessor controlled powertrain systems. The data analyzed in those studies consisted of the full spectrum of signals passing between the engine and the real-time microprocessor controller. The specific task of the classification system was to classify system operation as nominal or abnormal and to identify the fault present. The primary concern in earlier work was the identification of faults, in sensors or actuators in the powertrain system as it was exercised over its full operating range. The use of data from a variety of sources, each contributing some potentially useful information to the classification task, is commonly referred to as sensor fusion and typifies the type of problems successfully addressed using neural networks. In this work we explore the application of neural networks to a different diagnostic problem, the diagnosis of faults in newly manufactured engines and the utility of neural networks for process control.

Description: The Intelligent Satellite Control Software (ISACS) for the geoMagnetic tail observation satellite named GEOTAIL (launched in July 1992) has been successfully developed. ISACS has made it possible by applying Artificial Intelligence (AI) technology including an expert system to autonomously generate a tracking schedule, which originally used to be conducted manually. Using ISACS, a satellite operator can generate a maximum four day period of stored command stream autonomously and can easily confirm its safety. The ISACS system has another function -- to diagnose satellite troubles and to suggest necessary remedies. The workload of satellite operators has drastically been reduced since ISACS has been introduced into the operations of GEOTAIL.

Description: Time delay and small capacity of communication are the primary constraint in super long distance telerobotic systems such as astronautical robotic tasks. Intelligent telerobotics is thought to break this constraint. We aim to realize this super long distance telerobotic system with object handling knowledge base and intelligent monitoring. We will discuss physical and technical factors for this purpose.

Description: This paper reports experiments involving the handling of flexible parts (e.g. wires) when using a teleoperated system with time delay. The task is principally a peg-in-hole task involving the wrapping of a wire around two posts on the task-board. It is difficult to estimate the effects of the flexible parts; therefore, on-line teleoperation is indispensable for this class of unpredictable task. We first propose a teleoperation system based on the predictive image display, then describe an experimental teleoperation testbed with a four second transmission time delay. Finally, we report on wire handling operations that were performed to evaluate the performance of this system. Those experiments will contribute to future advanced experiments for the MITI ETS-7 mission.

Description: The focus of this paper is on the detection of natural, as opposed to human-made, objects. The distinction is important because, in the context of image analysis, natural objects tend to possess much greater variability in appearance than human-made objects. Hence, we shall focus primarily on the use of algorithms that 'learn by example' as the basis for image exploration. The 'learn by example' approach is potentially more generally applicable compared to model-based vision methods since domain scientists find it relatively easier to provide examples of what they are searching for versus describing a model.

Description: Remote terrestrial sensing (RTS) data is constantly being collected from a variety of space-based and earth-based sensors. The collected data, and especially 'value-added' analyses of the data, are finding growing application for commercial, government, and scientific purposes. The scale of this data collection and analysis is truly enormous; e.g., by 1995, the amount of data available in just one sector, NASA space science, will reach 5 petabytes. Moreover, the amount of data, and the value of analyzing the data, are expected to increase dramatically as new satellites and sensors become available (e.g., NASA's Earth Observing System satellites). Lockheed and other companies are beginning to provide data and analysis commercially. A critical issue for the exploitation of collected data is the dissemination of data and value-added analyses to a diverse and widely distributed customer base. Customers must be able to use their computational environment (eventually the National Information Infrastructure) to obtain timely and complete information, without having to know the details of where the relevant data resides and how it is accessed. Customers must be able to routinely use standard, widely available (and, therefore, low cost) analyses, while also being able to readily create on demand highly customized analyses to make crucial decisions. The diversity of user needs creates a difficult software problem: how can users easily state their needs, while the computational environment assumes the responsibility of finding (or creating) relevant information, and then delivering the results in a form that users understand? A software agent is a self-contained, active software module that contains an explicit representation of its operational knowledge. This explicit representation allows agents to examine their own capabilities in order to modify their goals to meet changing needs and to take advantage of dynamic opportunities. In addition, the explicit representation allows agents to advertize their capabilities and results to other agents, thereby allowing the collection of agents to reuse each others work.

Description: The optical rain gauge (ORG) by ScTI is designed to be tolerant to reduction of infrared beam intensity due to a dirty lens. Recently there is interest in long term use of optical gauges onboard buoys at sea. Because of logistics, these systems are serviced infrequently, i.e., every several months. Due to the proximity of the gauges to the sea surface, salt can be expected to be deposited on the lens. To obtain an indication of the potential for dirty lens to affect the ORG calibration, two simple experiments were conducted. In the control experiment, a sample ORG was compared to three other ORG's during natural rain events. Next a translucent mask was placed on the transmitter lens of the sample ORG and again data were collected under natural rain conditions. The mask reduced the gain of the perturbed ORG by about 30%. The perturbed ORG operated rather well in that the mask only causes a change in the gain and does not cause data drop out at low rain rates. However, the reduced gain would seriously impact an assessment of rain statistics.

Description: ScTI has performed a detailed analysis of four optical rain gauge ORG-105 sensors tested by Wallops Island on 8 May 1992. The four ORG's tested were S/N 2236, 2237, 2239, and 2241. Shown is a 30 minute time series of the individual ORG's, the ORG average, and the weighing gauge. The sensors tracked well with rainrates (RR) up to 45 mm/hr for the period. Also shown is a plot of accumulated rainfall over the same period. It can be seen that even though the ORG's tracked well, some ORG's tended to read higher and some read lower during the event.