ALBERT

Description: This document contains a description of a comprehensive database that is to be used for certification testing of airborne forward-look windshear detection systems. The database was developed by NASA Langley Research Center, at the request of the Federal Aviation Administration (FAA), to support the industry initiative to certify and produce forward-look windshear detection equipment. The database contains high resolution, three dimensional fields for meteorological variables that may be sensed by forward-looking systems. The database is made up of seven case studies which have been generated by the Terminal Area Simulation System, a state-of-the-art numerical system for the realistic modeling of windshear phenomena. The selected cases represent a wide spectrum of windshear events. General descriptions and figures from each of the case studies are included, as well as equations for F-factor, radar-reflectivity factor, and rainfall rate. The document also describes scenarios and paths through the data sets, jointly developed by NASA and the FAA, to meet FAA certification testing objectives. Instructions for reading and verifying the data from tape are included.

Description: The variabilities of the upper layer of the western Pacific warm pool (WPWP) were observed using satellite infrared data from 1982 to 1991 and altimeter data from November 1986 to September 1989. The warm pool was defined as the area where the sea surface temperatures are above 28 C. The eastern boundary oscillation, the centroid movement, and the upper-layer volume variation of the WPWP were intensively studied. Spectral analysis revealed that the eastern boundary oscillation of the WPWP was related to the El Nino event and the annual cycle. The centroid of the WPWP traced an ellipselike trajectory during a year and moved counterclockwise in most years. However, in 1982 and 1986, the years of the onset of El Nino events, the movements were clockwise. The upper-layer volume of the WPWP was divided latitudinally into three sections. The annual cycles in the northern (from 3 deg to 30 deg N) and southern (from 3 deg to 30 deg S) sections were dominant. No annual cycle was found in the equatorial section (from 3 deg S to 3 deg N), but the volume of warm water in the equatorial Pacific increased during the 1986/87 El Nino event. The equatorial section was further divided into the eastern and western sectors along 165 deg W. During the 1986/87 El Nino event, the volume of warm water increased in the eastern sector, but the variation was smaller in the western sector than that in the eastern sector. During the 1988 La Nina event, the warm water volumes decreased in both sectors.

Description: This paper compares the observed behavior of the (F2) layer of the ionosphere at Millstone Hill and Hobart with calculations from the field line interhemispheric plasma (FLIP) model for solar maximum, solstice conditions in 1990. During the study period the daily F(sub 10.7) index varied by more than a factor of 2 (123 to 280), but the 81-day mean F(sub 10.7) (F(sub 10.7 A)) was almost constant near 190. Calculations were performed with and without the effects of vibrationally excited N2 (N(sup *)(sub 2) which affects the loss rate of atomic oxygen ions. In the case without N(sup *)(sub 2) there is generally good agreement between the model and measurement for the daytime, peak density of the F region (NmF2). Both the model and the measurement show a strong seasonal anomaly with the winter noon densities a factor of 3 to 4 greater than the summer noon densities at Millstone Hill and a factor of 2 greater at Hobart. The seasonal anomaly in the model is caused by changes in the neutral composition as given by the mass spectrometer and incoherent scatter (MSIS) 86 neutral density model. There is generally little or no increase in the observed noon NmF2 as a function of daily F(sub 10.7) except at Millstone Hill in winter. In contrast to the generally good agreement between model and data at noon, the model badly underestimates the density at night at Millstone Hill at all seasons. At Hobart the model reproduces the nighttime density variations well in both winter and summer. The international reference ionosphere (IRI) model generally provides a good representation of the average behavior of noon NmF2 and hmF2 but because the data show a lot of day-to-day variability, there are often large differences. The FLIP model is able to reproduce this variability when hmF2 is specified. The IRI model peak densities are better than the FLIP densities at night, but the IRI model does not represent the Millstone Hill summer data very well at night in 1990.

Description: Formulas are presented that parameterize the heating rate and coefficient of turbulent heat conduction produced by saturated internal gravity waves (IGW) in the upper atmosphere. Estimates of these values are made using observational data. The parameterization of IGW influences are introduced into a one-dimensional model of global mean thermal and composition balances of the upper atmosphere. Computations are performed for different values of IGW energy fluxes entering into the upper atmosphere from below. It is shown that realistic vertical profiles of the global mean temperature can be obtained using different values of IGW energy flux into the upper atmosphere. Increasing the IGW intensity leads not only to an increase of the heating rate due to wave enery dissipation, but also to an increase of the heating rate due to wave energy dissipation, but also to an increase in the coefficient of turbulent heat conduction and cooling rate produced by turbulence generated by the wave. So, near an altitude of 100 km the main part of solar heating is compensated by infrared cooling on one hand, and the main part of wave dissipation heating is compensated by turbulent cooling on the other hand. These quasi-balances generally hold for different values of IGW intensity.

Description: The usefulness of the radiances measured by operational satellites in deriving radiation budgets is demonstrated by comparing the model calculations with the Earth Radiation Budget Experiment (ERBE) fluxes. The radiation budgets in the atmosphere and at the surface in the western tropical Pacific are computed by coupling radiative transfer models to satellite retrievals of cloud and the earth surface parameters for April 1985 and April 1987. The model-computed fluxes at the top of the atmosphere agree well with the ERBE fluxes in both the solar and IR spectral regions. The difference is greater than 10 W/sq m in the outgoing longwave flux and greater than 15 W/sq m in the net downward shortwave flux. The agreement indicates that long-term earth radiation budgets in the tropical Pacific can be computed using the radiances measured by operational satellites.

Description: Two independent datasets for the solar radiation at the surface derived from satellites are compared. The data derived from the Earth Radiation Budget Experiment (ERBE) is for the net solar radiation at the surface whereas the International Satellite Cloud Climatology Project (ISCCP) data is for the downward flux only and was corrected with a space- and time-varying albedo. The ISCCP net flux is at all times higher than the ERBE flux. The difference can be divided into an offset that decreases with latitude and another component that correlates with high tropical cloud cover. With this latter exception the two datasets provide spatial patterns of solar flux that are very similar. A tropical Pacific Ocean model is forced with these two datasets and observed climatological winds. The upward heat flux is parameterized taking into account separately the longwave radiative, latent, and sensible heat fluxes. Best fit values for the uncertain parameters are found using an optimization procedure that seeks to minimize the difference between model and observed SST by varying the parameters within a reasonable range of uncertainty. The SST field the model produces with the best fit parameters is the best the model can do. If the differences between the model and data are larger than can be accounted for by remaining uncertainties in the heat flux parameterization and forcing data then the ocean model must be held to be at fault. Using this method of analysis, a fundamental model fault is identified. Inadequate treatment of mixed layer/entrainment processes in upwelling regions of the eastern tropical Pacific leads to a large and seasonally varying error in the model SST. Elsewhere the model SST is insufficiently different from observed to be able to identify model errors.

Description: In this paper a scheme is proposed to use a point raingage to compare contemporaneous measurements of rain rate from a single-field-of-view (FOV) estimate based on a satellite remote sensor such as a microwave radiometer. Even in the ideal case the measurements are different because one is at a point and the other is an area average over the field of view. Also the point gage will be located randomly inside the field of view on different overpasses. A space-time spectral formalism is combined with a simple stochastic rain field to find the mean-square deviations between the two systems. It is found that by combining about 60 visits of the satellite to the ground-truth site, the expected error can be reduced to about 10% of the standard deviation of the fluctuations of the systems alone. This seems to be a useful level of tolerance in terms of isolating and evaluating typical biases that might be contaminating retrieval algorithms.

Description: Much of the new record of broadband earth radiation budget satellite measurements to be obtained during the late 1990s and early twenty-first century will come from the dual-radiometer Clouds and Earth's Radiant Energy System Instrument (CERES-I) flown aboard sun-synchronous polar orbiters. Simulation studies conducted in this work for an early afternoon satellite orbit indicate that spatial root-mean-square (rms) sampling errors of instantaneous CERES-I shortwave flux estimates will range from about 8.5 to 14.0 W/m on a 2.5 deg latitude and longitude grid resolution. Rms errors in longwave flux estimates are only about 20% as large and range from 1.5 to 3.5 W/sq m. These results are based on an optimal cross-track scanner design that includes 50% footprint overlap to eliminate gaps in the top-of-the-atmosphere coverage, and a 'smallest' footprint size to increase the ratio in the number of observations lying within to the number of observations lying on grid area boundaries. Total instantaneous measurement error also depends on the variability of anisotropic reflectance and emission patterns and on retrieval methods used to generate target area fluxes. Three retrieval procedures from both CERES-I scanners (cross-track and rotating azimuth plane) are used. (1) The baseline Earth Radiaton Budget Experiment (ERBE) procedure, which assumes that errors due to the use of mean angular dependence models (ADMs) in the radiance-to-flux inversion process nearly cancel when averaged over grid areas. (2) To estimate N, instantaneous ADMs are estimated from the multiangular, collocated observations of the two scanners. These observed models replace the mean models in computation of satellite flux estimates. (3) The scene flux approach, conducts separate target-area retrievals for each ERBE scene category and combines their results using area weighting by scene type. The ERBE retrieval performs best when the simulated radiance field departs from the ERBE mean models by less than 10%. For larger perturbations, both the scene flux and collocation methods produce less error than the ERBE retrieval. The scene flux technique is preferable, however, because it involves fewer restrictive assumptions.

Description: Vertical heat fluxes associated with mesoscale circulations generated by land-surface wetness discontinuities are often stronger than turbulent fluxes, especially in the upper part of the atmospheric planetary boundary layer. As a result, they contribute significantly to the subgrid-scale fluxes in large-scale atmospheric models. Yet they are not considered in these models. To provide some insights into the possible parameterization of these fluxes in large-scale models, a state-of-the-art mesoscale numerical model was used to investigate the relationships between mesoscale heat fluxes and atmospheric and land-surface characteristics that play a key role in the generation of mesoscale circulations. The distribution of land-surface wetness, the wavenumber and the wavelength of the land-surface discontinuities, and the large-scale wind speed have a significant impact on the mesoscale heat fluxes. Empirical functions were derived to characterize the relationships between mesoscale heat fluxes and the spatial distribution of land-surface wetness. The strongest mesoscale heat fluxes were obtained for a wavelength of forcing corresponding approximately to the local Rossby deformation radius. The mesoscale heat fluxes are weakened by large-scale background winds but remain significant even with moderate winds.

Description: A three-dimensional cloud model, radiative transfer model-based simulation system is tested and validated against the aircraft-based radiance observations of an intense convective system in southeastern Virginia on 29 June 1986 during the Cooperative Huntsville Meteorological Experiment. NASA's ER-2, a high-altitude research aircraft with a complement of radiometers operating at 11-micrometer infrared channel and 18-, 37-, 92-, and 183-GHz microwave channels provided data for this study. The cloud model successfully simulated the cloud system with regard to aircraft- and radar-observed cloud-top heights and diameters and with regard to radar-observed reflectivity structure. For the simulation time found to correspond best with the aircraft- and radar-observed structure, brightness temperatures T(sub b) are simulated and compared with observations for all the microwave frequencies along with the 11-micrometer infrared channel. Radiance calculations at the various frequencies correspond well with the aircraft observations in the areas of deep convection. The clustering of 37-147-GHz T(sub b) observations and the isolation of the 18-GHz values over the convective cores are well simulated by the model. The radiative transfer model, in general, is able to simulate the observations reasonably well from 18 GHz through 174 GHz within all convective areas of the cloud system. When the aircraft-observed 18- and 37-GHz, and 90- and 174-GHz T(sub b) are plotted against each other, the relationships have a gradual difference in the slope due to the differences in the ice particle size in the convective and more stratiform areas of the cloud. The model is able to capture these differences observed by the aircraft. Brightness temperature-rain rate relationships compare reasonably well with the aircraft observations in terms of the slope of the relationship. The model calculations are also extended to select high-frequency channels at 220, 340, and 400 GHz to simulate the Millimeter-wave Imaging Radiometer aircraft instrument to be flown in the near future. All three of these frequencies are able to discriminate the convective and anvil portions of the system, providing useful information similar to that from the frequencies below 183 GHz but with potentially enhanced spatial resolution from a satellite platform. In thin clouds, the dominant effect of water vapor is seen at 174, 340, and 400 GHz. In thick cloudy areas, the scattering effect is dominant at 90 and 220 GHz, while the overlaying water vapor can attenuate at 174, 340, and 400 GHz. All frequencies (90-400 GHz) show strong signatures in the core.

Description: The evolution of the stratopsheric flow during the major stratospheric sudden warming of February 1979 is studied using two primitive equation models of the stratosphere and mesosphere. The United Kingdom Meteorological Office Stratosphere-Mesosphere Model (SMM) uses log pressure as a vertical coordinate. A spectral, entropy coordinate version of the SMM (entropy coordinate model, or ECM) that has recently been developed is also used. Comparison of SMM simulations with forecasts performed using the University of California, Los Angeles general circulation model confirms the previously noted sensitivity of stratospheric forecasts to tropospheric forecasts and emphasizes the importance of adequate vertical resolution in modeling the stratosphere. The ECM simulations provide a schematic description of the three-dimensional evolution of the polar vortex and the motion of air through it. During the warming, the two cyclonic vortices tilt westward and equatorward with height. Strong upward velocities develop in the lower stratosphere on the west (cold) side of a baroclinic zone as it forms over Europe and Asia. Strong downward velocities appear in the upper stratosphere on the east (warm) side, strengthening the temperature gradients. After the peak of the warming, vertical velocities decrease, downward velocities move into the lower stratosphere, and upward velocities move into the upper stratosphere.

Description: Extensive observations indicate a distinct difference between maritime and continental effective drop size, r(sub e), for warm clouds. The latest version of the National Center for Atmospheric Research (NCAR) Community Climate Model (CCM2) is used to explore the sensitivity of differentiating between continental and maritime r(sub e) on the simulated climate. The results of this study indicate that the smaller drop size over continents leads to a reduction of surface-absorbed solar radiation from 20 to 60 W/sq m. This reduction in surface solar flux leads to a cooling of the continents by up to -3.5 K. The reduction in surface solar flux and temperature also leads to a reduction in latent heat flux and precipitation over land. In the January simulation, there is a significant shift in tropical precipitation associated with the Australian monsoon. This shift leads to a response in the extra tropical geopotential height field over the western United States. All of these changes reduce biases in the current version of CCM2.

Description: The author reexamines the Charney-Drazin problem with special attention to the concentration of potential vorticity gradient in the neighborhood of the tropopause. It is found that the degree of concentration has a profound effect on the response to stationary forcing, with greater concentration leading to greater response. Smoothing the concentration either analytically or numerically (by using coarser resolution) both lead to reduced responses, especially at higher wavenumbers. The results suggest a potentially important interaction between baroclinically unstable eddies and stationary waves. Insofar as the former act to mix potential vorticity in the troposphere while concentrating gradients at tropopause levels, they significantly condition the basic state for the latter.

Description: In a cloud formed during adiabatic expansion, the droplet size distribution will be systematically related to the critical supersaturation of the cloud condensation nuclei (CNN), but this relationship can be complicated in entraining clouds. Useful information about cloud processes, such as mixing, can be obtained from direct measurements of the CNN involved in droplet nucleation. This was accomplished by interfacing two instruments for a series of flights in maritime cumulus clouds. One instrument, the counterflow virtual impactor, collected cloud droplets, and the nonvolatile residual nuclei of the droplets was then passed to a CCN spectrometer, which measured the critical supersaturation (S(sub c)) spectrum of the droplet nuclei. The measured S(sub c) spectra of the droplet nuclei were compared with the S(sub c) spectra of ambient aerosol particles in order to identify which CCN were actually incorporated into droplets and to determine when mixing processes were active at different cloud levels. The droplet nuclei nearly always exhibited lower median S(sub c)'s than the ambient aerosol, as expected since droplets nucleate perferentially on particles with lower critical supersaturations. Critical supersaturation spectra from nuclei of droplets near cloud base were similar to those predicted for cloud regions formed adiabatically, but spectra of droplet nuclei from middle cloud levels showed some evidence that mixing had occurred. Near cloud top, the greatest variation in the spectra of the droplet nuclei was observed, and nuclei with high S(sub c)'s were sometimes present even within relatively large droplets. This suggests that the extent of mixing increases with height in cumulus clouds and that inhomogeneous mixing may be important near cloud top. These promising initial results suggest improvements to the experimental technique that will permit more quantitative results in future experiments.

Description: This paper discusses an application of polarimetric measurements at vertical incidence. In particular, the correlation coefficients between linear copolar components are examined, and measurements obtained with the National Severe Storms Laboratory (NSSL)'s and National Center for Atmospheric Research (NCAR)'s polarimetric radars are presented. The data are from two well-defined bright bands. A sharp decrease of the correlation coefficient, confined to a height interval of a few hundred meters, marks the bottom of the bright band.

Description: The classic Eady problem is modified to include beta does not equal 0, but with the basic distributions of temperature and zonal flow adjusted to preserve zero meridional gradients of basic-state potential vorticity in the fluid interior. Much of the mathematical simplicity of the classic problem is retained; however, the results differ in important ways. Specifically, the instability now has a long-wave cutoff in addition to the traditional short-wave cutoff. The former is associated with the fact that the phase speeds of the edge waves begin to differ so much as wavenumber is reduced that the two edge waves can no longer interact in order to form unstable modes. For the unstable modes, this manifests itself in that the steering level for unstable modes is always below the middle of the fluid and approaches the lower boundary near the long-wave cutoff. Relatedly, the amplitude of the unstable geopotential perturbations is larger at the upper boundary than at the lower boundary. Finally, below the long-wave cutoff, one of the neutral waves has a phase speed that becomes increasingly easterly as wavenumber decreases. This allows a resonant response to planetary-scale stationary forcing.

Description: Dynamical mechanisms underlying the equilibration of absolute instability are examined in a nonlinear, quasigeostrophic, two-layer model. The key to understanding the nonlinear equilibration is in recognizing that linear absolute instabilities can be stabilized both by a reduction of the vertical shear and by enhancement of the mean barotropic velocity. In a localized domain, the equilibration process proceeds with the creation of locally convectively unstable regions downstream, which encroach onto the locally absolutely unstable region until the local instability is suppressed. That local instabilities exist only if absolutely unstable regions span a minimum size is verified by eigenvalue calculations of three-dimensional flows. Numerical examples suggest that this critical size is at least 9000 km for a wide range of parameter values chosen to investigate the midlatitude storm tracks. Fluctuations arising from local absolute instability obtain maximum amplitude in the downstream convectively unstable regions rather than in the absolutely unstable regions themselves. Together, these results suggest that if an equilibrated absolute instability were to occur in midlatitudes, a zonal band of surface easterlies exceeding 900 km would be required and the associated enhanced variances would not be found coincident with the regions of absolute instability.

Description: Future climate change could have significant repercussions for lightning-caused wildfires. Two empirical fire models are presented relating the frequency of lightning fires and the area burned by these fires to the effective precipitation and the frequency of thunderstorm activity. One model deals with the seasonal variations in lightning fires, while the second model deals with the interannual variations of lightning fires. These fire models are then used with the Goddard Institute for Space Studies General Circulation Model to investigate possible changes in fire frequency and area burned in a 2 X CO2 climate. In the United States, the annual mean number of lightning fires increases by 44%, while the area burned increases by 78%. On a global scale, the largest increase in lightning fires can be expected in untouched tropical ecosystems where few natural fires occur today.

Description: Simulated climates using numerical atmospheric general circulation models (GCMs) have been shown to be highly sensitive to the fraction of GCM grid area assumed to be wetted during rain events. The model hydrologic cycle and land-surface water and energy balance are influenced by the parameter bar-kappa, which is the dimensionless fractional wetted area for GCM grids. Hourly precipitation records for over 1700 precipitation stations within the contiguous United States are used to obtain observation-based estimates of fractional wetting that exhibit regional and seasonal variations. The spatial parameter bar-kappa is estimated from the temporal raingauge data using conditional probability relations. Monthly bar-kappa values are estimated for rectangular grid areas over the contiguous United States as defined by the Goddard Institute for Space Studies 4 deg x 5 deg GCM. A bias in the estimates is evident due to the unavoidably sparse raingauge network density, which causes some storms to go undetected by the network. This bias is corrected by deriving the probability of a storm escaping detection by the network. A Monte Carlo simulation study is also conducted that consists of synthetically generated storm arrivals over an artificial grid area. It is used to confirm the bar-kappa estimation procedure and to test the nature of the bias and its correction. These monthly fractional wetting estimates, based on the analysis of station precipitation data, provide an observational basis for assigning the influential parameter bar-kappa in GCM land-surface hydrology parameterizations.

Description: Cloud temperature, liquid water content, and vertical air velocity are all considered in evaluating the microphysical growth state of ice phase precipitation particles in the atmosphere. The large-scale observations taken together with in situ measurements indicated that the most prevalent growth condition for large ice particles in active convection is sublimation during riming, whereas the most prevalent growth condition in stratiform precipitation is vapor deposition. The large-scale electrical observations lend further support to the idea that particles warmed by riming into sublimation charge negatively and particles in vapor deposition charge positively in collisions with small ice particles.

Description: The amount of upward current provided to the ionosphere by a thunderstorm that appeared over the Kennedy Space Center (KSC) on July 11, 1978, is reexamined using an analytic equation that describes a bipolar thunderstorm's current contribution to the global circuit in terms of its generator current, lightning currents, the altitudes of its charge centers, and the conductivity profile of the atmosphere. Ground-based measurements, which were obtained from a network of electric field mills positioned at various distances from the thunderstorm, were used to characterize the electrical activity inside the thundercloud. The location of the lightning discharges, the type of lightning, and the amount of charge neutralized during this thunderstorm were computed through a least squares inversion of the measured changes in the electric fields following each lightning discharge. These measurements provided the information necessary to implement the analytic equation, and consequently, a time-averaged estimate of this thunderstorm's current contribution to the global circuit was calculated. From these results the amount of conduction current supplied to the ionosphere by this small thunderstorm was computed to be less than 25% of the time-averaged generator current that flowed between the two vertically displaced charge centers.

Description: Radiative fluxes and cloud forcings for the ocean areas of the Arctic are computed from the monthly cloud product of the International Satellite Cloud Climatology Project (ISCCP) for 1983-90. Spatially averaged short-wave fluxes are compared well with climatological values, while downwelling longwave fluxes are significantly lower. This is probably due to the fact that the ISCCP cloud amounts are underestimates. Top-of-the-atmosphere radiative fluxes are in excellent agreement with measurements from the Earth Radiation Budget Experiment (ERBE). Computed cloud forcings indicate that clouds have a warming effect at the surface and at the top of the atmosphere during winter and a cooling effect during summer. The net radiative effect of clouds is larger at the surface during winter but greater at the top of the atmosphere during summer. Overall the net radiative effect of clouds at the top of the atmosphere is one of cooling. This is in contrast to a previous result from ERBE data showing arctic cloud forcings have a net warming effect. Sensitivities to errors in input parameters are generally greater during winter with cloud amount being the most important paarameter. During summer the surface radiation balance is most sensitive to errors in the measurements of surface reflectance. The results are encouraging, but the estimated error of 20 W/sq m in surface net radiative fluxes is too large, given that estimates of the net radiative warming effect due to a doubling of CO2 are on the order of 4 W/sq m. Because it is difficult to determine the accuracy of results with existing in situ observations, it is recommended that the development of improved algorithms for the retrieval of surface radiative properties be accompanied by the simultaneous assembly of validation datasets.

Description: Rawinsonde data from tropical Pacific stations were examined for westward-propagating 3-6-day meridional wind oscillations in the troposphere and lower stratosphere, 1973-1992. Four types were identified from cross-spectrum and principal component analysis. (1) The dominant oscillation, near 250 mb, had a period slightly greater than 5 days, zonal wavenumber 4-6, and modified Rossby-gravity structure near the date line. (2) In the western Pacific lower troposphere there was broadband activity with short zonal scale, coupled to upper-tropospheric waves in NH summer. (3) In the central Pacific, during NH autumn, there was a well-defined approximately 4 1/2-day oscillation with maximum amplitude in the lower troposphere and baroclinic phase tilt above. The vertical structure suggested coupling to deep tropical convection; this interpretation was supported by correlation of meridional wind with antisymmetric outgoing longwave radiation. (4) In the stratosphere, Rossby-gravity waves had periods less than or equal to 4 days and zonal wavenumber 3-4. Unlike tropospheric waves, these disturbances were coherent in a shallow layer, largest in west phase of quasi-biennial oscillation (QBO) and annual cycle (NH winter-spring).

Description: A set of continuous, high-resolution atmospheric radon (Rn-222) concentration time series and radon soil flux measurements were acquired during the summer of 1990 at a micrometeorological tower site 13 km northwest of Schefferville, Quebec, Canada. The tower was located in a dry upland, open-canopy lichen-spruce woodland. For the period July 23 to August 1, 1990, the mean radon soil flux was 41.1 +/- 4.8 Bq m(exp -2)/h. Radon surface flux from the two end-member forest floor cover types (lichen mat and bare soil) were 38.8 +/- 5.1 and 61.8 +/- 15.6 Bq m(exp -2)/h, respectively. Average total forest canopy resistances computed using a simple 'flux box' model for radon exchange between the forest canopy and the overlying atmosphere range from 0.47 +/- 0.24 s cm(exp -1) to 2.65 +/- 1.61 cm(exp -1) for daytime hours (0900-1700 LT) and from 3.44 +/- 0.91 s cm(exp -1) to 10.55 +/- 7.16 s cm(exp -1) for nighttime hours (2000-0600) for the period July 23 to August 6, 1990. Continuous radon profiling of canopy atmospheres is a suitable approach for determining rates of biosphere/atmosphere trace gas exchange for remote field sites where daily equipment maintenance is not possible. where daily equipment maintenance is not possible.

Description: Atmospheric aerosols were collected in the boundary layer and free troposphere over continental and coastal subarctic regions of Canada during the July - August 1990 joint U.S.-Canadian Arctic Boundary Layer Expedition (ABLE) 3B/Northern Wetlands Study (NOWES). The samples were analyzed for the following water soluble species: sulfate, nitrate, ammonium, potassium, sodium, chloride, oxalate, methylsulfonate, and total amine nitrogen. Ammonium and sulfate were the major water soluble components of these aerosols. The nearly neutral (overall) chemical composition of summertime aerosol particles contrasts their strongly acidic wintertime composition. Aerosol samples were separated into several air mass categories and characterized in terms of chemical composition, associated mixing ratios of gaseous compounds, and meteorological parameters. The fundamental category represented particles associated with 'background' air masses. The summertime atmospheric aerosols in background air over the North American subarctic and Arctic regions were characterized by relatively small and spatially uniform mixing ratios of the measured species. These aerosol particles were aged to the extent that they had lost their primary source signature. The chemical profile of the background air aerosols was frequently modified by additions from biomass fire plumes, aged tropical marine air, and intrusions of upper tropospheric/lower stratospheric air. Aerosols in boundary layer background air over the boreal forest region of Quebec-Labrador had significantly larger mixing ratios of ammonium and sulfate relative to the Hudson Bay region. This may reflect infiltration of anthropogenic pollution or be due to natural emissions from this region.

Description: Observations from a transect extending 100 km inland during the Northern Wetlands Study (NOWES) in 1990 show that the sea breeze develops on approximately 25% of days during summer and may penetrate up to 100 km inland on occasions. The sea breeze exhibits a marked diurnal clockwise rotation as a result of the Coriolis effect along the unobstructed coastline. The marine advective effect is shown to depend on gradient wind direction. With northwesterly upper level flow the sea breeze tends to be northeasterly in direction and is associated with decreased temperatures and vapor pressure deficits (VPD). With southwesterly upper level flow the sea breeze tends to have a southeasterly direction and less effect on temperatures and VPD. This is attributed to shorter residence times of air parcels over water. For two cases, Colorado State University mesoscale model simulations show good agreement with surface wind observations and suggest that under northwesterly gradient flow, Bowen ratios are increased in the onshore flow along western James Bay, while during southwesterly gradient flow these effects are negligible. These results have implications for the interpretation of local climate, ecology, and hydrology as well as land-based and airborne turbulent flux measurements made during NOWES.

Description: Net ecosystem CO2 exchange was measured during the 1990 growing season (June to August) along a transect starting 10 km inland from James Bay and extending 100 km interior to Kinosheo Lake, Ontario. Sites were chosen in three distinct areas: a coastal fen, an interior fen, and a bog. For the most productive sites in the bog, net daily uptake rates reached a maximum of 2.5 g C-CO2 m(exp -2)/d with an area-weighted exchange of 0.3 g C-CO2 m(exp -2)/d near midsummer. This site was estimated to be a net carbon source of 9 g C-CO2 m(exp -2) to the atmosphere over a 153-day growing season. The interior fen was less productive on a daily basis with a net maximum uptake of 0.5 g C-CO2 m(exp -2)/d and with corresponding area-weighted uptake of 0.1 g C-CO2 m(exp -2)/d during midsummer. Early and late season release of carbon to the atmosphere resulted in a net loss of 21 g C-CO2 m(exp -2) over the growing season from this site. The coastal fen was the most productive site with uptake rates peaking near 1.7 g C-CO2 m(exp -2)/d which corresponded to an area-weighted uptake of 0.8 g C-CO2 m(exp -2)/d during midsummer and an estimated net uptake of 6 g C-CO2 m(exp -2) for the growing season. Associated with net CO2 exchange measurements, multispectral reflectance properties of the sites were measured over the growing season using portable radiometers. These properties were related to exchange rates with the goal of examining the potential for satellite remote sensing to monitor biosphere/atmosphere CO2 exchange in this biome. The normalized difference vegetation index (NDVI) computed from surface reflectance was correlated with net CO2 exchange for all sites with the exception of areas with large proportions of Sphagnum moss cover. These mosses have greater near-infrared reflectance than typical surrounding vegetation and may require special adjustment for regional exchange/remote sensing applications.

Description: As part of the Canadian Northern Wetlands Study (NOWES) measurements of methane flux were made at the Kinosheo Lake tower site for a 1-month period during the 1990 summer intensive. The measurements were made with a diode-laser-based methane sensor using the eddy correlation technique. Measurements of the methane fluxes were made at two levels, 5 or 18 m. Approximately 900 half-hour average methane flux measurements were obtained. Weak temporal and diurnal trends were observed in the data. Fluxes averaged over the study period showed an overall methane emission of 16 mg CH4 m(exp -2)/d with a daytime average of 20 mg CH4 m(exp -2)/d and a nighttime average of 9 mg CH4 m(exp -2)/d. The effect of emission footprint was evident in the data. A strong relationship between the daily average methane flux and wet bog temperature at 20-cm depth was observed.

Description: The estimation of the rate of net CO2 uptake of vegetated land surfaces is essential for studies of global carbon cycle. The present paper demonstrates the use of spectral reflectance data from satellite remote sensing to model net CO2 flux (NCF) of a tallgrass canopy at the Konza prairie, Kansas. A bidirectional reflectance canopy model was used to estimate seasonal changes in canopy leaf area index (LAI) from surface reflectances remotely sensed by SPOT 1 and Landsat 5 satellites. The radiation model was also coupled with leaf conductance-photosynthesis models to scale up stomatal conductance and NCF from individual leaves to canopy level according to radiation distribution inside the canopy. The satellite-data-driven model was able to closely simulate the seasonal change in LAI as well as the short-term variation of canopy LAI caused by the dry period during late July and early August in the area. Modeled canopy stomatal conductance (g(sub c)) and NCF agree with measurements within 0.16 cm/s and 0.28 mg m(exp -2)/s, respectively, during the growth season from late May to late August. In October both measured and modeled NCF turned to small negative values as canopy photosynthesis diminished and predicted LAI approached zero. In addition to data scatter, some of the differences between modeled and measured g(sub c) and NCF may be attributed to uncertainties in seasonal changes of plant physiological status that were not detected by satellite data; some of the differences were caused by inadequate description of the dependence of nighttime CO2 flux of soil respiration on near-surface turbulent mixing.

Description: Investigations of the kinematics of the lower stratospheric Arctic vortex during the winter of 1991-1992 using the contour advection with surgery technique reveal three distinct events in which there was substantial intrusion of midlatitude air into the vortex, in apparent contradiction of the view that the polar vortex constitutes an isolated air mass. Two of these events, in late January and mid-February, were well documented. They were predicted in high-resolution forecasts by the European Centre for Medium-Range Weather Forecasts, most clearly in experimental forecasts with reduced diffusion. Direct confirmation of the presence of the intrusions and of their calculated locations was provided by aerosol observations from the airborne differential absorption laser lidar aboard the NASA DC-8, taken as part of the second Airborne Arctic Stratospheric Expedition campaign; aerosol-rich air of midlatitude origin was seen in the expected position of the intrusions. The reality of the February event was also confirmed by in situ measurements from the NASA ER-2. Such events may be significant for the chemical processes taking place within the winter vortex. The intrusions were evidently related to the meteorology of the northern stratosphere during this winter and in particular to persistent tropospheric blocking over the northeastern Atlantic Ocean and western Europe and concomitant ridging into the lower stratospheric vortex in this region. Nevertheless, preliminary investigations have indicated that such events are not uncommon in other northern hemisphere winters, although no such events were found in the southern hemisphere during the Antarctic winter of 1987.

Description: Based on a derivation of the two-stream daytime-mean equations of radiative flux transfer, a method for computing the daytime-mean actinic fluxes in the absorbing and scattering vertically inhomogeneous atmosphere is suggested. The method applies direct daytime integration of the particular solutions of the two-stream approximations or the source functions. It is valid for any duration of period of averaging. The merit of the method is that the multiple scattering computation is carried out only once for the whole averaging period. It can be implemented with a number of widely used two-stream approximations. The method agrees with the results obtained with 200-point multiple scattering calculations. The method was also tested in runs with a 1-km cloud layer with optical depth of 10, as well as with aerosol background. Comparison of the results obtained for a cloud subdivided into 20 layers with those obtained for a one-layer cloud with the same optical parameters showed that direct integration of particular solutions possesses an 'analytical' accuracy. In the case of the source function interpolation, the actinic fluxes calculated above the one-layer and 20-layer clouds agreed within 1%-1.5%, while below the cloud they may differ up to 5% (in the worst case). The ways of enhancing the accuracy (in a 'two-stream sense') and computational efficiency of the method are discussed.

Description: Monthly mean 2.5 deg x 2.5 deg resolution 10-m height wind speeds from the Special Sensor Microwave/Imager (SSM/I) instrument and the European Centre for Medium-Range Weather Forecasts (ECMWF) forecast-analysis system are compared between 60 deg S and 60 deg N during 1988-91. The SSM/I data were uniformly processed while numerous changes were made to the ECMWF forecast-analysis system. The SSM/I measurements, which were compared with moored-buoy wind observations, were used as a reference dataset to evaluate the influence of the changes made to the ECMWF system upon the ECMWF surface wind speed over the ocean. A demonstrable yearly decrease of the difference between SSM/I and ECMWF wind speeds occurred in the 10 deg S-10 deg N region, including the 5 deg S-5 deg N zone of the Pacific Ocean, where nearly all of the variations occurred in the 160 deg E-160 deg W region. The apparent improvement of the ECMWF wind speed occurred at the same time as the yearly decrease of the equatorial Pacific SSM/I wind speed, which was associated with the natural transition from La Nina to El Nino conditions. In the 10 deg S-10 deg N tropical Atlantic, the ECMWF wind speed had a 4-yr trend, which was not expected nor was it duplicated with the SSM/I data. No yearly trend was found in the difference between SSM/I and ECMWF surface wind speeds in middle latitudes of the Northern and Southern Hemispheres. The magnitude of the differences between SSM/I and ECMWF was 0.4 m/s or 100% larger in the Northern than in the Southern Hemisphere extratropics. In two areas (Arabian Sea and North Atlantic Ocean) where ECMWF and SSM/I wind speeds were compared to ship measurements, the ship data had much better agreement with the ECMWF analyses compared to SSM/I data. In the 10 deg S-10 deg N area the difference between monthly standard deviations of the daily wind speeds dropped significantly from 1988 to 1989 but remained constant at about 30% for the remaining years.

Description: In order to compare the clear-sky greenhouse effect and cloud-radiative forcing from general circulation models with Earth Radiation Budget Experiments (ERBE) data, it is necessary to calculate the general circulation model (GCM) clear-sky radiative fluxes in a way consistent with ERBE. This study discusses problems associated with the available methods for clear-sky radiative flux computations in GCMs and proposes a new approach, which uses a statistical relationship between the grid cloud cover and the availability of ERBE clear-sky measurement, established from ERBE pixel data, to sample the model radiative fluxes. Calculations with version 2 of the National Center for Atmospheric Research (NCAR) Community Climate Model using observed sea surface temperature (SST) show good agreement of clear-sky sampling from the proposed method with ERBE sampling. It is also shown that large improvements are achieved in the spatial variability of the model clear-sky radiative fluxes over ocean, with reference to ERBE, by using the new clear-sky sampling method.

Description: A rapid analytical radiative transfer model for reflection of solar radiation in plane-parallel atmospheres is developed based on the Sobolev approach and the delta function transformation technique. A distinct advantage of this model over alternative two-stream solutions is that in addition to yielding the irradiance components, which turn out to be mathematically equivalent to the delta-Eddington approximation, the radiance field can also be expanded in a mathematically consistent fashion. Tests with the model against a more precise multistream discrete ordinate model over a wide range of input parameters demonstrate that the new approximate method typically produces average radiance differences of less than 5%, with worst average differences of approximately 10%-15%. By the same token, the computational speed of the new model is some tens to thousands times faster than that of the more precise model when its stream resolution is set to generate precise calculations.

Description: The release of latent heat of condensation is the largest internal energy source of the atmosphere. Latent heating is most significant during the precipitation process. Our knowledge of the distribution of precipitation is poor. It is only well observed within limited areas of the globe. Over the oceans, for example, it is known only to within a factor of two. Thus, there are strong scientific requirements for observations of precipitation from instruments on board a satellite. The Tropical Rainfall Measuring Mission (TRMM) will be the first satellite to measure rainfall with adequate accuracy and provide information about the vertical distribution of precipitation, not only in tropics and subtropics, but to plus or minus 35 degrees of latitude. There is a need for the continuity of rain observations for climate modeling purposes and to expand observations to cover more of the globe than simply plus or minus 35 degrees latitude.

Description: Many polar lows are generated at the boundary between sea ice and the ocean, in regions of large temperature gradients, where in situ observations are rare or nonexistent. Since satellite observations are frequent in high-latitude regions, they can be used to detect polar lows and track their propagation and evolution. The Special Sensor Microwave/Imager (SSM/I) providing estimates of surface wind speed, integrated cloud liquid water content, water vapor content, and precipitation size ice-scattering signal over the ocean; the Geosat radar altimeter measuring surface wind speed and significant wave height; and the TIROS-N Operational Vertical Sounder (TOVS) allowing the determination of temperature and humidity profiles in the atmosphere have been used in synergy for a specific case which occurred in the Norwegian Sea on January, 23-24 1988. All three instruments show sharp atmospheric gradients associated with the propagation of this low across the ocean, which permit the detection of the polar low at a very early stage and tracking it during its development, propagation, and decay. The wind speed gradients are measured with good qualitative agreement between the altimeter and SSM/I. TOVS retrieved fields prior to the formation of the low confirm the presence of an upper level trough, while during the mature phase baroclinicity can be observed in the 1000-500 hPa geopotential thicknesses.

Description: A summary of the results of a detailed study of the vertical structure of mesoscale gravity waves conducted during the Cooperative Convective Precipitation Experiment (CCOPE) is presented. Pressure perturbation fields derived from the Doppler wind fields are compared with the vertical structure of eigenfunctions resulting from a solution to the Taylor-Goldstein linear wave equation for an atmosphere whose mean state is described by vertical profiles obtained from a representative CCOPE sounding. An analysis of the potential for shear instability is also performed on all of the soundings taken on this day to assess the representativeness of the one chosen for the linear theoretical analysis.

Description: The seasonal cycle of low stratiform clouds is studied using data from surface-based cloud climatologies. The impact of low clouds on the radiation budget is illustrated by comparison of data from the Earth Radiation Budget Experiment with the cloud climatologies. Ten regions of active stratocumulus convection are identified. These regions fall into four categories: subtropical marine, midlatitude marine, Arctic stratus, and Chinese stratus. With the exception of the Chinese region, all the regions with high amounts of stratus clouds are over the oceans. In all regions except the Arctic, the season of maximum stratus corresponds to the season of greatest lower-troposphere static stability. Interannual variations in stratus cloud amount also are related to changes in static stability. A linear analysis indicates that a 6 percent increase in stratus fractional area coverage is associated with each 1 C increase in static stability. Over midlatitude oceans, sky-obscuring fog is a large component of the summertime stratus amount. The amount of fog appears to be related to warm advection across sharp gradients of SST.

Description: This study investigates ice nucleation mechanisms in cold lenticular wave clouds, a cloud type characterized by quasi-steady-state air motions and microphysical properties. It is concluded that homogeneous ice nucleation is responsible for the ice production in these clouds at temperatures below about -33 C. The lack of ice nucleation observed above -33 C indicates a dearth of ice-forming nuclei, and hence heterogeneous ice nucleation, in these clouds. Aircraft measurements in the temperature range -31 to -41 C show the following complement of simultaneous and abrupt changes in cloud properties that indicate a transition from the liquid phase to ice: disappearance of liquid water; decrease in relative humidity from near water saturation to ice saturation; increase in mean particle size; change in particle concentration; and change in temperature due to the release of latent heat. A numerical model of cloud particle growth and homogeneous ice nucleation is used to aid in interpretation of our in situ measurements. The abrupt changes in observed cloud properties compare favorably, both qualitatively and quantitatively, with results from the homogeneous ice nucleation model. It is shown that the homogeneous ice nucleation rates from the measurements are consistent with the temperature-dependent rates employed by the model (within a factor of 100, corresponding to about 1 C in temperature) in the temperature range -35 deg to -38 C. Given the theoretical basis of the modeled rates, it may be reasonable to apply them throughout the -30 to -50 C temperature range considered by the theory.

Description: Current research on the effects of snow covers on the atmosphere over short term, interannual, and decadal timescales is reviewed. Data obtained show that anomalies of snow cover are clearly associated with significant local anomalies of air temperature, at least in the lowest 100-200 mb. The duration of the local anomalies generally ranges from several days to several weeks and is often limited by the weekly and monthly fluctuation of the snow anomalies. The snow-hydrology-soil moisture feedback is capable of producing a meaningful response in the atmosphere during spring and summer.

Description: The effects of sea ice on short-term climate variability are studied by comparing the speculated and model-derived climatic roles of sea ice with the observational evidence. It is concluded that sea ice variations have large local impacts on the lower atmosphere which can be useful in local forecasting at the 30-90 day range. Observational studies are found to be severely constrained by the fact that atmospheric forcing plays a major role in the development of sea ice anomalies. Global climate models exhibit significant sensitivities to large changes of sea ice coverage. While the responses of different models are similar in many respects, there are inconsistencies among the models in the changes of some fundamental variables.

Description: The Goddard Cumulus Ensemble Model (GCEM) has been used to demonstrate that cumulus-scale dynamics and microphysics play a major role in determining the vertical distribution of water vapor and clouds in the tropical atmosphere. The GCEM is described and is the basic structure of cumulus convection. The long-term equilibrium response to tropical convection to surface warming is examined. A picture of the water cycle within tropical cumulus clusters is developed.

Description: The Project for Intercomparison of Land-surface Parameterization Schemes (PILPS) is described and the first stage science plan outlined. PILPS is a project designed to improve the parameterization of the continental surface, especially the hydrological, energy, momentum, and carbon exchanges with the atmosphere. The PILPS Science Plan incorporates enhanced documentation, comparison, and validation of continental surface parameterization schemes by community participation. Potential participants include code developers, code users, and those who can provide datasets for validation and who have expertise of value in this exercise. PILPS is an important activity because existing intercomparisons, although piecemeal, demonstrate that there are significant differences in the formulation of individual processes in the available land surface schemes. These differences are comparable to other recognized differences among current global climate models such as cloud and convection parameterizations. It is also clear that too few sensitivity studies have been undertaken with the result that there is not yet enough information to indicate which simplifications or omissions are important for the near-surface continental climate, hydrology, and biogeochemistry. PILPS emphasizes sensitivity studies with and intercomparisons of existing land surface codes and the development of areally extensive datasets for their testing and validation.

Description: We investigate the relationship between mixed Rossby-gravity waves (MRGWs) and convection in a general circulation model. The experiments described are performed in a general circulation model with the lower boundary set to that of an ocean surface everywhere. Several experiments are run varying the convective parameterization scheme (using either a modified Kuo scheme or a moist convective adjustment scheme) and varying the tropical sea surface temperatures (specified to be zonally symmetric in all cases), thereby changing the location of the modeled intertropical convergence zones (ITCZs). The appearance of a robust MRGW occurs when the sea surface temperature is such that two ITCZs straddle the equator. The particular sea surface temperature distribution used and the parameterization scheme for convection also affect the structure and strength of the modeled MRGW. The vertical structure of MRGWs is analyzed in the experiment in which this wave mode is the most energetic. We show that MRGWs of several different zonal length scales exist in the troposphere in association with convection; however, it is only the longer length scales which can be discerned in the upper troposphere and lower stratosphere.

Description: Results are presented of numerical experiments performed to determine what constitutes preconditioning for stratospheric sudden warmings. A number of 20-d runs were performed with a 3D primitive equation model. Runs varied in that they had different (zonally symmetric) initial fields and different lower boundary planetary wave fields, both based on NMC data for particular days. Initial and boundary data were used from the two minor and the sole major sudden warmings of January-February 1979 and some minor and the sole major sudden warmings of January-February 1989. To isolate the properties of the initial conditions from preconditioned cases and those with no warming, hybrid initial fields were created by combining data from two different days, one preconditioned and the other not. Runs with the hybrid initial conditions indicate that it is the lower stratospheric winds that are important in characterizing a preconditioned flow. It is suggested that it would be possible for a preconditioned state to exist for at least a short time without a warming occurring.

Description: WINDII, the Wind Imaging Interferometer on the Upper Atmosphere Research Satellite, began atmospheric observations on September 28, 1991 and since then has been collecting data on winds, temperatures and emissions rates from atomic, molecular and ionized oxygen species, as well as hydroxyl. The validation of winds and temperatures is not yet complete, and scientific interpretation has barely begun, but the dominant characteristic of these data so far is the remarkable structure in the emission rate from the excited species produced by the recombination of atomic oxygen. The latitudinal and temporal variability has been noted before by many others. In this preliminary report on WINDII results we draw attention to the dramatic longitudinal variations of planetary wave character in atomic oxygen concentration, as reflected in the OI 557.7 nm emission, and to similar variations seen in the Meine1 hydroxyl band emission.

Description: The maximum entropy production principle suggested by Paltridge (1975) is applied to separating the satellite-determined required total transports into atmospheric and oceanic components. Instead of using the excessively restrictive equal energy dissipation hypothesis as a deterministic tool for separating transports between the atmosphere and ocean fluids, the satellite-inferred required 2D energy transports are imposed on Paltridge's energy balance model, which is then solved as a variational problem using the equal energy dissipation hypothesis only to provide an initial guess field. It is suggested that Southern Ocean transports are weaker than previously reported. It is argued that a maximum entropy production principle can serve as a governing rule on macroscale global climate, and, in conjunction with conventional satellite measurements of the net radiation balance, provides a means to decompose atmosphere and ocean transports from the total transport field.

Description: The objective of the study was to evaluate the use of a satellite-derived vegetation index and surface temperature estimates for the assessment of the difference in urban and rural air temperature due to the urban heat island effect. The difference in the ND (normalized difference) index between urban and rural regions appears to be an indicator of the difference in surface properties (evaporation and heat storage capacity) between the two environments that are responsible for the urban heat island effect. The use of the approach proposed here may provide a globally consistent method for assessing this phenomenon.

Description: Estimates of snow grain size for the near-surface snow layer were calculated for the Tioga Pass region and Mammoth Mountain in the Sierra Nevada, California, using an inversion technique and data collected by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). The Tioga Pass and Mammoth Mountain single-band AVIRIS radiance images were atmospherically corrected to obtain surface reflectance. A discrete-ordinate model was used to calculate directional reflectance as a function of snowpack grain size for a wide range of snow grain radii. The resulting radius vs. reflectance curves were each fit using a nonlinear least squares technique which provided a means of transforming surface reflectance in each AVIRIS image to optically equivalent grain size on a per-pixel basis. The model results and grain size estimates derived from the AVIRIS data show that, for solar incidence angles between 0 and 30, the technique provides good estimates of grain size. This work provides the first quantitative estimates for grain size using data acquired from an airborne remote sensing instrument and is an important step in improving our ability to retrieve snow physical properties independent of field measurements.

Description: Three numerical experiments that investigate the scaling of land-surface processes - either of the inputs or parameters - are reported, and the aggregated processes are compared to the spatially variable case. The first is the aggregation of the hydrologic response in a catchment due to rainfall during a storm event and due to evaporative demands during interstorm periods. The second is the spatial and temporal aggregation of latent heat fluxes, as calculated from SiB. The third is the aggregation of remotely sensed land vegetation and latent and sensible heat fluxes using TM data from the FIFE experiment of 1987 in Kansas. In all three experiments it was found that the surface fluxes and land characteristics can be scaled, and that macroscale models based on effective parameters are sufficient to account for the small-scale heterogeneities investigated.

Description: A Doppler-radar-based case study of a monsoonal rainband (MRB) observed on the evening of January 12, 1990 near Darwin, Austrialia, is presented. Convective and mesoscale circulation features of this system are elucidated. The structure of the MRB is similar to the conceptual model of Houze (1989). The convergence profiles observed in the stratiform region differed considerably from the classical Houze model. Divergence and not convergence existed at heights of 3-6 km. It is suggested that in the MRB stratiform region, consistent with the observed rearward advection of decaying convective cells, the observed mesoscale updraft/downdraft configuration is a result of decaying convection and not a separately evolving mesoscale process.

Description: Simple numerical experiments are performed in order to determine the effects of inconsistent combinations of horizontal and vertical resolution in both atmospheric models and observing systems. In both cases, we find that inconsistent spatial resolution is associated with enhanced noise generation. A rather fine horizontal resolution in a satellite-data observing system seems to be excessive when combined with the usually available relatively coarse vertical resolution. Using horizontal filters of different strengths, adjusted in such a way as to render the effective horizontal resolution more consistent with vertical resolution for the observing system, may result in improvement of the analysis accuracy. The increase of vertical resolution for a satellite data observing system with better vertically resolved data, the results are different in that little or no horizontal filtering is needed to make spatial resolution more consistent for the system. The obtained experimental estimates of consistent vertical and effective horizontal resolution are in a general agreement with consistent resolution estimates previously derived theoretically by the authors.

Description: Airborne radiometric measurements at frequencies near 92 GHz and 183 GHz were conducted over two precipitation events near Wallops Island, Virginia during February, 1986. The measured brightness temperatures are compared with those from calculations to estimate the snowfall and rainfall rates for both events. The estimated rates over water surface are within a factor of two of those derived from the concurrent measurements by the SPANDAR radar at the NASA Wallops Flight Facility. These estimated rates, however, suggest that both snowfall events are light and close to the threshold of radiometric detection especially over land surface. Observations of additional snowfall events with lower frequency channels are needed to demonstrate the approach and to account for cloud effects.

Description: The International Satellite Cloud Climatology Project dataset is used to correlate variations of cloud optical thickness and cloud temperature in today's atmosphere. The analysis focuses on low clouds in order to limit the importance of changes in cloud vertical extent, particle size, and water phase. Coherent patterns of change are observed on several time and space scales. On the planetary scale, clouds in colder, higher latitudes are found to be optically thicker than clouds in warmer, lower latitudes. On the seasonal scale, winter clouds are, for the most part, optically thicker than summer clouds. The logarithmic derivative of cloud optical thickness with temperature is used to describe the sign and magnitude of the optical thickness-temperature correlation. The seasonal, latitudinal, and day-to-day variations of this relation are examined for Northern Hemisphere clouds in 1984. In cold continental clouds, optical thickness increases with temperature, consistent with the temperature variation of the adiabatic cloud water content. In warm continental and in almost all maritime clouds, however, optical thickness decreases with temperature.

Description: Airborne observations using a downward-looking, dual-frequency, near-infrared, differential absorption lidar system provide the first measurements of the height-dependent pressure-perturbation field associated with a strong mesoscale gravity wave. A pressure-perturbation amplitude of 3.5 mb was measured within the lowest 1.6 km of the atmosphere over a 52-km flight line. Corresponding vertical displacements of 250-500 m were inferred from lidar-observed displacement of aerosol layers. Accounting for probable wave orientation, a horizontal wavelength of about 40 km was estimated. Satellite observations reveal wave structure of a comparable scale in concurrent cirrus cloud fields over an extended area. Smaller-scale waves were also observed. Local meteorological soundings are analyzed to confirm the existence of a suitable wave duct. Potential wave-generation mechanisms are examined and discussed. The large pressure-perturbation wave is attributed to rapid amplification or possible wave breaking of a gravity wave as it propagated offshore and interacted with a very stable marine boundary layer capped by a strong shear layer.

Description: A 2D model of a mesoscale convective flow is used to simulate the excitation and vertical propagation of gravity waves. Data obtained show that, in the absence of storm-relative mean winds in the stratosphere, the primary mode of excitation of gravity waves is by mechanical forcing owing to oscillatory updrafts. The stratospheric response consists of waves whose periods match the primary periods of the forcing. Due to the tendency of the oscillating updrafts to propagate toward the rear of the storm, gravity wave propagation is limited primarily to the rearward direction. Results suggest that squall-line-generated gravity waves arise from mechanical forcing rather than thermal effects.

Description: Monthly mean residual circulations were calculated from eight years of satellite data. The diabatic circulation is usually found to give a good approximation to the residual circulation, but this is not always the case. In particular, an example is shown at 60 deg S and 30 mbar where the diabatic and residual circulations show very different annual variations. Correlations between the vertical component of the residual circulation and temperature and ozone were computed. The computations indicate that yearly variations of temperatures in the tropics are under radiative control, except during stratospheric warmings. Interannual variations in seasonal mean temperatures are shown to be under dynamical control everywhere. Correlations between seasonal means of the vertical component of the residual circulation and ozone mixing ratios are consistent with what would be expected from the ozone variations being due to differences in the ozone transport, although transport effects cannot easily be distinguished from photochemical effects above the altitude of the ozone mixing ratio peak. Finally, variations in total ozone are examined in comparison with residual circulation variations. A one to two month phase lag is seen in the annual variation in the total ozone at 60 deg N with respect to the maximum downward residual motions. This phase lag is greater at 60 deg N than at 60 deg S. There is evidence at 60 deg S of a greater downward trend in the mean zonal ozone maxima than there is in the minima. A decreasing trend in the maximum descending motion is seen to accompany the ozone trend at 60 deg S.

Description: This report investigates the impact of differential net radiative heating on 2D energy transports within the atmosphere ocean system and the role of clouds on this process. The 2D mean energy transports, in answer to zonal and meridional gradients in the net radiation field, show an east-west coupled dipole structure in which the Pacific acts as the major energy source and North Africa as the major energy sink. It is demonstrated that the dipole is embedded in the secondary energy transports arising mainly from the differential heating between land and oceans in the tropics in which the tropical east-west (zonal) transports are up to 30 percent of the tropical north-south (meridional) transports.

Description: Precipitation recycling is the contribution of evaporation within a region to precipitation in that same region. The recycling rate is a diagnostic measure of the potential for interactions between land surface hydrology and regional climate. In this paper we present a model for describing the seasonal and spatial variability of the recycling process. The precipitation recycling ratio, rho, is the basic variable in describing the recycling process. Rho is the fraction of precipitation at a certain location and time which is contributed by evaporation within the region under study. The recycling model is applied in studyiing the hydrologic cycle in the Amazon basin. It is estimated that about 25% of all the rain that falls in the Amazon basin is contributed by evaporation within the basin. This estimate is based on analysis of a data set supplied by the European Centre for Medium-range Weather Forecasts (ECMWF). The same analysis is repeated using a different data set from the Geophysical Fluid Dynamics Laboratory (GFDL). Based on this data set, the recycling ratio is estimated to be 35%. The seasonal variability of the recycling ratio is small compared with the yearly average. The new estimates of the recycling ratio are compared with results of previous studies, and the differences are explained.

Description: Satellite imagery datasets and regional climate model results are intercompared for evaluation of model accuracy in the simulation of cloud cover. Both monthly average individual simulation times are analyzed. To provide a consistent comparison, satellite data are first mapped into the model's geographic projection, grid domain, and resolution. It is found that September 1988 monthly average cloud fraction results from the modeled simulations correspond to observations, in both spatial pattern and magnitude, with bias less than +/- 20% cloud fraction over the entire inland West. Agreement in the pattern of cloud fraction also is evident for monthly average cloud fraction in July, but there is no negative bias of 10%-30% cloud fraction in the model diagnosis of cloud cover. Correlations between the spatial distributions of model-derived and observed cloud fractions are found to exceed 0.80 for certain geographic regions of the West, and these correlations are largest over mountainous areas during summer. Case studies of a series of daily cloud cover demonstrate the ability of the model to simulate the effects of frontal passage on cloud distribution. The ability of the RegCM1 to simulate daily cloud fraction and diurnal cloud evolution is somewhat weak for the summer convective season. It is anticipated that a more recent version of the regional climate model may improve the simulation of summer season cloud cover, through changes in cloud parameterization and improvements in model resolution.

Description: A mesoscale atmospheric model was used to evaluate the impact of subgrid-scale landscape discontinuities on the vertical profiles of resolved temperature, moisture, and moist static energy in the planetary boundary layer (PBL) of general circulation models (GCMs). These profiles were produced with a 3D version of the model (using a horizontal grid resolution of 7.5 km and 13 vertical layers in the PBL) by averaging horizontally the various atmospheric variables over a 180 x 180 sq km domain-about the size of the horizontal domain represented by a single grid element in a GCM. They were compared to corresponding vertical profiles produced with a 1D version of the model, which simulates the PBL, as in a GCM, over a single horizontal grid element. Differences obtained between the horizontally averaged atmospheric variables produced with the 3D simulations and the 1D simulations emphasize the impact of subgrid-scale landscape discontinuities on GCM-resolved variables. Various types of landscape discontinuities, characterized by horizontal contrasts of surface wetness and size of land patches, were simulated under various background-wind conditions. Differences of temperature, specific humidity, and moist static energy as large as 4 K, 6 g/kg, and 10 kJ/kg were obtained in some cases. These differences were not affected significantly by moderate winds but were sensitive to the spatial distribution of surface wetness. These results emphasize the need to parameterize mesoscale processes induced by landscape discontinuities in GCMss.

Description: The character of cloud-to-ground lightning is examined during the life cycle of a distinct mesoscale segment of the 10-11 June 1985 mesoscale convective system (MCS). Three phases of lightning activity are identified and related to both the radar-observed structure of the convection and to the severe weather produced by the MCS. Positive strikes to ground are dominant when the MCS is first developing. Negative strikes then dominate during a period of intense leading-line convective activity with high storm tops. Finally, a period of relatively frequent positive strikes within the trailing stratiform region occurs during the demise of the MCS. This last phase begins after the vertical extent of the leading convective line decreases rapidly and markedly, with moderate intensity echoes (i.e., 30-40 dBZ) occurring mostly below the freezing level. The first period of frequent positive flashes results from the lightning associated with a single severe thunderstorm in southwest Kansas; however, a second severe storm occurs nearby and produces mainly negative strikes. An extended period of strong surface winds does not appear to have any direct relationship with the observed character of the lightning activity.

Description: We have performed one dimensional radiative transfer calculations to evaluate the impact of cirrus clouds on the tropical radiation budget. We investigate the sensitivity of solar and infrared fluxes to cloud optical depth, particle size distributions, and cloud height. If the observed solar cloud forcing in excess of 100 W/sq m is to be attributed to cirrus anvils alone, then the optical depth of these anvils must be at least 5 (assuming 50% cloud cover and an ice crystal effective radius of 15 microns). The net radiative forcing of cirrus near the tropical tropopause is positive (heating) for cloud optical depths less than about 16 and negative (cooling) for larger optical depths. If cirrus clouds alone are responsible for the equal and opposite shortwave and longwave cloud forcing in excess of 100 W/sq m observed by Earth Radiation Budget Experiment (ERBE), then the cirrus must typically take the form of deep, optically thick clouds with relatively small particles (radii of 10-20 microns) and cloud-tops well below the tropopause. The maintenance of this balance on monthly time scales can be attributed to a variety of correlations: The cloud cover of optically thick cirrus or thin cirrus overlying low-level stratus clouds could vary; or cirrus anvil height cloud increase along with a decrease in the ice crystal effective radius and an increase in optical depth. It would be of great interest to determine observationally which of these correlations is responsible for the observed lack of variation in cloud forcing.

Description: We have investigated the sensitivity of the intensity of convective activity and atmospheric radiative cooling to radiatively thick upper-tropospheric clouds using a new version of the Colorado State University General Circulation Model (CSU GCM). The model includes a bulk cloud microphysics scheme to predict the formation of cloud water, cloud ice, rain, and snow. The cloud optical properties are interactive and dependent upon the cloud water and cloud ice paths. We find that the formation of a persistent upper tropospheric cloud ice shield leads to decreased atmospheric radiative cooling and increased static stability. Convective activity is then strongly suppressed. In this way, upper-tropospheric clouds act as regulators of the global hydrologic cycle, and provide a negative feedback between atmospheric radiative cooling and convective activity.

Description: We have investigated the processes that control ice crystal nucleation in the upper troposphere using a numerical model. Nucleation of ice resulting from cooling was simulated for a range of aerosol number densities, initial temperatures, and cooling rates. In contrast to observations of stratus clouds, we find that the number of ice crystals that nucleate in cirrus is relatively insensitive to the number of aerosols present. The ice crystal size distribution at the end of the nucleation process is unaffected by the assumed initial aerosol number density. Essentially, nucleation continues until enough ice crystals are present such that their deposition growth rapidly depletes the vapor and shuts off any further nucleation. However, the number of ice crystals nucleated increases rapidly with decreasing initial temperature and increasing cooling rate. This temperature dependence alone could explain the large ice crystal number density observed in very cold tropical cirrus.

Description: Monte Carlo radiative transfer methods are employed here to estimate the plane-parallel albedo bias for marine stratocumulus clouds. This is the bias in estimates of the mesoscale-average albedo, which arises from the assumption that cloud liquid water is uniformly distributed. The authors compare such estimates with those based on a more realistic distribution generated from a fractal model of marine stratocumulus clouds belonging to the class of 'bounded cascade' models. In this model the cloud top and base are fixed, so that all variations in cloud shape are ignored. The model generates random variations in liquid water along a single horizontal direction, forming fractal cloud streets while conserving the total liquid water in the cloud field. The model reproduces the mean, variance, and skewness of the vertically integrated cloud liquid water, as well as its observed wavenumber spectrum, which is approximately a power law. The Monte Carlo method keeps track of the three-dimensional paths solar photons take through the cloud field, using a vectorized implementation of a direct technique. The simplifications in the cloud field studied here allow the computations to be accelerated. The Monte Carlo results are compared to those of the independent pixel approximation, which neglects net horizontal photon transport. Differences between the Monte Carlo and independent pixel estimates of the mesoscale-average albedo are on the order of 1% for conservative scattering, while the plane-parallel bias itself is an order of magnitude larger. As cloud absorption increases, the independent pixel approximation agrees even more closely with the Monte Carlo estimates. This result holds for a wide range of sun angles and aspect ratios. Thus, horizontal photon transport can be safely neglected in estimates of the area-average flux for such cloud models. This result relies on the rapid falloff of the wavenumber spectrum of stratocumulus, which ensures that the smaller-scale variability, where the radiative transfer is more three-dimensional, contributes less to the plane-parallel albedo bias than the larger scales, which are more variable. The lack of significant three-dimensional effects also relies on the assumption of a relatively simple geometry. Even with these assumptions, the independent pixel approximation is accurate only for fluxes averaged over large horizontal areas, many photon mean free paths in diameter, and not for local radiance values, which depend strongly on the interaction between neighboring cloud elements.

Description: A systematic approach is suggested for modeling the probability distribution of rain rate. Rain rate, conditional on rain and averaged over a region, is modeled as a temporally homogeneous diffusion process with appropiate boundary conditions. The approach requires a drift coefficient-conditional average instantaneous rate of change of rain intensity-as well as a diffusion coefficient-the conditional average magnitude of the rate of growth and decay of rain rate about its drift. Under certain assumptions on the drift and diffusion coefficients compatible with rain rate, a new parametric family-containing the lognormal distribution-is obtained for the continuous part of the stationary limit probability distribution. The family is fitted to tropical rainfall from Darwin and Florida, and it is found that the lognormal distribution provides adequate fits as compared with other members of the family and also with the gamma distribution.

Description: A physically based algorithm sensitive to emission and scattering is used to estimate rainfall using the Special Sensor Microwave/Imager (SSM/I). The algorithm is derived from radiative transfer calculations through an atmospheric cloud model specifying vertical distributions of ice and liquid hydrometeors as a function of rain rate. The algorithm is structured in two parts: SSM/I brightness temperatures are screened to detect rainfall and are then used in rain-rate calculation. The screening process distinguishes between nonraining background conditions and emission and scattering associated with hydrometeors. Thermometric temperature and polarization thresholds determined from the radiative transfer calculations are used to detect rain, whereas the rain-rate calculation is based on a linear function fit to a linear combination of channels. Separate calculations for ocean and land account for different background conditions. The rain-rate calculation is constructed to respond to both emission and scattering, reduce extraneous atmospheric and surface effects, and to correct for beam filling. The resulting SSM/I rain-rate estimates are compared to three precipitation radars as well as to a dynamically simulated rainfall event. Global estimates from the SSM/I algorithm are also compared to continental and shipboard measurements over a 4-month period. The algorithm is found to accurately describe both localized instantaneous rainfall events and global monthly patterns over both land and ovean. Over land the 4-month mean difference between SSM/I and the Global Precipitation Climatology Center continental rain gauge database is less than 10%. Over the ocean, the mean difference between SSM/I and the Legates and Willmott global shipboard rain gauge climatology is less than 20%.

Description: Estimates of monthly rainfall have been computed over the tropical Pacific using passive microwave satellite observations from the Special Sensor Microwave/Imager (SSM/I) for the preiod from July 1987 through December 1991. The monthly estimates were calibrated using measurements from a network of Pacific atoll rain gauges and compared to other satellite-based rainfall estimation techniques. Based on these monthly estimates, an analysis of the variability of large-scale features over intraseasonal to interannual timescales has been performed. While the major precipitation features as well as the seasonal variability distributions show good agreement with expected values, the presence of a moderately intense El Nino during 1986-87 and an intense La Nina during 1988-89 highlights this time period.

Description: The temperature dependence of the ClONO2 absorption spectrum has been measured between 220 and 298 K and between 195 and 430 nm using a diode array spectrometer. The absorption cross sections were determined using both: (1) absolute pressure measurements at 296 K and (2) measurements at various temperatures relative to 296 K using a dual absorption cell arrangement. The temperature dependence of the ClONO2 absorption spectrum shows very broad structure. The amplitude of the temperature dependence relative to that at 296 K is weak at short wavelengths, less than 2% at 215 nm and 220 K, but significant at the wavelengths important in the stratosphere, about 30% at 325 nm and 220 K. Our ClONO2 absorption cross section data are in good general agreement with the previous measurements of Molina and Molina (1979).

Description: This paper presents the methodologies and results of the multivariate modeling and two-dimensional spectral and correlation analysis of PRE-STORM rainfall gauge data. Estimated parameters of the models for the specific spatial averages clearly indicate the eastward and southeastward wave propagation of rainfall fluctuations. A relationship between the coefficients of the diffusion equation and the parameters of the stochastic model of rainfall fluctuations is derived that leads directly to the exclusive use of rainfall data to estimate advection speed (about 12 m/s) as well as other coefficients of the diffusion equation of the corresponding fields. The statistical methodology developed here can be used for confirmation of physical models by comparison of the corresponding second-moment statistics of the observed and simulated data, for generating multiple samples of any size, for solving the inverse problem of the hydrodynamic equations, and for application in some other areas of meteorological and climatological data analysis and modeling.

Description: Satellite data from the microwave sounding unit (MSU) channel 4, when carefully merged, provide daily zonal anomalies of lower-stratosphere temperature with a level of precision between 0.01 and 0.08 C per 2.5 deg latitude band. Global averages of these daily zonal anomalies reveal the prominent warming events due to volcanic aerosol in 1982 (El Chichon) and 1991 (Mt. Pinatubo), which are on the order of 1 C. The quasibiennial oscillation (QBO) may be extracted from these zonal data by applying a spatial filter between 15 deg N and 15 deg S latitude, which resembles the meridional curvature. Previously published relationships between the QBO and the north polar stratospheric temperatures during northern winter are examined but were not found to be reproduced in the MSU4 data. Sudden stratospheric warmings in the north polar region are represented in the MSU4 data for latitudes poleward of 70 deg N. In the Southern Hemisphere, there appears to be a moderate relationship between total ozone concentration and MSU4 temperatures, though it has been less apparent in 1991 and 1992. In terms of empirical modes of variability, the authors find a strong tendency in EOF 1 (39.2% of the variance) for anomalies in the Northern Hemisphere polar regions to be counterbalanced by anomalies equatorward of 40 deg N and 40 deg S latitudes. In addition, most of the modes revealed significant power in the 15-20 day period band.

Description: A detailed ice-phase bulk microphysical scheme has been developed for simulating the hydrometeor distributions of convective and stratiform precipitation in different large-scale environmental conditions. The proposed scheme involves 90 distinct microphysical processes, which predict the mixing ratios and the number concentrations of small ice crystals, snow, graupel, and frozen drops/hail, as well as the mixing ratios of liquid water on wet precipitation ice (snow, graupel, frozen drops). The number of adjustable coefficients has been significantly reduced in comparison with other bulk schemes. Additional improvements have been made to the parameterization in the following areas: (1) representing small ice crystals with nonzero terminal fall velocities and dispersive size distributions, (2) accurate and computationally efficient calculations of precipitation collection processes, (3) reformulating the collection equation to prevent unrealistically large accretion rates, (4) more realistic conversion by riming between different classes of precipitation ice, (5) preventing unrealistically large rates of raindrop freezing and freezing of liquid water on ice, (6) detailed treatment of various rime-splintering ice multiplication mechanisms, (7) a simple representation of the Hobbs-Rangno ice enhancement process, (8) aggregation of small ice crystals and snow, and (9) allowing explicit competition between cloud water condensation and ice deposition rates rather than using saturation adjustment techniques. For the purposes of conserving the higher moments of the particle distributions, preserving the spectral widths (or slopes) of the particle spectra is shown to be more important than strict conservation of particle number concentration when parameterizing changes in ice-particle number concentrations due to melting, vapor transfer processes (sublimation of dry ice, evaporation from wet ice), and conversion between different hydrometeor species. The microphysical scheme is incorporated into a nonhydrostatic cloud model in Part 2 of this study. The model performed well in simulating the radar and microphysical structures of a midlatitude-continental squall lines and a tropical-maritime squall system with minimal tuning of the parameterization, even though the vertical profiles of radar reflectivity differed substantially between these storms.

Description: A model of net primary production (NPP), decomposition, and nitrogen cycling in tundra ecosystems has been developed. The adjoint technique is used to study the sensitivity of the computed annual net CO2 flux to perturbation in initial conditions, climatic inputs, and model's main parameters describing current seasonal CO2 exchange in wet sedge tundra at Barrow, Alaska. The results show that net CO2 flux is most sensitive to parameters characterizing litter chemical composition and more sensitive to decomposition parameters than to NPP parameters. This underlines the fact that in nutrient-limited ecosystems, decomposition drives net CO2 exchange by controlling mineralization of main nutrients. The results also indicate that the short-term (1 year) response of wet sedge tundra to CO2-induced warming is a significant increase in CO2 emission, creating a positive feedback to atmosphreic CO2 accumulation. However, a cloudiness increase during the same year can severely alter this response and lead to either a slight decrease or a strong increase in emitted CO2, depending on its exact timing. These results demonstrate that the adjoint method is well suited to study systems encountering regime changes, as a single run of the adjoint model provides sensitivities of the net CO2 flux to perturbations in all parameters and variables at any time of the year. Moreover, it is shown that large errors due to the presence of thresholds can be avoided by first delimiting the range of applicability of the adjoint results.

Description: The NASA/Goddard three-dimensional chemistry and transport model is driven by winds from a stratospheric data assimilation system. Synoptic- and planetary-scale patterns, apparent in satellite observations of trace constituents, are successfully reproduced for seasonal integrations. As model integrations proceed, however, the quality of simulations decreases, and systematic differences between calculation and measurement appear. The differences are explained by examining the zonal-mean residual circulation. The vertical velocity w-bar (sup star) is calculated two ways: (1) from the diabatic heating rates and temperature tendency and (2) from the Eulerian vertical velocity and the horizontal eddy heat flux convergence. The results from these calculations differ substantially. Periodic insertion of observational data during the assimilation process continually shocks the general circulation model and produces these differences, which leads to an overestimate of the mean vertical heat and cconstituent transport. Such differences are expected to be general to all data assimilation products. This interpretation is corroborated by two-dimensional (2D) model calculations. When w-bar(sup star) is calculated from (2), the 2D ozone evolution is unrealistic and qualitatively similar to the 3D model simulation. The 2D ozone evolution is reasonable when w-bar (sup star) from (1).

Description: Total Ozone Mapping Spectrometer (TOMS) images of the springtime Southern Hemisphere commonly show concentric layers in the total ozone field outside the Antarctic polar vortex. The layering appears to result from horizontal folding and stretching of regions on the equatorward flank of the polar vortex near the midlatitude ozone maximum. This folding and stretching interleaves low and high ozone air from the subtropics and midlatitudes, respectively. Occaisional large amplitude wave events can extract very low ozone air from the interior of the polar vortex (the Antarctic ozone hole), but the folding and stretching occurs in midlatitiudes even when wave amplitudes are not exceptionally large. The folding and stretching results in relatively rapid horizontal mixing of the atmosphere on the equatorward flank of the jet. This type of Lagrangian behavior may be common in the atmosphere, but is visible when local tracer gradients are large and observations with high spatial resolution are available. Also, experimentation has shown that gray-scale images of TOMS data show the details of the spatial distribution of ozone much more clearly than contour maps or false-color images.

Description: In the surface hydrologic parameterization of general circulation models (GCMs), it is commonly assumed that the precipitation processes are homogeneous over a GCM grid square and that the precipitation intensity is uniformly distributed. Based on evidence that the spatial distribution of precipitation within a GCM grid square is crucial for the land surface hydrology parameterization, a few researchers have explored the impacts of assuming that the precipitation is exponentially distributed. This paper explores the suitability of the aforementioned assumptions. First, a statistical analysis is conducted of historical precipitation data for three GCM grids in different regions of the United States. The analysis suggests that neither the uniform nor the exponential distribution assumption may be suitable at the GCM grid scale and, that instead, the spatial variability in precipitation is characterized by statistical patterns that are inhomogeneous. These patterns vary from grid to grid and are induced by the interaction between atmospheric conditions and various land surface characteristics, such as topographical features, surface properties, etc. Within the same grid square, however, the statistical patterns are generally constant from year to year. Based on this analysis, a computationally viable (i.e., usable with GCMs) stochastic precipitation disaggregation scheme that utilizes these stable statistical patterns is proposed. The method was used to generate spatially distributed hourly rainfall for a summer season in the southwestern region of the continental United States. Analysis of the results shows that the methodology preserves the seasonal characteristics of spatial variability in precipitation that is observed in the long-term historical data.

Description: Radiative transfer simulations are performed to determine how water vapor and nonprecipitating cloud liquid water and ice particles within typical midlatitude atmospheres affect brightness temperatures T(sub B)'s of moisture sounding channels used in the Advanced Microwave Sounding Unit (AMSU) and AMSU-like instruments. The purpose is to promote a general understanding of passive top-of-atmosphere T(sub B)'s for window frequencies at 23.8, 89.0, and 157.0 GHz, and water vapor frequencies at 176.31, 180.31, and 182.31 GHz by documenting specific examples. This is accomplished through detailed analyses of T(sub B)'s for idealized atmospheres, mostly representing temperate conditions over land. Cloud effects are considered in terms of five basic properties: droplet size distribution, phase, liquid or ice water content, altitude, and thickness. Effects on T(sub B) of changing surface emissivity also are addressed. The brightness temperature contribution functions are presented as an aid to physically interpreting AMSU T(sub B)'s. Both liquid and ice clouds impact the T(sub B)'s in a variety of ways. The T(sub B)'s at 23.8 and 89 GHz are more strongly affected by altostratus liquid clouds than by cirrus clouds for equivalent water paths. In contrast, channels near 157 and 183 GHz are more strongly affected by ice clouds. Higher clouds have a greater impact on 157- and 183-GHz T(sub B)'s than do lower clouds. Clouds depress T(sub B)'s of the higher-frequency channels by suppressing, but not necessarily obscuring, radiance contributions from below. Thus, T(sub B)'s are less closely associated with cloud-top temperatures than are IR radiometric temperatures. Water vapor alone accounts for up to 89% of the total attenuation by a midtropospheric liquid cloud for channels near 183 GHz. The Rayleigh approximation is found to be adequate for typical droplet size distributions; however, Mie scattering effects from liquid droplets become important for droplet size distribution functions with modal radii greater than 20 micrometers near 157 and 183 GHz, and greater than 30-40 micrometers at 89 GHz. This is due mainly to the relatively small concentrations of droplets much larger than the mode radius. Orographic clouds and tropical cumuli have been observed to contain droplet size distributions with mode radii in the 30-40 micrometers range. Thus, as new instruments bridge the gap between microwave and infrared to frequencies even higher than 183 GHz, radiative transfer modelers are cautioned to explicitly address scattering characteristics of such clouds.

Description: High-latitude environments are thought to play several critical roles in the global balance of radiatively active trace gases. Adequate documentation of the source and sink strengths for trace gases requires long time series of detailed measurements, including heat and moisture budgets. A fen near Schefferville, Quebec, was instrumented during the summer of 1990 for the measurement of the surface energy, radiation, and moisture balances as well as for eddy correlation estimates of ozone and methane flux. Despite the limited fetch at this site, analysis of the tower flux 'footprint' indicates that at least 80% of the flux observed originates from sources within the fen. Sensible heat fluxes averaged 25% of the daytime net radiation at the site, while the latent heat flux, determined from the energy balance, was 63%; the Bowen ratio varied from 0.2 to 0.8 from day to day, without a seasonal trend to the variation. The competing effects of rooted macrophyte development (with concomitant effects on roughness and transpiration) and the normal shift in synoptic pattern around day 200 to warm, dry conditions results in a lack of net seasonal effect on the energy partitioning. Over the period from days 170 to 230, the evaporation (167 mm) was double the rainfall, while the decline in water level was 107 mm, leaving a net runoff of 0.44 mm/d. The total hydrocarbon flux was 75-120 mg m(exp -2)/d, following a diurnal pattern similar to heat or moisture flux, while the daytime ozone flux was about -1.11 x 10(exp 11) molecules cm(exp -2)/s. A period near the end of the experiment, during week 30, produced the strongest total hydrocarbon flux, associated with warmer deep (1 m) soil temperatures, lower fen water levels, and the late summer shift in wind direction at that time. An early summer 'flush' of total hydrocarbon was not observed.

Description: Biomass-burning impacted air masses sampled over central and eastern Canada during the summer of 1990 as part of ABLE 3B contained enhanced mixing ratios of gaseous HNO3, HCOOH, CH3COOH, and what appears to be (COOH)2. These aircraft-based samples were collected from a variety of fresh burning plumes and more aged haze layers from different source regions. Values of the enhancement factor, delta X/delta CO, where X represents an acidic gas, for combustion-impacted air masses sampled both near and farther away from the fires, were relatively uniform. However, comparison of carboxylic acid emission ratios measured in laboratory fires to field plume enhancement factors indicates significant in-plume production of HCOOH. Biomass-burning appears to be an important source of HNO3, HCOOH, and CH3COOH to the troposphere over subarctic Canada.

Description: The Arctic Boundary Layer Expedition (ABLE) 3B used data from ground-based, aircraft, and satellite platforms to characterize the chemistry and dynamics of the troposphere in subarctic and Arctic regions of midcontinent and eastern Canada during July - August 1990. This paper reports the experimental design for ABLE 3B and a brief overview of results. The detailed results are presented in a series of papers in this issue. The chemical composition of the atmospheric mixed layer over remote tundra, boreal wetland, and forested environments was influenced by emissions of CH4 and nonmethane hydrocarbons from biogenic sources, emissions of gases and aerosols from local biomass burning, and transport of pollutants into the study areas from urban/industrial sources. Minimum concentrations of both trace gas and aerosol species in boundary layer air were associated with Arctic source areas. In the free troposphere the biospheric influence was undetectable, and major sources of chemical variability were related to long-range transport of pollutants into the study areas from biomass burning and industrial sources in Alaska and the Great Lakes regions, respectively. Minimum concentrations of both trace gas and aerosol species in the free troposphere were associated with a persistent, widespread air mass which both chemistry and air mass trajectory analyses suggested had originated in the tropical Pacific. Subsidence of air from the upper troposphere and lower stratosphere frequently enhanced ozone and influenced other trace gas and aerosol species at midtropospheric altitudes. The North American Arctic is a complex dynamical and chemical environment with considerable spatial and temporal variability in aerosol and trace gas concentrations. The use of atmospheric chemical indicators for climate change detection will require a much more comprehensive Arctic monitoring program than currently exists.

Description: Peat chemistry appears to exert primary control over methane production rates in the Canadian Northern Wetlands Study (NOWES) area. We determined laboratory methane production rate potentials in anaerobic slurries of samples collected from a transect of sites through the NOWES study area. We related methane production rates to indicators of resistance to microbial decay (peat C: N and lignin: N ratios) and experimentally manipulated substrate availability for methanogenesis using ethanol (EtOH) and plant litter. We also determined responses of methane production to pH and temperature. Methane production potentials declined along the gradient of sites from high rates in the coastal fens to low rates in the interior bogs and were generally highest in surface layers. Strong relationships between CH4 production potentials and peat chemistry suggested that methanogenesis was limited by fermentation rates. Methane production at ambient pH responded strongly to substrate additions in the circumneutral fens with narrow lignin: N and C: N ratios (delta CH4/delta EtOH = 0.9-2.3 mg/g) and weakly in the acidic bogs with wide C: N and lignin: N ratios (delta CH4/delta EtOH = -0.04-0.02 mg/g). Observed Q(sub 10) values ranged from 1.7 to 4.7 and generally increased with increasing substrate availability, suggesting that fermentation rates were limiting. Titration experiments generally demonstrated inhibition of methanogenesis by low pH. Our results suggest that the low rates of methane emission observed in interior bogs during NOWES likely resulted from pH and substrate quality limitation of the fermentation step in methane production and thus reflect intrinsically low methane production potentials. Low methane emission rates observed during NOWES will likely be observed in other northern wetland regions with similar vegetation chemistry.

Description: Ponds on peatlands of the Hudson Bay lowlands (HBLs) are complex ecosystems in which the fluxes to the atmosphere of CH4 and CO2 were controlled by interacting physical and biological factors. This resulted in strong diel variations of both dissolved gas concentrations and gas fluxes to the atmosphere, necessitating frequent sampling on a 24-hour schedule to enable accurate estimates of daily fluxes. Ponds at three sites on the HBL were constant net sources of CH4 and CO2 to the atmosphere at mean rates of 110-180 mg CH4 m(exp -2)/d and 3700-11,000 mg CO2 m(exp -2)/d. Rates peaked in August and September. For CH4 the pond fluxes were 3-30 times higher than adjacent vegetated surfaces. For CO2 the net pond fluxes were similar in magnitude to the vegetated fluxes but the direction of the flux was opposite, toward atmosphere. Even though ponds cover only 8-12% of the HBL area, they accounted for 30% of its total CH4 flux to the atmosphere. There is some circumstantial evidence that the ponds are being formed by decomposition of the underlying peat and that this decomposition is being stimulated by the activity of N2 fixing cyanobacteria that grow in mats at the peat-water interface. The fact that the gas fluxes from the ponds were so different from the surrounding vegetated surfaces means that any change in the ratio of pond to vegetated area, as may occur in response to climate change, would affect the total HBL fluxes.

Description: This study examines 1211 cases of coincident ozone profiles derived from 1164 Umkehrs and 928 Stratospheric Aerosol and Gas Experiment II (SAGE II) profiles within 1000 km and 12 hours between October 1984 and April 1989 to study the stratospheric-aerosol effect on Umkehr ozone profiles. Because of the close correspondence of stratospheric aerosol optical depth at the SAGE II-measured 0.525-micrometer wavelength and the extrapolated 0.32 Umkehr wavelength determined in this study we use the 0.525-micrometer data to determine the aerosol effect on Umkehr profiles. At the 95% confidence level, we find the following errors to the Umkehr ozone amounts: in Umkehr layer 9 (-2.9 +/- 2.1), layer 8 (-2.3 +/- 1.1), layer 7 (0.1 +/- 1.1), layer 6 (2.2 +/- 1.0), layer 5 (-1.5 +/- 0.8), and layer 4 (-2.4 +/- 1.7) in percent ozone amount per 0.01 stratospheric aerosol optical depth. These results agree with previous theoretical and empirical studies within their respective error bounds in layers 9, 8, and 7. The results in layers 6, 5, and 4 differ significantly from those in previous works. Using only those eight stations with more than 47 coincidences results in mean aerosol effects that are not significantly different from the 14-station results. Because SAGE II and Umkehr produce different ozone retrievals in layer 9 and because the intralayer correlation of SAGE II ozone and aerosol in layer 9 is nonzero, one must exercise some caution in attributing the entire SAGE II-Umkehr difference in this layer to an aerosol effect.

Description: Advanced data assimilation methods are applied to simple but highly nonlinear problems. The dynamical systems studied here are the stochastically forced double well and the Lorenz model. In both systems, linear approximation of the dynamics about the critical points near which regime transitions occur is not always sufficient to track their occurrence or nonoccurrence. Straightforward application of the extended Kalman filter yields mixed results. The ability of the extended Kalman filter to track transitions of the double-well system from one stable critical point to the other depends on the frequency and accuracy of the observations relative to the mean-square amplitude of the stochastic forcing. The ability of the filter to track the chaotic trajectories of the Lorenz model is limited to short times, as is the ability of strong-constraint variational methods. Examples are given to illustrate the difficulties involved, and qualitative explanations for these difficulties are provided. Three generalizations of the extended Kalman filter are described. The first is based on inspection of the innovation sequence, that is, the successive differences between observations and forecasts; it works very well for the double-well problem. The second, an extension to fourth-order moments, yields excellent results for the Lorenz model but will be unwieldy when applied to models with high-dimensional state spaces. A third, more practical method--based on an empirical statistical model derived from a Monte Carlo simulation--is formulated, and shown to work very well. Weak-constraint methods can be made to perform satisfactorily in the context of these simple models, but such methods do not seem to generalize easily to practical models of the atmosphere and ocean. In particular, it is shown that the equations derived in the weak variational formulation are difficult to solve conveniently for large systems.

Description: In this paper, fundamental multiscale circulation modes in the global atmosphere are identified with the objective of providing better understanding of atmospheric low-frequency variabilities over a wide range of spatial and temporal scales. With the use of a combination of rotated principal component technique, singular spectrum analysis, and phase space portraits, three categories of basic multiscale modes in the atmosphere are found. The first is the interannual-mode (IAM), which is dominated by time scales longer than a year and can be attributed to heating and circulation anomalies associated with the coupled tropical ocean-atmosphere, in particular the El Nino-Southern Oscillation. The second is a set of tropical intraseasonal modes consisting of three separate multiscale patterns (ISO-1, -2, -3) related to tropical heating that can be identified with the different phases of the Madden-Julian Oscillation (MJO), including its teleconnection to the extratropics. The ISO spatial and temporal patterns suggest that the extratropical wave train in the North Pacific and North America is related to heating over the Maritime Continent and that the evolution of the MJO around the equator may require forcing from the extratropics spawning convection over the Indian Ocean. The third category represents extratropical intraseasonal oscillations arising from internal dynamics of the basic-state circulation. In the Northern Hemisphere, there are two distinct circulation modes with multiple frequencies in this category: the Pacific/North America (PNA) and the North Atlantic/Eurasia (NAE). In the Southern Hemisphere, two phase-locked modes (PSA-1 and PSA-2) are found depicting an eastward propagating wave train from eastern Australia, via the Pacific South America to the South Atlantic. The extratropical modes exhibit temporal characteristics such as phase locking and harmonic oscillations possibly associated with quadratically nonlinear dynamical systems. Additionally, the observed monthly and seasonal anomalies arise from a complex interplay of the various multiscale low-frequency modes. The relative dominance of the different modes varies widely from month to month and from year to year. On the monthly time scale, while one or two mechanisms may dominate in one year, no single mechanism seems to dominate for all years. There are indications that when the IAM, that is, ENSO heating patterns are strong, the extratropical modes may be suppressed and vice versa. For the seasonal mean, the interannual mode tends to dominate and the contribution from the PNA remains quite significant.

Description: Results are presented of experiments performed with the 3D primitive equation model of the middle atmosphere to elucidate the connections between the behavior of planetary waves in the stratosphere and in the upper troposphere during the seasonal cycle. In these experiments the geopotential height fields at the model's lower boundary consisted of a time-varying zonally symmetric component plus a time-independent wave 1. Experiments were also performed with the same time-varying zonally symmetric component plus a time-independent combination of wave 1 and 2. In each of the wave 1-only experiments, the wave 1 amplitude in the middle stratosphere has two maxima, one in fall and the other in spring, separated by a midwinter minimum. It is suggested that some important transient features in the seasonal evolution of stratospheric planetary waves do not depend on time variations in the amplitude or phase of the corresponding tropospheric waves.

Description: General circulation model (GCM) simulations of low-frequency variability with time scales of 20 to 70 days are analyzed for the Pacific sector during boreal winter. The GCM's leading mode in the upper-tropospheric zonal wind is associated with fluctuations of the East Asian jet; this mode resembles, in both structure and amplitude, the Pacific/North American (PNA) pattern found in the observations on these time scales. In both the model and observations, the PNA anomaly is characterized by: (1) a linear balance in the upper-tropospheric vorticity budget with no significant Rossby wave source in the tropics, (2) a barotropic conversion of kinetic energy from the time mean Pacific jet, and (3) a north/south displacement of the Pacific storm track. In the GCM, the latter is associated with synoptic eddy heat flux and latent heat anomalies that appear to contribute to a strong lower-tropospheric source of wave activity over the North Pacific. This is in contrast to the observations, which show only a weak source of wave activity in this region.

Description: The estimates of rainfall rate derived from a multiparameter radar based on reflectivity factor (R sub ZH), differential reflectivity (R sub DR), and specific differential propagation phase (R sub DP) have widely varying accuracies over the dynamic range of the natural occurrence of rainfall. This paper presents a framework to optimally combine the three estimates, R sub zH, R sub DR, and R sub DP, to derive the best estimate of rainfall using coherent multiparameter radars. The optimization procedure is demonstrated for application to multiparameter radar measurements at C band.

Description: The formal retrieval error analysis of Rodgers (1990) allows the quantitative determination of such retrieval properties as measurement error sensitivity, resolution, and inversion bias. This technique was applied to five numerical inversion techniques and two nonlinear iterative techniques used for the retrieval of middle atmospheric constituent concentrations from limb-scanning millimeter-wave spectroscopic measurements. It is found that the iterative methods have better vertical resolution, but are slightly more sensitive to measurement error than constrained matrix methods. The iterative methods converge to the exact solution, whereas two of the matrix methods under consideration have an explicit constraint, the sensitivity of the solution to the a priori profile. Tradeoffs of these retrieval characteristics are presented.

Description: The accuracy of scatterometer-derived winds using a frequency-domain analysis and simulation is studied. The wavenumber spectra of the Seasat-A-scatterometer (SASS)-derived wind fields have been observed to be accentuated relative to the input wind field. The results of extensive simulations designed to test this observation are reported. Actual SASS measurements of the normalized radar backscatter (NRB) over an orbit (rev) are used as a template to generate simulated NRB measurements. Pointwise estimation of winds from the simulated NRB measurements is accomplished with the ambiguity closest to the true wind selected as the unique wind vector estimate. For comparison, winds are also estimated using a new model-based approach. After wind retrieval, the spectra of the estimated wind fields are computed and compared to the input wind field. The high wavenumber portion of the spectra of the pointwise estimated winds was higher than the spectra of the true winds by an amount which depends on the wind speed variance.

Description: Recent studies on the Amazon deforestation problem predict that removal of the forest will result in a higher surface temperature, a significant reduction in evaporation and precipitation, and possibly significant changes in the tropical circulation. Here, we discuss the basic mechanisms contributing to the response of the tropical atmosphere to deforestation. A simple linear model of the tropical atmosphere is used in studying the effects of deforestation on climate. It is suggested that the impact of large-scale deforestation on the circulation of the tropical atmosphere consists of two components: the response of the tropical circulation to the negative change in precipitation (heating), and the response of the same circulation to the positive change in surface temperature. Owing to their different signs, the changes in predicted temperature and precipitation excite competing responses working in opposite directions. The predicted change in tropical circulation determines the change, if any, in atmospheric moisture convergence, which is equivalent to the change in run-off. The dependence of run-off predictions on the relative magnitudes of the predicted changes in precipitation and surface temperature implies that the predictions about run-off are highly sensitive, which explains, at least partly, the disagreement between the different models concerning the sign of the predicted change in Amazonian run-off.

Description: The vertical structure of monsoon thermal forcing by precipitating convection is diagnosed in terms of horizontal divergence. Airborne Doppler-radar divergence profiles from nine diverse mesoscale convective systems (MCSs) are presented. The MCSs consisted of multicellular convective elements which in time gave rise to areas of stratiform precipitation. Each of the three basic building blocks of the MCSs - convective, intermediary, and stratiform precipitation areas - has a consistent, characteristic divergence profile. Convective areas have low-level convergence, with its peak at 2-4 km altitude, and divergence above 6 km. Intermediary areas have convergence aloft, peaked near 10 km, feeding into mean ascent high in the upper troposphere. Stratiform areas have mid-level convergence, indicating a mesoscale downdraught below the melting level, and a mesoscale updraught aloft. Rawinsonde composite divergence profiles agree with the Doppler data in at least one important respect: the lower-tropospheric convergence into the MCSs peaks 2-4-km above the surface. Rawinsonde vorticity profiles show that monsoonal tropical cyclones spin-up at these elevated levels first, then later descend to the surface. Rawinsonde observations on a larger, continental scale demonstrate that at large horizontal scales only the 'gravest vertical mode' of MCS heating is felt, while the effects of shallower components of the heating (or divergence) profiles are trapped near the heating, as predicted by geostrophic adjustment theory.

Description: Aircraft temperature and pressure measurements as well as satellite imagery are used to establish the amplitudes and the space and time scale of potential temperature disturbances over convective systems. A conceptual model is proposed for the generation of mesoscale gravity waves by convection. The momentum forcing that a reasonable distribution of convection might exert on the tropical stratosphere through convectively excited mesoscale gravity waves of the observed amplitudes is estimated. Aircraft measurements show that presence of mesoscale disturbances in the lower stratospheric temperature, disturbances that appear to be associated with underlying convection. If the disturbances are convectively excited mesoscale gravity waves, their amplitude is sufficient that their breakdown in the upper stratosphere will exert a zonal force comparable to but probably smaller than the planetary-scale Kelvin waves.

Description: A pattern recognition scheme that is used to generate initial seeds or 'first guess' fields for clustering the discrete lightning discharges into storm cells is presented. The clustering process is critically dependent on the prior accuracy of the lightning location estimates. The generation of subsequent storm life-cycle time series (for extrapolation forecasting) also relies on the cluster analysis procedure assigning the correct number of lightning discharges (objects) to the proper storms (groups). The advantages and limitations of different clustering strategies for storm identification and tracking are examined. Storm identification with lightning data alone is compared to storm identification with radar alone, and some synergies for sensor fusion are explored.

Description: Monthly rain falling on the Indian Ocean is mapped for the period 1979 through 1981 by means of observations of the Nimbus-7 Scanning Multichannel Microwave Radiometer. Both stationary and mobile parts were found in the pattern of rain. The stationary part consisted of three zonal and two meridional bands. Only one, the band along and south of the equator, maintained a strong presence through all seasons. A north equatorial counterpart to this south equatorial band also was persistent, but weak. The mobile part of the pattern took the form of a wave. The locus of this wave was an eastward-tilted figure eight, which straddled the equator. The wave moved clockwise along the north loop of the figure eight, counterclockwise along the south loop. The crest of the wave crossed the equator from south to north in May or June and crossed the equator from north to south between August and October. Along its path the equatorial bands were alternately amplified and damped, and the transient bands were activated and suppressed. The effect of the bands and wave was to produce a strong 'monsoon' (annual cycle, summer peak) signature in rain falling over both the northeastern and southwestern reaches of the Indian Ocean.

Description: We examine integrated water vapor fields and rain intensity patterns derived from the Scanning Multichannel Microwave Radiometer (SMMR) and Special Sensor Microwave/Imager (SSM/I) for several rapidly deepening and non-rapidly deepening midlatitude cyclones in the North Atlantic. Our goal is to identify features in the satellite data unique to the rapidly deepening cases, and to explore how these data can potentially be used in the analysis and forecasting of these events.

Description: The use of the Advanced Microwave Sounding Unit (AMSU) moisture imagery for the dynamic, kinematic, and synoptic interpretations requires the availability of idealized 'soundings' of temperature, water vapor, and cloudiness as input to a radiative transfer model for determining how atmospheric conditions control the top-of-atmosphere brightness temperatures. This paper uses idealized vertical configurations of atmospheric water vapor and liquid cloud water to illustrate their effects on the AMSU moisture channels. Calculations with various cloud-drop size distribution functions (SDFs) showed that, at least for the specific SDFs used in this work, the Rayleigh approximation is invalid for mode radii greater than 30-40 microns, for frequencies above 89.0 GHz. It is also invalid above about 20 microns for frequencies above 157 GHz. It is concluded that AMSU moisture channels might not be suitable for detecting convective instability over land using image-based approaches.

Description: Rainfall over tropical ocean regions, particularly in the tropical Pacific, is estimated using Special Sensor Microwave/Imager (SSM/I) data. Instantaneous rainfall estimates are derived from brightness temperature values obtained from the satellite data using the Hughes D-Matrix algorithm. Comparisons with other satellite techniques are made to validate the SSM/I results for the tropical Pacific. The correlation coefficients are relatively high for the three data sets investigated, especially for the annual case.

Description: Full-resolution GOES data archived over North America at hourly time intervals are used to examine the seasonal climatology of deep convective cloudiness for two warm seasons. The analyses discussed here concern convective activity across subtropical latitudes from east of Florida to westward of Baja California. It is demonstrated that high-resolution GOES data can be processed on long-time scales to elucidate aspects of regional convective activity and its monthly, seasonal, and interannual variability. The results can be directly compared with corresponding analyses of mean tropospheric features to identify key forcing mechanisms for regional convection.