ALBERT

Description: In June 1978, the Seasat satellite was launched carrying, among other instruments, the Seasat-A scatterometer system (SASS), which produced ambiguous wind speed and direction data at the ocean surface. A fifteen day subset of dealiased wind vector data with the inherent ambiguities removed was produced for the period of September 6-20, 1978. On September 8, SASS began to observe a development of frontal cyclogenesis in the South Pacific off the east coast of New Zealand, in an area of few surface observations. A large mature cyclone contained weak warm and cold fronts and an occlusion with a strong horizontal wind shear. Satellite imagery shows that a strong upper-level jet streak was moving rapidly over the area of the surface frontal occlusion and as the jet passed over this area a new vortex formed. This cyclogenesis event was studied using 50-km resolution scatterometer surface wind data. High-resolution fields of wind vectors, divergence and vorticity are computed and plotted from the scatterometer data to study the structure and development of the newly formed cyclonic vortex, not otherwise possible using conventional observations.

Description: The NASA scatterometer (NSCAT) is a spaceborne scatterometer scheduled to be deployed in the mid-1990s. An analysis of the wind retrieval error distribution for wind estimates based on backscatter measurements made by the NSCAT instrument is presented. The results are based on an end-to-end simulation of the scatterometer instrument and data processing. In general, the distribution of the wind speed error, when normalized, is independent of the true wind speed and direction. The wind speed error can be characterized by a normal distribution. The wind direction error is independent of the true wind speed, but depends on the true wind direction. Details for wind vectors with true wind speeds from 3 m/s to 33 m/s and true wind directions from 0 to 360 deg are presented.

Description: Modulation of a rain wave pattern by longer waves has been studied. An analytical model taking into account capillarity effects and obliquity of short waves has been developed. Modulation rates in wave number and amplitude have been computed. Experiments were carried out in a wave tank. First results agree with theoretical models, but higher values of modulation rates are measured. These results could be taken into account for understanding the radar response from the sea surface during rain.

Description: Rainfall modification of directional scatterometer response from the sea surface was simulated in wind-wave tank experiments. Data show that for the range of conditions in laboratory experiments, rain enhances radar cross section for all azimuthal angles relative to wind direction. This result broadens previous measurements, which showed that scatterometer response increases with increasing rainfall for radars pointing upwind. But more to the point, the data also show that the directional dynamic-range of scatterometry diminishes rapidly as rainfall rate increases. Thus, while it may be possible to determine wind speed and direction during rain, it will require adequate system sensitivity.

Description: An account of the construction of surface pressure fields from Seasat-A satellite scatterometer (SASS) winds as carried out by different methods, and the comparison of these pressure fields with those derived from in situ ship observations is presented. On the assumption that the pressure adjusts itself instantaneously to the motion field, it may be computed by various methods. One of these makes use of planetary boundary theory, and of the possible techniques in this category a two-layer iterative scheme admitting of the parametrization of diabatic and baroclinic effects and of secondary flow was chosen. A second method involves the assumption of zero two-dimensional divergence, leading to a Laplace's equation (the balance equation) in pressure, with the wind field serving as a forcing function. This method does not accommodate adiabatic or baroclinic effects, and requires a knowledge of the pressure at all boundary points. Two comparison fields are used for validation: the conventional operational analyses of the US National Meteorological Center (NMC), and the special analyses of the Gulf of Alaska Experiment (GOASEX), which were done by hand. The results of the computations were as follows: (1) The pressure fields, as computed from the SASS winds alone, closely approximated the NMC fields in regions where reasonable in situ coverage was available (typically, one or two mb differences over most of the chart, three to four mb in extreme cases); (2) In some cases the SASS-derived pressure fields displayed high-resolution phenomena not detected by the NMC fields, but evident in the GOASEX data; and, (3) As expected, the pressure fields derived from the balance equation were much smoother and less well resolved than the SASS-derived or NMC fields. The divergence as measured from the SASS winds is smaller than, but of the same order of magnitude as, the vorticity.

Description: The accuracy of temperature, pressure, potential temperature, and horizontal wind measurements is discussed in connection with the use of Meteorological Measurement System data in the AAOE. The vertical distribution of temperature measurements and latitudinal variations of the zonal wind for 12 flights over Antarctica during the 1987 AAOE campaign are summarized. Model atmospheres from 0 to 32 km at 70 deg and 55 deg S for the August-September period are constructed. Above the 420 K isentropic surface, the polar vortex remains strong throughout August and September of 1987.

Description: The formulation and evolution of polar stratospheric ice clouds are simulated using a one-dimensional model of cloud microphysics. It is found that the optical thickness and particle size of ice clouds depend on the cooling rate of the air in which the cloud formed. It is necessary that there be an energy barrier to ice nucleation upon the preexisting aerosols in order to account for the cooling rate dependence of the cloud properties.

Description: The history of minimum temperatures at 50 and 70 mb is examined from NMC, UK Met O and ECMWF analyses. MSU channel 24 data are similarly inspected. South Pole sonde data are used to calculate saturation humidity mixing ratio as a function of altitude and time throughout 1987. Saturation with respect to ice could be maintained for water mixing ratios of 3.5 ppmv for a period of about 80 days from mid-June to mid-September. Dehydration to mixing ratios of 1 ppmv or less was possible sporadically. Data from the ER-2 flights between 53 S and 72 S are used in conjunction with particle size measurements and air parcel trajectories to demonstrate the dehydration occurring over Antarctica. Water mixing ratios at the latitude of Punta Arens (53 S), in conjunction with tracer measurements and trajectory analysis, show that at potential temperatures from about 325 to 400 K, the dryness (less than 3 ppmv) had its origin over Antarctica rather than in the tropics. Water mixing ratios within the Antarctic vortex varied from 1.5 to 3.8 ppmv, with a strong isentropic gradient being evident in the region of high potential vorticity gradients.

Description: Measurements of total reactive nitrogen, NOy, total water vapor, and aerosols were made as part of the Airborne Antarctic Ozone Experiment. The measurements were made using instruments located onboard the NASA ER-2 aircrafts which conducted twelve flights over the Antarctic continent reaching altitudes of 18 km at 72 S latitude. Each instrument utilized an ambient air sample and provided a measurement up to 1 Hz or every 200 m of flight path. The data presented focus on the flights of Aug. 17th and 18th during which Polar Stratospheric Clouds (PSCs) were encountered containing concentrations of 0.5 to 1.0 micron diameter aerosols greater than 1 cm/cu. The temperature pressure during these events ranged as low as 184 K near 75 mb pressure, with water values near 3.5 ppm by volume (ppmv). With the exception of two short periods, the PSC activity was observed at temperatures above the frost point of water over ice. The data gathered during these flights are analyzed and presented.

Description: A multichannel statistical approach is used to retrieve rainfall rates from the brightness temperature T(B) observed by passive microwave radiometers flown on a high-altitude NASA aircraft. T(B) statistics are based upon data generated by a cloud radiative model. This model simulates variabilities in the underlying geophysical parameters of interest, and computes their associated T(B) in each of the available channels. By further imposing the requirement that the observed T(B) agree with the T(B) values corresponding to the retrieved parameters through the cloud radiative transfer model, the results can be made to agree quite well with coincident radar-derived rainfall rates. Some information regarding the cloud vertical structure is also obtained by such an added requirement. The applicability of this technique to satellite retrievals is also investigated. Data which might be observed by satellite-borne radiometers, including the effects of nonuniformly filled footprints, are simulated by the cloud radiative model for this purpose.

Description: This paper presents the results of a field program using a ground-based Raman lidar system to observe changes in moisture profiles as a cold and a warm front passed over the NASA/Goddard Space Flight Center in Greenbelt, Maryland. The lidar operating only during darkness is capable of providing continuous high vertical resolution profiles of water vapor mixing ratio and aerosol scattering ratio from near the surface to about 7 km altitude. The lidar data acquired on three consecutive nights from shortly after sunset to shortly before sunrise, along with upper air data from specially launched rawinsondes, have provided a unique visualization of the detailed structure of the two fronts.

Description: An improved version of the GISS Model II cumulus parameterization designed for long-term climate integrations is used to study the effects of entrainment and multiple cloud types on the January climate simulation. Instead of prescribing convective mass as a fixed fraction of the cloud base grid-box mass, it is calculated based on the closure assumption that the cumulus convection restores the atmosphere to a neutral moist convective state at cloud base. This change alone significantly improves the distribution of precipitation, convective mass exchanges, and frequencies in the January climate. The vertical structure of the tropical atmosphere exhibits quasi-equilibrium behavior when this closure is used, even though there is no explicit constraint applied above cloud base.

Description: The joint frequency distribution technique was used to analyze buoyancy fluxes in the marine atmospheric boundary layer (MABL) for the cloud street regime noted during the Genesis of Atlantic Lows Experiment. It is found that for the lower half of the MABL, the buoyancy flux is mainly generated by the rising thermals and the sinking compensating ambient air, and is mainly consumed by the entrainment and detrainment of thermals, penetrative convection, and the entrainment from the MABL top. If the buoyancy flux is primarily driven by the temperature flux, these buoyancy-flux generating processes should be the same for the lower boundary layers over land and ocean. The results of the scale analysis of the buoyancy flux agree well with those obtained for mesoscale cellular convection during the Air-Mass Transformation Experiment.

Description: Results are presented on numerical experiments that were carried out to investigate the mechanisms of the observed variabilities in wind and convection associated with supercloud clusters (SCCs), westerly wind bursts, and 30-60 day oscillations in the western Pacific region. It is shown that the generation of a 30-60 day eastward propagating precipitation pattern in the Lau and Peng (1987) model, which can be identified as SCC, is accompanied by convective clusters coming in opposite direction to that of the SCC itself. The results suggest that the westward propagating cloud clusters are produced at the initial stage of the 30-60 day disturbance due to mutual adjustment of the large-scale flow and heating.

Description: Multidecadal time series of surface winds from central tropical Pacific islands are used to compute trends in the trade winds between the end of WWII and 1985. Over this period, averaged over the whole region, there is no statistically significant trend in speed or zonal or meridional wind (or pseudostress). However, there is some tendency, within a few degrees of the equator, toward weakening of the easterlies and increased meridional flow toward the equator. Anomalous conditions subsequent to the 1972-73 ENSO event make a considerable contribution to the long-term trends. The period 1974-80 has been noted previously to have been anomalous, and trends over that period are sharply greater than those over the longer records.

Description: This paper presents an analysis of the uncertainties expected in vertical velocities using a vertically pointing airborne Doppler radar which has a nadir or zenith-pointing beam. To examine the expected uncertainty, the Doppler velocity equation for a moving platform is derived and it is applied to cases of nadir-fixed and stabilized beams. The main emphasis of the paper is on the effect of platform stability on the deduced vertical air motions and it is shown that the antenna must be stabilized to obtain desired accuracy in the vertical velocity measurements.

Description: Monthly fields of shortwave radiation (SR) and latent heat flux (LE) over the central and eastern tropical Pacific between 1980 and 1983 have been computed using satellite data. They are the dominant variable components of surface thermal forcing on the ocean in this time scale. During the 1982-1983 ENSO episode, surface-wind convergence and cloudiness associated with the displacement of equatorial organized convection caused a reduction in both the SR into the ocean and the LE out of the ocean. The lag-correlation coefficients between the forcing (SR-LE) and the sea surface temperature are found to be significantly high outside the equatorial region, showing that surface thermal forcing is the dominant factor in sea surface temperature change. In the narrow equatorial wave guide, ocean dynamics play a more important role, and surface heat flux is a consequence rather than the cause of sea surface temperature change.

Description: An aircraft experiment has been conducted with a dual-frequency (10 GHz and 35 GHz) radar/radiometer system and an 18-GHz radiometer to test various rain-rate retrieval algorithms from space. In the experiment, which took place in the fall of 1988 at the NASA Wallops Flight Facility, VA, both stratiform and convective storms were observed. A ground-based radar and rain gauges were also used to obtain truth data. An external radar calibration is made with rain gauge data, thereby enabling quantitative reflectivity measurements. Comparisons between path attenuations derived from the surface return and from the radar reflectivity profile are made to test the feasibility of a technique to estimate the raindrop size distribution from simultaneous radar and path-attenuation measurements.

Description: Azimuthal response of a scatterometer to radiation scattered by the sea surface was studied in a wind-wave tank. The variation of the normalized radar cross section with the azimuth angle is fitted by a three-term series. Results show that the upwind-downwind asymmetry decreases as the wind speed increases. The crosswind modulation depends on the wind velocity. The results show that the evolution of the long-wind-crosswind ratio evolves with wind speed in a manner similar to the evolution of the isotropy of short capillary-gravity waves. The maximum of the isotropy of the short wind waves is obtained for wind velocities close to 4 m/s. For the same value of the velocity, the variations of radar response between long-wind and crosswind directions is minimum. For lower or higher values of wind velocities the directional accuracy of the radar increases, since the wind-wave field tends to align in the wind direction.

Description: The present interpretation of the radar cross section sigma exp 0 measured by satellite altimeters implies that the rms wave slope gamma is controlled solely by the local wind. However, parameters of wave spectra, including the exponent in the power law for the equilibrium range, depend on sea maturity. The latter is characterized by the nondimensional fetch, x = gX/U-squared. Consequently, gamma and sigma exp 0 are controlled by both U and the wind fetch X. Geosat data for one year are used jointly with in-situ wind and wave observations to assess the fetch-related error trend in altimeter wind speeds. This trend results in overestimated winds in the regions and seasons characterized by a high x, and vice versa. A procedure for wind speed retrieval based on processing sigma exp 0 jointly with the significant wave height information contained in the altimeter wave forms is proposed.

Description: Seasat scatterometer data over the Arabian Sea are used to build wind-stress fields during July and August 1978. They are first compared with 3-day wind analyses from ship data along the Somali coast. Seasat scatterometer specifications of 2-m/s and 20-deg accuracy are fulfilled in almost all cases. The exceptions are for winds stronger than 14 m/s, which are underestimated by the scatterometer by 15 percent. Wind stress is derived from these wind data using a bulk formula with a drag coefficient depending on the wind intensity. A successive-correction objective analysis is used to build the wind-stress field over the Arabian Sea with 2 x 2-deg and 6-day resolution. The final wind-stress fields are not significantly dependent on the objective analysis because of the dense coverage of the scatterometer. The combination of scatterometer and coastal ship data gives the best coverage to resolve monsoon wind structures even close to the coast. The final wind stress fields show wind features consistent with other monthly mean wind stress field. However, a high variability is observed on the 6-day time scale.

Description: Cumulus processes involved in the interaction and merging of clouds under the influence of different imposed conditions (including large-scale lifting forcing, environmental wind shear, and cloud microphysical processes) were studied using simulations with a three-dimensional model. The design of the study was to generate several convective clouds randomly inside the model domain, and then to observe and analyze the interactions and merging between the simulated clouds. Ten merged clouds were identified. Seven of these, each involving two previously separated clouds, generally lie along a line parallel to the initial environmental wind shear vector, while one (also a two-cloud system) lies along a line perpendicular to the wind shear vector prior to merging. The remaining two merging systems involve three parent clouds each; they are a combination of parallel and perpendicular cells. The merging mechanisms associated with three-cloud merging cases are studied by examining the temperature, pressure, and wind fields prior to, during, and following the merging.

Description: The possibility of initiating the growth of ice sheets by solar insolation variations is examined. The study is conducted using a climate model with three different orbital configurations corresponding to 116,000 and 106,000 yr before the present and a modified insolation field with greater reductions in summer insolation at high northern latitudes. Despite the reduced summer and fall insolation, the model fails to maintain snow cover through the summer at locations of suspected ice sheet initiation. The results suggest that there is a discrepancy between the model's response to Milankovitch perturbations and the geophysical evidence of ice sheet initiation. If the model results are correct, the growth of ice shown by geophysical evidence would have occurred in an extremely ablative environment, demanding a complicated strategy.

Description: Thirty-three Landsat TM scenes of California stratocumulus cloud fields were acquired as part of the FIRE Marine Stratocumulus Intensive Field Observations in July 1987. They exhibit a wide variety of stratocumulus structures. Analysis has so far focused upon the July 7 scene, in which aircraft from NASA, NCAR, and the British Meteorological Office repeatedly gathered data across a stratocumulus-fair weather cumulus transition. The aircraft soundings validate the cloud base temperature threshold determined by spatial coherence analysis of the TM thermal band. Brightness variations in the stratocumulus region exhibit a -5/3 power-law decrease of the wavenumber spectra for scales larger than the cloud thickness, about 200 m, changing to a -3 power at smaller scales. Observations by an upward-looking three-channel microwave radiometer on San Nicolas Island also show the -5/3 power-law in total integrated liquid water, suggesting that the largest-scale TM brightness variations are primarily due to variations in the liquid water. The Kolmogorov 5/3 power suggests that for some purposes liquid water in turbulent stratocumulus clouds may be treated as a passive scalar, simply reflecting variations in vertical velocity. This may be tested using the velocities measured by the aircraft.

Description: New nonlinear partial differential equations containing terrain curvature and its rate of change are derived that describe the flow of an atmospheric density current. Unlike the classical hydraulic-type equations for density currents, the new equations are valid for two-dimensional, gradually varied flow over highly curved terrain, hence suitable for computing unsteady (or steady) flows over arbitrary mountain/valley profiles. The model assumes the atmosphere above the density current exerts a known arbitrary variable pressure upon the unknown interface. Later this is specialized to the varying hydrostatic pressure of the atmosphere above. The new equations yield the variable velocity distribution, the interface position, and the pressure distribution that contains a centrifugal component, often significantly larger than its hydrostatic component. These partial differential equations are hyperbolic, and the characteristic equations and characteristic directions are derived. Using these to form a characteristic mesh, a hypothetical unsteady curved-flow problem is calculated, not based upon observed data, merely as an example to illustrate the simplicity of their application to unsteady flows over mountains.

Description: The sources of sub-Saharan precipitation are studied using diagnostic procedures integrated into the code of the GISS climate model. Water vapor evaporating from defined source regions is 'tagged', allowing the determination of the relative contributions of each evaporative source to the simulated July rainfall in the Sahel. Two June-July simulations are studied to compare the moisture sources, moisture convergence patterns and the spatial variations of precipitation for rainy and drought conditions. Results for this case study indicate that patterns of moisture convergence and divergence over northern Africa had a stronger influence on model rainfall over the sub-Sahara than did evaporation rates over the adjacent oceans or moisture advection from ocean to continent. While local continental evaporation contributed significant amounts of water to Sahelian precipitation in the'rainy' simulation, moisture from the Indian Ocean did not precipitate over the Sahel in either case.

Description: The long-term climatic evolution of the earth has been studied on the basis of one-dimensional, globally-averaged climate models yielding only a qualitative understanding of climatic history, and in any case proceeding under a series of potentially invalid assumptions. One such major assumption, which invites comparison with a three-dimensional GCM, is that of fixed relative humidity. A GCM may also be used to study the problem of water loss from both the earth and Venus.

Description: An analysis of the relationship between the IMF section boundary crossing, solar flares, the sunspot 11 year cycle variation and the thunderstorm index is given, using the superposition epoch method, for data from more than 13,000 thunderstorms from 10 meteorological stations in the Beijing area and the Northeast region during 1957 to 1978. The results show that for some years a correlation exists between the thunderstorm index and the positive IMF section boundary crossing. The thunderstorm index increases obviously within three days near the crossing and on the seventh day after the crossing. The influence of the crossing on thunderstorms is stronger in the first half year than the latter half year. For different classes of solar flares, the influences are not equally obvious. The solar flares which appeared on the west side, especially in the western region (from 0 to 30 deg) have the most obvious influence. There is no discernible correlation between the thunderstorm index and the sunspot eleven-year cycle.

Description: The possibility that solar variations associated with the 11-year solar cycle may be the cause of the changes in tropospheric weather and climate has been the subject to scientific investigation for several decades. Meteorologists are greatly concerned with the changes in tropospheric phenomena. An attempt was made to find solar activity related changes in tropospheric weather, by the modulation of the quasi-biennial oscillation (QBO) of zonal wind at 50 mb. Rainfall and surface temperature data for a period of about three solar cycles, 1953 to 1988, from various stations in the Indian subcontinent were utilized. By extension, a possible teleconnection was looked for between the temperature changes in middle atmospheric levels and surface temperature when the data are stratified according to east or west phase of the QBO. The temperature data were averaged for January and February to represent the winter temperature and for July and August to represent the summer temperature.

Description: The part of energy of the planetary waves which enters the stratosphere depends on conditions of planetary wave generation and propagation through the tropopause, and the part of planetary wave energy which enters the mesosphere depends on conditions of planetary wave propagation through the stratopause. An attempt is made to estimate connections between extratropical middle atmosphere temperature long term variations and portions of energy of planetary waves which enter the mesosphere and stratosphere during winter seasons in Northern and Southern Hemispheres. Interannual variations of temperatures at the 30 km and 70 km levels are investigated for the central winter months of the period 1970 to 1986. This period includes the descending branch of the 20th solar cycle and the whole 21st cycle. Calculations are made on the basis of measurements at Heiss Island and Molodezhnaya.

Description: The direct modulation of temperature of the mid-latitude mesosphere by the solar-cycle EUV variation, which leads to greater heat input at higher solar activity, is well established. Middle atmosphere temperature modulation by the solar cycle is independently confirmed by the variation of reflection heights of low frequency radio waves in the lower ionosphere, which are regularly monitored over about 30 years. As explained elsewhere in detail, these reflection heights depend on the geometric altitude of a certain isobaric surface (near 80 k), and on the solar ionizing Lyman-alpha radiation flux. Knowing the solar cycle variation of Lyman-alpha how much the measured reflection heights would be lowered with the transition from solar minimum to maximum can be calculated, if the vertical baric structure of the neutral atmosphere would remain unchanged. An discrepancy between expected and observed height change must be explained by an uplifting of the isobaric level from solar minimum to maximum, caused by the temperature rise in the mesosphere. By integrating the solar cycle temperature changes over the height region of the middle atmosphere, and assuming that the lower boundary (tropopause) has no solar cycle variation, the magnitude of this uplifting can be estimated. It is given for the Lidar-derived and for the rocket-measured temperature variations. Comparison suggests that the real amplitude of the solar cycle temperature variation in the mesosphere is underestimated when using the rocket data, but probably overestimated with the Lidar data.

Description: Before the introduction of the Quasi Biennial Oscillation (Q.B.O.) in the study of the solar atmosphere relationship by Labitzke (1987) and Labitzke and Van Loon (1988), the only region of the atmosphere where an effect of a change in solar activity was generally admitted was the mesosphere. The response of the mesosphere, in phase with the solar activity, was found to be about one order of magnitude above model expectancy (around 10 to 20 Kelvin). It was observed independently of the season and maximized around 70 km (Chanin et al. 1987). However, from the same study, it was shown that the response of the stratosphere of opposite sign, clearly seen during winter and autumn, was at the threshold of detection in spring and summer. In the stratosphere, it was shown later that the separation of the data taking into account the sign of the Q.B.O. amplifies the negative correlation of the stratospheric temperature with solar activity in winter; it then becomes more significantly negative for the East phase of the Q.B.O. than when the data are all mixed (Labitzke and Chanin 1988). The studies of the seasonal response of the atmosphere to solar effect is crucial to understand the possible mechanism responsible of such a solar activity Q.B.O. relationship, knowing that the global dynamic circulation is quite different according to the seasons. The question is examined as to whether such separation of the data according to the phase of the Q.B.O. has any impact on the solar response of the middle atmosphere for seasons other than winter.

Description: There is no doubt that the antropogenic effect play an important role in the effects of corpuscular radiation on weather and climate. The task, however, is to distinguish between antropogenic effect in the atmosphere due to human activities and natural climatic fluctuations influencing biological systems. The increase in global temperature during the past 100 years is in relatively good coincidence with the increase in geomagnetic (corpuscular) activity. It is concluded that it could have been the increase in temperature on the Northern Hemisphere, due to the processes occurring in the auroral oval under enhanced corpuscular radiation which led to an increased atmospheric concentration of CO2 in the past. Both processes, i.e., antropogenic and solar activity effects, should be therefore intensively studied due to their important role for elucidating the past and present global change mainly in temperature, climate and biological systems.

Description: A blocking pattern which formed over eastern North America following the landfall of Hurricane Juan during November 1985 was investigated. It is hypothesize that latent heat released in the Hurricane's rainfall was either directly or indirectly responsible for the large observed 500 mb height rises over eastern Canada during the formation of this block. This idea is evaluated with a diagnostic model for the height tendency field which includes latent heat release as a forcing function. The total column heating is calculated using satellite-derived precipitation estimates. These estimates are qualitatively congruent with observations, but overestimate light rainfall and underestimate heavy rainfall. The calculations reveal that the direct contribution of the heating to the 500 mb height tendency field is small relative to the quasigeostrophic forcing. However, maxima in heating coincide with regions where anticyclonic potential vorticity is generated. Once such region is just upstream of the location of large 500 mb height rises in the incipient block. An indirect role is proposed for the heating in this case. Specifically, anticyclonic potential vorticity is generated near the heating maxima; this vorticity is then advected downstream, forcing the 500 mb heights to rise and the block to develop.

Description: Simple models are being developed to simulate interaction of planetary and synoptic-scale waves incorporating the effects of large-scale topography; eddy heat and momentum fluxes (or nonlinear dynamics); radiative heating/cooling; and latent heat release (precipitation) in synoptic-scale waves. The importance of latent heat release is determined in oceanic storm tracks for temporal variability and time-mean behavior of planetary waves. The model results were compared with available observations of planetary and synoptic-scale wave variability and time-mean circulation. The usefulness of monitoring precipitation in oceanic storm tracks by satellite observing systems was ascertained. The modeling effort includes two different low-order quasi-geostrophic models-time-dependent version and climatological mean version. The modeling also includes a low-order primitive equation model. A time-dependent, multi-level version will be used to validate the two-level Q-G models and examine effects of spherical geometry.

Description: Applications of mesosphere stratosphere troposphere radar to mesoscale meteorology are discussed. The applications include using the radar either as a research tool to improve our understanding of certain dynamical systems or as part of a network used to provide input data for weather forecasting. The workhorse of the operational observing network is the radiosonde balloon which provides measurements of pressure, temperature, humidity, and winds up to heights of 16 to 20 km. Horizontal and vertical measurement capabilities, reflectivity data, derivable quantities and parameters, and special operational requirements are surveyed.

Description: An understanding is being developed for processes which may be important in the atmosphere, and the definition and analysis of baroclinic experiments utilizing the geophysical fluid flow cells (GFFC) apparatus in microgravity space flights. Included are studies using numerical codes, theoretical models, and terrestrial laboratory experiments. The numerical modeling is performed in three stages: calculation of steady axisymmetric flow, calculation of fastest-growing linear eigenmodes, and nonlinear effects (first, wave-mean flow interactions, then wave-wave interactions). The code can accommodate cylindrical, spherical, or channel geometry. It uses finite differences in the vertical and meridional directions, and is spectral in the azimuthal. The theoretical work was mostly in the area of effects of topography upon the baroclinic instability problem. The laboratory experiments are performed in a cylindrical annulus which has a temperture gradient imposed upon the lower surface and an approximately isothermal outer wall, with the upper and inner surfaces being nominally thermally insulating.

Description: Advanced methods were developed to determine time varying winds and turbulence from digital flight data recorders carried aboard modern airliners. Analysis of several cases involving severe clear air turbulence encounters at cruise altitudes has shown that the aircraft encountered vortex arrays generated by destabilized wind shear layers above mountains or thunderstorms. A model was developed to identify the strength, size, and spacing of vortex arrays. This model is used to study the effects of severe wind hazards on operational safety for different types of aircraft. It is demonstrated that small remotely piloted vehicles and executive aircraft exhibit more violent behavior than do large airliners during encounters with high-altitude vortices. Analysis of digital flight data from the accident at Dallas/Ft. Worth in 1985 indicates that the aircraft encountered a microburst with rapidly changing winds embedded in a strong outflow near the ground. A multiple-vortex-ring model was developed to represent the microburst wind pattern. This model can be used in flight simulators to better understand the control problems in severe microburst encounters.

Description: An essential component in global climate research is accurate cloud cover and type determination. Of the two approaches to texture-based classification (statistical and textural), only the former is effective in the classification of natural scenes such as land, ocean, and atmosphere. In the statistical approach that was adopted, parameters characterizing the stochastic properties of the spatial distribution of grey levels in an image are estimated and then used as features for cloud classification. Two types of textural measures were used. One is based on the distribution of the grey level difference vector (GLDV), and the other on a set of textural features derived from the MaxMin cooccurrence matrix (MMCM). The GLDV method looks at the difference D of grey levels at pixels separated by a horizontal distance d and computes several statistics based on this distribution. These are then used as features in subsequent classification. The MaxMin tectural features on the other hand are based on the MMCM, a matrix whose (I,J)th entry give the relative frequency of occurrences of the grey level pair (I,J) that are consecutive and thresholded local extremes separated by a given pixel distance d. Textural measures are then computed based on this matrix in much the same manner as is done in texture computation using the grey level cooccurrence matrix. The database consists of 37 cloud field scenes from LANDSAT imagery using a near IR visible channel. The classification algorithm used is the well known Stepwise Discriminant Analysis. The overall accuracy was estimated by the percentage or correct classifications in each case. It turns out that both types of classifiers, at their best combination of features, and at any given spatial resolution give approximately the same classification accuracy. A neural network based classifier with a feed forward architecture and a back propagation training algorithm is used to increase the classification accuracy, using these two classes of features. Preliminary results based on the GLDV textural features alone look promising.

Description: Low-altitude microburst wind shear encounters can significantly affect aircraft performance during approach or takeoff. Over the past 25 years, hazardous wind shear was a contributing factor in over two dozen commercial airline accidents in which there were over 500 fatalities. NASA, the FAA, and the National Center for Atmospheric Research were involved in the design and testing of various sensors to detect the hazard. Among the sensors being tested are the ground-based Terminal Doppler Weather Radar (TDWR) and airborne Doppler radar and LIDAR systems. While these sensor systems do measure horizontal wind shear, they do not adequately account for the vertical wind, which is a key component of the microburst hazard to the aircraft. A technique is defined to estimate aircraft hazard from the combined effects of horizontal and vertical winds, given only horizontal wind information.

Description: In general, work has continued on developing and evaluating algorithms designed to manage the Laser Atmospheric Wind Sounder (LAWS) lidar pulses and to compute the horizontal wind vectors from the line-of-sight (LOS) measurements. These efforts fall into three categories: Improvements to the shot management and multi-pair algorithms (SMA/MPA); observing system simulation experiments; and ground-based simulations of LAWS.

Description: Optimized Royal Signals and Radar Establishment (RSRE), Laser True Airspeed System (LATAS) algorithm for low backscatter conditions was developed. The algorithm converts backscatter intensity measurements from focused continuous-wave (CW) airborne Doppler lidar into backscatter coefficients. The performance of optimized algorithm under marginal backscatter signal conditions was evaluated. The 10.6 micron CO2 aerosol backscatter climatologies were statistically analyzed. Climatologies reveal clean background aerosol mode near 10(exp -10)/kg/sq m/sr (mixing ratio units) through middle and upper troposhere, convective mode associated with planetary boundary layer convective activity, and stratospheric mode associated with volcanically-generated aerosols. Properties of clean background mode are critical to design and simulation studies of Laser Atmospheric Wind Sounder (LAWS), a MSFC facility Instrument on the Earth Observing System (Eos). Previous intercomparisons suggested correlation between aerosol backscatter at CO2 wavelength and water vapor. Field measurements of backscatter profiles with MSFC ground-based Doppler lidar system (GBDLS) were initiated in late FY-88 to coincide with independent program of local rawinsonde releases and overflights by Multi-spectral Atmospheric Mapping Sensor (MAMS), a multi-channel infrared radiometer capable of measuring horizontal and vertical moisture distributions. Design and performance simulation studies for LAWS would benefit from the existence of a relationship between backscatter and water vapor.

Description: The orographic flow data set was obtained from a flight program to measure the influence of orographic features on turbulence momentum, heat, and moisture fluxes. The NASA B-57 aircraft instrumented with probes for measuring the three fluctuating wind speed components, temperature, and humidity was the primary measuring vehicle. Ancillary measurements were made with several ground-based sensors. These include the NOAA radar wind profilers, the Boulder wind network, the PROFS mesoscale surface network, the Boulder Atmospheric Observatory 300 m tower, special rawinsonde observations, and the NOAA/WPL Doppler lidar. The major objective of the flight program was to provide planetary boundary layer parameter information for new and current general circulation computer models. A numerical code, WINDER, based on a discrete element technique was run to numerically model the water tunnel simulated flow. Comparison of the analytical model with the experimental results is very good. Physical fluid dynamic principles embedded in the computational model and visual and hot wire anemometer measurements from the simulation are being rationalized to develop a physical understanding of the vortex flow. The results will be interpreted as they pertain to full scale atmospheric flows.

Description: Significant accomplishments during 1987 to 1988 are summarized with regard to each of the major project components. Model 1 requires satisfaction of two nonlinear horizontal momentum equations, the integrated continuity equation, and the hydrostatic equation. Model 2 requires satisfaction of model 1 plus the thermodynamic equation for a dry atmosphere. Model 3 requires satisfaction of model 2 plus the radiative transfer equation. Model 4 requires satisfaction of model 3 plus a moisture conservation equation and a parameterization for moist processes.

Description: Theoretical studies of low frequency and seasonal weather variability; dynamical properties of observational and general circulation model (GCM)-generated records; effects of the hydrologic cycle and latent heat release on extratropical weather; and Earth-system science studies are summarized.

Description: Research efforts during FY-88 have focused on completion of several projects relating to analysis of FGGE data during SOP-1 and on expanded studies of global atmospheric moisture. In particular, a revised paper on the relationship between diabatic heating and baroclinicity in the South Pacific Convergence Zone (SPCZ) was submitted. A summary of completed studies on diagnostic convective parameterization was presented at the Satellite Meteorology and Oceanography Convergence last February. These investigations of diabatic heating in the SPCZ have demonstrated the requirement for a more quantitative description of atmospheric moisture. As a result, efforts were directed toward use of passive remote microwave measurements from the Nimbus-7 SMMR and the DOD's Special Sensor Microwave Imager (SSMI/I) as critical sources of moisture data. Activities this year are summarized.

Description: The focus was on the development of a blocking anticyclone that formed over the North Atlantic in January 1979 and a marine cyclone that deepened explosively prior to the onset of the block. The extended height tendency equation was used as the primary diagnostic tool. Focusing on the domain encompassing the migrating ridge that eventually formed the block, it was found that vorticity advection played the dominant role in the development of the ridge and the formation of the block. Also of interest was an attempt to evaluate the relative importance of synoptic-scale, planetary-scale, and synoptic/planetary-scale interactions as the block developed. To accomplish this, all data fields were partitioned into synoptic and planetary-scale components using a Barnes-type filter. Finally, the cyclone was diagnosed by examining the low level static stability fields associated with the cyclone's development.

Description: A description of each of the criteria used to decide where to locate the Atmospheric Visibility Monitoring (AVM) telescope systems is given, along with a weighting factor for each of them. These criteria include low probability of clouds, fog, smog, haze, low scattering, low turbulence, availability of security and maintenance, and suitability of a site for a potential optical reception station. They will be used to determine which three of several sites under consideration will be used for monitoring visibility through the atmosphere as it applies to an optical ground-based receiving network as may be used in NASA space missions in decades to come.

Description: Background information on the rocket-triggered lightning project at Kennedy Space Center (KSC), a summary of the forecasting problem there, the facilities and equipment available for undertaking field experiments at KSC, previous research activity performed, a description of the atmospheric science field laboratory near Mosquito Lagoon on the KSC complex, methods of data acquisition, and present results are discussed. New sources of data for the 1989 field experiment include measuring the electric field in the lower few thousand feet of the atmosphere by suspending field measuring devices below a tethered balloon. Problems encountered during the 1989 field experiment are discussed. Future prospects for both triggered lightning and lightning-kindled forest fire research at KSC are listed.

Description: The research carried out in the past year dealt with nonlinear baroclinic wave dynamics. The model consisted of an Eady baroclinic basic state and uneven Elkman dissipation at the top and bottom boundaries with/without slopes. The method of solution used a truncated spectral expansion with three zonal waves and one or two meridional modes. Numerical experiments were performed on synoptic scale waves or planetary scale waves with/without wave-wave interaction.

Description: Research objectives are to: (1) describe the synoptic scale variability of moisture over the tropical Pacific Ocean; (2) describe the systems leading to this variability; and (3) develop and implement satellite analysis procedures to facilitate (1) and (2) over the data sparse Pacific.

Description: Numerical modeling simulations of tropical squall lines were begun to determine the role of large scale terrain features over Arizona and Mexico in their initiation and propagation. Installation was completed for a short-base, high resolution lightning location and detection network in and around Tucson. Data from a Doppler wind profiler is being analyzed to determine the role of large scale heating over the inter-mountain plateau region in governing local diurnal wind variations and possible relationships to the monsoon flow. The portable solar photometer for determining high temporal resolution values of the local precipitable water vapor was completed and calibrated. The assembly is nearly completed for a multi-channel microwave passive radiometer to determine local temperature and water vapor profiles.

Description: Global archives were established for ECMWF 12-hour, multilevel analysis beginning 1 January 1985; day and night IR temperatures, and solar incoming and solar absorbed. Routines were written to access these data conveniently from NASA/MSFC MASSTOR facility for diagnostic analysis. Calculations of diabatic heating rates were performed from the ECMWF data using 4-day intervals. Calculations of precipitable water (W) from 1 May 1985 were carried out using the ECMWF data. Because a major operational change on 1 May 1985 had a significant impact on the moisture field, values prior to that date are incompatible with subsequent analyses.

Description: Mesoscale phenomena such as thunderstorm and sea breeze frontal circulations are being investigated using a 50 MHz Doppler wind profiler at the Kennedy Space Center. The profiler installation will begin October 1, 1988 and will be completed by February 17, 1989. The focus of current research and plans for next year include: examination of vertical velocities associated with local thunderstorm activity and sea breeze frontal circulations and compare the vertical velocities to conceptual mesoscale models; implementation of space-time conversion analysis techniques to blend profiler data with National Meteorological Center's model output and other wind data such as jimsphere, windsonde and rawinsonde for mesoscale analysis; development of suggestions for use of wind profiler data in mesoscale analysis and forecasting at Kennedy Space Center; and problems detection in the quality of the profiler data during this research project. Researchers will work closely with MSFC to identify and solve the data quality problems.

Description: A numerical advection model which can be run on a local computer in a real-time forecast environment is described. This isentropic forecast model provides the local forecast office with easy access to temporally and spatially detailed estimates of atmospheric temperature, moisture, and wind field changes between 12-h rawinsonde observations. Case studies are presented based on the use of the model to predict the preconvective environment in spring and summer situations. Short-term forecasts of midtropospheric static stability patterns and stability index changes are traced for several severe storm events with and without the inclusion of surface data observed during the day. Forecast images of VAS low- and midlevel moisture fields and vertical moisture gradients are compared with the observations to determine the utility of the combined model/VAS imagery as a nowcasting guide.

Description: Results are presented on the analysis of the Nimbus-7 satellite data set obtained on regional and seasonal variations in global cloud cover. Four midseason months (April, July, and October 1979 and January 1980) were analyzed for the total cloud amount, the cloud amounts at high, middle, and low altitudes, the cirrus and deep convective clouds, and the cloud and clear-sky 11.5 micron-derived radiances; in addition, noon versus midnight cloud amounts were examined. The Nimbus-7 data are compared with three previously published cloud climatologies, and the differences among these data sets are discussed.

Description: The propagation of orographic gravity waves into an atmosphere with exponentially decreasing density is simulated with a two-dimensional, nonlinear, time-dependent numerical model. After the stationary wave is established over the mountain, the model predicts that wave breaking causes a large reduction of the vertical momentum flux in the flow, not only at levels where wave breaking is present, but also far below the lowest occurrence of overturning. More than half of the decrease in momentum flux is explained by the presence of large amplitude, downward propagating waves, which are generated in regions of wave breaking. The downward propagating waves appear almost simultaneously with overturning, and have nonzero phase speeds, suggesting a strongly nonlinear generation mechanism that depends on local wave properties. The generation of these downward propagating waves is a robust process, insensitive to mountain height, mountain width, or density scale height. These results have important implications for observational studies of orographically generated waves as well as for schemes that seek to parameterize the effects of orogoraphy in large-scale models.

Description: A GCM with new parameterizations of solar and terrestrial radiation, parameterized cloud optical properties, and a simple representation of the cloud liquid water feedback is used with several observational data sets to analyze the effects of cloudiness on the earth's radiation budget. The January and July results from the model are in reasonable agreement with data from Nimbus-7. It is found that the simulated cloudiness overpredicts subtropical and midlatitude cloudiness. The simulated atmospheric cloud radiative forcing is examined. The clear-sky radiation fields obtained by two methods of Cess and Potter (1987) are compared. Also, a numerical experiment was performed to determine the effects of the water vapor continuum on the model results.

Description: A parameterization is presented for computing downward, upward, and net longwave radiation at the earth's surface using data from NOAA sun-synchronous satellites. The parameterization is applied to satellite soundings for April, 1982 over a large region of the tropical Pacific Ocean. Sensitivity studies were used to estimate the random and systematic errors in computed fluxes due to probable errors in TOVS-derived parameters. It is suggested that large biases in the results due to errors in TOVS-derived parameters may be corrected with data from the International Satellite Cloud Climatology Project.

Description: Global, daily, visible and IR radiance measurements from the NOAA-5 Scanning Radiometer (SR) are analyzed for the months of January, April, July and October, 1977 to infer cloud and surface radiative properties. In this paper, the data and analysis method are described. A unique feature of the method is that it utilizes radiative transfer models that simulate the SR measurements using explicit parameters representing the properties of the surface, atmosphere, and clouds. The accuracy of all the results depends primarily on the proper separation of the total radiance distribution into those parts representing clear and cloudy scenes.

Description: Marine stratus clouds were simultaneously observed by nadir Nd:YAG lidar measurements and in situ cloud physics measurements. A procedure was applied to derive the two-dimensional vertical cross section of the liquid water from within the cloud top lidar observations. A comparison to direct in-cloud liquid water observations gave good results. The liquid water retrieval was limited to an effective optical depth of 1.5. The true cloud optical thickness was also obtained from the retrieval procedure to a corresponding limit of 3.8. The optical thickness of the observed marine stratus clouds was predominantly below 3.0.

Description: Consideration is given to the use of the Defense Meteorological Satellite Program's Special Sensor Microwave/Imager (SSM/I) to identify precipitation in warm and cold land and ocean environments. It is shown that the polarization diversity of the SSM/I operating at 85.5 GHz makes it possible to discriminate between low brightness temperatures due to surface water bodies and those due to precipitation. The theoretical sensitivity of SSM/I between 19.35 and 8.5. GHz and the polarization correction for water surfaces and the effects of cloud water are discussed. Examples are presented of observational studies using SSM/I for the delineation of precipitation over land and oceans.

Description: Global, daily, visible, and infrared radiance measurements from the NOAA-5 Scanning Radiometer (SR) are analyzed for the months of January, April, July, and October 1977 to infer surface radiative properties. A radiative transfer model that simulates the spectral and angular characteristics of the NOAA-5 SR measurements is used to retrieve monthly mean surface visible reflectances and temperature at 25 km resolution. These surface properties were found sufficiently accurate for simulation of clear sky radiances to determine global, seasonal variations in cloudiness. Further comparisons of these results with other data highlight the analysis difficulties and radiative model shortcomings that must be overcome to monitor regional and seasonal variations of earth's surface. These preliminary results also provide an estimate of the magnitude of these variations.

Description: The dynamic and microphysical processes associated with convective and stratiform regions in a tropical squall-type convective line are studied with a multidimensional, time-dependent cloud scale model. It is found that a two-dimensional model reproduces many features of tropical squall systems, including the fractional portion of stratiform rain and the stratiform fractional area coverage. A three-dimensional model did not capture stratiform features, but did simulate explosive growth and a convex-leading edge associated with the convective region. Results are presented from sensitivity tests on ice-phase microphysical processes and mesoscale middle and upper level ascent.

Description: Nineteen years of wind data over the North Atlantic are used to calculate a field of wind stress curl. An empirical orthogonal function (EOF) analysis is performed on this field, resulting in spatial patterns of wind stress curl and associated time series. A Monte Carlo technique is used to establish the statistical significance of each spatial pattern, and the associated time series are spectrally analyzed. The first four statistically significant EOF modes represent more than 50 percent of the curl variance, and the spatial patterns of curl associated with these modes exhibit the major elements of North Atlantic climatology. Most of the time series spectral variance is contained in annual and semiannual frequencies. The features observed include the individual annual variation of the subtropical high and the subpolar low, the annual oscillation of intensity between pressure centers, the influence of localized strong SST gradients and associated cyclogenesis regions, and the constant nature of the trades.

Description: The stratocumulus (SC) entrainment and thermodynamics are studied using data from the Dynamics and Chemistry of Marine Stratocumulus (DYCOMS) experiment. The rate of entrainment of air from the free troposphere into the cloud-topped PBL is estimated using a technique based on the measurement of ozone flux and mean distribution. The average measured value of entrainment rate was found to be 3.0 mm/s with a range of 1.0 to 5.0 mm/s for cloudy cases. Thermodynamic budgets are constructed for eight DYCOMS cases. It was found that the divergence of the solar radiative flux is an important component of the boundary-layer energetics during midday, the factor which must be accounted for in mixed-layer models. The net longwave radiative flux profiles show good agreement with theoretical models, but an unambiguous partitioning of the flux divergence between inversion and mixed layers could not be established.

Description: REMOTE sensing using an airborne infrared lidar has shown an unexpected capability to detect open leads in Arctic sea ice and their associated meteorology in winter. It is shown here that vertical profiles of backscattered radiation demonstrate strong returns from hydrometeor plumes originating from leads having a surface water temperature near -1.8 C. Recently refrozen leads are also distinguishable by the lidar backscatter from adjacent thicker, older sea ice. Wide leads release enough energy to create buoyant plumes which penetrate the Arctic boundary layer inversion, transporting heat and moisture into the troposphere. These results show that the role of the Arctic as a global heat sink may need to be reevaluated, and that lead plumes have a significant effect on the radiation budget.

Description: Five-day National Meteorological Center wind and outgoing longwave radiation (OLR) data are used here to study the seasonal, intraseasonal variation of angular momentum and its relationship with tropical convection. It is found that the zonally averaged u component and the zonally averaged OLR show pronounced periodicities around 30-60 days. The globally averaged zonal wind at 250 mbar exhibits strong seasonal variation in connection with the expansion and contraction of the easterly and westerly belt. At upper levels this is manifested as the poleward and equatorward migration of the subtropical jet streams. Results show that the 30- to 60-day oscillation in 250 mbar zonal mean wind and OLR can be represented by three basic meridional structures: tropical modes, tropical-subtropical modes, and extratropical modes.

Description: The characteristics of linear and finite amplitude baroclinic waves during Southern Hemisphere winter were investigated using a three-dimensional global primitive-equation spectral mode, with emphasis on the dynamics of planetary scale modes which propagate into the stratosphere. Linear stability computations showed that the winter zonally averaged basic state is baroclinically unstable, with the most unstable modes being Charney-type modes; however, the results showed that Green-type modes also exist; these give much larger relative amplitude in the stratosphere although they have smaller growth rates. The results suggest that the existence of both Charney and Green modes might explain certain observed properties of the eastward traveling waves, such as the lack of coherence between troposphere and stratosphere at zonal wavenumbers 1 and 2.

Description: NOAA-7 Advanced Very High Resolution Radiometer (AVHRR) Global Vegetation Indices (GVI) were used during the 1986 rainy season (June-September) over Senegal to monitor rainfall. The satellite data were used in conjunction with ground-based measurements so as to derive empirical relationships between rainfall and GVI. The regression obtained was then used to map the total rainfall corresponding to the growing season, yielding good results. Normalized Difference Vegetation Indices (NDVI) derived from High Resolution Picture Transmission (HRPT) data were also compared with actual evapotranspiration (ET) data and proved to be closely correlated with it with a time lapse of 20 days.

Description: An engineering estimate of the spectrum of atmospheric microburst noise radiation in the range 2-20 Hz is developed. This prediction is obtained via a marriage of standard aeroacoustic theory with a numerical computation of the relevant fluid dynamics. The 'computational aeroacoustics' technique applied here to the interpretation of atmospheric noise measurements is illustrative of a methodology that can now be employed in a wide class of problems.

Description: A pulsed CO2 lidar with coherent detection has been used to measure the correlation time of backscatter from an ensemble of atmospheric aerosol particles which are illuminated by the pulsed radiation. The correlation time of the backscatter of the return signal, which is directly related to the velocity spectral width, can be used to study the velocity structure constant of atmospheric turbulence and wind shear. Various techniques for correlation time measurement are discussed, and several measurement results are presented for the technique using the information contained in the statistical distribution of a set of lidar return signal intensities.

Description: A brief review of observations of clouds using satellites highlights open issues and directions for future studies. The key one is improved treatment of the effects of small-scale spatial inhomogeneity in remote sensing data analyses and in the treatment of radiation in climate models, though studies and observations of the spectral dependence of cloud-radiation interactions are also limited. Significant progress in understanding the role of clouds in climate, especially regarding cloud-radiation budget relationships, is expected in the next several years because of an unprecedent suite of global and regional observation and analysis programs.

Description: A maximum likelihood estimation (MLE) technique to the problem of cloud identification using coarse resolution broadband satellite data is developed and tested using simulated satellite observations. The results suggest that, in the determination of cloud conditions for the inversion of satellite-measured radiances to fluxes, the MLE method is an improvement over a Lambertian earth assumption and the clear/cloud threshold used in the inversion of Nimbus 3 and Nimbus 7 data. The use of the MLE method in the operational processing of Earth Radiation Budget Experiment scanner data is considered.

Description: Most empirical studies estimating the climate sensitivity due to cloud amount are based on calculating a cloud sensitivity parameter by means of satellite measurements. A different way to estimate the interaction between clouds and radiation is to determine the effect of radiative heating and cooling by clouds on the generation of available potential energy (GAPE). This provides a measure of the influence of clouds on the general circulation. The radiative contribution of clouds to the GAPE is given by the net cloud generated radiative heating (CGRH) and an efficiency factor. The CGRH, parameterized in terms of cloud type, height, and optical depth, can be estimated from satellite measurements. The efficiency factor is only a function of the temperature distribution of the atmosphere. Thus, taking satellite measurements as for instance those of the ISCCP, the radiative impact of clouds on the general circulation can be inferred.

Description: The Earth Radiation Budget Experiment (ERBE) has been observing the earth during the past 4 years from three satellites. Numerous validation procedures have been applied to the data. Particularly important have been intercomparisons between three channels of scanning radiometers and wide and medium field-of-view radiometers. These intercomparisons and onboard calibration targets have provided assurance of high data quality. In addition to the classic radiation budget parameters: global absorbed and emitted energy, ERBE is producing fluxes on scales of 250 km that can significantly increase understanding of the earth's climate. Of particular interest are ERBE measurements of clear-sky albedos and longwave fluxes, which are part of the ERBE data products. Use of data from the precessing ERBE satellite together with data from the sun-synchronous NOAA-9 and NOAA-10 satellites also provides information on the diurnal cycle of shortwave and longwave fluxes.

Description: This paper examines a parameterization of a quasi-geostrophic eddy transport that takes into account the time variation of eddy transfer coefficients according to Green's (1970) theory. It was found that, in the original eddy transfer relationship of Green, connecting the integral of the northward eddy entropy flux through midlatitudes with the second power of the difference in 500-mb entropy across the region of baroclinic activity, a value of 4 for the exponent is obtained when the temperature gradients at 500 mb are used. When the gradients at 1000 mb are used, an exponent of 1.5 is obtained. The differences in the powers in the eddy transfer relation were explored in a two-level zonally averaged model. It was found that an appropriate choice of power may be of special importance if the model is devised to simulate the seasonal climate cycle or to test astronomical changes inducing different seasonalities.

Description: An airborne lidar and a scanning radiometer aboard an ER-2 aircraft were used to observe a cold cirrus cloud, and a Forward Scattering Spectrometer Probe (FSSP) was used to obtain simultaneous in situ microphysical observations at two altitudes within the cloud. Lidar depolarization ratio data show that the clouds were composed predominantly of ice crystals. At an altitude where the temperature was -62.7 C, the lidar and radiometer analysis gave a visible extinction to infrared absorption ratio (alpha) of 2.3, while the cloud microphysics data provided an alpha value of 3.77. The discrepancy is attributed to undersizing of particles by the FSSP. Direct and remote measurements showed better agreement for a lower layer where the temperature was -47.3 C.

Description: Several attenuation-based methods for estimating the rainfall rate were applied to measurements made by an airborne dual-wavelength radar operating at 0.87 cm, the K(a)-band, and at 3 cm, the X-band. These methods included the traditional Z-R methods, designated Z(X)-R and Z(K)-R for the X- and K(a) band wavelengths, respectively; single- and dual-wavelength surface reference techniques (SRT and DSRT, respectively); and standard dual-wavelength methods with and without range-averaging. As the primary sources of error for these methods are nearly independent, agreement among the rain rates obtained with these methods would lend confidence in the results. Correlation coefficients obtained between the rainfall rates with the Z(X)-R and DSRT methods were generally between 0.7 and 0.9. Good agreement among the methods occurred most often in stratiform rain for rain rates betwen a few mm/hr to about 15 mm/hr, i.e., where attenuation at the shorter wavelength is significant but not so severe as to result in a loss of signal.

Description: This paper describes a method for the retrieval of total precipitable water (W) in dry atmospheres, which relies on the strong water vapor absorption line at 183 GHz (for W less than 0.6 g/sq cm) and on the absorption near 90 GHz (for W above 0.6 g/sq cm). The method is very sensitive to the variations of W less than 0.5 g/sq cm and is complementary to the established methods that use the weak 22 GHz water vapor absorption line to retrieve W in the 1-6 g/sq cm. The technique was demonstrated by the analysis of two Advanced Microwave Moisture Sounder observations of dry atmospheres following cold-air outbreaks on March 13, 1983, and February 23, 1986.

Description: This paper examines the feasibility of deriving the zonal wind at the equator using mean temperature and geopotential height data from satellite and radiosonde/rocketsonde measurements. Using climatological data of the stratosphere and mesosphere based on monthly mean Nimbus 5 Selective Chopper Radiometer and Nimbus 6 Pressure Modulator Radiometer measurements, and stratospheric monthly mean data from the National Meteorological Center for 1979-1986, meridional cross sections and time profiles of the zonal wind in the tropical middle atmosphere are presented. The derived zonal wind at the equator reproduces the mean climatology of the tropical middle atmosphere as well as the general characteristics of the equatorial quasi-biennial (QBO) and semiannual oscillations (SAO) observed in monthly mean radiosonde and rocketsonde data. Although the amplitude of the derived wind QBO is for the most part underestimated relative to direct-wind measurements, the amplitude of the derived wind SAO compares fairly well with rocketsonde observations.

Description: The effect of the viewing zenith angle (VZA) on the cloudiness values observed by a satellite was investigated using a combination of two cloud-amount data sets derived from nearly simultaneous collocated GOES-E and GOES-W radiance measurements over the Pacific Ocean during May 1979 and July 1983. A hybrid bispectral threshold method was used to analyze data for single-layer and total cloudiness. It was found that the cloud fraction values increased with increasing VZA for almost all cases. Low clouds exhibited the greatest increases with a VZA increase for cloud amounts in the 0.1 range, whereas high clouds showed greatest increases for cloud amounts around 0.5. Midlevel clouds showed only a slight dependence on VZA. Total cloudiness increased the most, reflecting its predominantly low-cloud composition.

Description: A reasonably accurate and noniterative saturation adjustment scheme is proposed to calculate: (1) the amount of condensation and/or deposition necessary to remove any supersaturated vapor, or (2) the amount of evaporation and/or sublimation necessary to remove any subsaturation in the presence of cloud droplets and/or cloud ice. This proposed scheme can be implemented for a nonhydrostatic cloud model. The derivation of the scheme, an evaluation of its performance, and tests for sensitivity to variations in a few key parameters are presented.

Description: The classical picture of the influence of midlatitude troughs on cloud patterns is studied in a GCM and in satellite and NMC data, by comparing the cross correlation of the poleward component of the wind and the outgoing longwave radiation (OLR). Over storm tracks, the variance of the OLR is explained by a correlation with the poleward component of the horizontal wind, forced by stronger correlations with the vertical velocity through the cloud and humidity fields. It is found that stratiform clouds have a dominant impact on the model OLR fluctuations, even over much of the tropics.

Description: This paper presents a series of meteorological analyses used to aid the interpretation of the in situ Airborne Antarctic Ozone Experiment (AAOE) observations obtained aboard the ER-2 and DC-8 aircraft and examines the basis and accuracy of the analytical procedure. Maps and sections of meteorological variables derived from the UK Meteorological Office Global Model are presented for ER-2 and DC-8 flight days. It is found that analyzed temperatures and winds are generally in good agreement with AAOE observations at all levels; minor discrepancies were evident only at DC-8 altitudes. Maps of potential vorticity presented on the 428-K potential temperature surface show that the vortex is essentially circumpolar, although there are periods when major distortions are apparent.

Description: The convergence-lightning nowcasting method that began to be operationally tested and evaluated by USAF personnel at the KSC in 1986 uses surface convergence, especially the areal average, to identify the potential for new thunderstorm growth. When this areal method fails, normalized cell convergence, which is defined as the integral of cellular convergence with respect to its area, can be used as a complementary technique in association with the areal one. A detailed examination is conducted for one day, and the relationship of normalized cell convergence to lightning flash counts is examined in a 68-day ensemble.

Description: Techniques for extraction boundary layer parameters from measurements of a short-pulse CO2 Doppler lidar are described. The measurements are those collected during the First International Satellites Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE). By continuously operating the lidar for about an hour, stable statistics of the radial velocities can be extracted. Assuming that the turbulence is horizontally homogeneous, the mean wind, its standard deviations, and the momentum fluxes were estimated. Spectral analysis of the radial velocities is also performed from which, by examining the amplitude of the power spectrum at the inertial range, the kinetic energy dissipation was deduced. Finally, using the statistical form of the Navier-Stokes equations, the surface heat flux is derived as the residual balance between the vertical gradient of the third moment of the vertical velocity and the kinetic energy dissipation. Combining many measurements would normally reduce the error provided that, it is unbiased and uncorrelated. The nature of some of the algorithms however, is such that, biased and correlated errors may be generated even though the raw measurements are not. Data processing procedures were developed that eliminate bias and minimize error correlation. Once bias and error correlations are accounted for, the large sample size is shown to reduce the errors substantially. The principal features of the derived turbulence statistics for two case studied are presented.

Description: A novel coronal current-determining instrument is being used at NASA-KSC which overcomes previous difficulties with wind sensitivity and a voltage-threshold 'deadband'. The mounting of the corona needle at an elevated location reduces coronal and electrode layer space-charge influences on electric fields, rendering the measurement of space charge density possible. In conjunction with a space-charge compensation model, these features allow a more realistic estimation of cloud base electric fields and the potential for lightning strike than has previously been possible with ground-based sensors.

Description: Cirrus cloud radiative and physical characteristics are determined using a combination of ground-based, aircraft, and satellite measurements taken as part of the First ISCCP Regional Experiment (FIRE) Cirrus Intensive Field Observations (IFO) during October and November 1986. Lidar backscatter data are used to define cloud base, center, and top heights and the corresponding temperatures. Coincident GOES 4 km visible (0.65 microns) and 8 km infrared window (11.5 microns) radiances are analyzed to determine cloud emittances and reflectances. Infrared optical depth is computed from the emittance results. Visible optical depth is derived from reflectance using a theoretical ice crystal scattering model and an empirical bidirectional reflectance mode. No clouds with visible optical depths greater than 5 or infrared optical depths less than 0.1 were used in the analysis. Average cloud thickness ranged from 0.5 km to 8 km for the 71 scenes. An average visible scattering efficiency of 2.1 was found for this data set. The results reveal a significant dependence of scattering efficiency on cloud temperature.

Description: Observations of cirrus and altocumulus clouds during the First International Satellite Cloud Climatology Project Regional Experiment (FIRE) are compared to theoretical models of cloud radiative properties. Three tests are performed. First, LANDSAT radiances are used to compare the relationship between nadir reflectance ot 0.83 micron and beam emittance at 11.5 microns with that predicted for model calculations using spherical and nonspherical phase functions. Good agreement is found between observations and theory when water droplets dominate. Poor agreement is found when ice particles dominate, especially using scattering phase functions for spherical particles. Even when compared to a laboratory measured ice particle phase function, the observations show increased side scattered radiation relative to the theoretical calculations. Second, the anisotropy of conservatively scattered radiation is examined using simultaneous multiple angle views of the cirrus from LANDSAT and ER-2 aircraft radiometers. Observed anisotropy gives good agreement with theoretical calculations using the laboratory measured ice particle phase function and poor agreement with a spherical particle phase function. Third, Landsat radiances at 0.83, 1.65, and 2.21 microns are used to infer particle phase and particle size. For water droplets, good agreement is found with King Air FSSP particle probe measurements in the cloud. For ice particles, the LANDSAT radiance observations predict an effective radius of 60 microns versus aircraft observations of about 200 microns. It is suggested that this descrepancy may be explained by uncertainty in the imaginary index of ice and by inadequate measurements of small ice particles by microphysical probes.

Description: It has been known for many years that turbulence significant for aviation often occurs near thunderstorm tops. However, that turbulence is not well-predicted because of an incomplete understanding of the processes which generate it and because of inadequate observations. The current study seeks to alleviate these problems via: (1) a comprehensive review of recent theoretical and experimental studies related to turbulence near thunderstorm tops (TNTT), and (2) three case studies designed to examine the feasibility of using data derived from commercial aircraft to study TNTT. The literature review revealed extensive evidence which showed that convection often produces significant barrier effects; several mesoscale phenomena capable of producing turbulence may occur, depending on wind and stability conditions near the thunderstorm tops. These include two- and three-dimensional lee waves, rotors, Kelvin-Helmholtz instabilities, and Karman vortices. Conventional meteorological data were combined with data derived from the aircraft flight tapes to produce quantitative descriptions of the turbulence and its mesoscale environment for the three cases.

Description: Often in the short-range limited-area numerical weather prediction (NWP) of extratropical weather systems the effects of planetary boundary layer (PBL) processes are considered secondarily important. However, it may not be the case for the regional NWP of mesoscale convective systems over the arid and semi-arid highlands of the southwestern and south-central United States in late spring and summer. Over these dry regions, the PBL can grow quite high up into the lower middle troposphere (600 mb) due to very effective solar heating and hence a vigorous air-land thermal interaction can occur. The interaction representing a major heat source for regional dynamical systems can not be ignored. A one-dimensional PBL-radiation model was developed. The model PBL consists of a constant-flux surface layer superposed with a well-mixed (Ekman) layer. The vertical eddy mixing coefficients for heat and momentum in the surface layer are determined according to the surface similarity theory, while their vertical profiles in the Ekman layer are specified with a cubic polynomial. Prognostic equations are used for predicting the height of the nonneutral PBL. The atmospheric radiation is parameterized to define the surface heat source/sink for the growth and decay of the PBL. A series of real-data numerical experiments has been carried out to obtain a physical understanding how the model performs under various atmospheric and surface conditions. This one-dimensional model will eventually be incorporated into a mesoscale prediction system. The ultimate goal of this research is to improve the NWP of mesoscale convective storms over land.

Description: Numerical simulations were performed with a very high resolution (7.25 km) version of the MASS model (Version 4.0) in an effort to diagnose the vertical wind shear and static stability structure during the Shuttle Challenger disaster which occurred on 28 January 1986. These meso-beta scale simulations reveal that the strongest vertical wind shears were concentrated in the 200 to 150 mb layer at 1630 GMT, i.e., at about the time of the disaster. These simulated vertical shears were the result of two primary dynamical processes. The juxtaposition of both of these processes produced a shallow (30 mb deep) region of strong vertical wind shear, and hence, low Richardson number values during the launch time period. Comparisons with the Cape Canaveral (XMR) rawinsonde indicates that the high resolution MASS 4.0 simulation more closely emulated nature than did previous simulations of the same event with the GMASS model.

Description: ITCZ intensity across the Pacific for the 1975-1976 season is studied. Analysis of the data reveal: zonal bands of convective development; isolated patches of stronger convection; and several zonal bands which merge as tropical intraseasonal oscillations. The ITCZ intensity at 160 deg W is compared with the intensity averaged all the way across the Pacific. Composite outgoing longwave radiation anomaly maps for 22 days of local active ITCZ at 160 deg W, 26 days of zonally averaged active ITCZ, and 16 days of both local and zonal active ITCZ are presented and examined. It is noted that: systematic, yet different, patterns, movement, and temporal behavior correspond to local and zonally extensive convection; intensifying zonal convection moves southward; local convection occurs at two preferential latitudes; and local convection does not mean there is a generally active ITCZ.

Description: The application of satellite data to the study of synoptic characteristics is analyzed. A radiative transfer model is used to generate a database of satellite information on synoptic features. Canonical discriminant analysis is employed to reveal the differences among synoptic sounding classes; and the fine structures within each sounding class are examined with rotated factor analysis. Diagrams of wind and frontal inversions are presented. It is noted that the applicability of satellite data depends on the method used to analyze it and multivariable statistical techniques may be useful for deriving additional information from satellite data.

Description: Nitrous oxide (N2O) flux from vented chambers was measured over intensively and minimally managed bermuda grass hay meadows in a humid, subtropical ecosystem for several years during scheduled sampling protocol following harvest, fertilization and rainfall events while measuring diel N2O emissions once during each of 5 seasonal day growth cycles which divided each calendar year. Soil pools of nitrite NO2(-), nitrate NO3(-), and ammonia (NH3) were measured in soil samples taken at 2 and 10 cm depths during each emission collection to determine transformations of the nitrogen pools coupled with N2O emissions. The highest diel N2O emission occur midday in the Spring cycle, measuring 9.0 g N/ha/d only for several weeks, while emissions dropped to less than 1.0 g N/ha/day during hot, dry, and colder months. Intensively managed meadows (4 fertilizations and harvests per year plus pest management) induced higher seasonal N2O emissions than minimal treatment (1 fertilization and harvest) averaging 2.75 and 5.97 g N/ha/day. Nitrous oxide emission data as responses to soil parameters and environmental parameters were also measured where air temperature, soil moisture, and fertilization were the most powerful factors.

Description: Two nonlinear triggered lightning models have been developed for use in finite difference calculations. Both are based on three species of air chemistry physics and couple nonlinearly calculated air conductivity to Maxwell's equations. The first model is suitable for use in three-dimensional modeling and has been applied to the analysis of triggered lightning on the NASA F106B Thunderstorm Research Aircraft. The model calculates number densities of positive ions, negative ions, and electrons as a function of time and space through continuity equations, including convective derivative terms. The set of equations is closed by using experimentally determined mobilities, and the mobilities are also used to determine the air conductivity. Results from the model's application to the F106B are shown. The second model is two-dimensional and incorporates an enhanced air chemistry formulation. Momentum conservation equations replace the mobility assumption of the first model. Energy conservation equations for neutrals, heavy ions, and electrons are also used. Energy transfer into molecular vibrational modes is accounted for. The purpose for the enhanced model is to include the effects of temperature into the air breakdown, a necessary step if the model is to simulate more than the very earliest stages of breakdown. Therefore, the model also incorporates a temperature-dependent electron avalanche rate. Results from the model's application to breakdown around a conducting ellipsoid placed in an electric field are shown.

Description: This video tape addresses ongoing lightning research and how data is valuable to upcoming projects.