ALBERT

Description: Anthropogenic aerosols are intricately linked to the climate system and to the hydrologic cycle. The net effect of aerosols is to cool the climate system by reflecting sunlight. Depending on their composition, aerosols can also absorb sunlight in the atmosphere, further cooling the surface but warming the atmosphere in the process. These effects of aerosols on the temperature profile, along with the role of aerosols as cloud condensation nuclei, impact the hydrologic cycle, through changes in cloud cover, cloud properties and precipitation. Unravelling these feedbacks is particularly difficult because aerosols take a multitude of shapes and forms, ranging from desert dust to urban pollution, and because aerosol concentrations vary strongly over time and space. To accurately study aerosol distribution and composition therefore requires continuous observations from satellites, networks of ground-based instruments and dedicated field experiments. Increases in aerosol concentration and changes in their composition, driven by industrialization and an expanding population, may adversely affect the Earth's climate and water supply.

Description: Pulsed thermoelectrically cooled QC-DFB lasers operating at 15.6 micrometers were characterized for spectroscopic gas sensing applications. A new method for wavelength scanning based on repetition rate modulation was developed. A non-wavelength-selective pyroelectric detector was incorporated in the sensor configuration giving the advantage of room-temperature operation and low cost. Absorption lines of CO2 and H2O were observed in ambient air, providing information about the concentration of these species.

Description: A new method for identifying the structure and other characteristics of extreme weather events is introduced and applied to both model simulations and observations. The approach is based on a linear regression model that links daily extreme precipitation amounts for a particular point on the globe to precipitation and related quantities at all other points. We present here some initial results of our analysis of extreme precipitation events over the United States, including how they are influenced by ENSO and various large-scale teleconnection patterns such as the PNA. The results are based on simulations made with the NASA/NCAR AGCM (Lin and Rood 1996). The quality of the simulated climate for the NASA/NCAR AGCM forced with observed SSTs is described in Chang et al. (2001). The runs analyzed here consist of three 20-year runs forced with idealized cold, neutral and warm ENSO SST anomalies (superimposed on the mean seasonal cycle of SST). The idealized warm or cold SST anomalies are fixed throughout each 20- year simulation and consist of the first EOF (+/- 3 standard deviations) of monthly SST data. Comparisons are made with the results obtained from a similar analysis that uses daily NOAA precipitation observations (Higgins et al. 1996) over the United States and NCEP/NCAR reanalysis data for the period 1949-1998.

Description: The tropics and extratropics are two dynamically distinct regimes. The coupling between these two regimes often defies simple analytical treatment. Progress in understanding of the dynamical interaction between the tropics and extratropics relies on better observational descriptions to guide theoretical development. However, global analyses currently contain significant errors in primary hydrological variables such as precipitation, evaporation, moisture, and clouds, especially in the tropics. Tropical analyses have been shown to be sensitive to parameterized precipitation processes, which are less than perfect, leading to order-one discrepancies between estimates produced by different data assimilation systems. One strategy for improvement is to assimilate rainfall observations to constrain the analysis and reduce uncertainties in variables physically linked to precipitation. At the Data Assimilation Office at the NASA Goddard Space Flight Center, we have been exploring the use of tropical rain rates derived from the TRMM Microwave Imager (TMI) and the Special Sensor Microwave/ Imager (SSM/I) instruments in global data assimilation. Results show that assimilating these data improves not only rainfall and moisture fields but also related climate parameters such as clouds and radiation, as well as the large-scale circulation and short-range forecasts. These studies suggest that assimilation of microwave rainfall observations from space has the potential to significantly improve the quality of 4-D assimilated datasets for climate investigations (Hou et al. 2001). In the next few years, there will be a gradual increase in microwave rain products available from operational and research satellites, culminating to a target constellation of 9 satellites to provide global rain measurements every 3 hours with the proposed Global Precipitation Measurement (GPM) mission in 2007. Continued improvements in assimilation methodology, rainfall error estimates, and model parameterizations are needed to ensure that we derive maximum benefits from these observations.

Description: Data assimilation brings together atmospheric observations and atmospheric models-what we can measure of the atmosphere with how we expect it to behave. NASA's Data Assimilation Office (DAO) sponsors research projects in data reanalysis, which take several years of observational data and analyze them with a fixed assimilation system, to create an improved data set for use in atmospheric studies. Using NCCS computers, one group of NASA researchers employs reanalysis to examine the role of summertime low-level jet (LLJ) winds in regional seasonal climate. Prevailing winds that blow strongly in a fixed direction within a vertically and horizontally confined region of the atmosphere are known as jets. Jets can dominate circulation and have an enormous impact on the weather in a region. Some jets are as famous as they are influential. The jet stream over North America, for instance, is the wind that blows eastward across the continent, bringing weather from the west coast and increasing the speed of airplanes flying to the east coast. The jet stream, while varying in intensity and location, is present in all seasons at the very high altitude of 200-300 millibars - more than 6 miles above Earth's surface.

Description: The intraseasonal variation (ISV) in the 30-60 day band, also known as Madden-Julian oscillation (MJO), has been studied for decades. Madden and Julian showed that the oscillation originated from the western Indian Ocean, propagated eastward, got enhanced over the maritime continent and weakened after passing over the dateline. Composite studies showed evidences of a signal in upper and lower level zonal wind propagating around the globe during an oscillation. Theoretical studies pointed out that the interaction with the warm ocean surface and the coupling with the convective and radiative processes in the atmosphere could manifest the oscillation, which propagates eastward via mutual feedbacks between the wave motions and the cumulus heating. Over tropical South America, no independent 30-60 day oscillation has been reported so far, despite that Amazon is the most distinct tropical convection center over the western hemisphere and the fluxes from its surface of tropical rainforests are close to that from the warm tropical ocean. Liebmann et al. showed a distinct spectral peak of 40-50 day oscillation in outgoing longwave radiation (OLR) over tropical South America and considered that was manifested by the MJO propagation. Nogues-Paegle et al. (2000) focused on a dipole pattern of the OLR anomaly with centers of action over the South Atlantic Convergence Zone (SACZ) and the subtropical plain. They used the regional 10-90 day filtered data and demonstrated this pattern could be represented by the fifth mode of the rotated empirical orthogonal function. Its principal component was further analyzed using the singular spectrum analysis. Their result showed two oscillatory modes with periods of 36-40 days and 22-28 days, of which the former was related to the MJO influence and the latter linked to the remote forcing over southwest of Australia, which produced a wave train propagating southeastward, rounding the southern tip of South America and returning back toward the northeast. The 22-28 day mode has distinct impact on SACZ, responsible for the regional seesaw pattern of alternating dry and wet conditions. In this study we will focus on the 30-60-day spectral band and investigate whether the independent oscillation source over tropical South America is existed. First, we will show the seasonal dependence of the tropical South American ISV in Section 3. Then, the leading principal modes of 30-60 day bandpass filtered 850-hPa velocity potential (VP850) will be computed to distinguish the stationary ISV over tropical South America (SISA) from the propagating MJO in the austral summertime in Section 4. The importance of SISA in representing the regional ISV over South America will be discussed. In Section 5, we will demonstrate the mass oscillation regime of SISA, which is well separated from that of MJO by the Andes, and the convective coupling with rainfall. The dynamical response of SISA and the impact on the South American summer monsoon (SASM) will be presented. Finally, we will give the concluding remarks.

Description: Hydrologists have long speculated that soil moisture information can be used to increase skill in monthly to seasonal forecast systems. For this to be true, though, three conditions must be satisfied: (1) an imposed initial soil moisture anomaly in the forecast system must have some memory, so that it persists into the forecast period; (2) the modeled atmosphere must respond in a predictable way to the persisted anomaly; and (3) the forecast model must correctly represent both the soil moisture memory and the atmospheric response as they occur in nature. In this short paper, we review some recent work at NSIPP (NASA Seasonal-to-Interannual Prediction Project) that addresses all three conditions.

Description: Noting the similarities among the spatial patterns of outgoing longwave radiation among MJO and ENSO, Lau and Chan speculated a possible relationship between the two phenomena. This speculation received a substantial boost in credibility after the 1997-98 El Nino, when MJO activities were found to be substantially enhanced prior to the onset of the warm phase, and clear signals of oceanic Kelvin waves forced by MJO induced anomalous surface wind were detected as possible triggers of ENSO. Yet statistical and modeling studies have so far yielded either nil or at best, very weak relationship between MJO activities and SST. Recently Kessler suggested using an MJO index which includes convective variability in the equatorial central Pacific lead to a more robust MJO-ENSO relationship. Clearly, while MJO might have been instrumental in triggering some El Nino, there are other events that can occur without any MJO trigger.

Description: In this study, we have applied GCM water vapor tracers (WVT) to simulate the North American water cycle. WVTs allow quantitative computation of the geographical source of water for precipitation that occurs anywhere in the model simulation. This can be used to isolate the impact that local surface evaporation has on precipitation, compared to advection and convection. A 15 year 1 deg, 1.25 deg. simulation has been performed with 11 global and 11 North American regional WVTs. Figure 1 shows the source regions of the North American WVTs. When water evaporates from one of these predefined regions, its mass is used as the source for a distinct prognostic variable in the model. This prognostic variable allows the water to be transported and removed (precipitated) from the system in an identical way that occurs to the prognostic specific humidity. Details of the model are outlined by Bosilovich and Schubert (2002) and Bosilovich (2002). Here, we present results pertaining to the onset of the simulated North American monsoon.

Description: The objectives of this study are to (1) develop a better understanding of how observations constrain/impact the MJO in a data assimilation system with the aim of improving the representation of the MJO, and (2) to carry out AGCM predictability/forecast experiments under various observational constraints to assess model errors and sensitivity to initial conditions. Our current focus is on the second objective.

Description: Some 250,000 weather reports are collected by the National Weather Service (NWS) every day. Important measurements are taken by satellites, weather balloons, ground weather stations, airplanes, oceangoing ships, and tethered ocean buoys. Local or global weather models rely on these reports to provide the raw data used as initial conditions for the models to produce a weather prediction.

Description: A portable modular gas sensor for measuring the 13C/12C isotopic ratio in CO2 with a precision of 0.8%(+/-1 sigma) was developed for volcanic gas emission studies. This sensor employed a difference frequency generation (DFG)-based spectroscopic source operating at 4.35 micrometers (approximately 2300 cm-1) in combination with a dual-chamber gas absorption cell. Direct absorption spectroscopy using this specially designed cell permitted rapid comparisons of isotopic ratios of a gas sample and a reference standard for appropriately selected CO2 absorption lines. Special attention was given to minimizing undesirable precision degrading effects, in particular temperature and pressure fluctuations.

Description: Considerable uncertainty surrounds the issue of whether precipitation over the tropical oceans (30 deg N/S) systematically changes with interannual sea-surface temperature (SST) anomalies that accompany El Nino (warm) and La Nina (cold) events. Time series of rainfall estimates from the Tropical Rainfall Measuring Mission (TRMM Precipitation Radar (PR) over the tropical oceans show marked differences with estimates from two TRMM Microwave Imager (TMI) passive microwave algorithms. We show that path-integrated attenuation derived from the effects of precipitation on the radar return from the ocean surface exhibits interannual variability that agrees closely with the TMI time series. Further analysis of the frequency distribution of PR (2A25 product) rain rates suggests that the algorithm incorporates the attenuation measurement in a very conservative fashion so as to optimize the instantaneous rain rates. Such an optimization appears to come at the expense of monitoring interannual climate variability.

Description: A workshop on cumulus parameterization took place at the NASA Goddard Space Flight Center from December 3-5, 2001. The major objectives of this workshop were (1) to review the problem of representation of moist processes in large-scale models (mesoscale models, Numerical Weather Prediction models and Atmospheric General Circulation Models), (2) to review the state-of-the-art in cumulus parameterization schemes, and (3) to discuss the need for future research and applications. There were a total of 31 presentations and about 100 participants from the United States, Japan, the United Kingdom, France and South Korea. The specific presentations and discussions during the workshop are summarized in this paper.

Description: The results presented here show that tropical convection plays a role in each of the three primary processes involved in the in situ formation of tropopause cirrus. First, tropical convection transports moisture from the surface into the upper troposphere. Second, tropical convection excites Rossby waves that transport zonal momentum toward the ITCZ, thereby generating rising motion near the equator. This rising motion helps transport moisture from where it is detrained from convection to the cold-point tropopause. Finally, tropical convection excites vertically propagating tropical waves (e.g. Kelvin waves) that provide one source of large-scale cooling near the cold-point tropopause, leading to tropopause cirrus formation.

Description: We show the comparisons between ground-based measurements of spectrally integrated (300 nm to 380 nm) ultraviolet (UV) irradiance with satellite estimates from the Total Ozone Mapping Spectrometer (TOMS) total ozone and reflectivity data for the whole period of TOMS measurements (1979-2000) over the Meteorological Observatory of Moscow State University (MO MSU), Moscow, Russia. Several aspects of the comparisons are analyzed, including effects of cloudiness, aerosol, and snow cover. Special emphasis is given to the effect of different spatial and temporal averaging of ground-based data when comparing with low-resolution satellite measurements (TOMS footprint area 50-200 sq km). The comparisons in cloudless scenes with different aerosol loading have revealed TOMS irradiance overestimates from +5% to +20%. A-posteriori correction of the TOMS data accounting for boundary layer aerosol absorption (single scattering albedo of 0.92) eliminates the bias for cloud-free conditions. The single scattering albedo was independently verified using CIMEL sun and sky-radiance measurements at MO MSU in September 2001. The mean relative difference between TOMS UV estimates and ground UV measurements mainly lies within 1 10% for both snow-free and snow period with a tendency to TOMS overestimation in snow-free period especially at overcast conditions when the positive bias reaches 15-17%. The analysis of interannual UV variations shows quite similar behavior for both TOMS and ground measurements (correlation coefficient r=0.8). No long-term trend in the annual mean bias was found for both clear-sky and all-sky conditions with snow and without snow. Both TOMS and ground data show positive trend in UV irradiance between 1979 and 2000. The UV trend is attributed to decreases in both cloudiness and aerosol optical thickness during the late 1990's over Moscow region. However, if the analyzed period is extended to include pre-TOMS era (1968-2000 period), no trend in ground UV irradiance is detected.

Description: In a companion paper, the temperature dependence of Raman scattering and its influence on the Raman water vapor signal and the lidar equations was examined. New forms of the lidar equation were developed to account for this temperature sensitivity. Here we use those results to derive the temperature dependent forms of the equations for the aerosol scattering ratio, aerosol backscatter coefficient, extinction to backscatter ratio and water vapor mixing ratio. Pertinent analysis examples are presented to illustrate each calculation.

Description: Based on the single-scattering optical properties that are pre-computed using an improve geometric optics method, the bulk mass absorption coefficient, single-scattering albedo, and asymmetry factor of ice particles have been parameterized as a function of the mean effective particle size of a mixture of ice habits. The parameterization has been applied to compute fluxes for sample clouds with various particle size distributions and assumed mixtures of particle habits. Compared to the parameterization for a single habit of hexagonal column, the solar heating of clouds computed with the parameterization for a mixture of habits is smaller due to a smaller cosingle-scattering albedo. Whereas the net downward fluxes at the TOA and surface are larger due to a larger asymmetry factor. The maximum difference in the cloud heating rate is approx. 0.2 C per day, which occurs in clouds with an optical thickness greater than 3 and the solar zenith angle less than 45 degrees. Flux difference is less than 10 W per square meters for the optical thickness ranging from 0.6 to 10 and the entire range of the solar zenith angle. The maximum flux difference is approximately 3%, which occurs around an optical thickness of 1 and at high solar zenith angles.

Description: A case of torrential precipitation process in the Mei-yu front, an Asian monsoon system east to the Tibetan Plateau, is studied with the coupled Penn State University/NCAR MM5 and NASA/GSFC PLACE (Parameterization for Land - Atmosphere - Cloud Exchange) models. Remote and local impacts of water vapor on the location and intensity of Mei-yu precipitation are studied by numerical experiments. Results demonstrate that the water vapor source for this heavy precipitation case in Yangtze river basin is derived mostly from the Bay of Bengal, transported by the southwesterly low-level Jet (LLJ) southeast to the Tibetan Plateau. The moist convection is a critical process in the development and maintenance of the front. The meridional and zonal secondary circulations resulted from Mei-yu condensation heating both act to increase the wind speed in the LLJ. The condensation induced local circulation strengthens the moisture transport in the LLJ, providing a positive feedback to sustain the Mei-yu precipitation system. It is found that local precipitation recycling shifts heavy rain toward the warm side of the Mei-yu front. This shift of rainfall location is due to the pronounced increase of atmospheric moisture and decrease of surface temperature over the warm side of the front.

Description: The aim of this paper is to report extreme winter/early-spring air temperature (hereinafter temperature) anomalies in mid-latitude Europe, and to discuss the underlying forcing to these interannual fluctuations. Warm advection from the North Atlantic in late winter controls the surface-air temperature, as indicated by the substantial correlation between the speed of the surface southwesterlies over the eastern North Atlantic (quantified by a specific Index Ina) and the 2-meter level air temperatures (hereinafter Ts) over Europe, 45-60 deg N, in winter. In mid-March and subsequently, the correlation drops drastically (quite often it is negative). This change in the relationship between Ts and Ina marks a transition in the control of the surface-air temperature: absorption of insolation replaces the warm advection as the dominant control. This forcing by maritime-air advection in winter was demonstrated in a previous publication, and is re-examined here in conjunction with extreme fluctuations of temperatures in Europe. We analyze here the interannual variability at its extreme by comparing warm-winter/early-spring of 1989/90 with the opposite scenario in 1995/96. For these two December-to-March periods the differences in the monthly mean temperature in Warsaw and Torun, Poland, range above 10 C. Short-term (shorter than a month) fluctuations of the temperature are likewise very strong. We conduct pentad-by-pentad analysis of the surface-maximum air temperature (hereinafter Tmax), in a selected location, examining the dependence on Ina. The increased cloudiness and higher amounts of total precipitable water, corollary effects to the warm low-level advection. in the 1989/90 winter, enhance the positive temperature anomalies. The analysis of the ocean surface winds is based on the Special Sensor Microwave/Imager (SSM/I) dataset; ascent rates, and over land wind data are from the European Centre for Medium-Range Weather Forecasts (ECMWF); maps of 2-m temperature, cloud cover and precipitable water are from the National Centers for Environmental Prediction (NCEP) Reanalysis.

Description: The potential role of soil moisture initialization in seasonal forecasting is illustrated through ensembles of simulations with the NASA Seasonal-to-Interannual Prediction Project (NSIPP) model. For each boreal summer during 1997-2001, we generated two 16-member ensembles of 3-month simulations. The first, "AMIP-style" ensemble establishes the degree to which a perfect prediction of SSTs would contribute to the seasonal prediction of precipitation and temperature over continents. The second ensemble is identical to the first, except that the land surface is also initialized with "realistic" soil moisture contents through the continuous prior application (within GCM simulations leading up to the start of the forecast period) of a daily observational precipitation data set and the associated avoidance of model drift through the scaling of all surface prognostic variables. A comparison of the two ensembles shows that soil moisture initialization has a statistically significant impact on summertime precipitation and temperature over only a handful of continental regions. These regions agree, to first order, with regions that satisfy three conditions: (1) a tendency toward large initial soil moisture anomalies, (2) a strong sensitivity of evaporation to soil moisture, and (3) a strong sensitivity of precipitation to evaporation. The degree to which the initialization improves forecasts relative to observations is mixed, reflecting a critical need for the continued development of model parameterizations and data analysis strategies.

Description: Data from three cloudy days (March 3, 21, 29, 2000) of the ARM Enhanced Shortwave Experiment II (ARESE II) were analyzed. Grand averages of broadband absorptance among three sets of instruments were compared. Fractional solar absorptances were approx. 0.21-0.22 with the exception of March 3 when two sets of instruments gave values smaller by approx. 0.03-0.04. The robustness of these values was investigated by looking into possible sampling problems with the aid of 500 nm spectral fluxes. Grand averages of 500 nm apparent absorptance cover a wide range of values for these three days, namely from a large positive (approx. 0.011) average for March 3, to a small negative (approximately -0.03) for March 21, to near zero (approx. 0.01) for March 29. We present evidence suggesting that a large part of the discrepancies among the three days is due to the different nature of clouds and their non-uniform sampling. Hence, corrections to the grand average broadband absorptance values may be necessary. However, application of the known correction techniques may be precarious due to the sparsity of collocated flux measurements above and below the clouds. Our analysis leads to the conclusion that only March 29 fulfills all requirements for reliable estimates of cloud absorption, that is, the presence of thick, overcast, homogeneous clouds.

Description: This study uses a twenty-three year (1979-2001) satellite-gauge merged community data set to further describe the relationship between El Nino Southern Oscillation (ENSO) and precipitation. The globally complete precipitation fields reveal coherent bands of anomalies that extend from the tropics to the polar regions. Also, ENSO-precipitation relationships were analyzed during the six strongest El Ninos from 1979 to 2001. Seasons of evolution, Pre-onset, Onset, Peak, Decay, and Post-decay, were identified based on the strength of the El Nino. Then two simple and independent models, first order harmonic and linear, were fit to the monthly time series of normalized precipitation anomalies for each grid block. The sinusoidal model represents a three-phase evolution of precipitation, either dry-wet-dry or wet-dry-wet. This model is also highly correlated with the evolution of sea surface temperatures in the equatorial Pacific. The linear model represents a two-phase evolution of precipitation, either dry-wet or wet-dry. These models combine to account for over 50% of the precipitation variability for over half the globe during El Nino. Most regions, especially away from the Equator, favor the linear model. Areas that show the largest trend from dry to wet are southeastern Australia, eastern Indian Ocean, southern Japan, and off the coast of Peru. The northern tropical Pacific and Southeast Asia show the opposite trend.

Description: The fixed-lag Kalman smoother (FLKS) has been proposed as a framework to construct data assimilation procedures capable of producing high-quality climate research datasets. Fixed-lag Kalman smoother-based systems, referred to as retrospective data assimilation systems, are an extension to three-dimensional filtering procedures with the added capability of incorporating observations not only in the past and present time of the estimate, but also at future times. A variety of simplifications are necessary to render retrospective assimilation procedures practical. In this article, we present an FLKS-based retrospective data assimilation system implementation for the Goddard Earth Observing System (GOES) Data Assimilation System (DAS). The practicality of this implementation comes from the practicality of its underlying (filter) analysis system, i.e., the physical-space statistical analysis system (PSAS). The behavior of two schemes is studied here. The first retrospective analysis (RA) scheme is designed simply to update the regular PSAS analyses with observations available at times ahead of the regular analysis times. Although our GEOS DAS implementation is general, results are only presented for when observations 6-hours ahead of the analysis time are used to update the PSAS analyses and thereby to calculate the so-called lag-1 retrospective analyses. Consistency tests for this RA scheme show that the lag-1 retrospective analyses indeed have better 6-hour predictive skills than the predictions from the regular analyses. This motivates the introduction of the second retrospective analysis scheme which, at each analysis time, uses the 6-hour retrospective analysis to replace the first-guess normally used in the PSAS analysis, and therefore allows the calculation of a revised (filter) PSAS analysis. Since in this scheme the lag-1 retrospective analyses influence the filter results, this procedure is referred to as the retrospective-based iterative analysis (RIA) scheme. Results from the RIA scheme indicate its potential for improving the overall quality of the assimilation.

Description: The implications of using different control variables for the analysis of moisture observations in a global atmospheric data assimilation system are investigated. A moisture analysis based on either mixing ratio or specific humidity is prone to large extrapolation errors, due to the high variability in space and time of these parameters and to the difficulties in modeling their error covariances. Using the logarithm of specific humidity does not alleviate these problems, and has the further disadvantage that very dry background estimates cannot be effectively corrected by observations. Relative humidity is a better choice from a statistical point of view, because this field is spatially and temporally more coherent and error statistics are therefore easier to obtain. If, however, the analysis is designed to preserve relative humidity in the absence of moisture observations, then the analyzed specific humidity field depends entirely on analyzed temperature changes. If the model has a cool bias in the stratosphere this will lead to an unstable accumulation of excess moisture there. A pseudo-relative humidity can be defined by scaling the mixing ratio by the background saturation mixing ratio. A univariate pseudo-relative humidity analysis will preserve the specific humidity field in the absence of moisture observations. A pseudorelative humidity analysis is shown to be equivalent to a mixing ratio analysis with flow-dependent covariances. In the presence of multivariate (temperature-moisture) observations it produces analyzed relative humidity values that are nearly identical to those produced by a relative humidity analysis. Based on a time series analysis of radiosonde observed-minus-background differences it appears to be more justifiable to neglect specific humidity-temperature correlations (in a univariate pseudo-relative humidity analysis) than to neglect relative humidity-temperature correlations (in a univariate relative humidity analysis). A pseudo-relative humidity analysis is easily implemented in an existing moisture analysis system, by simply scaling observed-minus background moisture residuals prior to solving the analysis equation, and rescaling the analyzed increments afterward.

Description: Surface-air temperatures measured in winter at 3 meteorological stations in central Europe rise substantially for most of the second-half of the 20th century. This means shorter winter, and longer growing season, which has positive implications for regional agriculture. However, these positive trends stopped in winter of 1996, and for the recent 7 years no further climatic amelioration is reported.

Description: In this study, a technique for estimating vertical profiles of precipitation from multifrequency, multiresolution active and passive microwave observations is investigated using both simulated and airborne data. The technique is applicable to the Tropical Rainfall Measuring Mission (TRMM) satellite multi-frequency active and passive observations. These observations are characterized by various spatial and sampling resolutions. This makes the retrieval problem mathematically more difficult and ill-determined because the quality of information decreases with decreasing resolution. A model that, given reflectivity profiles and a small set of parameters (including the cloud water content, the intercept drop size distribution, and a variable describing the frozen hydrometeor properties), simulates high-resolution brightness temperatures is used. The high-resolution simulated brightness temperatures are convolved at the real sensor resolution. An optimal estimation procedure is used to minimize the differences between simulated and observed brightness temperatures. The retrieval technique is investigated using cloud model synthetic and airborne data from the Fourth Convection And Moisture Experiment. Simulated high-resolution brightness temperatures and reflectivities and airborne observation strong are convolved at the resolution of the TRMM instruments and retrievals are performed and analyzed relative to the reference data used in observations synthesis. An illustration of the possible use of the technique in satellite rainfall estimation is presented through an application to TRMM data. The study suggests improvements in combined active and passive retrievals even when the instruments resolutions are significantly different. Future work needs to better quantify the retrievals performance, especially in connection with satellite applications, and the uncertainty of the models used in retrieval.

Description: The GEWEX Cloud System Study (GCSS; GEWEX is the Global Energy and Water Cycle Experiment) was organized to promote development of improved parameterizations of cloud systems for use in climate and numerical weather prediction models, with an emphasis on the climate applications. The strategy of GCSS is to use two distinct kinds of models to analyze and understand observations of the behavior of several different types of clouds systems. Cloud-system-resolving models (CSRMs) have high enough spatial and temporal resolutions to represent individual cloud elements, but cover a wide enough range of space and time scales to permit statistical analysis of simulated cloud systems. Results from CSRMs are compared with detailed observations, representing specific cases based on field experiments, and also with statistical composites obtained from satellite and meteorological analyses. Single-column models (SCMs) are the surgically extracted column physics of atmospheric general circulation models. SCMs are used to test cloud parameterizations in an un-coupled mode, by comparison with field data and statistical composites. In the original GCSS strategy, data is collected in various field programs and provided to the CSRM Community, which uses the data to "certify" the CSRMs as reliable tools for the simulation of particular cloud regimes, and then uses the CSRMs to develop parameterizations, which are provided to the GCM Community. We report here the results of a re-thinking of the scientific strategy of GCSS, which takes into account the practical issues that arise in confronting models with data. The main elements of the proposed new strategy are a more active role for the large-scale modeling community, and an explicit recognition of the importance of data integration.

Description: New state of the art methodology is described to analyze AIRS/AMSU/HSB data in the presence of multiple cloud formations. The methodology forms the basis for the AIRS Science Team algorithm which will be used to analyze AIRS/AMSU/HSB data on EOS Aqua. The cloud clearing methodology requires no knowledge of the spectral properties of the clouds. The basic retrieval methodology is general and extracts the maximum information from the radiances, consistent with the channel noise covariance matrix. The retrieval methodology minimizes the dependence of the solution on the first guess field and the first guess error characteristics. Results are shown for AIRS Science Team simulation studies with multiple cloud formations. These simulation studies imply that clear column radiances can be reconstructed under partial cloud cover with an accuracy comparable to single spot channel noise in the temperature and water vapor sounding regions, temperature soundings can be produced under partial cloud cover with RMS errors on the order of, or better than, 1deg K in 1 km thick layers from the surface to 700 mb, 1 km layers from 700 mb to 300 mb, 3 km layers from 300 mb to 30 mb, and 5 km layers from 30 mb to 1 mb, and moisture profiles can be obtained with an accuracy better than 20% absolute errors in 1 km layers from the surface to nearly 200 mb.

Description: Over tropical land regions, rain rate maxima in mesoscale convective systems revealed by the Precipitation Radar (PR) flown on the Tropical Rainfall Measuring Mission (TRMM) satellite are found to correspond to thunderstorms, i.e., Cbs. These Cbs are reflected as minima in the 85 GHz brightness temperature, T85, observed by the TRMM Microwave Imager (TMI) radiometer. Because the magnitude of TMI observations do not discriminate satisfactorily convective and stratiform rain, we developed here a different TMI discrimination method. In this method, two types of Cbs, strong and weak, are inferred from the Laplacian of T85 at minima. Then, to retrieve rain rate, where T85 is less than 270 K, a weak (background) rain rate is deduced using T85 observations. Furthermore, over a circular area of 10 km radius centered at the location of each T85 minimum, an additional Cb component of rain rate is added to the background rain rate. This Cb component of rain rate is estimated with the help of (T19-T37) and T85 observations. Initially, our algorithm is calibrated with the PR rain rate measurements from 20 MCS rain events. After calibration, this method is applied to TMI data taken from several tropical land regions. With the help of the PR observations, we show that the spatial distribution and intensity of rain rate over land estimated from our algorithm are better than those given by the current TMI-Version-5 Algorithm. For this reason, our algorithm may be used to improve the current state of rain retrievals on land.

Description: The Project for Intercomparison of Land Surface Parameterization Schemes (PILPS) Phase 2(e) showed that in cold regions the annual runoff production in Land Surface Schemes (LSSs) is closely related to the maximum snow accumulation, which in turn is controlled in large part by winter sublimation. To help further explain the relationship between snow cover, turbulent exchanges and runoff production, a simple equivalent model-(SEM) was devised to reproduce the seasonal and annual fluxes simulated by 13 LSSs that participated in PILPS Phase 2(e). The design of the SEM relates the annual partitioning of precipitation and energy in the LSSs to three primary parameters: snow albedo, effective aerodynamic resistance and evaporation efficiency. Isolation of each of the parameters showed that the annual runoff production was most sensitive to the aerodynamic resistance. The SEM was somewhat successful in reproducing the observed LSS response to a decrease in shortwave radiation and changes in wind speed forcings. SEM parameters derived from the reduced shortwave forcings suggested that increased winter stability suppressed turbulent heat fluxes over snow. Because winter sensible heat fluxes were largely negative, reductions in winter shortwave radiation imply an increase in annual average sensible heat.

Description: Urban heat islands (UHIs) are caused by the heat-retaining properties of surfaces usually found in urban cities like asphalt and concrete. The UHI can typically be observed on the evening TV weather map as warmer temperatures over the downtown of major cities and cooler temperatures in the suburbs and surrounding rural areas. The UHI has now become a widely acknowledged, observed, and researched phenomenon because of its broad environmental and societal implications. Interest in the UHI will intensify in the future as existing urban areas expand and rural areas urbanize. By the year 2025, more than 60% of the world's population will live in cities, with higher percentages expected in developed nations. The urban growth rate in the United States, for example, is estimated to be 12.5%, and the recent 2000 Census found that more than 80% of the population currently lives in urban areas. Furthermore, the U.S. population is not only growing but is tending to concentrate more in urban areas within the environmentally sensitive coastal zones. Urban growth creates unique and often contentious issues for policymakers related to land use zoning, transportation planning, agricultural production, housing and development, pollution, and natural resources protection. Urban expansion and its associated UHIs also have measurable impacts on weather and climate processes. The UHI has been documented to affect local and regional temperature, wind patterns, and air quality.

Description: The simulations of climatology and response of the South American summer monsoon (SASM) to the 1997/98 El Nino are investigated using six atmospheric general circulation models. Results show all models simulate the large-scale features of the SASM reasonably well. However, both stationary and seasonal components of the surface pressure are overestimated, resulting in an excessively strong SASM in the model climatology. The low-level northwesterly jet over eastern foothills of the Andes is not well resolved because of the coarse resolution of the models. Large rainfall simulation biases are found in association with the Andes and the Atlantic ITCZ, indicating model problems in handling steep mountains and parameterization of convective processes. The simulation of the 1997/98 El Nino impact on SASM is examined based on an ensemble of ten two-year (September 1996 - August 1998) integration. Results show that most models can simulate the large-scale tropospheric warming response over the tropical central Pacific, including the dynamic response of Rossby wave propagation of the Pacific-South America (PSA) pattern that influences remote areas. Deficiencies are found in simulating the regional impacts over South America. Model simulation fails to capture the southeastward expansion of anomalously warm tropospheric air. As a result, the upper tropospheric anomalous high over the subtropical Andes is less pronounced, and the enhancement of subtropical westerly jet is displaced 5deg-10deg equatorward compared to the observed. Over the Amazon basin, the shift of Walker cell induced by El Nino is not well represented, showing anomalous easterlies in both upper and lower troposphere.

Description: When from the southwest, North Atlantic ocean surface winds are known to bring warm and moist airmasses into central Europe in winter. By tracing backward trajectories from western Europe, we establish that these airmasses originate in the southwestern North Atlantic, in the very warm regions of the Gulf Stream. Over the eastern North Atlantic, Lt the gateway to Europe, the ocean-surface winds changed directions in the second half of the XXth century, those from the northwest and from the southeast becoming so infrequent, that the direction from the southwest became even more dominant. For the January-to-March period, the strength of south-westerlies in this region, as well as in the source region, shows in the years 1948-1995 a significant increase, above 0.2 m/sec/ decade. Based on the sensitivity of the surface temperature in Europe, slightly more than 1 C for a 1m/sec increase in the southwesterly wind, found in the previous studies, the trend in the warm advection accounts for a large part of the warming in Europe established for this period in several reports. However, for the most recent years, 1996-2001, the positive trend in the southwesterly advection appears to be is broken, which is consistent with unseasonally cold events reported in Europe in those winters. This study had, some bearing on evaluating the respective roles of the North Atlantic Oscillation and the Greenhouse Gas Global warming, GGG, in the strong winter warming observed for about half a century over the northern-latitude continents. Changes in the ocean-surface temperatures induced by GGG may have produced the dominant southwesterly direction of the North Atlantic winds. However, this implies a monotonically (apart from inherent interannual variability) increasing advection, and if the break in the trend which we observe after 1995 persists, this mechanism is counter-indicated. The 1948-1995 trend in the south-westerlies could then be considered to a large degree attributable to the North Atlantic Oscillation.

Description: The 17-year (1982-1998) trend in surface temperature shows a general cooling over the Antarctic continent, warming of the sea ice zone, with moderate changes over the oceans. Warming of the peripheral seas is associated with negative trends in the regional sea ice extent. Effects of the Southern Hemisphere Annular Mode (SAM) and the extrapolar Southern Oscillation (SO) on surface temperature are quantified through regression analysis. Positive polarities of the SAM are associated with cold anomalies over most of Antarctica, with the most notable exception of the Antarctic Peninsula. Positive temperature anomalies and ice edge retreat in the Pacific sector are associated with El Nino episodes. Over the past two decades, the drift towards high polarity in the SAM and negative polarity in the SO indices couple to produce a spatial pattern with warmer temperatures in the Antarctic Peninsula and peripheral seas, and cooler temperatures over much of East Antarctica.

Description: Validation of satellite remote-sensing methods for estimating rainfall against rain-gauge data is attractive because of the direct nature of the rain-gauge measurements. Comparisons of satellite estimates to rain-gauge data are difficult, however, because of the extreme variability of rain and the fact that satellites view large areas over a short time while rain gauges monitor small areas continuously. In this paper, a statistical model of rainfall variability developed for studies of sampling error in averages of satellite data is used to examine the impact of spatial and temporal averaging of satellite and gauge data on intercomparison results. The model parameters were derived from radar observations of rain, but the model appears to capture many of the characteristics of rain-gauge data as well. The model predicts that many months of data from areas containing a few gauges are required to validate satellite estimates over the areas, and that the areas should be of the order of several hundred km in diameter. Over gauge arrays of sufficiently high density, the optimal areas and averaging times are reduced. The possibility of using time-weighted averages of gauge data is explored.

Description: A NOAA P-3 instrumented aircraft observed an intense, fast-moving narrow cold frontal Farmhand as it approached the Pacific Northwest coast on 19 February 2001 during the Pacific Coastal Jets Experiment. Pseudo-dual-Doppler analyses performed on the airborne Doppler radar data while the frontal system was well offshore indicated that a narrow ribbon of very high radar reflectively convective cores characterized the Farmhand at low levels with echo tops to approximately 4-5 km. The NCFR exhibited gaps in its narrow ribbon of high reflectively, probably as a result of hydrodynamic instability all no its advancing cold pool leading edge. In contrast to some earlier studies of cold frontal rainbands, density current theory described well the motion of the overall front. The character of the updraft structure associated with the heavy rainfall at its leading edge varied across the gap region. The vertical shear of the cross-frontal low-level ambient flow exerted a strong influence on the updraft character, consistent with theoretical arguments developed for squall lines describing the balance of vorticity at the leading edge. In short regions south of the gaps the vertical wind shear was strongest with the updrafts and rain shafts more intense, narrower, and more erect or even downshear tilted. North of the gaps the wind shear weakened with less intense Dihedrals which tilted upshear with a broader band of rainfall. Simulations using a nonhydrostatic mesoscale nested grid model are used to investigate the gap regions, particularly the balance of cold pool induced to pre-frontal ambient shears at the leading edge. Observations confirm the model results that the updraft character depends on the balance of vorticity at the leading edge. Downshear-tilted updrafts imply that convection south of the gap regions would weaken with time relative to the frontal segments north of the gaps since inflow air would be affected by passage through the heavy rain region before ascent, suggesting a mechanism for gap filling.

Description: The sensitivities to surface friction and the Coriolis parameter in tropical cyclogenesis are studied using an axisymmetric version of the Goddard cloud ensemble model. Our experiments demonstrate that tropical cyclogenesis can still occur without surface friction. However, the resulting tropical cyclone has very unrealistic structure. Surface friction plays an important role of giving the tropical cyclones their observed smaller size and diminished intensity. Sensitivity of the cyclogenesis process to surface friction. in terms of kinetic energy growth, has different signs in different phases of the tropical cyclone. Contrary to the notion of Ekman pumping efficiency, which implies a preference for the highest Coriolis parameter in the growth rate if all other parameters are unchanged, our experiments show no such preference.

Description: The summer climate of southern Mexico and Central America is characterized by a mid summer drought (MSD), where rainfall is reduced by 40% in July as compared to June and September. A mid-summer reduction in the climatological number of eastern Pacific tropical cyclones has also been noted. Little is understood about the climatology and interannual variability of these minima. The present study uses a novel approach to quantify the bimodal distribution of summertime rainfall for the globe and finds that this feature of the annual cycle is most extreme over Pan America and adjacent oceans. One dominant interannual signal in this region occurs the summer before a strong winter El Nino/Southern Oscillation ENSO. Before El Nino events the region is dry, the MSD is strong and centered over the ocean, and the mid-summer minimum in tropical cyclone frequency is most pronounced. This is significantly different from Neutral cases (non-El Nino and non-La Nina) when the MSD is weak and positioned over the land bridge. The MSD is highly variable for La Nina years, and there is not an obvious mid-summer minimum in the number of tropical cyclones.

Description: The autoconversion rate is a key process for the formation of precipitation in warm clouds. In climate models, physical processes such as autoconversion rate, which are calculated from grid mean values, are biased, because they do not take subgrid variability into account. Recently, statistical cloud schemes have been introduced in large-scale models to account for partially cloud-covered grid boxes. However, these schemes do not include the in-cloud variability in their parameterizations. In this paper, a new statistically based autoconversion rate considering the in-cloud variability is introduced and tested in three cases using the Canadian Single Column Model (SCM) of the global climate model. The results show that the new autoconversion rate improves the model simulation, especially in terms of liquid water path in all three case studies.

Description: A recently developed technique called cloud slicing used for deriving upper tropospheric ozone from the Nimbus 7 Total Ozone Mapping Spectrometer (TOMS) instrument combined together with temperature-humidity and infrared radiometer (THIR) is no longer applicable to the Earth Probe TOMS (EPTOMS) because EPTOMS does not have an instrument to measure cloud top temperatures. For continuing monitoring of tropospheric ozone between 200-500hPa and testing the feasibility of this technique across spacecrafts, EPTOMS data are co-located in time and space with the Geostationary Operational Environmental Satellite (GOES)-8 infrared data for 2001 and early 2002, covering most of North and South America (45S-45N and 120W-30W). The maximum column amounts for the mid-latitudinal sites of the northern hemisphere are found in the March-May season. For the mid-latitudinal sites of the southern hemisphere, the highest column amounts are found in the September-November season, although overall seasonal variability is smaller than those of the northern hemisphere. The tropical sites show the weakest seasonal variability compared to higher latitudes. The derived results for selected sites are cross validated qualitatively with the seasonality of ozonesonde observations and the results from THIR analyses over the 1979-1984 time period due to the lack of available ozonesonde measurements to study sites for 2001. These comparisons show a reasonably good agreement among THIR, ozonesonde observations, and cloud slicing-derived column ozone. With very limited co-located EPTOMS/GOES data sets, the cloud slicing technique is still viable to derive the upper tropospheric column ozone. Two new variant approaches, High-Low (HL) cloud slicing and ozone profile derivation from cloud slicing are introduced to estimate column ozone amounts using the entire cloud information in the troposphere.

Description: The Global Precipitation Measurement (GPM) program is a mission to measure precipitation from space, and is a similar but much expanded mission of the Tropical Rainfall Measuring Mission. Its scope is not limited to scientific research, but includes practical and operational applications such as weather forecasting and water resource management. To meet the requirements of operational use, the GPM uses multiple low-orbiting satellites to increase the sampling frequency and to create three-hourly global rain maps that will be delivered to the world in quasi-real time. A dual-frequency radar (DPR) will be installed on the primary satellite that plays an important role in the whole mission. The DPR will realize measurement of precipitation with high sensitivity, high precision and high resolutions. This paper describes an outline of the GPM program, its issues and the roles and development of the DPR.

Description: ENSCO, Inc., developed the Meteorological and Atmospheric Real-time Safety Support (MARSS) system for real-time assessment of meteorological data displays and toxic material spills. MARSS also provides mock scenarios to guide preparations for emergencies involving meteorological hazards and toxic substances. Developed under a Small Business Innovation Research (SBIR) contract with Kennedy Space Center, MARSS was designed to measure how safe NASA and Air Force range safety personnel are while performing weather sensitive operations around launch pads. The system augments a ground operations safety plan that limits certain work operations to very specific weather conditions. It also provides toxic hazard prediction models to assist safety managers in planning for and reacting to releases of hazardous materials. MARSS can be used in agricultural, industrial, and scientific applications that require weather forecasts and predictions of toxic smoke movement. MARSS is also designed to protect urban areas, seaports, rail facilities, and airports from airborne releases of hazardous chemical substances. The system can integrate with local facility protection units and provide instant threat detection and assessment data that is reportable for local and national distribution.

Description: Through a cooperative venture with NASA's Stennis Space Center, WorldWinds, Inc., developed a unique weather and wave vector map using space-based radar satellite information and traditional weather observations. Called WorldWinds, the product provides accurate, near real-time, high-resolution weather forecasts. It was developed for commercial and scientific users. In addition to weather forecasting, the product's applications include maritime and terrestrial transportation, aviation operations, precision farming, offshore oil and gas operations, and coastal hazard response support. Target commercial markets include the operational maritime and aviation communities, oil and gas providers, and recreational yachting interests. Science applications include global long-term prediction and climate change, land-cover and land-use change, and natural hazard issues. Commercial airlines have expressed interest in the product, as it can provide forecasts over remote areas. WorldWinds, Inc., is currently providing its product to commercial weather outlets.

Description: Pollution is often considered a localized phenomenon, but it is now clear that it travels from region-to-region, country to country, and even continent to continent. In addition to urban pollution in developed countries, large emissions from developing nations and large-scale biomass fires add to the global pollution burden. Ozone and aerosols are two components of pollution that contribute to radiative forcing of the earth s climate. In turn, as climate changes, rates of chemical and microphysical reactions may be perturbed. Considering the earth as a coupled chemical-microphysical-climate system poses challenges for models and observations alike. These issues were the topic of a Workshop held in May 2002 at NASA GSFC s Laboratory for Atmospheres. Highlights of the Workshop are summarized in this article.

Description: Hundreds of Vaisala sondes with a RS80-H Humicap thin-film capacitor humidity sensor were launched during the Tropical Rainfall Measuring Mission (TRMM) field campaigns in Large Scale Biosphere-Atmosphere held in Brazil (LBA) and in Kwajalein experiment (KWAJEX) held in the Republic of Marshall Islands. Using Six humidity error correction algorithms by Wang et al., these sondes were corrected for significant dry bias in the RS80-H data. It is further shown that sonde surface temperature error must be corrected for a better representation of the relative humidity. This error becomes prominent due to sensor arm-heating in the first 50-s data.

Description: The primary cause of the midwestern North American drought in the summer of 1988 has been identified to be the La Nina SST anomalies. Yet with the SST anomalies prescribed, this drought has not been simulated satisfactorily by any general circulation model. Seven simulation-experiments, each containing an ensemble of 4-sets of simulations, were conducted with the GEOS GCM for both 1987 and 1988. All simulations started from January 1 and continued through the end of August. In the first baseline case, Case 1, only the SST anomalies and some vegetation parameters were prescribed, while everything else (such as soil moisture, snow-cover, and clouds) was interactive. The GCM did produce some of the circulation features of a drought over North America, but they could only be identified on the planetary scales. The 1988 minus 1987 precipitation differences show that the GCM was successful in simulating reduced precipitation in the mid-west, but the accompanying circulation anomalies were not well simulated, leading one to infer that the GCM has simulated the drought for the wrong reason. To isolate the causes for this unremarkable circulation, analyzed winds and soil moisture were prescribed in Case 2 and Case 3 as continuous updates by direct replacement of the GCM-predicted fields. These cases show that a large number of simulation biases emanate from wind biases that are carried into the North American region from surroundings regions. Inclusion of soil moisture also helps to ameliorate the strong feedback, perhaps even stronger than that of the real atmosphere, between soil moisture and precipitation. Case 2 simulated one type of surface temperature anomaly pattern, whereas Case 3 with the prescribed soil moisture produced another.

Description: In this paper we describe our plans to better understand this phenomenon using TEC measurements from ground and space-borne GPS receivers.

Description: The objectives of this work are twofold. First, to provide real-time meteorological satellite guidance to airborne field missions for NASA's Upper Atmosphere Research Program, the Global Tropospheric Experiment, and the Atmospheric Effects of Aviation Project. Extensive meteorological satellite datasets were provided for use by the mission scientist and by the science team. These same data were then archived for postdeployment data analysis by the science team. Second, to provide scientific analysis of the data from the airborne field missions supported. The results of these analyses were made public through peer-reviewed publications.

Description: A new web browser for viewing and manipulating meteorological data sets is located on a web server at NASA, Langley Research Center. The browser uses a live access server (LAS) developed by the Thermal Modeling and Analysis Project at NOAA's Pacific Marine Environmental Laboratory. LAS allows researchers to interact directly with the data to view, select, and subset the data in terms of location (latitude, longitude) and time such as day, month, or year. In addition, LAS can compare two data sets and can perform averages and variances, LAS is used here to show how it functions as an internet/web browser for use by the scientific and educational community. In particular its versatility in displaying and manipulating data sets of atmospheric measurements in the earth s radiation budget (ERB) or energy balance, which includes measurements of absorbed solar radiation, reflected shortwave radiation (RSW), thermal outgoing longwave radiation (OLR), and net radiation is demonstrated. These measurements are from the Clouds and the Earth s Radiant Energy System (CERES) experiment and the surface radiation budget (SRB) experiment.

Description: Urban heat islands (UHIs) are caused by the heat-retaining properties of surfaces usually found in urban cities like asphalt and concrete. The UHI can typically be observed on the evening TV weather map as warmer temperatures over the downtown of major cities and cooler temperatures in the suburbs and surrounding rural areas. The UHI has now become a widely acknowledged, observed, and researched phenomenon because of its broad environmental and societal implications. Interest in the UHI will intensify in the future as existing urban areas expand and rural areas urbanize. By the year 2025, more than 60% of the world s population will live in cities, with higher percentages expected in developed nations. The urban growth rate in the United States, for example, is estimated to be 12.5%, and the recent 2000 Census found that more than 80% of the population currently lives in urban areas. Furthermore, the U.S. population is not only growing but is tending to concentrate more in urban areas within the environmentally sensitive coastal zones. Urban growth creates unique and often contentious issues for policymakers related to land use zoning, transportation planning, agricultural production, housing and development, pollution, and natural resources protection. Urban expansion and its associated TJHIs also have measurable impacts on weather and climate processes. The UHI has been documented to affect local and regional temperature, wind patterns, and air quality

Description: A general circulation model (GCM) that includes water vapor tracer (WVT) diagnostics is used to delineate the dominant sources of water vapor for precipitation during the North American monsoon. A 15-year model simulation carried out with one-degree horizontal resolution and time varying sea surface temperature is able to produce reasonable large-scale features of the monsoon precipitation. Within the core of the Mexican monsoon, continental sources provide much of the water for precipitation. Away from the Mexican monsoon (eastern Mexico and Texas), continental sources generally decrease with monsoon onset. Tropical Atlantic Ocean sources of water gain influence in the southern Great Plains states where the total precipitation decreases during the monsoon onset. Pacific ocean sources do contribute to the monsoon, but tend to be weaker after onset. Evaluating the development of the monsoons, soil water and surface evaporation prior to monsoon onset do not correlate with the eventual monsoon intensity. However, the most intense monsoons do use more local sources of water than the least intense monsoons, but only after the onset. This suggests that precipitation recycling is an important factor in monsoon intensity.

Description: Atmospheric data assimilation is the name scientists give to the techniques of blending atmospheric observations with atmospheric model results to obtain an accurate idea of what the atmosphere looks like at any given time. Because two pieces of information are used, observations and model results, the outcomes of data assimilation procedure should be better than what one would get by using one of these two pieces of information alone. There is a number of different mathematical techniques that fall under the data assimilation jargon. In theory most these techniques accomplish about the same thing. In practice, however, slight differences in the approaches amount to faster algorithms in some cases, more economical algorithms in other cases, and even give better overall results in yet some other cases because of practical uncertainties not accounted for by theory. Therefore, the key is to find the most adequate data assimilation procedure for the problem in hand. In our Data Assimilation group we have been doing extensive research to try and find just such data assimilation procedure. One promising possibility is what we call retrospective iterated analysis (RIA) scheme. This procedure has recently been implemented and studied in the context of a very large data assimilation system built to help predict and study weather and climate. Although the results from that study suggest that the RIA scheme produces quite reasonable results, a complete evaluation of the scheme is very difficult due to the complexity of that problem. The present work steps back a little bit and studies the behavior of the RIA scheme in the context of a small problem. The problem is small enough to allow full assessment of the quality of the RIA scheme, but it still has some of the complexity found in nature, namely, its chaotic-type behavior. We find that the RIA performs very well for this small but still complex problem which is a result that seconds the results of our early studies.

Description: The current status of monitoring global precipitation amounts and patterns is described using data sets from the Global Precipitation Climatology Project (GPCP) of the World Climate Research Program (WCRP) and from recent research satellites, especially the Tropical Rainfall Measuring Mission (TRMM). The GPCP monthly (and pentad) data set is a 23-year, globally complete precipitation analysis that is used to explore global and regional variations and trends. The data set is a blend of data mainly from low-orbit microwave satellites and geosynchronous infrared satellites, with additional input from satellite sounder data, Outgoing Longwave Radiation (OLR) data and raingauges. The monthly GPCP data set shows no significant global trend in precipitation over the twenty years, unlike the positive trend in global surface temperatures over the past century. Regional trends are also analyzed. A trend pattern that is a combination of both El Nino and La Nina precipitation features is evident in the 23-year data set. This pattern is related to an increase with time in the number of combined months of El Nino and La Nina during the 23-year period. This apparent trend may be a short-term variation, but also might be related to the increase with time of extreme precipitation events reported elsewhere. Patterns of precipitation variation related to ENSO and other phenomena are shown with clear signals extending from the Tropics into middle and high latitudes of both hemispheres. Also shown, as an example of higher time resolution data is the GPCP daily analysis, which is available for the last six years. A second focus of the talk is on TRMM precipitation data and how these newer data sets incorporating information from the first space-borne meteorological radar compare with the established GPCP data sets.

Description: The United States Great Plains (USGP) experienced a number of multi-year droughts during the last century, most notably the droughts of the 1930s and 1950s. This study examines the causes of such droughts using ensembles of long term (1930-1999) simulations carried out with the NASA Seasonal-to-Interannual Prediction Project (NSIPP-1) atmospheric general circulation model (AGCM) forced with observed sea surface temperatures (SSTs). The results show that the model produces long-term (multi-year) variations in the USGP precipitation that are similar to those observed. A correlative analysis suggests that the ensemble mean low frequency (time scales longer than about 6 years) rainfall variations in the USGP are linked to a pan-Pacific pattern of SST variability that is the leading empirical orthogonal function (EOF) in the low frequency SST data. The link between the SST and the Great Plains precipitation is confirmed in idealized AGCM simulations, in which the model is forced by the 2 polarities of the pan-Pacific SST pattern. The idealized simulations further show that it is primarily the tropical part of the SST anomalies that influence the USGP. As such, the USGP tend to have above normal precipitation when the tropical Pacific SSTs are above normal, while there is a tendency for drought when the tropical SSTs are cold. The upper tropospheric response to the pan-Pacific SST EOF shows a global-scale pattern with a strong wave response in the Pacific and a substantial zonally-symmetric component in which USGP pluvial (drought) conditions are associated with reduced (enhanced) heights throughout the extra-tropics. The potential predictability of rainfall in the USGP associated with SSTs is rather modest, with on average about 1/3 of the total low frequency rainfall variance forced by SST anomalies. Further idealized experiments with climatological SST, suggest that the remaining low frequency variance in the USGP precipitation is the result of interactions with soil moisture. In particular, simulations with soil moisture feedback show a six-fold increase in the variance in annual USGP precipitation compared with simulations in which the soil feedback is excluded. In addition to increasing variance, the interactions with the soil introduce year-to-year memory in the hydrological cycle that is consistent with a red noise process, in which the low frequencies in the deep soil are the result of integrating a net forcing (precipitation-evaporation-runoff) that is white noise on interannual time scales. As such, the role of low frequency SST variability is to introduce a bias to the net forcing on the soil moisture that drives the random process preferentially to either wet or dry conditions.

Description: Information on the turbulent fluxes of momentum, moisture, and heat at the air-sea interface is essential in improving model simulations of climate variations and in climate studies. We have derived a 13.5-year (July 1987-December 2000) dataset of daily surface turbulent fluxes over global oceans from the Special Sensor Mcrowave/Imager (SSM/I) radiance measurements. This dataset, version 2 Goddard Satellite-based Surface Turbulent Fluxes (GSSTF2), has a spatial resolution of 1 degree x 1 degree latitude-longitude and a temporal resolution of 1 day. Turbulent fluxes are derived from the SSM/I surface winds and surface air humidity, as well as the 2-m air and sea surface temperatures (SST) of the NCEP/NCAR reanalysis, using a bulk aerodynamic algorithm based on the surface layer similarity theory.

Description: This is an invited review paper intended to be published as a Chapter in a book entitled "The Global Climate System: Patterns, Processes and Teleconnections" Cambridge University Press. The author begins with an introduction followed by a primer of climate models, including a description of various modeling strategies and methodologies used for climate diagnostics and predictability studies. Results from the CLIVAR Monsoon Model Intercomparison Project (MMIP) were used to illustrate the application of the strategies to modeling the Asian monsoon. It is shown that state-of-the art atmospheric GCMs have reasonable capability in simulating the seasonal mean large scale monsoon circulation, and response to El Nino. However, most models fail to capture the climatological as well as interannual anomalies of regional scale features of the Asian monsoon. These include in general over-estimating the intensity and/or misplacing the locations of the monsoon convection over the Bay of Bengal, and the zones of heavy rainfall near steep topography of the Indian subcontinent, Indonesia, and Indo-China and the Philippines. The intensity of convection in the equatorial Indian Ocean is generally weaker in models compared to observations. Most important, an endemic problem in all models is the weakness and the lack of definition of the Mei-yu rainbelt of the East Asia, in particular the part of the Mei-yu rainbelt over the East China Sea and southern Japan are under-represented. All models seem to possess certain amount of intraseasonal variability, but the monsoon transitions, such as the onset and breaks are less defined compared with the observed. Evidences are provided that a better simulation of the annual cycle and intraseasonal variability is a pre-requisite for better simulation and better prediction of interannual anomalies.

Description: Reliable cloud pressure estimates are needed for accurate retrieval of ozone and other trace gases using satellite-borne backscatter ultraviolet (buv) instruments such as the global ozone monitoring experiment (GOME). Cloud pressure can be derived from buv instruments by utilizing the properties of rotational-Raman scattering (RRS) and absorption by O2-O2. In this paper we estimate cloud pressure from GOME observations in the 355-400 nm spectral range using the concept of a Lambertian-equivalent reflectivity (LER) surface. GOME has full spectral coverage in this range at relatively high spectral resolution with a very high signal-to-noise ratio. This allows for much more accurate estimates of cloud pressure than were possible with its predecessors SBUV and TOMS. We also demonstrate the potential capability to retrieve chlorophyll content with full-spectral buv instruments. We compare our retrieved LER cloud pressure with cloud top pressures derived from the infrared ATSR instrument on the same satellite. The findings confirm results from previous studies that showed retrieved LER cloud pressures from buv observations are systematically higher than IR-derived cloud-top pressure. Simulations using Mie-scattering radiative transfer algorithms that include O2-O2 absorption and RRS show that these differences can be explained by increased photon path length within and below cloud.

Description: Weatherwise is probably the most popular newstand magazine focusing on the subject of weather. It is published six times per year and includes features on weather, climate, and technology. This article (to appear in the January/February Issue) provides a comprehensive review of NASA s past, present, and future contributions in satellite remote sensing for weather and climate processes. The article spans the historical strides of the TIROS program through the scientific and technological innovation of Earth Observer-3 and Global Precipitation Measurement (GPM). It is one of the most thorough reviews of NASA s weather and climate satellite efforts to appear in the popular literature.

Description: Global distribution of aeolian dust is simulated from 1981 to 1996 with the Goddard Ozone Chemistry Aerosol Radiation and Transport (GOCART) model. The results are assessed with in-situ measurements and the Total Ozone Mapping Spectrometer (TOMS) aerosol products. The annual budget over the different continents and oceans are analyzed. It is found that there is a maximum of 25% difference of global annual emission from the minimum in 1996 to the maximum in 1988. There is a downward trend of dust emission over Africa and East Asia, of 6 and 2 Tg/yr, respectively. The inter-annual variability of dust distribution is analyzed over the North Atlantic and Africa. It is found that in winter most of the North Atlantic and Africa dust loading is correlated with the North Atlantic Oscillation. The GOCART model indicates that a controlling factor of such correlation can be attributed to dust emission from the Sahel. The Bodele depression is the major dust source in winter and its inter-annual variability is highly correlated with the NAO. However, it is not possible to conclude without further analysis that the North Atlantic Oscillation is forcing the inter-annual variability of dust emission and in-turn dust concentration over the North Atlantic.

Description: This atlas presents detailed incomparisons of several climatological wind and temperature data sets which cover the middle atmosphere (over altitudes approx. 10-80 km). A number of middle atmosphere climatologies have been developed in the research community based on a variety of meteorological analyses and satellite data sets. Here we present comparisons between these climatological data sets for a number of basic circulation statistics, such as zonal mean temperature, winds and eddy flux statistics. Special attention is focused on tropical winds and temperatures, where large differences exist among separate analyses. We also include comparisons between the global climatologies and historical rocketsonde wind and temperature measurements, and also with more recent lidar temperature data. These comparisons highlight differences and uncertainties in contemporary middle atmosphere data sets, and allow biases in particular analyses to be isolated. In addition, a brief atlas of zonal mean temperature and wind statistics is provided to highlight data availability and as a quick-look reference. This technical report is intended as a companion to the climatological data sets held in archive at the SPARC Data Center (http://www.sparc.sunysb.edu).

Description: We present the basic ideas of the dynamics system of the finite-volume General Circulation Model developed at NASA Goddard Space Flight Center for climate simulations and other applications in meteorology. The dynamics of this model is designed with emphases on conservative and monotonic transport, where the property of Lagrangian conservation is used to maintain the physical consistency of the computational fluid for long-term simulations. As the model benefits from the noise-free solutions of monotonic finite-volume transport schemes, the property of Lagrangian conservation also partly compensates the accuracy of transport for the diffusion effects due to the treatment of monotonicity. By faithfully maintaining the fundamental laws of physics during the computation, this model is able to achieve sufficient accuracy for the global consistency of climate processes. Because the computing algorithms are based on local memory, this model has the advantage of efficiency in parallel computation with distributed memory. Further research is yet desirable to reduce the diffusion effects of monotonic transport for better accuracy, and to mitigate the limitation due to fast-moving gravity waves for better efficiency.

Description: Global maps of rainfall are of great importance in connection with modeling of the earth s climate. Comparison between the maps of rainfall predicted by computer-generated climate models with observation provides a sensitive test for these models. To make such a comparison, one typically needs the total precipitation amount over a large area, which could be hundreds of kilometers in size over extended periods of time of order days or months. This presents a difficult problem since rain varies greatly from place to place as well as in time. Remote sensing methods using ground radar or satellites detect rain over a large area by essentially taking a series of snapshots at infrequent intervals and indirectly deriving the average rain intensity within a collection of pixels , usually several kilometers in size. They measure area average of rain at a particular instant. Rain gauges, on the other hand, record rain accumulation continuously in time but only over a very small area tens of centimeters across, say, the size of a dinner plate. They measure only a time average at a single location. In making use of either method one needs to fill in the gaps in the observation - either the gaps in the area covered or the gaps in time of observation. This involves using statistical models to obtain information about the rain that is missed from what is actually detected. This paper investigates such a statistical model and validates it with rain data collected over the tropical Western Pacific from ship borne radars during TOGA COARE (Tropical Oceans Global Atmosphere Coupled Ocean-Atmosphere Response Experiment). The model incorporates a number of commonly observed features of rain. While rain varies rapidly with location and time, the variability diminishes when averaged over larger areas or longer periods of time. Moreover, rain is patchy in nature - at any instant on the average only a certain fraction of the observed pixels contain rain. The fraction of area covered by rain decreases, as the size of a pixel becomes smaller. This means that within what looks like a patch of rainy area in a coarse resolution view with larger pixel size, one finds clusters of rainy and dry patches when viewed on a finer scale. The model makes definite predictions about how these and other related statistics depend on the pixel size. These predictions were found to agree well with data. In a subsequent second part of the work we plan to test the model with rain gauge data collected during the TRMM (Tropical Rainfall Measuring Mission) ground validation campaign.

Description: This paper describes the LASE system and presents measurements acquired during CAMEX-3 and CAM EX-4.

Description: The 1997/98 is a strong El Nino warm event, while the 1998/99 is a moderate La Nina cold event. We have investigated surface heat budgets and sea surface temperature (SST) tendency for these two events in the tropical western Pacific and eastern Indian Oceans using satellite-retrieved surface radiative and turbulent fluxes. The radiative fluxes are taken from the Goddard Satellite-retrieved Surface Radiation Budget (GSSRB), derived from radiance measurements of the Japanese Geostationary Meteorological Satellite 5. The GSSRB covers the domain 40 deg S - 4 deg N, 90 deg E-17 deg W and a period from October 1997 to December 2000. The spatial resolution is 0.5 deg x 0.5 deg lat-long and the temporal resolution is 1 day. The turbulent fluxes are taken from Version 2 of the Goddard Satellite-based Surface Turbulent Fluxes (GSSTF-2). The GSSTF-2 has a spatial resolution of 1 deg x 1 deg lat-long over global Oceans and a temporal resolution of 1 day covering the period July 1987-December 2000. Daily turbulent fluxes are derived from the S S M (Special Sensor Microwave/Imager) surface wind and surface air humidity, and the SST and 2-m air temperature of the NCEP/NCAR reanalysis, using a stability-dependent bulk flux algorithm. The changes of surface heat budgets, SST and tendency, cloudiness, wind speed, and zonal wind stress of the 1997/98 El Nino relative to the1998/99 La Nina for the northern winter and spring seasons are analyzed. The relative changes of surface heat budgets and SST tendency of the two events are quite different between the tropical eastern Indian and western Pacific Oceans. For the tropical western Pacific, reduced solar heating (more clouds) is generally associated with decreased evaporative cooling (weaker winds), and vise versa. The changes in evaporative cooling over-compensate that of solar heating and dominate the spatial variability of the changes in net surface heating. Both solar heating and evaporative cooling offset each other to reduce interannual variability of net surface heating. In addition, the area of increased SST tendency is larger than that of increased net surface heating, due to less solar radiation penetration through the bottom of deeper ocean mixed layer (stronger winds). For the tropical eastern Indian Ocean, enhanced solar heating (less clouds) is generally associated with reduced evaporative cooling (weaker winds). Both solar heating and evaporative cooling reinforce each other to increase interannual variability of net surface heating. In addition, the area of increased SST tendency is smaller than that of increased net surface heating in the southern domain. The relative changes in wind and zonal wind stress indicate more solar radiation penetration through the ocean mixed layer and more northward heat transport by Ocean current from the south to the north Indian Ocean for the El Nino than for the La Nina.

Description: Progress in understanding of the role of water in global weather and climate is currently limited by our knowledge of the spatial and temporal variability of primary hydrological fields such as precipitation and evaporation. The Tropical Rainfall Measuring Mission (TRMM) has recently demonstrated that use of microwave-based rainfall observations from space in data assimilation can provide better climate data sets and improve short-range weather forecasting. At NASA, we have been exploring non-traditional approaches to assimilating TRMM Microwave Imager (TMI) and Special Sensor Microwavehager (SSM/I) surface rain rate and latent heating profile information in global systems. In this talk we show that assimilating microwave rain rates using a continuous variational assimilation scheme based on moisture tendency corrections improves quantitative precipitation estimates (QPE) and related clouds, radiation energy fluxes, and large-scale circulations in the Goddard Earth Observing System (GEOS) reanalyses. Short-range forecasts initialized with these improved analyses also yield better QPE scores and storm track predictions for Hurricanes Bonnie and Floyd. We present a status report on current efforts to assimilate convective and stratiform latent heating profile information within the general variational framework of model parameter estimation to seek further improvements. Within the next 5 years, there will be a gradual increase in microwave rain products available from operational and research satellites, culminating to a target constellation of 9 satellites to provide global rain measurements every 3 hours with the proposed Global Precipitation Measurement (GPM) mission in 2007/2008. Based on what has been learned from TRMM, there is a high degree of confidence that these observations can play a'major role in improving weather forecasts and producing better global datasets for understanding the Earth's water and energy cycle. The key to success is to adopt an integrated approach to retrieval, validation, modeling, and data assimilation in a coordinated end-to-end observation-application program.

Description: The main scientific goal of the GPM mission, currently planned for start in the 2007 time frame, is to investigate important scientific problems arising within the context of global and regional water cycles. These problems cut across a hierarchy of scales and include climate-water cycle interactions, techniques for improving weather and climate predictions, and better methods for combining observed precipitation with hydrometeorological prediction models for applications to hazardous flood-producing storms, seasonal flood/draught conditions, and fresh water resource assessments. The GPM mission will expand the scope of precipitation measurement through the use of a constellation of some 9 satellites, one of which will be an advanced TRMM-like "core" satellite carrying a dual-frequency Ku-Ka band precipitation radar and an advanced, multifrequency passive microwave radiometer with vertical-horizontal polarization discrimination. The other constellation members will include new dedicated satellites and co-existing Operational/research satellites carrying similar (but not identical) passive microwave radiometers. The goal of the constellation is to achieve approximately 3-hour sampling at any spot on the globe. The constellation's orbit architecture will consist of a mix of sun-synchronous and non-sun-synchronous satellites with the core satellite providing measurements of cloud-precipitation microphysical processes plus calibration-quality rainrate retrievals to be used with the other retrieval information to ensure bias-free constellation coverage. GPM is organized internationally, currently involving a partnership between NASA in the US and the National Space Development Agency in Japan. Additionally, the program is actively pursuing agreements with other international partners and domestic scientific agencies and institutions, as well as participation by individual scientists from academia, government, and the private sector to fulfill mission goals and to pave the way for what ultimately is expected to become an internationally-organized operational global precipitation observing system. Notably, the broad societal applications of GPM are reflected in the United Nation s identification of this mission as a foremost candidate for its Peaceful Uses of Space Program. In this presentation, an overview of the GPM mission design will be presented, followed by an explanation of its scientific agenda as an outgrowth of making improvements in rain retrieval accuracy, microphysics dexterity, sampling frequency, and global coverage. All of these improvements offer new means to observe variability in precipitation and water cycle fluxes and to achieve improved predictability of weather, climate, and hydrometeorology. Specifically, the scientific agenda of GPM has been designed to leverage the measurement improvements to improve prognostic model performance, particularly quantitative precipitation forecasting and its linked phenomena at short, intermediate, and extended time scales. The talk will address how GPM measurements will enable better detection of accelerations and decelerations in regional and global water cycle processes and their relationship to climate variability, better impacts of precipitation data assimilation on numerical weather prediction and global climate reanalysis, and better performance from basin scale hydrometeorological models for short and long term flood-drought forecasting and seasonal fresh water resource assessment. Improved hydrometeorological forecasting will be possible by using continuous global precipitation observations to obtain better closure in water budgets and to generate more realistic forcing of the models themselves to achieve more accurate estimates of interception, infiltration, evaporation/transpiration fluxes, storage, and runoff.

Description: In the talk I shall review the MODIS use of spectral information to derive aerosol size distribution, optical thickness and reflected spectral flux. The accuracy and validation of the MODIS products will be discussed. A few applications will be shown: inversion of combined MODIS+lidar data, aerosol Anthropogenic direct forcing, and dust deposition in the Atlantic Ocean. I shall also discuss the aerosol information that MODIS is measuring: real ref index, single scattering albedo, size of fine and coarse modes, and describe the AEROSAT concept that uses bright desert and glint to derive aerosol absorption.

Description: A recent publication by Shepherd et al. (2002) demonstrated the feasibility of using TRMM precipitation radar (PR) estimates to identify precipitation anomalies caused by urbanization. The approach is particularly useful for investigating this global process because TRMM data span large portions of the globe and comprise an extended temporal dataset. Recent literature suggests that urbanized regions of Houston, Texas may be influencing lightning and precipitation formation over and downwind of the city. Possible mechanisms include: (1) enhanced convergence through interactions between the sea breeze, Galveston bay breeze, and urban heat island circulations, (2) enhanced convergence due to increased surface roughness over the city and/or destabilization of the boundary layer by the UHI, or (3) enhanced cloud condensation nuclei due to urban and industrial aerosol sources. In this study, a downscaling analysis of spatial and temporal trends in rainfall around the Houston Area is being conducted. The downscaling analysis concept involves identifying and quantifying urban rainfall anomalies at progressively smaller spatial and temporal scales using the TRMM satellite, ground-based radar, and a dense network of rain gauges. The goal is to test the hypothesis that the Houston urban district and regions in the climatological downwind region of the city exhibit enhanced rainfall amounts relative to the climatological upwind regions. TRMM was launched in 1997 and currently operates in a low-inclination (35 deg), non-sun-synchronous orbit at an altitude of 402 km (350 km prior to August 2001). The satellite analysis follows the methodologies described in Shepherd et al. (2002). Nearly five years of TRMM PR-derived mean monthly rainfall estimates are utilized to produce annual and warm season isohyetal analyses around Houston. Early results indicate that rainfall rates (mm/h) for the entire period are largest within 100 km northeast and east of Houston (e.g. the "hypothesized downwind region"). The mean rainfall rate over the Houston urban center is 30.5% larger than the upwind control region. The mean rainfall rate in the downwind region is 34.4% larger than the upwind region. An analysis of a parameter called the urban rainfall ratio (URR) illustrates that 65% (88%) of the satellite-derived rainfall rates in the downwind (upwind control) region are greater (less) than the mean background rainfall rate of the entire study region. When the data is stratified by summer months from 1998 to 2001 (June-August), even greater influence over and downwind of the urban area is observed in the statistics. This result is consistent with published reports of urban-generated rainfall being more prevalent in the warm season. The research demonstrates that the evolving TRMM satellite climatology is a credible way to detect mesoscale precipitation signatures that may be linked to urbanization. Early results also corroborate recent findings on Houston-induced convection/drainfall anomalies. Burian and Shepherd will report on other aspects of the downscaling analysis in future forums, but early rain gauge results are consistent with the satellite-based observations.

Description: Coordinated ground, aircraft, and satellite observations are analyzed from the 1999 TRMM Kwajalein Atoll field experiment (KWAJEX) to better understand the relationships between cloud microphysical processes and microwave radiation intensities in the context of physical evaluation of the Level 2 TRMM radiometer rain profile algorithm and uncertainties with its assumed microphysics-radiation relationships. This talk focuses on the results of a multi-dataset analysis based on measurements from KWAJEX surface, air, and satellite platforms to test the hypothesis that uncertainties in the passive microwave radiometer algorithm (TMI 2a12 in the nomenclature of TRMM) are systematically coupled and correlated with the magnitudes of deviation of the assumed 3-dimensional microphysical properties from observed microphysical properties. Re-stated, this study focuses on identifying the weaknesses in the operational TRMM 2a12 radiometer algorithm based on observed microphysics and radiation data in terms of over-simplifications used in its theoretical microphysical underpinnings. The analysis makes use of a common transform coordinate system derived from the measuring capabilities of the aircraft radiometer used to survey the experimental study area, i.e., the 4-channel AMPR radiometer flown on the NASA DC-8 aircraft. Normalized emission and scattering indices derived from radiometer brightness temperatures at the four measuring frequencies enable a 2-dimensional coordinate system that facilities compositing of Kwajalein S-band ground radar reflectivities, ARMAR Ku-band aircraft radar reflectivities, TMI spacecraft radiometer brightness temperatures, PR Ku-band spacecraft radar reflectivities, bulk microphysical parameters derived from the aircraft-mounted cloud microphysics laser probes (including liquid/ice water contents, effective liquid/ice hydrometeor radii, and effective liquid/ice hydrometeor variances), and rainrates derived from any of the individual ground, aircraft, or satellite algorithms applied to the radar or radiometer measurements, or their combination. The results support the study's underlying hypothesis, particularly in context of ice phase processes, in that the cloud regions where the 2a12 algorithm's microphysical database most misrepresents the microphysical conditions as determined by the laser probes, are where retrieved surface rainrates are most erroneous relative to other reference rainrates as determined by ground and aircraft radar. In reaching these conclusions, TMI and PR brightness temperatures and reflectivities have been synthesized from the aircraft AMPR and ARMAR measurements with the analysis conducted in a composite framework to eliminate measurement noise associated with the case study approach and single element volumes obfuscated by heterogeneous beam filling effects. In diagnosing the performance of the 2a12 algorithm, weaknesses have been found in the cloud-radiation database used to provide microphysical guidance to the algorithm for upper cloud ice microphysics. It is also necessary to adjust a fractional convective rainfall factor within the algorithm somewhat arbitrarily to achieve satisfactory algorithm accuracy.

Description: The North Alabama Lightning Mapping Array became operational in November 2001 as a principal component of a severe weather test bed to infuse new science and technologies into the short-term forecasting of severe and hazardous weather and the warning decision-making process. The LMA project is a collaboration among NASA scientists, National Weather Service (NWS) weather forecast offices (WFOs), emergency managers, and other partners. The time rate-of-change of storm characteristics and life-cycle trending are accomplished in real-time through the second generation Lightning Imaging Sensor Data Applications Display (LISDAD II) system, initially developed in T997 through a collaboration among NASA/MSFC, MIT/Lincoln Lab and the Melbourne, FL WFO. LISDAD II is now a distributed decision support system with a JAVA-based display application that allows anyone, anywhere to track individual storm histories within the Tennessee Valley region of the southeastern U.S. Since the inauguration of the LMA there has been an abundance of severe weather. During 23-24 November 2001, a major tornado outbreak was monitored by LMA in its first data acquisition effort (36 tornadoes in Alabama). Since that time the LMA has collected a vast amount of data on hailstorms and damaging wind events, non-tornadic supercells, and ordinary non-severe thunderstorms. In this paper we provide an overview of LMA observations and discuss future prospects for improving the short-term forecasting of convective weather.

Description: A general circulation model (GCM) relies on various physical parameterizations and provides a solution to the atmospheric equations of motion. A data assimilation system (DAS) combines information from observations with a GCM forecast and produces analyzed meteorological fields that represent the observed atmospheric state. An off-line chemistry and transport model (CTM) can use winds and temperatures from a either a GCM or a DAS. The latter application is in common usage for interpretation of observations from various platforms under the assumption that the DAS transport represents the actual atmospheric transport. Here we compare the transport produced by a DAS with that produced by the particular GCM that is combined with observations to produce the analyzed fields. We focus on transport in the tropics and middle latitudes by comparing the age-of-air inferred from observations of SF6 and CO2 with the age-of-air calculated using GCM fields and DAS fields. We also compare observations of ozone, total reactive nitrogen, and methane with results from the two simulations. These comparisons show that DAS fields produce rapid upward tropical transport and excessive mixing between the tropics and middle latitudes. The unrealistic transport produced by the DAS fields may be due to implicit forcing that is required by the assimilation process when there is bias between the GCM forecast and observations that are combined to produce the analyzed fields. For example, the GCM does not produce a quasi-biennial oscillation (QBO). The QBO is present in the analyzed fields because it is present in the observations, and systematic implicit forcing is required by the DAS. Any systematic bias between observations and the GCM forecast used to produce the DAS analysis is likely to corrupt the transport produced by the analyzed fields. Evaluation of transport in the lower tropical stratosphere in a global chemistry and transport model.

Description: There are two important but observationally uncertain parameters in the grid averaged surface energy budgets of mesoscale models - surface moisture availability and thermal heat capacity. A technique has been successfully developed for assimilating Geostationary Operational Environmental Satellite (GOES) skin temperature tendencies during the mid-morning time frame to improve specification of surface moisture. In a new application of the technique, the use of satellite skin temperature tendencies in early evening is explored to improve specification of the surface thermal heat capacity. Together, these two satellite assimilation constraints have been shown to significantly improve the characterization of the surface energy budget of a mesoscale model on fine spatial scales. The GOES assimilation without the adjusted heat capacity was run operationally during the International H2O Project on a 12-km grid. This paper presents the results obtained when using both the moisture availability and heat capacity retrievals in concert. Preliminary results indicate that retrieved moisture availability alone improved the verification statistics of 2-meter temperature and dew point forecasts. Results from the 1.5 month long study period using the bulk heat capacity will be presented at the meeting.

Description: The TRMM Precipitation Radar (PR) measures the vertical profile of reflectivity from which the surface rain rate is estimated after attenuation corrections in the 2A21 algorithm. Characteristics of the vertical reflectivity profile is important for various reasons ranging from scientific to instrument algorithms. It is well known that different types of precipitation such as stratiform or convection, have different heating profiles. The vertical profile of reflectivity can provide information on precipitation classification. The vertical reflectivity structure also provides information on precipitation processes such as growth and aggregation. In terms of TRMM algorithms, an independent estimate of the vertical profiles are also extremely important since the PR returns can be attenuated in the rain layer near the surface. Corrections for attenuation are required in the lowest few kilometers, necessitating some assumptions about the rain size distributions and the reflectivity profile below the lowest measurement unaffected by the surface return. Furthermore, some assumptions about the vertical reflectivity profile are required for Ground Validation (GV) radars, since their lowest scan may be 1 or more kilometers above the surface. Statistics on the vertical reflectivity and Doppler structure are presented from the ER-2 Doppler Radar (EDOP) which participated in several TRMM field campaigns (TEFLUN-A, TEFLUN-B, and LBA) and CAMEX-3. The ER-2 aircraft overflew diverse precipitation types during these campaigns. EDOP is an X-band (9.6 GHz) radar for which returns are less attenuated than at the TRMM PR frequency. The EDOP profiles are first corrected for attenuation using the SRT method. The data from all the ER-2 campaigns are then classified by type (convection, stratiform, and other) and then statistics were performed on the vertical reflectivity and Doppler profiles in the form of CFAD's. These CFADs are compared and discussed. The computed CFAD's indicate significant differences as a function of precipitation type and location (hurricane versus non-hurricane, Brazil versus Florida). The implications of these profiles will be discussed.

Description: The NASA Global Modeling Initiative (GMI) has completed two 35-year simulations with WMO future baseline boundary conditions that simulate increasing N2O and CH4 emissions and decreasing organic chlorine and bromine emissions. Simulations were done with the GMI offline chemistry and transport model using 1) 1 year of winds from the Finite-Volume General Circulation Model (FV-GCM), repeated for the 35 years, and 2) 1 year of winds from the Finite-Volume Data Assimilation System (FV-DAS), repeated for 35-years. The simulations have full stratospheric chemistry. To understand differences in simulated ozone recoveries, basic transport and circulation differences between these models are evaluated. The distribution of mean age of stratospheric air in the FV-GCM run agrees well with observations in the lower stratosphere but the FV-DAS ages are generally too low. This implies circulation and mixing differences that will affect the distributions of other trace species such as CH4, NO, and the organic halogens, all of which are responding to changing boundary conditions and are involved in ozone loss. Realism of model transport is evaluated, with particular attention given to regions and seasons where ozone recovery is expected. Preliminary results indicate increasing ozone trends in the lowermost stratosphere in summer and in the Antarctic and Arctic lower stratosphere in winter and spring.

Description: Direct current measurements are available near the attachment point from both natural cloud-to-ground lightning and rocket-triggered lightning, but little is known about the rise time and peak amplitude of return-stroke currents aloft. We present, as functions of height, current amplitudes, rise times, and effective propagation velocities that have been estimated with a novel remote-sensing technique from data on 24 subsequent return strokes in six different lightning flashes that were triggering at the NASA Kennedy Space Center, FL, during 1987. The unique feature of this data set is the stereo pairs of still photographs, from which three-dimensional channel geometries were determined previously. This has permitted us to calculate the fine structure of the electric-field-change (E) waveforms produced by these strokes, using the current waveforms measured at the channel base together with physically reasonable assumptions about the current distributions aloft. The computed waveforms have been compared with observed E waveforms from the same strokes, and our assumptions have been adjusted to maximize agreement. In spite of the non-uniqueness of solutions derived by this technique, several conclusions seem inescapable: 1) The effective propagation speed of the current up the channel is usually significantly (but not unreasonably) faster than the two-dimensional velocity measured by a streak camera for 14 of these strokes. 2) Given the deduced propagation speed, the peak amplitude of the current waveform often must decrease dramatically with height to prevent the electric field from being over-predicted. 3) The rise time of the current wave front must always increase rapidly with height in order to keep the fine structure of the calculated field consistent with the observations.

Description: Wintertime observations of MCSs (Mesoscale Convective Systems) over the Sea of Japan - 2001 (WMO-01) were collected from January 12 to February 1, 2001. One of the major objectives is to better understand and forecast snow systems and accompanying disturbances and the associated key physical processes involved in the formation and development of these disturbances. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, wind profilers, radiometers, etc.) during WMO-01 provided a first attempt at investigating the detailed characteristics of convective storms and air pattern changes associated with winter storms over the Sea of Japan region. WMO-01 also provided estimates of the apparent heat source (Q1) and apparent moisture sink (Q2). The vertical integrals of Q1 and Q2 are equal to the surface precipitation rates. The horizontal and vertical adjective components of Q1 and Q2 can be used as large-scale forcing for the Cloud Resolving Models (CRMs). The Goddard Cumulus Ensemble (GCE) model is a CRM (typically run with a 1-km grid size). The GCE model has sophisticated microphysics and allows explicit interactions between clouds, radiation, and surface processes. It will be used to understand and quantify precipitation processes associated with wintertime convective systems over the Sea of Japan (using data collected during the WMO-01). This is the first cloud-resolving model used to simulate precipitation processes in this particular region. The GCE model-simulated WMO-01 results will also be compared to other GCE model-simulated weather systems that developed during other field campaigns (i.e., South China Sea, west Pacific warm pool region, eastern Atlantic region and central USA).

Description: We use kinematic and diabatic back trajectory calculations, driven by winds from a general circulation model (GCM) and two different data assimilation systems (DAS), to compute the age spectrum at three latitudes in the lower stratosphere. The age-spectra are compared to chemical transport model (CTM) calculations, and the mean ages from all of these studies are compared to observations. The age spectra computed using the GCM winds show a reasonably isolated tropics in good agreement with observations; however, the age spectra determined from the DAS differ from the GCM spectra. For the DAS diabatic trajectory calculations there is too much exchange between the tropics and mid-latitudes. The age spectrum is thus too broad and the tropical mean age is too old as a result of mixing older mid latitude air with tropical air. Likewise the mid latitude mean age is too young due to the in mixing of tropical air. The DAS kinematic trajectory calculations show excessive vertical dispersion of parcels in addition to excessive exchange between the tropics and mid latitudes. Because air is moved rapidly to the troposphere from the vertical dispersion, the age spectrum is shifted toward the young side. The excessive vertical and meridional dispersion compensate in the kinematic case giving a reasonable tropical mean age. The CTM calculation of the age spectrum using the DAS winds shows the same vertical and meridional dispersive characteristics of the kinematic trajectory calculation. These results suggest that the current DAS products will not give realistic trace gas distributions for long integrations; they also help explain why the extra tropical mean ages determined in a number of previous DAS driven CTM s are too young compared with observations. Finally, we note trajectory-generated age spectra . show significant age anomalies correlated with the seasonal cycles. These anomalies can be linked to year-to-year variations in the tropical heating rate. The anomalies are suppressed in the CTM spectra suggesting that the CTM transport scheme is too diffusive.

Description: The development of a satellite infrared (IR) technique for estimating convective and stratiform rainfall and its application in studying the diurnal variability of rainfall on a global scale is presented. The Convective-Stratiform Technique (CST), calibrated by coincident, physically retrieved rain rates from the Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR), is applied over the global tropics during 2001. The technique is calibrated separately over land and ocean, making ingenious use of the IR data from the TRMM Visible/Infrared Scanner (VIRS) before application to global geosynchronous satellite data. The low sampling rate of TRMM PR imposes limitations on calibrating IR-based techniques; however, our research shows that PR observations can be applied to improve IR-based techniques significantly by selecting adequate calibration areas and calibration length. The diurnal cycle of rainfall, as well as the division between convective and stratiform rainfall will be presented. The technique is validated using available data sets and compared to other global rainfall products such as Global Precipitation Climatology Project (GPCP) IR product, calibrated with TRMM Microwave Imager (TMI) data. The calibrated CST technique has the advantages of high spatial resolution (4 km), filtering of non-raining cirrus clouds, and the stratification of the rainfall into its convective and stratiform components, the latter being important for the calculation of vertical profiles of latent heating.

Description: Global precipitation analysis covering the last few decades and the impact of the new TRMM precipitation observations are discussed. The 20+ year, monthly, globally complete precipitation analysis of the World Climate Research Program's (WCRP/GEWEX) Global Precipitation Climatology Project (GPCP) is used to explore global and regional variations and trends and is compared to the much shorter TRMM(Tropica1 Rainfall Measuring Mission) tropical data set. A trend pattern that is a combination of both El Nino and La Nina precipitation features is evident in the 20-year data set. This pattern is related to an increase with time in the number of combined months of El Nino and La Nina during the 20 year period. Monthly anomalies of precipitation are related to ENSO variations with clear signals extending into middle and high latitudes of both hemispheres. The GPCP daily, 1deg latitude-longitude analysis, which is available from January 1997 to the present is described and the evolution of precipitation patterns on this time scale related to El Nino and La Nina is described. Finally, a TRMM-based 3-hr analysis is described that uses TRMM to calibrate polar-orbit microwave observations from SSM/I and geosynchronous IR observations and merges the various calibrated observations into a final, 3-hr resolution map. This TRMM standard product will soon be available for the entire TRMM period (January 1998- present). A real-time version of this merged product is being produced and is available at 0.25deg latitude-longitude resolution over the latitude range from 50degN-50degS. Images from this data set can be seen at the U.S. TRMM web site (trmm.gsfc.nasa.gov). Examples will be shown, including its use in monitoring flood conditions and relating weather-scale events to climate variations.

Description: Analysis will be presented which explores the impact of land conditions on monthly to seasonal climate simulations in a variety of atmospheric general circulation models (AGCMs). In one set of experiments, the Geophysical Fluid Dynamics Laboratory (GDFL) AGCM is used to explore the nature of soil-moisture predictability and associated climate predictability as an initial value problem. For another set of experiments, the Center for Ocean Land Atmosphere (COLA) and the Goddard Earth Observing System 2 (GEOS-2) AGCMs are used to investigate the impact of realistic snow initialization and assimilation in retrospective climate forecasts for the northern hemisphere spring (March-June).

Description: Despite the fact that the low-level jet of the southern Great Plains (the GPLLJ) of the U.S. is primarily a nocturnal phenomenon that virtually vanishes during the daylight hours, it is one of the most persistent and stable climatological features of the low-level continental flow during the warm-season months, May through August. We have used significant-level data to validate the skill of the GEOS-1 Data Assimilation System (DAS) in realistically detecting this jet and inferring its structure and evolution. We have then carried out a 15-year reanalysis with the GEOS-1 DAS to determine its climatology and mean diurnal cycle and to study its interannual variability. Interannual anomalies of the meridional flow associated with the GPLLJ are much smaller than the mean diurnal fluctuations, than random intraseasonal anomalies, and than the mean wind itself. There are three maxima of low-level meridional flow variance over the Great Plains and the Gulf of Mexico: a 1.2 m2 s-2 peak over the southeast Texas, to the east and south of the mean velocity peak, a 1.0 m2 s-2 peak over the western Gulf of Mexico, and a .8 m2 s-2 peak over the upper Great Plains (UGP), near the Nebraska/South Dakota border. Each of the three variance maxima corresponds to a spatially coherent, jet-like pattern of low-level flow interannual variability. There are also three dominant modes of interannual variability corresponding to the three variance maxima, but not in a simple one-to-one relationship. Cross-sectional profiles of mean southerly wind over Texas remain relatively stable and recognizable from year to year with only its eastward flank showing significant variability. This variability, however, exhibits a distinct, biennial oscillation during the first six to seven years of the reanalysis period and only then. This intermittent biennial oscillation (IBO, one of the three modes discussed in the previous paragraph) in the lowlevel flow is restricted to the region surrounding eastern Texas and is also evident in the NCEP/NCAR reanalysis data set from about 1978 to 1985 or 1986 and again from 1995 to 2000. It is evident as well in surface pressure in both the GEOS-1 and NCEP/NCAR sets. The interannual anomalies do not necessarily persist uniformly throughout an entire season, but can fluctuate from one part of the season to the next. To estimate the characteristic sub-seasonal time scales for coherence of these fluctuations, we have taken the weekly anomaly of low-level wind at each point of the domain from the climatological average for that given point and that given week of the season and computed the covariance of its fluctuations over all weeks and over all years with the weekly climatological anomaly of the meridional wind at each of the three reference points discussed above. The typical duration of a coherent interannual anomaly within a given warm season increases with decreasing latitude from 2 to 3 weeks over the UGP, to 6 to 7 weeks over eastern Texas. Coherence over the western Gulf of Mexico is intermediate between the two with a typical duration of 4 to 5 weeks. There appears to be evidence that the interannual anomalies over Texas the Gulf propagate to the UGP after a week and those over the Gulf propagate there after 2 to 3 weeks. There also appears to be some reverse propagation of interannual anomalies over the UGP to Texas and to the Gulf after a period of about one week. The interannual anomalies in southerly flow over eastern Texas seem to correlate well with interannual anomalies of surface temperature and (negative) ground wetness and over western Texas.

Description: In the past few years, the capabilities of NOAA Geostationary Operational Environmental Satellites (GOES) have increased dramatically. Hourly vertical sounder data is now generally available, but may be unavailable depending upon cloud conditions, satellite operations, and computer system problems at NOAA#s National Environmental Satellite Display and Information Service (NESDIS). Meteorologists at NESDIS have used vertical sounder data to develop experimental products for forecasting the probability of convective downbursts. The two products of interest are the Microburst Day Predictive Index (MDPI), which provides an indication of microburst potential and the WINDEX which is a forecast of maximum winds assuming a microburst does occur. Data analyses were made for the central Florida convective season, that is, the period beginning May 1 and ending September 30. The MDPI showed significant potential as an aid in forecasting convective downbursts. MDPI calculated from GOES soundings were well correlated with those calculated from Cape Canaveral RAOBs.

Description: The Severe Thunderstorm Observations, Research, and Monitoring network (STORMnet) became operational in 2001 as a test bed to infuse new science and technologies into the severe and hazardous weather forecasting and warning process. STORMnet is collaboration among NASA scientists, National Weather Service (NWS) forecasters, emergency managers and other partners. STORMnet integrates total lightning observations from a ten-station 3-D VHF regional lightning mapping array, the National Lightning Detection Network (NLDN), real-time regional NEXRAD Doppler radar, satellite visible and infrared imagers, and a mobile atmospheric profiling system to characterize storms and their evolution. The storm characteristics and life-cycle trending are accomplished in real-time through the second generation Lightning Imaging Sensor Demonstration and Display (LISDAD II), a distributed processing system with a JAVA-based display application that allows anyone, anywhere to track individual storm histories within the Tennessee Valley region of north Alabama and Tennessee, a region of the southeastern U.S. well known for abundant severe weather.

Description: The National Weather Service Office (WFO) in Huntsville, Alabama (HUN) is slated to begin full-time operations in early 2003. With the opening of the Huntsville WFO, a unique opportunity has arisen for close and productive collaboration with scientists at NASA Marshall Space Flight Center (MSFC) and the University of Alabama Huntsville (UAH). As a part of the collaboration effort, NASA has developed the Short-term Prediction Research and Transition (SPoRT) Center. The mission of the SPoRT center is to incorporate NASA earth science technology and research into the NWS operational environment. Emphasis will be on improving mesoscale and short-term forecasting in the first 24 hours of the forecast period. As part of the collaboration effort, the NWS and NASA will develop an implementation and evaluation plan to streamline the integration of the latest technologies and techniques into the operational forecasting environment. The desire of WFO HUN, NASA, and UAH is to provide a model for future collaborative activities between research and operational communities across the country.

Description: The ThOR mission uses a lightning mapping sensor in geostationary Earth orbit to provide continuous observations of thunderstorm activity over the Americas and nearby oceans. The link between lightning activity and cloud updrafts is the basis for total lightning observations indicating the evolving convective intensification and decay of storms. ThOR offers a national operational demonstration of the utility of real-time total lightning mapping for earlier and more reliable identification of potentially severe and hazardous storms. Regional pilot projects have already demonstrated that the dominance in-cloud lightning and increasing in-cloud lash rates are known to precede severe weather at the surface by tens of minutes. ThOR is currently planned for launch in 2005 on a commercial or research satellite. Real-time data will be provided to selected NWS Weather Forecast Offices and National Centers (EMC/AWC/SPC) for evaluation.

Description: This section on Urban Climates provides a basic understanding of what comprises the urban climate and what factors control the overall development of the urban climate. We also discuss in this section, methods for evaluating urban climate characteristics and forcing functions as well as how the urban heat island effect comes into play as a dynamic influence on urban climatology. Additionally, we examine and discuss the major radiation and energy balance of city (i.e., shortwave and longwave radiation, albedo, net all-wave radiation, total energy balance, and sensible latent, and storage heat) and the interactions of these energy balances with the lower atmosphere. The use of remote sensing to measure urban surface temperatures as a driving force in the development of the urban heat island effect is presented. We also discuss how the overall moisture, precipitation, humidity, and air movement in cities (i,e,, wind speeds and wind direction) and wind environment of the city affects urban climatology.

Description: This paper describes the procedures and algorithms for the laboratory calibration and the field data retrieval of the NASA Langley / Ames Diode Laser Hygrometer as implemented during the NASA Trace-P mission during February to April 2000. The calibration is based on a NIST traceable dewpoint hygrometer using relatively high humidity and short pathlength. Two water lines of widely different strengths are used to increase the dynamic range of the instrument in the course of a flight. The laboratory results are incorporated into a numerical model of the second harmonic spectrum for each of the two spectral window regions using spectroscopic parameters from the HITRAN database and other sources, allowing water vapor retrieval at upper tropospheric and lower stratospheric temperatures and humidity levels. The data retrieval algorithm is simple, numerically stable, and accurate. A comparison with other water vapor instruments on board the NASA DC-8 and ER-2 aircraft is presented.

Description: This study presents a comprehensive examination of the spectrum of storm types and their attributes worldwide (between 35N and 35S latitude), and as a function of season, location, and convective regime using the observed lightning, microwave scattering, and reflectivity signatures from NASA's Tropical Rainfall Measuring Mission (TRMM) low-Earth orbiting observatory. A global, multi-year data set (1998-2000) indicates that the deepest thunderstorms (having reflectivity in excess of 50 dBZ at 9 km altitude) occur in all the sub-tropical continents and occasionally over the open ocean, but are most common over the Americas. The most intense storms have the greatest lightning rates, lowest brightness temperatures and greatest depth of reflectivity- all indicative of strong updrafts and a well-developed volume of precipitation-sized ice particles. Mesoscale convective systems occurring within or in association with forcing from the sub-tropical continents are the most prolific lightning producers. The greatest flash rate to date of 993 flashes per minute was observed by NASA's Lightning Imaging Sensor on May 6, 1999 during an overpass of a pre-frontal squall line extending from Tennessee to Louisiana. The global distribution and frequency of thunderstorms, and the most recent summary of the extreme storms observed from space, in particular, will be discussed in greater detail.

Description: The bi-spectral threshold (BTH) for cloud detection and height assignment is now operational at NASA's Global Hydrology and Climate Center (GHCC). This new approach is similar in principle to the bi-spectral spatial coherence (BSC) method with improvements made to produce a more robust cloud-filtering algorithm for nighttime cloud detection and subsequent 24-hour operational cloud top pressure assignment. The method capitalizes on cloud and surface emissivity differences from the GOES 3.9 and 10.7-micrometer channels to distinguish cloudy from clear pixels. Separate threshold values are determined for day and nighttime detection, and applied to a 20-day minimum composite difference image to better filter background effects and enhance differences in cloud properties. A cloud top pressure is assigned to each cloudy pixel by referencing the 10.7-micrometer channel temperature to a thermodynamic profile from a locally -run regional forecast model. This paper and supplemental poster will present an objective validation of nighttime cloud detection by the BTH approach in comparison with previous methods. The cloud top pressure will be evaluated by comparing to the NESDIS operational CO2 slicing approach.

Description: The process by which liquid cloud droplets homogeneously crystallize into ice is still not well-understood. The ice nucleation process based on the standard and classical theory of homogeneous freezing, initiates within the interior volume of a cloud droplet. Current experimental data on homogeneous freezing rates of ice in droplets of supercooled water, both in air and emulsion oil samples, show considerable scatter. For example, at -33 C, the reported volume-based freezing rates of ice in supercooled water vary by as much as 5 orders of magnitude, which is well outside the range of measurement uncertainties. Here, we show that the process of ice nucleus formation at the air (or oil)-liquid water interface may help to explain why experimental results on ice nucleation rates yield different results in different ambient phases. Our results also suggest that surface crystallization of ice in cloud droplets can explain why low amounts of supercooled water have been observed in the atmosphere near -40 C.