ALBERT

Description: The goal of this study was to determine, through modeling, the impact of aircraft emissions on regional air quality, especially in regard to fine particulate matter (PM(2.5)) as well as ozone and other pollutants. For this, we focused on Hartsfield-Jackson Atlanta International Airport which is the busiest airport in the world based on passenger traffic (AIC, 2003). Hartsfield-Jackson serves the metropolitan Atlanta area where air quality does not meet national standards. Emissions from mobile and industrial sources (including several large electric power generating utilities) are the major contributors to the area's air pollution. In this study, we assessed the impact of Hartsfield-Jackson Airport on air quality around Atlanta, Georgia, and compared it to the impacts of other emission sources in the area. The assessment was built upon other, related air quality studies involving both field and modeling components. To achieve the objectives, first a detailed inventory was developed for aircraft and other emissions at Hartsfield-Jackson Atlanta International Airport. Then, air quality simulations were performed to relate these emissions to regional air quality around Atlanta. The Community Multiscale Air Quality Model (CMAQ) was used as the modeling platform. The period of August 11-20 2000 was selected as the episode to be modeled in this study. Prior modeling of this episode during the Fall Line Air Quality Study (FAQS) and availability of additional PM(2.5) measurements for evaluation played a major role in this selection. Meteorological data for this episode as well as emission data for sources other than aircrafts were already available from FAQS.

Description: The project goals are: Make data analysis faster and cheaper. Increase use of NASA data by removing barriers to data access. Cope with data heterogeneity. Support code reuse and rapid application development. Support multiple applications, users. Including fire and health domains. Improve QOS. Always provide an answer. Tell user how good it is, where it come from.

Description: Several meteorological datasets, including U.K. Met Office (MetO), European Centre for Medium-Range Weather Forecasts (ECMWF), National Centers for Environmental Prediction (NCEP), and NASA's Goddard Earth Observation System (GEOS-4) analyses, are being used in studies of the 2002 Southern Hemisphere (SH) stratospheric winter and Antarctic major warming. Diagnostics are compared to assess how these studies may be affected by the meteorological data used. While the overall structure and evolution of temperatures, winds, and wave diagnostics in the different analyses provide a consistent picture of the large-scale dynamics of the SH 2002 winter, several significant differences may affect detailed studies. The NCEP-NCAR reanalysis (REAN) and NCEP-Department of Energy (DOE) reanalysis-2 (REAN-2) datasets are not recommended for detailed studies, especially those related to polar processing, because of lower-stratospheric temperature biases that result in underestimates of polar processing potential, and because their winds and wave diagnostics show increasing differences from other analyses between similar to 30 and 10 hPa (their top level). Southern Hemisphere polar stratospheric temperatures in the ECMWF 40-Yr Re-analysis (ERA-40) show unrealistic vertical structure, so this long-term reanalysis is also unsuited for quantitative studies. The NCEP/Climate Prediction Center (CPC) objective analyses give an inferior representation of the upper-stratospheric vortex. Polar vortex transport barriers are similar in all analyses, but there is large variation in the amount, patterns, and timing of mixing, even among the operational assimilated datasets (ECMWF, MetO, and GEOS-4). The higher-resolution GEOS-4 and ECMWF assimilations provide significantly better representation of filamentation and small-scale structure than the other analyses, even when fields gridded at reduced resolution are studied. The choice of which analysis to use is most critical for detailed transport studies (including polar process modeling) and studies involving synoptic evolution in the upper stratosphere. The operational assimilated datasets are better suited for most applications than the NCEP/CPC objective analyses and the reanalysis datasets.

Description: The general climate of the island of Puerto Rico is dominated by the easterly trade winds from the Atlantic Ocean, and during synoptically calm days by the topographic and local land surface characteristics [1]. The urban canopy of the metropolitan area of San Juan, capital city of the Island, may introduce a new microclimate that changes the characteristics of the low atmosphere and interacts with the other microclimates already present in the island. The primitive land cover and land use (LCLU) of the metropolitan area of San Juan was composed by broadleaf trees, moist soils, and very dense vegetation in general. The urban LCLU changes the balance for the mass, momentum and energy between the bottom boundary and the lower atmosphere, creating different climate conditions over urban and rural regions. Some of these differences are low relative humidity and high temperatures observed in urban areas when compared to rural areas. These in turn produces a convective circulation over the urban areas, a phenomenon compared to the sea and land breezes, commonly known as heat islands (UHI). Factors that contribute to the formation of the UHI are anthropogenic heat sources, aerosols from pollutants, fast water canalization due to the presence of buildings and streets, among others. The comparison between urban and rural climates is the most common approach to analyze the UHI. These contrasts are larger in clear and calm conditions and tend to disappear in cloudy and windy weather. The UHI was recognized in the early 1950 s as closed isotherms that separates the city from the general temperature field [2]. The impact of the urban LCLU in San Juan, Puerto Rico, was quantified calculating the difference between historical data sets for the air temperature over an identified urban area and a rural area dT(U-R). The analysis of the climatological data revealed that a UHI exists in the metropolitan area of San Juan, Puerto Rico. The data reveals a permanent urban heat island effect present in the SJMA during the year, which is increasing at a rate of 0.41oC/decade. These findings encouraged the planning and execution of an intense field campaign in February 2004 referred as the ATLAS San Juan mission. The focus of the remaining of this report is the analysis of the data for this field campaign.

Description: During its first three years, the Tropical Rainfall Measuring Mission (TRMM) satellite observed nearly six million precipitation features. The population of precipitation features is sorted by lightning flash rate, minimum brightness temperature, maximum radar reflectivity. areal extent, and volumetric rainfall. For each of these characteristics, essentially describing the convective intensity or the size of the features, the population is broken into categories consisting of the top 0.001%, top 0.01%, top 0.1%, top 1%, top 2.4%. and remaining 97.6%. The set of weakest/smallest features composes 97.6% of the population because that fraction does not have detected lightning, with a minimum detectable flash rate of 0.7 flashes (fl) per minute. The greatest observed flash rate is 1351 fl per minute; the lowest brightness temperatures are 42 K (85 GHz) and 69 K (37 GHz). The largest precipitation feature covers 335 000 square kilometers and the greatest rainfall from an individual precipitation feature exceeds 2 x 10 kg per hour of water. There is considerable overlap between the greatest storms according to different measures of convective intensity. The largest storms are mostly independent of the most intense storms. The set of storms producing the most rainfall is a convolution of the largest and the most intense storms. This analysis is a composite of the global Tropics and subtropics. Significant variability is known to exist between locations. seasons, and meteorological regimes. Such variability will be examined in Part II. In Part I, only a crude land-ocean separation is made. The known differences in bulk lightning flash rates over land and ocean result from at least two differences in the precipitation feature population: the frequency of occurrence of intense storms and the magnitude of those intense storms that do occur. Even when restricted to storms with the same brightness temperature, same size, or same radar reflectivity aloft, the storms over water are considerably less likely to produce lightning than are comparable storms over land.

Description: In this paper, GA-based techniques are used to invert the equations of an electromagnetic model based on Dense Medium Radiative Transfer Theory (DMRT) under the Quasi Crystalline Approximation with Coherent Potential to retrieve snow depth, mean grain size and fractional volume from microwave brightness temperatures. The technique is initially tested on both noisy and not-noisy simulated data. During this phase, different configurations of genetic algorithm parameters are considered to quantify how their change can affect the algorithm performance. A configuration of GA parameters is then selected and the algorithm is applied to experimental data acquired during the NASA Cold Land Process Experiment. Snow parameters retrieved with the GA-DMRT technique are then compared with snow parameters measured on field.

Description: Use of dual-wavelength radar, with properly chosen wavelengths, will significantly lessen the ambiguities in the retrieval of microphysical properties of hydrometeors. In this paper, a dual-wavelength algorithm is described to estimate the characteristic parameters of the snow size distributions. An analysis of the computational results, made at X and Ka bands (T-39 airborne radar) and at S and X bands (CP-2 ground-based radar), indicates that valid estimates of the median volume diameter of snow particles, D(sub 0), should be possible if one of the two wavelengths of the radar operates in the non-Rayleigh scattering region. However, the accuracy may be affected to some extent if the shape factors of the Gamma function used for describing the particle distribution are chosen far from the true values or if cloud water attenuation is significant. To examine the validity and accuracy of the dual-wavelength radar algorithms, the algorithms are applied to the data taken from the Convective and Precipitation-Electrification Experiment (CaPE) in 1991, in which the dual-wavelength airborne radar was coordinated with in situ aircraft particle observations and ground-based radar measurements. Having carefully co-registered the data obtained from the different platforms, the airborne radar-derived size distributions are then compared with the in-situ measurements and ground-based radar. Good agreement is found for these comparisons despite the uncertainties resulting from mismatches of the sample volumes among the different sensors as well as spatial and temporal offsets.

Description: Rainfall variability on seasonal and interannual-to-interdecadal time scales in the tropical Atlantic is quantified using a 25-year (1979-2003) monthly rainfall dataset from the Global Precipitation Climatology Project (GPCP). The ITCZ measured by monthly rainfall between 15-37.5 deg W attains its peak as moving to the northernmost latitude (4-10 deg N) during July-September in which the most total rainfall is observed in the tropical Atlantic basin (17.5 deg S-22.5 deg N, 15 deg-37.5 deg W); the ITCZ becomes weakest during January-February with the least total rainfall as it moves to the south. In contrast, rainfall variability on interannual to interdecadal time scales shows a quite different seasonal preference. The most intense interannual variability occurs during March-May when the ITCZ tends to be near the equator and becomes weaker. Significant, negative correlations between the ITCZ strength and latitude anomalies are observed during boreal spring and early summer. The ITCZ strength and total rainfall amount in the tropical Atlantic basin are significantly modulated by the Pacific El Nino and the Atlantic equatorial mode (or Atlantic Nino) particularly during boreal spring and summer; whereas the impact of the Atlantic interhemispheric mode is considerably weaker. Regarding the anomalous latitudes of the ITCZ, the influence can come from both local, i.e., the Atlantic interhemispheric and equatorial modes, and remote forcings, i. e., El Nino; however, a direct impact of El Nino on the latitudes of the ITCZ can only be found during April-July, not in winter and early spring in which the warmest SST anomalies are usually observed in the equatorial Pacific.

Description: An atmospheric general circulation model simulation for 1948-1997 of the water budgets for the MacKenzie, Mississippi and Amazon River basins is presented. In addition to the water budget, we include passive tracers to identify the geographic sources of water for the basins, and the analysis focuses on the mechanisms contributing to precipitation recycling in each basin. While each basin s precipitation recycling has a strong dependency on evaporation during the mean annual cycle, the interannual variability of the recycling shows important relationships with the atmospheric circulation. The MacKenzie River basin has only a weak interannual dependency on evaporation, where the variations in zonal moisture transport from the Pacific Ocean can affect the basin water cycle. On the other hand, the Mississippi River basin has strong interannual dependencies on evaporation. While the precipitation recycling weakens with increased low level jet intensity, the evaporation variations exert stronger influence in providing water vapor for convective precipitation at the convective cloud base. High precipitation recycling is also found to be partly connected to warm SSTs in the tropical Pacific Ocean. The Amazon River basin evaporation exhibits small interannual variations, so that the interannual variations of precipitation recycling are related to atmospheric moisture transport from the tropical south Atlantic Ocean. Increasing SSTs over the 50-year period are causing increased easterly transport across the basin. As moisture transport increases, the Amazon precipitation recycling decreases (without real time varying vegetation changes). In addition, precipitation recycling from a bulk diagnostic method is compared to the passive tracer method used in the analysis. While the mean values are different, the interannual variations are comparable between each method. The methods also exhibit similar relationships to the terms of the basin scale water budgets.

Description: Comparison of polar sea ice results from 11 major global climate models and satellite-derived observations for 1979-2004 reveals that each of the models is simulating seasonal cycles that are phased at least approximately correctly in both hemispheres. Each is also simulating various key aspects of the observed ice cover distributions, such as winter ice not only throughout the central Arctic basin but also throughout Hudson Bay, despite its relatively low latitudes. However, some of the models simulate too much ice, others too little ice (in some cases varying depending on hemisphere and/or season), and some match the observations better in one season versus another. Several models do noticeably better in the Northern Hemisphere than in the Southern Hemisphere, and one does noticeably better in the Southern Hemisphere. In the Northern Hemisphere all simulate monthly average ice extents to within +/-5.1 x 10(exp 6)sq km of the observed ice extent throughout the year; and in the Southern Hemisphere all except one simulate the monthly averages to within +/-6.3 x 10(exp 6) sq km of the observed values. All the models properly simulate a lack of winter ice to the west of Norway; however, most do not obtain as much absence of ice immediately north of Norway as the observations show, suggesting an under simulation of the North Atlantic Current. The spread in monthly averaged ice extents amongst the 11 model simulations is greater in the Southern Hemisphere than in the Northern Hemisphere and greatest in the Southern Hemisphere winter and spring.

Description: Can cloud screening of an aerosol data set, affect the aerosol optical thickness (AOT) climatology? Aerosols, humidity and clouds are correlated. Therefore, rigorous cloud screening can systematically bias towards less cloudy conditions, underestimating the average AOT. Here, using AERONET data we show that systematic rejection of variable atmospheric optical conditions can generate such bias in the average AOT. Therefore we recommend (1) to introduce more powerful spectral variability cloud screening and (2) to change the philosophy behind present aerosol climatologies: Instead of systematically rejecting all cloud contaminations, we suggest to intentionally allow the presence of cloud contamination, estimate the statistical impact of the contamination and correct for it. The analysis, applied to 10 AERONET stations with approx. 4 years of data, shows almost no change for Rome (Italy), but up to a change in AOT of 0.12 in Beijing (PRC). Similar technique may be explored for satellite analysis, e.g. MODIS.

Description: This paper presents a comparison of cloud-top altitude retrieval methods applied to S-HIS (Scanning High Resolution Interferometer Sounder) measurements. Included in this comparison is an improvement to the traditional CO2 Slicing method. The new method, CO2 Sorting, determines optimal channel pairs to apply the CO2 Slicing. Measurements from collocated samples of the Cloud Physics Lidar (CPL) and Modis Airborne Simulator (MAS) instruments assist in the comparison. For optically thick clouds good correlation between the S-HIS and lidar cloud-top retrievals are found. For tenuous ice clouds there can be large differences between lidar (CPL) and S-HIS retrieved cloud-tops. It is found that CO2 Sorting significantly reduces the cloud height biases for the optically thin cloud (total optical depths less then 1.0). For geometrically thick but optically thin cirrus clouds large differences between the S-HIS infrared cloud top retrievals and the CPL detected cloud top where found. For these cases the cloud height retrieved by the S-HIS cloud retrievals correlated closely with the level the CPL integrated cloud optical depth was approximately 1.0.

Description: Numerical model studies have suggested that the ongoing global warming will likely affect tropical cyclone activity. But so far little observed evidence has been detected to support the projected future changes. Using satellite-supported best-track data from 1965 to 2003, we show for the first time that over the past four decades the two prevailing typhoon tracks in the western North Pacific (WNP) have shifted westward significantly; the typhoon activity over the South China Sea has considerably decreased; and East Asia has experienced increasing typhoon influence. Our trajectory model simulation indicates that the long-term shifts in the typhoon tracks result primarily from the changes in the mean translation velocity of typhoons or the large-scale steering flow, which is associated with the westward expansion and strengthening of the WNP subtropical high.

Description: Microwave brightness temperatures of snow covered terrains can be modeled by means of the Dense Radiative Transfer Medium Theory (DMRT). In a dense medium, such as snow, the assumption of independent scattering is no longer valid and the scattering of correlated scatterers must be considered. In the DMRT, this is done considering a pair distribution function of the particles position. In the electromagnetic model, the snowpack is simulated as a homogeneous layer having effective permittivity and albedo calculated through the DMRT. In order to account for clustering of snow crystals, a model of cohesive particles can be applied, where the cohesion between the particles is described by means of a dimensionless parameters called stickiness (z), representing a measure of the inversion of the attraction of the particles. The lower the z the higher the stickiness. In this study, microwave signatures of melting and refreezing cycles of seasonal snowpacks at high altitudes are studied by means of both experimental and modeling tools. Radiometric data were collected 24 hours per day by the University of Michigan Tower Mounted Radiometer System (TMRS). The brightness temperatures collected by means of the TMRS are simulated by means of a multi-layer electromagnetic model based on the dense medium theory with the inputs to the model derived from the data collected at the snow pits and from the meteorological station. The paper is structured as follows: in the first Section the temperature profiles recorded by the meteorological station and the snow pit data are presented and analyzed; in the second Section, the characteristics of the radiometric system used to collect the brightness temperatures are reported together with the temporal behavior of the recorded brightness temperatures; in the successive Section the multi-layer DMRT-based electromagnetic model is described; in the fourth Section the comparison between modeled and measured brightness temperatures is discussed. We dedicate the last Section to the conclusions and future works.

Description: The direct effects of aerosols on global and regional climate during boreal spring are investigated based on simulations using the NASA Global Modeling and Assimilation Office (GMAO) finite-volume general circulation model (fvGCM) with Microphyics of clouds in Relaxed Arakawa Schubert Scheme (McRAS). The aerosol loading are prescribed from three-dimensional monthly distribution of tropospheric aerosols viz., sulfate, black carbon, organic carbon, soil dust, and sea salt from output of the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The aerosol extinction coefficient, single scattering albedo, and asymmetric factor are computed as wavelength-dependent radiative forcing in the radiative transfer scheme of the fvGCM, and as a function of the aerosol loading and ambient relative humidity. We find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excites a planetary scale teleconnection pattern in sea level pressure, temperature and geopotential height spanning North Africa through Eurasia to the North Pacific. Surface cooling due to direct effects of aerosols is found in the vicinity and downstream of the aerosol source regions, i.e., South Asia, East Asia, and northern and western Africa. Additionally, atmospheric heating is found in regions with large loading of dust (over Northern Africa, and Middle East), and black carbon (over South-East Asia). Paradoxically, the most pronounced feature in aerosol-induced surface temperature is an east-west dipole anomaly with strong cooling over the Caspian Sea, and warming over central and northeastern Asia, where aerosol concentration are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection driven by atmospheric heating anomalies induced by absorbing aerosols in the source regions, but the influence was conveyed globally through barotropic energy dispersion and sustained by feedback processes associated with the regional circulations.

Description: The solar cycle (SC) effect in the lower atmosphere has been linked observationally to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves. Salby and Callaghan analyzed the QBO observations covering more than 40 years and found that it contains a relatively large SC signature at 20 km. Following up on a 2D study with our Numerical Spectral Model (NSM), we discuss here a 3D study in which we simulated the QBO under the influence of the SC. For a SC period of 10 years, the amplitude of the relative variations of radiative forcing is taken to vary from 0.2% at the surface to 2% at 50 km to 20% at 100 km and above. This model produces in the lower stratosphere a relatively large modulation of the QBO, which appears to be related to the SC and is in qualitative agreement with the observations. Further studies are needed, (1) to determine whether the effect is real and the results are robust and (2) to explore the mechanism(s) that may amplify the SC effect. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with or aid the SC effect.

Description: The twelve weather and climate models participating in the Global Land-Atmosphere Coupling Experiment (GLACE) show both a wide variation in the strength of land-atmosphere coupling and some intriguing commonalities. In this paper, we address the causes of variations in coupling strength - both the geographic variations within a given model and the model-to-model differences. The ability of soil moisture to affect precipitation is examined in two stages, namely, the ability of the soil moisture to affect evaporation, and the ability of evaporation to affect precipitation. Most of the differences between the models and within a given model are found to be associated with the first stage - an evaporation rate that varies strongly and consistently with soil moisture tends to lead to a higher coupling strength. The first stage differences reflect identifiable differences in model parameterization and model climate. Intermodel differences in the evaporation-precipitation connection, however, also play a key role.

Description: Radiative transfer calculations for nadir-viewing satellites normally assume the atmosphere to be horizontally homogeneous. Yet it has been shown recently that horizontal gradients can lead to significant errors in satellite infrared and microwave soundings. We extend the methodology to backscatter ultra-violet observations of ozone, and present a first estimate of the effect s magnitude. The Solar Backscatter Ultra-Violet/2 (SBUV/2) instrument, a pure nadir sounder, serves as our test bed. Our results indicate that in a vast majority of cases the abovementioned errors can be neglected. However, occurrence of higher errors, particularly at wavelengths longer than 300 nm, coincides with some of the most interesting atmospheric phenomena like tropopause folds and the South polar ozone hole. This leads to a seasonal variation of the magnitude of the effect. Due to the mostly zonal geometry of the ozone distribution, there is also the possibility that biases may be introduced, which is particularly critical if the data are to be assimilated or used to determine trends. The results presented are tested for robustness using different model atmospheres. The influence of horizontal inhomogeneities will be even more pronounced for cross-track sounders and limb viewers, and easier to detect once higher resolution atmospheric models are available. This will be investigated in future studies.

Description: An ideal and simple formulation is successfully derived that well represents a quasi-linear relationship found between the domain-averaged water vapor, q (mm), and temperature, T (K), fields obtained from a series of quasi-equilibrium (long-term) simulations for the Tropics using the two-dimensional Goddard Cumulus Ensemble (GCE) model. Earlier model work showed that the forced maintenance of two different wind profiles in the Tropics leads to two different equilibrium states. Investigating this finding required investigation of the slope of the moisture-temperature relations, which turns out to be linear in the Tropics. The extra-tropical climate equilibriums become more complex, but insight on modeling sensitivity can be obtained by linear stepwise regression of the integrated temperature and humidity. A globally curvilinear moisture-temperature distribution, similar to the famous Clausius-Clapeyron curve (i.e., saturated water vapor pressure versus temperature), is then found in this study. Such a genuine finding clarifies that the dynamics are crucial to the climate (shown in the earlier work) but the thermodynamics adjust. The range of validity of this result is further examined herein. The GCE-modeled tropical domain-averaged q and T fields form a linearly-regressed "q-T" slope that genuinely resides within an ideal range of slopes obtained from the aforementioned formulation. A quantity (denoted as dC2/dC1) representing the derivative between the static energy densities due to temperature (C2) and water vapor (C1) for various quasi-equilibrium states can also be obtained. A dC2/dC1 value near unity obtained for the GCE-modeled tropical simulations implies that the static energy densities due to moisture and temperature only differ by a pure constant for various equilibrium states. An overall q-T relation also including extra-tropical regions is, however, found to have a curvilinear relationship. Accordingly, warm/moist regions favor change in water vapor faster than temperature, while cold/dry regions favor an increase in temperature quicker than water vapor.

Description: Land surface models (LSMs) are computer programs, similar to weather and climate prediction models, which simulate the stocks and fluxes of water (including soil moisture, snow, evaporation, and runoff) and energy (including the temperature of and sensible heat released from the soil) after they arrive on the land surface as precipitation and sunlight. It is not currently possible to measure all of the variables of interest everywhere on Earth with sufficient accuracy and space-time resolution. Hence LSMs have been developed to integrate the available observations with our understanding of the physical processes involved, using powerful computers, in order to map these stocks and fluxes as they change in time. The maps are used to improve weather forecasts, support water resources and agricultural applications, and study the Earth's water cycle and climate variability. NASA's Global Land Data Assimilation System (GLDAS) project facilitates testing of several different LSMs with a variety of input datasets (e.g., precipitation, plant type). Precipitation is arguably the most important input to LSMs. Many precipitation datasets have been produced using satellite and rain gauge observations and weather forecast models. In this study, seven different global precipitation datasets were evaluated over the United States, where dense rain gauge networks contribute to reliable precipitation maps. We then used the seven datasets as inputs to GLDAS simulations, so that we could diagnose their impacts on output stocks and fluxes of water. In terms of totals, the Climate Prediction Center (CPC) Merged Analysis of Precipitation (CMAP) had the closest agreement with the US rain gauge dataset for all seasons except winter. The CMAP precipitation was also the most closely correlated in time with the rain gauge data during spring, fall, and winter, while the satellitebased estimates performed best in summer. The GLDAS simulations revealed that modeled soil moisture is highly sensitive to precipitation, with differences in spring and summer as large as 45% depending on the choice of precipitation input.

Description: Cloud development, the onset of precipitation and the effect of aerosol on clouds depend on the structure of the cloud profiles of droplet size and phase. Aircraft measurements of cloud profiles are limited in their temporal and spatial extent. Satellites were used to observe cloud tops not cloud profiles with vertical profiles of precipitation-sized droplets anticipated from Cloudsat. The recently proposed CLAIM-3D satellite mission (cloud aerosol interaction mission in 3D) suggests to measure profiles of cloud microphysical properties by retrieving them from the solar and infrared radiation reflected or emitted from cloud sides. Inversion of measurements from the cloud sides requires rigorous understanding of the 3-dimensional (3D) properties of clouds. Here we discuss the reflected sunlight from the cloud sides and top at two wavelengths: one nonabsorbing to solar radiation (0.67 micrometers) and one with liquid water efficient absorption of solar radiation (2.1 micrometers). In contrast to the plane-parallel approximation, a conventional approach to all current operational retrievals, 3D radiative transfer is used for interpreting the observed reflectances. General properties of the radiation reflected from the sides of an isolated cloud are discussed. As a proof of concept, the paper shows a few examples of radiation reflected from cloud fields generated by a simple stochastic cloud model with the prescribed vertically resolved microphysics. To retrieve the information about droplet sizes, we propose to use the probability density function of the droplet size distribution and its first two moments instead of the assumption about fixed values of the droplet effective radius. The retrieval algorithm is based on the Bayesian theorem that combines prior information about cloud structure and microphysics with radiative transfer calculations.

Description: This paper presents a comparison of volume extinction coefficients in tropical ice clouds retrieved from two instruments : the 532-nm Cloud Physics Lidar (CPL), and the in-situ probe Cloud Integrating Nephelometer (CIN). Both instruments were mounted on airborne platforms during the CRYSTAL-FACE campaign and took measurements in ice clouds up to 17km. Coincident observations from three cloud cases are compared : one synoptically-generated cirrus cloud of low optical depth, and two ice clouds located on top of convective systems. Emphasis is put on the vertical variability of the extinction coefficient. Results show small differences on small spatial scales (approx. 100m) in retrievals from both instruments. Lidar retrievals also show higher extinction coefficients in the synoptic cirrus case, while the opposite tendency is observed in convective cloud systems. These differences are generally variations around the average profile given by the CPL though, and general trends on larger spatial scales are usually well reproduced. A good agreement exists between the two instruments, with an average difference of less than 16% on optical depth retrievals.

Description: How do the continents affect large-scale hydrological cycles? How important can one continent be to the climate system? To address these questions, 4-years of National Aeronautics and Space Administration (NASA) Terra Moderate Resolution Imaging Spectroradiometer (MODIS) observations, Tropical Rainfall Measuring Mission (TRMM) observations, and the Global Precipitation Climatology Project (GPCP) global precipitation analysis, were used to assess the land impacts on clouds, rainfall, and water vapor at continental scales. At these scales, the observations illustrate that continents are integrated regions that enhance the seasonality of atmospheric and surface hydrological parameters. Specifically, the continents of Eurasia and North America enhance the seasonality of cloud optical thickness, cirrus fraction, rainfall, and water vapor. Over land, both liquid water and ice cloud effective radii are smaller than over oceans primarily because land has more aerosol particles. In addition, different continents have similar impacts on hydrological variables in terms of seasonality, but differ in magnitude. For example, in winter, North America and Eurasia increase cloud optical thickness to 17.5 and 16, respectively, while in summer, Eurasia has much smaller cloud optical thicknesses than North America. Such different land impacts are determined by each continent s geographical condition, land cover, and land use. These new understandings help further address the land-ocean contrasts on global climate, help validate global climate model simulated land-atmosphere interactions, and help interpret climate change over land.

Description: The radar return powers from a three-frequency radar, with center frequency at 22.235 GHz and upper and lower frequencies chosen with equal water vapor absorption coefficients, can be used to estimate water vapor density and parameters of the precipitation. A linear combination of differential measurements between the center and lower frequencies on one hand and the upper and lower frequencies on the other provide an estimate of differential water vapor absorption. The coupling between the precipitation and water vapor estimates is generally weak but increases with bandwidth and the amount of non-Rayleigh scattering of the hydrometeors. The coupling leads to biases in the estimates of water vapor absorption that are related primarily to the phase state and the median mass diameter of the hydrometeors. For a down-looking radar, path-averaged estimates of water vapor absorption are possible under rain-free as well as raining conditions by using the surface returns at the three frequencies. Simulations of the water vapor attenuation retrieval show that the largest source of error typically arises from the variance in the measured radar return powers. Although the error can be mitigated by a combination of a high pulse repetition frequency, pulse compression, and averaging in range and time, the radar receiver must be stable over the averaging period. For fractional bandwidths of 20% or less, the potential exists for simultaneous measurements at the three frequencies with a single antenna and transceiver, thereby significantly reducing the cost and mass of the system.

Description: High-resolution infrared radiance spectra obtained from near nadir observations provide atmospheric, surface, and cloud property information. A fast radiative transfer model, including cloud effects, is used for atmospheric profile and cloud parameter retrieval. The retrieval algorithm is presented along with its application to recent field experiment data from the NPOESS Airborne Sounding Testbed - Interferometer (NAST-I). The retrieval accuracy dependence on cloud properties is discussed. It is shown that relatively accurate temperature and moisture retrievals can be achieved below optically thin clouds. For optically thick clouds, accurate temperature and moisture profiles down to cloud top level are obtained. For both optically thin and thick cloud situations, the cloud top height can be retrieved with an accuracy of approximately 1.0 km. Preliminary NAST-I retrieval results from the recent Atlantic-THORPEX Regional Campaign (ATReC) are presented and compared with coincident observations obtained from dropsondes and the nadir-pointing Cloud Physics Lidar (CPL).

Description: This paper presents a new Bayesian algorithm for retrieving surface rain rate from Tropical Rainfall Measurements Mission (TRMM) Microwave Imager (TMI) over the ocean, along with validations against estimates from the TRMM Precipitation Radar (PR). The Bayesian approach offers a rigorous basis for optimally combining multichannel observations with prior knowledge. While other rain rate algorithms have been published that are based at least partly on Bayesian reasoning, this is believed to be the first self-contained algorithm that fully exploits Bayes Theorem to yield not just a single rain rate, but rather a continuous posterior probability distribution of rain rate. To advance our understanding of theoretical benefits of the Bayesian approach, we have conducted sensitivity analyses based on two synthetic datasets for which the true conditional and prior distribution are known. Results demonstrate that even when the prior and conditional likelihoods are specified perfectly, biased retrievals may occur at high rain rates. This bias is not the result of a defect of the Bayesian formalism but rather represents the expected outcome when the physical constraint imposed by the radiometric observations is weak, due to saturation effects. It is also suggested that the choice of the estimators and the prior information are both crucial to the retrieval. In addition, the performance of our Bayesian algorithm is found to be comparable to that of other benchmark algorithms in real-world applications, while having the additional advantage of providing a complete continuous posterior probability distribution of surface rain rate.

Description: The recent Department of Energy Atmospheric Radiation Measurement (ARM) Aerosol Intensive Operations Period (MOP, May 2003) yielded one of the best measurement sets obtained to-date to assess our ability to measure the vertical profile of ambient aerosol extinction sigma(sub ep)(lambda) in the lower troposphere. During one month, a heavily instrumented aircraft with well characterized aerosol sampling ability carrying well proven and new aerosol instrumentation, devoted most of the 60 available flight hours to flying vertical profiles over the heavily instrumented ARM Southern Great Plains (SGP) Climate Research Facility (CRF). This allowed us to compare vertical extinction profiles obtained from 6 different instuments: airborne Sun photometer (AATS-14), airborne nephelometer/absorption photometer, airborne cavity ring-down system, ground-based Raman lidar and 2 ground-based elastic backscatter lidars. We find the in-situ measured sigma(sub ep)(lambda) to be lower than the AATS-14 derived values. Bias differences are 0.002 - 0.004 K/m equivalent to 12-17% in the visible, or 45% in the near-infrared. On the other hand, we find that with respect to AATS-14, the lidar sigma(sub ep)(lambda) are higher. An unnoticed loss of sensitivity of the Raman lidar had occurred leading up to AIOP and we expect better agreement from the recently restored system looking at the collective results from 6 field campaigns conducted since 1996, airborne in situ measurements of sigma(sub ep)(lambda) tend to be biased slightly low (17% at visible wavelengths) when compared to airborne Sun photometer sigma(sub ep)(lambda). On the other hand, sigma(sub ep)(lambda) values derived from lidars tend to have no or positive biases. From the bias differences we conclude that the typical systematic error associated with measuring the tropospheric vertical profile of the ambient aerosol extinction with current state of-the art instrumentation is 15-20% at visible wavelengths and potentially larger in the UV and near-infrared.

Description: How do the continents affect large-scale hydrological cycles? How important can one continent be to the climate system? To address these questions, five years of National Aeronautics and Space Administration (NASA) Terra Moderate Resolution Imaging Spectroradiometer (MODIS) observations, Tropical Rainfall Measuring Mission (TRMM) observations, and the Global Precipitation Climatology Project (GPCP) global precipitation analysis, were used to assess the land impacts on clouds, ~ainfall, and water vapor at continental scales. At these scales, Empirical Orthogonal Function (EOF) and continentally averaged analyses illustrate that continents as integrated regions enhance the seasonality of atmospheric and surface hydrological parameters. Specifically, the continents of Eurasia and North America enhance the seasonality of cloud optical thickness, cirrus fraction, rainfall, and water vapor. Over land, both liquid water and ice cloud effective radii are smaller than over oceans, primarily because land has more aerosol particles. In addition, different continents have similar impacts on hydrological variables in terms of seasonality, but differ in magnitude. For example, in winter, North America and Eurasia increase cloud optical thickness to 17.5 and 16, respectively, while in summer, Eurasia has much smaller cloud optical thicknesses than North An1erica. Such different land impacts are determined by each continent's geographical condition, land cover, and land use. These new understandings help further address the land-ocean contrasts on global climate, help validate global climate model simulated land-atmosphere interactions, and shed light on interpreting climate change over land.

Description: Measurements of BrO suggest that inorganic bromine (Br(sub y)) at and above the tropopause is 4 to 8 ppt greater than assumed in models used in past ozone trend assessment studies. This additional bromine is likely carried to the stratosphere by short-lived biogenic compounds and their decomposition products, including tropospheric BrO. Including this additional bromine in an ozone trend simulation increases the computed ozone depletion over the past approx.25 years, leading to better agreement between measured and modeled ozone trends. This additional Br(sub y) (assumed constant over time) causes more ozone depletion because associated BrO provides a reaction partner for ClO, which increases due to anthropogenic sources. Enhanced Br(sub y) causes photochemical loss of ozone below approx.14 km to change from being controlled by HO(sub x) catalytic cycles (primarily HO2+O3) to a situation where loss by the BrO+HO2 cycle is also important.

Description: Two techniques for retrieving the slope and intercept parameters of an assumed exponential raindrop size distribution (RSD), vertical air velocity, and attenuation by precipitation and water vapor in light stratiform rain using observations by airborne, nadir looking dual-wavelength (X-band, 3.2 cm and W-band, 3.2 mm) radars are presented. In both techniques, the slope parameter of the RSD and the vertical air velocity are retrieved using only the mean Doppler velocities at the two wavelengths. In the first method, the intercept of the RSD is estimated from the observed reflectivity at the longer wavelength assuming no attenuation at that wavelength. The attenuation of the shorter wavelength radiation by precipitation and water vapor are retrieved using the observed reflectivity at the shorter wavelength. In the second technique, it is assumed that the longer wavelength suffers attenuation only in the melting band. Then, assuming a distribution of water vapor, the melting band attenuation at both wavelengths and the rain attenuation at the shorter wavelength are retrieved. Results of the retrievals are discussed and several physically meaningful results are presented.

Description: The bright band is the enhanced radar echo associated with the melting of hydrometeors in stratiform rain where the melting process usually occurs below 0 C isotherm over a distance of about 500m. To simulate this radar signature, a scattering model of melting snow is proposed in which the fractional water content is prescribed as a function of the radius of a spherical mixed- phase particle consisting of air, ice and water. The model is based on the observation that melting starts at the surface of the particle and then gradually develops towards the center. To compute the scattering parameters of a non-uniform melting particle, the particle is modeled as a sphere represented by a collection of 64(exp 3) cubic cells of identical size where the probability of water at any cell is prescribed as a function of the radius. The internal field of the particle, used for deriving the effective dielectric constant, is computed by the Conjugate Gradient and Fast Fourier Transform (CGFFT) numerical methods. To make computations of the scattering parameters more efficient, a multi-layer stratified-sphere scattering model is introduced after demonstrating that the scattering parameters of the non-uniformly melting particle can be accurately reproduced by the stratified sphere. In conjunction with a melting layer model that describes the melting fractions and fall velocities of hydrometeors as a function of the distance from the 0 C isotherm, the stratified-sphere model is used to simulate the radar bright band profiles. These simulated profiles are shown to compare well with measurements from the Precipitation Radar (PR) aboard the Tropical Rainfall Measuring Mission (TRMM) satellite and a dual-wavelength airborne radar. The results suggest that the proposed model of a melting snow particle may be useful in studying the characteristics of the bright-band in particular and mixed- phase hydrometeors in general.

Description: The NASA/NOAA Electronic Theater presents Earth science observations and visualizations from space in a historical perspective. Fly in from outer space to Athens and site of the 2004 Summer Olympics and the Far East using 1 m IKONOS "Spy Satellite" data. Contrast the 1972 Apollo 17 "Blue Marble" image of the Earth with the latest US and International global satellite images that allow us to view our Planet from any vantage point. See the latest spectacular images from NASA/NOAA/Commercial remote sensing missions like Terra, GOES, TRMM, SeaWiFS, & Landsat 7, QuickBird of the SE Asia Tsunami, devastation of Hurricane Katrina this year in New Orleans, and the LA/San Diego Fires of 2003. See how High Definition Television (HDTV) is revolutionizing the way we do science communication. Take the pulse of the planet on a daily, annual and 30-year time scale. See daily thunderstorms, the annual blooming of the northern hemisphere land masses and oceans, fires in Africa, dust storms in Iraq, and carbon monoxide exhaust from global burning. See visualizations featured on Newsweek, TIME, National Geographic, Popular Science covers & National & International Network TV. Spectacular new global visualizations of the observed and simulated atmosphere & oceans are shown. See the currents and vortexes in the oceans that bring up the nutrients to feed tiny plankton and draw the fish, whales and fishermen. See the how the ocean blooms in response to El Nino/La Nina climate changes. The Etheater will be presented using the latest High Definition TV (HDTV) and video projection technology on a large screen. See city lights around the globe and in your area observed by the "night-vision" DMSP satellite, Also see how Keyhole and Google Maps are using satellite and aerial photography to help you find your house and plan your vacation.

Description: By 2025, 60 to 80 percent of the world s population will live in urban environments. Additionally, the following facts published by the United Nations further illustrates how cities will evolve in the future. Urban areas in the developing world are growing very rapidly. The urban growth rate will continue to be particularly rapid in the urban areas of less developed regions, averaging 2.4 per cent per year during 2000-2030, consistent with a doubling time of 29 years. The urbanization process will continue worldwide. The concentration of population in cities is expected to continue so that, by 2030, 84 percent of the inhabitants of more developed countries will be urban dwellers. Urbanization impacts the whole hierarchy of human settlements. In 2000,24.8 per cent of the world population lived in urban settlements with fewer than 500,000 inhabitants and by 2015 that proportion will likely rise to 27.1 per cent.

Description: Attitude measurements using gyros and magnetometers placed on a stratospheric balloon during a non-pointed test flight were used to observe the natural azimuth and elevation motions of a balloon/load train/gondola at an altitude of 36 km over a total flight time of 400 minutes. Time traces of the entire flight are presented. This flight, conducted under nominal atmospheric conditions, had significant motion about the azimuth. Some discussion on balloon disturbances is also included.

Description: Seeing clouds in polar regions has been a problem for the imagers used on satellites. Both clouds and snow and ice are white, which makes clouds over snow hard to see. And for thermal infrared imaging both the surface and the clouds cold. The Geoscience Laser Altimeter System (GLAS) launched in 2003 gives an entirely new way to see clouds from space. Pulses of laser light scatter from clouds giving a signal that is separated in time from the signal from the surface. The scattering from clouds is thus a sensitive and direct measure of the presence and height of clouds. The GLAS instrument orbits over Antarctica 16 times a day. All of the cloud observations for October 2003 were summarized and compared to the results from the MODIS imager for the same month. There are two basic cloud types that are observed, low stratus with tops below 3 km and high cirrus form clouds with cloud top altitude and thickness tending at 12 km and 1.3 km respectively. The average cloud cover varies from over 93 % for ocean and coastal regions to an average of 40% over the East Antarctic plateau and 60-90% over West Antarctica. When the GLAS monthly average cloud fractions are compared to the MODIS cloud fraction data product, differences in the amount of cloud cover are as much as 40% over the continent. The results will be used to improve the way clouds are detected from the imager observations. These measurements give a much improved understanding of distribution of clouds over Antarctica and may show how they are changing as a result of global warming.

Description: CEOP has contributed and will continue to provide integrated data sets from diverse platforms for better understanding of the water and energy cycles, and for validating models. In this talk, I will show examples of how CEOP has contributed to the formulation of a strategy for the study of the monsoon as a system. The CEOP data concept has led to the development of the CEOP Inter-Monsoon Studies (CIMS), which focuses on the identification of model bias, and improvement of model physics such as the diurnal and annual cycles. A multi-model validation project focusing on diurnal variability of the East Asian monsoon, and using CEOP reference site data, as well as CEOP integrated satellite data is now ongoing. Similar validation projects in other monsoon regions are being started. Preliminary studies show that climate models have difficulties in simulating the diurnal signals of total rainfall, rainfall intensity and frequency of occurrence, which have different peak hours, depending on locations. Further more model diurnal cycle of rainfall in monsoon regions tend to lead the observed by about 2-3 hours. These model bias offer insight into lack of, or poor representation of key components of the convective,and stratiform rainfall. The CEOP data also stimulated studies to compare and contrasts monsoon variability in different parts of the world. It was found that seasonal wind reversal, orographic effects, monsoon depressions, meso-scale convective complexes, SST and land surface land influences are common features in all monsoon regions. Strong intraseasonal variability is present in all monsoon regions. While there is a clear demarcation of onset, breaks and withdrawal in the Asian and Australian monsoon region associated with climatological intraseasonal variability, it is less clear in the American and Africa monsoon regions. The examination of satellite and reference site data in monsoon has led to preliminary model experiments to study the impact of aerosol on monsoon variability. I will show examples of how the study of the dynamics of aerosol-water cycle interactions in the monsoon region, can be best achieved using the CEOP data and modeling strategy.

Description: Ozone data from the Improved Limb Atmospheric Spectrometer-II (ILAS-II) were included in addition to other satellite observations in the ozone assimilation system at the Global Modeling and Assimilation Office (GMAO) of NASA/Goddard. The control run assimilated data from NOAA 16 Solar Backscatter Ultraviolet/2 (SBUV/2) and Polar Ozone and Aerosol Measurement III (POAM III) instruments. Persistent impacts over Antarctica and transient impacts over northern middle and high latitudes are seen from April to October 2003, when ILAS-II provided good coverage. The largest improvements with respect to independent ozone sonde data are seen over the South Pole station. Ozone analyses and forecasts from the assimilation of SBUV/2, POAM III and ILAS-II data b e used to investigate the transport of ozone to southern middle latitudes following the breakup of the Antarctic vortex. The quality of analyses and forecasts is evaluated by comparison with independent Stratospheric Aerosol and Gas Experiment III (SAGE III) ozone data near 46degs. Anomaly correlations between SAGE III data and forecasts'exceed 0.6 for up to five to seven days at 30,50, and 70 ma. The loss of skill with advancing forecast length is related to dynamical errors due to an excessively persistent vortex in longer forecasts, which hampers the transport of low ozone air into middle latitudes.

Description: Water is essential to life in the Earth system. The water cycle components that sustain life are becoming more scarce and polluted. The most recent (1999-2004) drought experienced in the southwestern United States is the seventh worst in the approximately 500-year proxy tree-ring record. As a result, many regions contemplated drought emergencies in which severe water restrictions are implemented. Though larger weather and climate processes likely control drought processes, there is increasing evidence that anthropogenic or human-related activities can significantly alter precipitation processes. Urbanization is an example of anthropogenic forcing. Recent studies continue to provide evidence that urban environments can modify or induce precipitation under a specific set of conditions. Arid and semi-arid regions of the southwestern United States and other parts of the world are rapidly developing and placing greater demands on the environmental system. In the past fifty years, Phoenix has expanded from a predominantly agricultural center to an urbanized region with extent 700 percent larger than its size in the middle of the twentieth century. Riyadh's population grew from about a half million people in 1972 to almost two million by 2000. Saudi Arabia experienced urbanization later than many other countries; in the early 1970s its urban-rural ratio was still about 1:3. By 1990 the ratio had reversed to about 3:l. In the mid-1970s Riyadh's population was increasing by about 10 percent a year. Irrigation also significantly increased between 1972 and 1990 southeast of Riyadh. The study employs a 108-year precipitation historical data record, global climate observing network observations and satellite data to identify possible anomalies in rainfall in and around two major arid urban areas, Phoenix, Arizona and Riyadh, Saudi Arabia. It provides statistically sound evidence that rainfall distribution and magnitude is statistically different in post-urban than in pre-urban (1900-1950) Phoenix. The study hypothesis that a complex interaction between the city landscape, irrigated lands, and nearby mountains have created preferred regions for rainfall development. The study also provides early evidence that rapidly urbanizing parts of the arid Middle East may also be experiencing different precipitation regimes in response to urbanization and irrigation.

Description: Monthly mean fields of select geophysical parameters derived from analysis of AIRS/AMSU data, and their interannual differences, are shown and compared with analogous fields derived from other sources. All AIRS fields are derived using the AIRS Science Team Version 4 algorithm. Monthly mean results are shown for January 2004, as are interannual differences between January 2004 and January 2003. AIRS temperature and water vapor profile fields are compared with monthly mean collocated ECMWF 3 hour forecast and monthly mean TOVS Pathfinder Path A data. AIRS Tropospheric and Stratospheric coarse climate indicators are compared with analogous MSU products derived by Spencer and christy and found in the TOVS Pathfinder Path A data set. Total ozone is compared with results produced by TOMS. OLR is compared with OLR derived using CERES data and found in the TOVS Pathfinder Path A data set. AIRS results agree well in all cases, especially in the interannual difference sense.

Description: Vertical wind measurements are presented of an atmospheric thermal in which ospreys (Pandion haliaetus) were soaring. These observations were made with a Doppler lidar, allowing high spatial and high temporal resolution wind profiles in clear air. The thermal was generated at the onset of a cloud bank, producing a rolling eddy upon which ospreys were seen to be riding. A determination is made on the size and shape of the thermal, wind speeds involved, and the altitude to which the birds could have ridden the thermal.

Description: An algorithm is proposed for the computation of streamfunction and velocity potential from given horizontal velocity vectors based on solving a minimization problem. To guarantee the uniqueness of the solution and computational reliability of the algorithm, a Tikhonov regularization is applied. The solution implies that the obtained streamfunction and velocity potential have minimal magnitude, while the given velocity vectors can be accurately reconstructed from the computed streamfunction and velocity potential. Because the formulation of the minimization problem allows for circumventing the explicit specification of separate boundary conditions on the streamfunction and velocity potential, the algorithm is easily applicable to irregular domains. By using an advanced minimization algorithm with the use of adjoint techniques, the method is computationally efficient and suitable for problems with large dimensions. An example is presented for coastal oceans to illustrate the practical application of the algorithm.

Description: A mechanistic model simulation initialized on 14 September 2002, forced by 100-hPa geopotential heights from Met Office analyses, reproduced the dynamical features of the 2002 Antarctic major warming. The vortex split on approx.25 September; recovery after the warming, westward and equatorward tilting vortices, and strong baroclinic zones in temperature associated with a dipole pattern of upward and downward vertical velocities were all captured in the simulation. Model results and analyses show a pattern of strong upward wave propagation throughout the warming, with zonal wind deceleration throughout the stratosphere at high latitudes before the vortex split, continuing in the middle and upper stratosphere and spreading to lower latitudes after the split. Three-dimensional Eliassen-Palm fluxes show the largest upward and poleward wave propagation in the 0(deg)-90(deg)E sector prior to the vortex split (coincident with the location of strongest cyclogenesis at the model's lower boundary), with an additional region of strong upward propagation developing near 180(deg)-270(deg)E. These characteristics are similar to those of Arctic wave-2 major warmings, except that during this warming, the vortex did not split below approx.600 K. The effects of poleward transport and mixing dominate modeled trace gas evolution through most of the mid- to high-latitude stratosphere, with a core region in the lower-stratospheric vortex where enhanced descent dominates and the vortex remains isolated. Strongly tilted vortices led to low-latitude air overlying vortex air, resulting in highly unusual trace gas profiles. Simulations driven with several meteorological datasets reproduced the major warming, but in others, stronger latitudinal gradients at high latitudes at the model boundary resulted in simulations without a complete vortex split in the midstratosphere. Numerous tests indicate very high sensitivity to the boundary fields, especially the wave-2 amplitude. Major warmings occurred for initial fields with stronger winds and larger vortices, but not smaller vortices, consistent with the initiation of wind-deceleration by upward-propagating waves near the poleward edge of the region where wave 2 can propagate above the jet core. Thus, given the observed 100-hPa boundary forcing, stratospheric preconditioning is not needed to reproduce a major warming similar to that observed. The anomalously strong forcing in the lower stratosphere can be viewed as the primary direct cause of the major warming.

Description: The 2003-2004 Arctic winter was remarkable in the approximately 50-year record of meteorological analyses. A major warming beginning in early January 2004 led to nearly 2 months of vortex disruption with high-latitude easterlies in the middle to lower stratosphere. The upper stratospheric vortex broke up in late December, but began to recover by early January, and in February and March was the strongest since regular observations began in 1979. The lower stratospheric vortex broke up in late January. Comparison with 2 previous years, 1984-1985 and 1986-1987, with prolonged midwinter warming periods shows unique characteristics of the 2003-2004 warming period: The length of the vortex disruption, the strong and rapid recovery in the upper stratosphere, and the slow progression of the warming from upper to lower stratosphere. January 2004 zonal mean winds in the middle and lower stratosphere were over 2 standard deviations below average. Examination of past variability shows that the recent frequency of major stratospheric warmings (7 in the past 6 years) is unprecedented. Lower stratospheric temperatures were unusually high during 6 of the past 7 years, with 5 having much lower than usual potential for polar stratospheric cloud (PSC) formation and ozone loss (nearly none in 1998-1999, 2001-2002, and 2003-2004, and very little in 1997-1998 and 2000-2001). Middle and upper stratospheric temperatures, however, were unusually low during and after February. The pattern of 5 of the last 7 years with very low PSC potential would be expected to occur randomly once every 850 years. This cluster of warm winters, immediately following a period of unusually cold winters, may have important implications for possible changes in interannual variability and for determination and attribution of trends in stratospheric temperatures and ozone.

Description: Estimates of ocean circulation constrained by in situ and remotely sensed observations have become routinely available during the past five years, and they are being applied to myriad scientific and operational problems [Stammer et al.,2002]. Under the Global Ocean Data Assimilation Experiment (GODAE), several regional and global estimates have evolved for applications in climate research, seasonal forecasting, naval operations, marine safety, fisheries,the offshore oil industry, coastal management, and other areas. This article reports on recent progress by one effort, the consortium for Estimating the Circulation and Climate of the Ocean (ECCO), toward a next-generation synthesis of ocean and sea-ice data that is global, that covers the full ocean depth, and that permits eddies.

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: The Moderate Resolution Imaging Spectroradiometer (MODIS) was developed by NASA and launched onboard the Terra spacecraft on December 18, 1999 and Aqua spacecraft on May 4, 2002. It achieved its final orbit and began Earth observations on February 24, 2000 for Terra and June 24, 2002 for Aqua. A comprehensive set of remote sensing algorithms for cloud masking and the retrieval of cloud physical and optical properties has been developed by members of the MODIS atmosphere science team. The archived products from these algorithms have applications in climate change studies, climate modeling, numerical weather prediction, as well as fundamental atmospheric research. In addition to an extensive cloud mask, products include cloud-top properties (temperature, pressure, effective emissivity), cloud thermodynamic phase, cloud optical and microphysical parameters (optical thickness, effective particle radius, water path), as well as derived statistics. We will describe the various cloud properties being analyzed on a global basis from both Terra and Aqua. These include the latitudinal distribution of cloud optical and radiative properties of both liquid water and ice clouds, as well as joint histograms of cloud optical thickness and effective radius for selected geographical locations around the world.

Description: The Upper Atmosphere Research Satellite (UARS) was deployed from the Space Shuttle in September 1991 and continued to collect observations to September 2005. The ten instruments aboard the UARS made measurements on atmospheric constituents, atmospheric winds and external forcings. These measurements led to incredible advancements in our understanding of the upper atmosphere. Some of the scientific achievements, which resulted from UARS measurements will be discussed in this review talk.

Description: Chemistry-climate models attempt to bring together our best knowledge of the key processes that govern the composition of the atmosphere and its response to changes in forcing. We test these models on a process by process basis by comparing model results to data from many sources. A more difficult task is testing the model response to changes. One way to do this is to use the natural and anthropogenic experiments that have been done on the atmosphere and are continuing to be done. These include the volcanic eruptions of El Chichon and Pinatubo, the solar cycle, and the injection of chlorine and bromine from CFCs and methyl bromide. The test of the model's response to these experiments is their ability to produce the long-term variations in ozone and the trace gases that affect ozone. We now have more than 25 years of satellite ozone data. We have more than 15 years of satellite and ground-based data of HC1, HN03, and many other gases. I will discuss the testing of models using long-term satellite data sets, long-term measurements from the Network for Detection of Stratospheric Change (NDSC) , long-term ground-based measurements of ozone.

Description: Cloud microphysics are inevitable affected by the smoke particle (CCN, cloud condensation nuclei) size distributions below the clouds, Therefore, size distributions parameterized as spectral bin microphysics are needed to explicitly study the effect of atmospheric aerosol concentration on cloud development, rainfall production, and rainfall rates for convective clouds. Recently, a detailed spectral-bin microphysical scheme was implemented into the the Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bim microphysical processes is based on solving stochastic kinetic equations for the size distribution functions of water droplets (i.e., cloud droplets and raindrops), and several types of ice particles [i.e., pristine ice crystals (columnar and plate-like), snow (dendrites and aggregates), graupel and frozen drops/hail]. Each type is described by a special size distribution function containing many categories (i.e., 33 bins). Atmospheric aerosols are also described using number density size-distribution functions.

Description: We have produced a profile merged ozone data set (MOD) based on the SBUV/SBUV2 series of nadir-viewing satellite backscatter instruments, covering the period from November 1978 - December 2003. In 2004, data from the Nimbus 7 SBUV and NOAA 9, ll, and 16 SBUV/2 instruments were reprocessed using the Version 8 (V8) algorithm and most recent calibrations. More recently, data from the Nimbus 4 BUT instrument, which was operational from 1970 - 1977, were also reprocessed using the V8 algorithm. As part of the V8 profile calibration, the Nimbus 7 and NOAA 9 (1993-1997 only) instrument calibrations have been adjusted to match the NOAA 11 calibration, which was established based on comparisons with SSBUV shuttle flight data. Differences between NOAA 11, Nimbus 7 and NOAA 9 profile zonal means are within plus or minus 5% at all levels when averaged over the respective periods of data overlap. NOAA 16 SBUV/2 data have insufficient overlap with NOAA 11, so its calibration is based on pre-flight information. Mean differences over 4 months of overlap are within plus or minus 7%. Given the level of agreement between the data sets, we simply average the ozone values during periods of instrument overlap to produce the MOD profile data set. Initial comparisons of coincident matches of N4 BUV and Arosa Umkehr data show mean differences of 0.5 (0.5)% at 30km; 7.5 (0.5)% at 35 km; and 11 (0.7)% at 40 km, where the number in parentheses is the standard error of the mean. In this study, we use the MOD profile data set (1978-2003) to estimate the change in profile ozone due to changing stratospheric chlorine levels. We use a standard linear regression model with proxies for the seasonal cycle, solar cycle, QBO, and ozone trend. To account for the non-linearity of stratospheric chlorine levels since the late 1990s, we use a time series of Effective Chlorine, defined as the global average of Chlorine + 50 * Bromine at 1 hPa, as the trend proxy. The Effective Chlorine data are taken from the 3-D Goddard CTM. We will show the latest trend results using this statistical model. In addition, the Nimbus 4 BUV data offer an opportunity to test the physical properties of our statistical model. From ground-based comparisons we will establish an uncertainty range for the Nimbus 4 data. We then extrapolate our statistical model fit backwards in time and compare to the Nimbus 4 data. We compare the characteristics of the residual, defined as the difference between the data and statistical regression fit, during the Nimbus 4 time period and the 1978-2003 period over which the statistical model coefficients were estimated, and present these results.

Description: About 60% of world population reside in Asia, in term of which sheer population density presents a major environmental stress. Economic expansion in this region is, in fact, accompanied by increases in bio-fuel burning, industrial pollution, and land cover and land use changes. With a growth rate of approx. 8%/yr for Indian economy, more than 600 million people from Lahore, Pakistan to Calcutta, India over the Indo-Gangetic Basin have particularly witnessed increased frequencies of floods and droughts as well as a dramatic increase in atmospheric loading of aerosols (i.e., anthropogenic and natural aerosol) in recent decades. This regional change (e.g., aerosol, cloud, precipitation, etc.) will constitute a vital part of the global change in the 21st century. Better understanding of the impacts of aerosols in affecting monsoon climate and water cycles is crucial in providing the physical basis to improve monsoon climate prediction and for disaster mitigation. Based on climate model simulations, absorbing aerosols (dust and black carbon) play a critical role in affecting interannual and intraseasonal variability of the Indian monsoon. An initiative on the integrated (aerosols, clouds, and precipitation) measurements approach over the Indo-Gangetic Basin will be discussed. An array of ground-based (e.g., AERONET, MPLNET, SMART-COMMIT, etc.) and satellite (e.g., Terra, A-Train, etc.) sensors will be utilized to acquire aerosol characteristics, sources/sinks, and transport processes during the pre-monsoon (April-May, aerosol forcing) season, and to obtain cloud and precipitation properties during the monsoon (May-June, water cycle response) season. Close collaboration with other international programs, such as ABC, CLIVAR, GEWEX, and CEOP in the region is anticipated.

Description: The Cirrus Regional Study of Tropical Anvils and Cirrus Layers - Florida Area Cirrus Experiment (CRYSTAL-FACE) was a measurement campaign designed to investigate tropical cirrus cloud physical properties and formation processes. Understanding the production of upper tropospheric cirrus clouds is essential for the successful modeling of the Earth's climate. The mission was staged in July 2002 with flights of 6 aircraft from Key West, Florida. Several aircraft were used for in situ and remote sensing of aerosols, ice crystals, meteorological fields, radiative fluxes, and gas concentrations. The NASA ER-2 and WB-57, the Proteus aircraft, owned by Northrop Grumman and operated by Scaled Composites, CIRPAS provided the DeHavilland UV-l8A, "Twin Otter" aircraft, the University of North Dakota provided a Cessna Citation II aircraft , and NSF supported the ELDORA radar onboard the Naval Research Laboratory P-3 aircraft. In this presentation, I will describe some of the flights, the conditions, and some of the results from the mission.

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. The technique makes use of the IR data from the TRMM Visible/Infrared Scanner (VIRS) before application to global geosynchronous satellite data. The calibrated CST technique has the advantages of high spatial resolution (4 km), filtering of nonraining cirrus clouds, and the stratification of the rainfall into its convective and stratiform components, the last being important for the calculation of vertical profiles of latent heating. 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. Results from five years of PR data will show the global-tropical partitioning of convective and stratiform rainfall.

Description: The damage surveys conducted by the NWS in the aftermath of a reported tornadic event are used to document the location of the tornado ground damage track (path length and width) and an estimation of the tornado intensity. This study explored the possibility of using near real-time medium and high-resolution satellite imagery from the NASA EOS satellites to provide additional information for the surveys. MODIS and ASTER data were used to study the damage tracks from three tornadic storms; the La Plata, Maryland storm of 28 April 2002 and the Carter-Butler Counties and Madison County Missouri storms of 24 April 2002. These storms varied in intensity (from F0-F4) and occurred over regions with different land use. It was found that, depending on the nature of land use, tornado damage tracks from intense storms (F2 or greater) may be evident in both ASTER and MODIS satellite imagery. In areas of dense vegetation the scar patterns show up very clearly, while in areas of grassland and regions with few trees, scar patterns are not at all obvious in the satellite imagery. The detection of previously unidentified segments of a damage track caused by the 24 April 2004 Madison County, Missouri tornado demonstrates the utility of satellite imagery for damage surveys. However, the capability to detect tornado tracks in satellite imagery appears to be as much dependent on the nature of the underlying surface and land use as on the severity of the tornadic storm. The imaging sensors on the NPOESS operational satellites to be launched in 2006 will continue the unique observing capabilities of the EOS instruments.

Description: The North Alabama Lightning Mapping Array (LMA) is s a 3-D VHF regional lightning detection system that provides on-orbit algorithm validation and instrument performance assessments for the NASA Lightning Imaging Sensor, as well as information on storm kinematics and updraft evolution that offers the potential to improve severe storm warning lead time by up t o 50% and decrease te false alarm r a t e ( for non-tornado producing storms). In support of this latter function, the LMA serves as a principal component of a severe weather test bed to infuse new science and technology into the short-term forecasting of severe and hazardous weather, principally within nearby National Weather Service forecast offices. The LMA, which became operational i n November 2001, consists of VHF receivers deployed across northern Alabama and a base station located at the National Space Science and Technology Center (NSSTC), which is on t h e campus of the University of Alabama in Huntsville. The LMA system locates the sources of impulsive VHF radio signals s from lightning by accurately measuring the time that the signals aririve at the different receiving stations. Each station's records the magnitude and time of the peak lightning radiation signal in successive 80 ms intervals within a local unused television channel (channel 5, 76-82 MHz in our case ) . Typically hundreds of sources per flash can be reconstructed, which i n t u r n produces accurate 3-dimensional lightning image maps (nominally 〈50 m error within 150 la. range). The data are transmitted back t o a base station using 2.4 GHz wireless Ethernet data links and directional parabolic grid antennas. There are four repeaters in the network topology and the links have an effective data throughput rate ranging from 600 kbits s -1 t o 1.5 %its s -1. This presentation provides an overview of t h e North Alabama network, the data processing (both real-time and post processing) and network statistics.

Description: We ask how lightning measurements might best be used to improve short-term (0-24 hr) weather forecasting. We examine recently developed strategies for the integration of lightning data into Short-term forecasts (nowcasts) of convective and severe weather hazards and the assimilation of lightning data into numerical weather prediction models. In each strategy we define specific metrics of forecast improvement and a progress assessment. We also address the conventional observing system deficiencies and potential gap-filling information that can be addressed through the use of the lightning measurement.

Description: Activation of some fraction of increased concentrations of sub-micron soluble aerosol particles lead to enhanced cloud droplet concentrations and hence smaller droplets, increasing their total cross sectional area and thus reflecting solar radiation more efficiently (the Twomey, or first indirect, effect). However, because of competition during condensational growth, droplet distributions tend to broaden as numbers increase, reducing the sensitivity of cloud albedo to droplet concentration on the order of 10%. Also, smaller droplets less effectively produce drizzle through collisions and coalescence, and it is widely expected (and found in large-scale models) that decreased precipitation leads to clouds with more cloud water on average (the so-called cloud lifetime, or second indirect, effect). Much of the uncertainty regarding the overall indirect aerosol effect stems from inadequate understanding of such changes in cloud water. Detailed simulations based on FIRE-I, ASTEX, and DYCOMS-II conditions show that suppression of precipitation from increased droplet concentrations leads to increased cloud water only when sufficient precipitation reaches the surface, a condition favored when the overlying air is-humid or droplet concentrations are very low. Otherwise, aerosol induced suppression of precipitation enhances entrainment of overlying dry air, thereby reducing cloud water and diminishing the indirect climate forcing.

Description: Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud- resolving models (CRMs) agree with observations better than traditional single-column models in simulating various types of clouds and cloud systems from different geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a super-parameterization or multi-scale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameterization. NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. A seed fund is available at NASA Goddard to build a MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM). A prototype MMF in being developed and production runs will be conducted at the beginning of 2005. In this talk, I will present: (1) A brief review on GCE model and its applications on precipitation processes, ( 2 ) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), (3) A cloud library generated by Goddard MMF, and 3D GCE model, and (4) A brief discussion on the GCE model on developing a global cloud simulator.

Description: We report on a NASA supported modeling effort to reduce uncertainty in carbon cycle processes that create the so-called missing sink of atmospheric CO2. Our overall objective is to improve characterization of CO2 source/sink processes globally with improved formulations for atmospheric transport, terrestrial uptake and release, biomass and fossil fuel burning, and observational data analysis. The motivation for this study follows from the perspective that progress in determining CO2 sources and sinks beyond the current state of the art will rely on utilization of more extensive and intensive CO2 and related observations including those from satellite remote sensing. The major components of this effort are: 1) Continued development of the chemistry and transport model using analyzed meteorological fields from the Goddard Global Modeling and Assimilation Office, with comparison to real time data in both forward and inverse modes; 2) An advanced biosphere model, constrained by remote sensing data, coupled to the global transport model to produce distributions of CO2 fluxes and concentrations that are consistent with actual meteorological variability; 3) Improved remote sensing estimates for biomass burning emission fluxes to better characterize interannual variability in the atmospheric CO2 budget and to better constrain the land use change source; 4) Evaluating the impact of temporally resolved fossil fuel emission distributions on atmospheric CO2 gradients and variability. 5) Testing the impact of existing and planned remote sensing data sources (e.g., AIRS, MODIS, OCO) on inference of CO2 sources and sinks, and use the model to help establish measurement requirements for future remote sensing instruments. The results will help to prepare for the use of OCO and other satellite data in a multi-disciplinary carbon data assimilation system for analysis and prediction of carbon cycle changes and carbodclimate interactions.

Description: The geomagnetic storms of April 17-21,2002 and May 29-30,2003 caused large decreases in the O/N2 column density ratio in the thermosphere. For these storms, ON2 column density decreases of greater than 50% were observed to extend to mid-to-low latitudes with the FUV sensitive Earth Camera of the Visible Imaging System (VIS) on the Polar spacecraft. Simultaneously in these same regions, the ground-based GPS network observed approximately 80% reductions in the Total Electron Content (TEC) of the ionosphere. The reduction in the Om2 column density ratio is due mainly to increases in the molecular species that have welled-up into the thermosphere from the lower levels of the atmosphere due to auroral heating. The geomagnetic-storm driven increase in molecular densities at typical ionospheric heights rapidly charge exchange with the ambient ionized atoms and subsequently dissociatively recombine with the ionospheric electrons leading to a reduction in the total charge density. The transition boundaries between high and low regions of O/N2 as well as TEC can be tracked in the images and the thermospheric winds can be determined from the motion of the boundaries. The motion of these boundaries during the development of the geomagnetic storm will be discussed.

Description: Ozone data from Aura Microwave Limb Sounder (MLS) and Ozone Monitoring Instrument (OMI) were assimilated into the ozone model at NASA's Global Modeling and Assimilation Office (GMAO). This assimilation produces ozone fields that are superior to those from the operational GMAO assimilation of Solar Backscatter Ultraviolet (SBUV/2) instrument data. Assimilation of Aura data improves the representation of the "ozone hole" and the agreement with independent Stratospheric Aerosol and Gas Experiment (SAGE) III and ozone sonde data. Ozone in the lower stratosphere is captured better: mean state, vertical gradients, spatial and temporal variability are all improved. Inclusion of OMI and MLS data together, or separately, in the assimilation system provides a way of checking how consistent OMI and MLS data are with each other, and with the ozone model. We found that differences between OMI total ozone column data and model forecasts decrease after MLS data are assimilated. This indicates that MLS stratospheric ozone profiles are consistent with OMI total ozone columns. The evaluation of error characteristics of OMI and MLS ozone will continue as data from newer versions of retrievals becomes available. We report on the initial step in obtaining global assimilated ozone fields that combine measurements from different Aura instruments, the ozone model at the GMAO, and their respective error characteristics. We plan to use assimilated ozone fields in estimation of tropospheric ozone. We also plan to investigate impacts of assimilated ozone fields on numerical weather prediction through their use in radiative models and in the assimilation of infrared nadir radiance data from NASA's Advanced Infrared Sounder (AIRS).

Description: Recently Grabowski (2001) and Khairoutdinov and Randall (2001) have proposed the use of 2D CFWs as a "super parameterization" [or multi-scale modeling framework (MMF)] to represent cloud processes within atmospheric general circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization used in conventional GCMs. A prototype Goddard MMF based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM) is now being developed. The prototype includes the fvGCM run at 2.50 x 20 horizontal resolution with 32 vertical layers from the surface to 1 mb and the 2D (x-z) GCE using 64 horizontal and 32 vertical grid points with 4 km horizontal resolution and a cyclic lateral boundary. The time step for the 2D GCE would be 15 seconds, and the fvGCM-GCE coupling frequency would be 30 minutes (i.e. the fvGCM physical time step). We have successfully developed an fvGCM-GCE coupler for this prototype. Because the vertical coordinate of the fvGCM (a terrain-following floating Lagrangian coordinate) is different from that of the GCE (a z coordinate), vertical interpolations between the two coordinates are needed in the coupler. In interpolating fields from the GCE to fvGCM, we use an existing fvGCM finite- volume piecewise parabolic mapping (PPM) algorithm, which conserves the mass, momentum, and total energy. A new finite-volume PPM algorithm, which conserves the mass, momentum and moist static energy in the z coordinate, is being developed for interpolating fields from the fvGCM to the GCE. In the meeting, we will discuss the major differences between the two MMFs (i.e., the CSU MMF and the Goddard MMF). We will also present performance and critical issues related to the MMFs. In addition, we will present multi-dimensional cloud datasets (i.e., a cloud data library) generated by the Goddard MMF that will be provided to the global modeling community to help improve the representation and performance of moist processes in climate models and to improve our understanding of cloud processes globally (the software tools needed to produce cloud statistics and to identify various types of clouds and cloud systems from both high-resolution satellite and model data will be also presented).

Description: CEOP has contributed and will continue to provide integrated data sets from diverse platforms for better understanding of the water and energy cycles, and for validaintg models. In this talk, I will show examples of how CEOP has contributed to the formulation of a strategy for the study of the monsoon as a system. The CEOP data concept has led to the development of the CEOP Inter-Monsoon Studies (CIMS), which focuses on the identification of model bias, and improvement of model physics such as the diurnal and annual cycles. A multi-model validation project focusing on diurnal variability of the East Asian monsoon, and using CEOP reference site data, as well as CEOP integrated satellite data is now ongoing. Preliminary studies show that climate models have difficulties in simulating the diurnal signals of total rainfall, rainfall intensity and frequency of occurrence, which have different peak hours, depending on locations. Further more model diurnal cycle of rainfall in monsoon regions tend to lead the observed by about 2-3 hours. These model bias offer insight into lack of, or poor representation of, key components of the convective and stratiform rainfall. The CEOP data also stimulated studies to compare and contrasts monsoon variability in different parts of the world. It was found that seasonal wind reversal, orographic effects, monsoon depressions, meso-scale convective complexes, SST and land surface land influences are common features in all monsoon regions. Strong intraseasonal variability is present in all monsoon regions. While there is a clear demarcation of onset, breaks and withdrawal in the Asian and Australian monsoon region associated with climatological intraseasonal variabillity, it is less clear in the American and Africa monsoon regions. The examination of satellite and reference site data in monsoon has led to preliminary model experiments to study the impact of aerosol on monsoon variability. I will show examples of how the study of the dynamics of aerosol-water cycle interactions in the monsoon region, can be best achieved using the CEOP data and modeling strategy.

Description: Trace gases measurements are used to diagnose both the chemistry and transport of the atmosphere. These lectures emphasize the interpretation of trace gases measurements and techniques used to untangle chemistry and transport effects. I will discuss PV transform, trajectory techniques, and age-of-air as far as the circulation of the stratosphere.

Description: The ODASI consortium is focused activity of the NOAA/OGP/Climate Diagnostics and Experimental Prediction Program with the goal of improving ocean data assimilation methods and their implementations in support of seasonal forecasts with coupled general circulation models. The consortium is undertaking coordinated assimilation experiments, with common forcing data sets and common input data streams. With different assimilation systems and different models, we aim to understand what approach works best in improving forecast skill in the equatorial Pacific. The presentation will provide an overview of the consortium goals and plans and recent results focused towards evaluating data impacts.

Description: The Geoscience Laser Altimeter System (GLAS) launched in 2003 is the first polar orbiting satellite lidar. The instrument was designed for high performance observations of the distribution and optical scattering cross sections of clouds and aerosol. The backscatter lidar operates at two wavelengths, 532 and 1064 nm. Both receiver channels meet and exceed their design goals, and beginning with a two month period through October and November 2003, an excellent global lidar data set now exists. The data products for atmospheric observations include the calibrated, attenuated backscatter cross section for cloud and aerosol; height detection for multiple cloud layers; planetary boundary layer height; cirrus and aerosol optical depth and the height distribution of aerosol and cloud scattering cross section profiles. The data sets are now in open release through the NASA data distribution system. The initial results on global statistics for cloud and aerosol distribution has been produced and in some cases compared to other satellite observations. The sensitivity of the cloud measurements is such that the 70% global cloud coverage result should be the most accurate to date. Results on the global distribution of aerosol are the first that produce the true height distribution for model inter-comparison.

Description: Clouds have a powerful influence on atmospheric radiative transfer and hence are crucial to understanding and interpreting the exchange of radiation between the Earth's surface, the atmosphere, and space. Because clouds are highly variable in space, time and physical makeup, it is important to be able to observe them in three dimensions (3-D) with sufficient resolution that the data can be used to generate and validate parameterizations of cloud fields at the resolution scale of global climate models (GCMs). Simulation of photon transport in three dimensionally inhomogeneous cloud fields show that spatial inhomogeneities tend to decrease cloud reflection and absorption and increase direct and diffuse transmission, Therefore it is an important task to characterize cloud spatial structures in three dimensions on the scale of GCM grid elements. In order to validate cloud parameterizations that represent the ensemble, or mean and variance of cloud properties within a GCM grid element, measurements of the parameters must be obtained on a much finer scale so that the statistics on those measurements are truly representative. High spatial sampling resolution is required, on the order of 1 km or less. Since the radiation fields respond almost instantaneously to changes in the cloud field, and clouds changes occur on scales of seconds and less when viewed on scales of approximately 100m, the temporal resolution of cloud properties should be measured and characterized on second time scales. GCM time steps are typically on the order of an hour, but in order to obtain sufficient statistical representations of cloud properties in the parameterizations that are used as model inputs, averaged values of cloud properties should be calculated on time scales on the order of 10-100 s. The Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE) provides exceptional temporal (100 ms) and spatial (30 m) resolution measurements of aerosol and cloud backscatter in three dimensions. HARLIE was used in a ground-based configuration in several recent field campaigns. Principal data products include aerosol backscatter profiles, boundary layer heights, entrainment zone thickness, cloud fraction as a function of altitude and horizontal wind vector profiles based on correlating the motions of clouds and aerosol structures across portions of the scan. Comparisons will be made between various cloud detecting instruments to develop a baseline performance metric.

Description: Numerical cloud models, which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Because cloud-scale dynamics are treated explicitly, uncertainties stemming from convection that have to be parameterized in (hydrostatic) large-scale models are obviated, or at least mitigated, in cloud models. Global models will use the non-hydrostatic framework when their horizontal resolution becomes about 10 kilometers, the theoretical limit for the hydrostatic approximation. This juncture will be reached one to two decades from now. Over the past generation, voluminous datasets on atmospheric convection have been accumulated from radar, instrumented aircraft, satellites, and rawinsonde measurements in field campaigns, enabling the detailed evaluation of models. Improved numerical methods have resulted in more accurate and efficient dynamical cores in models. Improvements have been made in the parameterizations of microphysical processes, radiation, boundary-layer effects, and turbulence; however, microphysical parameterizations remain a major source of uncertainty in all classes of atmospheric models. In recent years, exponentially increasing computer power has extended cloud-resolving-model integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-kilometer scales are resolved in horizontal domains as large as 10,000 kilometers in two dimensions, and 1,000 x 1,000 square kilometers in three-dimensions. Cloud models now provide statistical information useful for developing more realistic physically-based parameterizations for climate models and numerical weather prediction models. A review of developments and applications of cloud models in the past, present and future will be presented in this talk. In particular, a new approach to using cloud-resolving models to study the interactions between clouds, precipitation and aerosols will be presented.

Description: Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud-resolving models (CRMs) agree with observations better than traditional singlecolumn models in simulating various types of clouds and cloud systems from Merent geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloudscale model (termed a super-parameterization or multiscale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and more sophisticated physical parameteridon NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. A seed fund is available at NASA Goddard to build a MMF based on the 2D Goddard cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM). A prototype MMF in being developed and production nms will be conducted at the beginning of 2005. In this talk, I will present: (1) A brief review on GCE model and its applications on precipitation processes, (2) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), (3) A cloud library generated by Goddard MMF, and 3D GCE model, and (4) A brief discussion on the GCE model on developing a global cloud simulator.

Description: The 45th Weather Squadron (45 WS) at Cape Canaveral Air-Force Station (CCAFS)ln Florida issues a probability of lightning occurrence in their daily 24-hour and weekly planning forecasts. This information is used for general planning of operations at CCAFS and Kennedy Space Center (KSC). These facilities are located in east-central Florida at the east end of a corridor known as 'Lightning Alley', an indication that lightning has a large impact on space-lift operations. Much of the current lightning probability forecast is based on a subjective analysis of model and observational data and an objective forecast tool developed over 30 years ago. The 45 WS requested that a new lightning probability forecast tool based on statistical analysis of more recent historical warm season (May-September) data be developed in order to increase the objectivity of the daily thunderstorm probability forecast. The resulting tool is a set of statistical lightning forecast equations, one for each month of the warm season, that provide a lightning occurrence probability for the day by 1100 UTC (0700 EDT) during the warm season.

Description: Each morning, the forecasters at the National Weather Service in Melbourn, FL (NWS MLB) produce an experimental cloud-to-ground (CG) lightning threat index map for their county warning area (CWA) that is posted to their web site (http://www.srh.weather.gov/mlb/ghwo/lightning.shtml) . Given the hazardous nature of lightning in central Florida, especially during the warm season months of May-September, these maps help users factor the threat of lightning, relative to their location, into their daily plans. The maps are color-coded in five levels from Very Low to Extreme, with threat level definitions based on the probability of lightning occurrence and the expected amount of CG activity. On a day in which thunderstorms are expected, there are typically two or more threat levels depicted spatially across the CWA. The locations of relative lightning threat maxima and minima often depend on the position and orientation of the low-level ridge axis, forecast propagation and interaction of sea/lake/outflow boundaries, expected evolution of moisture and stability fields, and other factors that can influence the spatial distribution of thunderstorms over the CWA. The lightning threat index maps are issued for the 24-hour period beginning at 1200 UTC (0700 AM EST) each day with a grid resolution of 5 km x 5 km. Product preparation is performed on the AWIPS Graphical Forecast Editor (GFE), which is the standard NWS platform for graphical editing. Currently, the forecasters create each map manually, starting with a blank map. To improve efficiency of the forecast process, NWS MLB requested that the Applied Meteorology Unit (AMU) create gridded warm season lightning climatologies that could be used as first-guess inputs to initialize lightning threat index maps. The gridded values requested included CG strike densities and frequency of occurrence stratified by synoptic-scale flow regime. The intent is to increase consistency between forecasters while enabling them to focus on the mesoscale detail of the forecast, ultimately benefiting the end-users of the product. Several studies took place at the Florida State University (FSU) and NWS Tallahassee (TAE) for which they created daily flow regimes using Florida 1200 UTC synoptic soundings and CG strike densities from National Lightning Detection Network (NLDN) data. The densities were created on a 2.5 km x 2.5 km grid for every hour of every day during the warm seasons in the years 1989-2004. The grids encompass an area that includes the entire state of Florida and adjacent Atlantic and Gulf of Mexico waters. Personnel at the two organizations provided this data and supporting software for the work performed by the AMU. The densities were first stratified by flow regime, then by time in 1-, 3-, 6-, 12-, and 24-hour increments while maintaining the 2.5 km x 2.5 km grid resolution. A CG frequency of occurrence was calculated for each stratification and grid box by counting the number of days with lightning and dividing by the total number of days in the data set. New CG strike densities were calculated for each stratification and grid box by summing the strike number values over all warm seasons, then normalized by dividing the summed values by the number of lightning days. This makes the densities conditional on whether lightning occurred. The frequency climatology values will be used by forecasters as proxy inputs for lightning prObability, while the density climatology values will be used for CG amount. In addition to the benefits outlined above, these climatologies will provide improved temporal and spatial resolution, expansion of the lightning threat area to include adjacent coastal waters, and potential to extend the forecast to include the day-2 period. This presentation will describe the lightning threat index map, discuss the work done to create the maps initialized with climatological guidance, and show examples of the climatological CG lightning densities and frequencies of occurren based on flow regime.

Description: This report summarizes the Applied Meteorology Unit (AMU) activities for the first quarter of Fiscal Year 2005 (October - December 2005). Tasks reviewed include: (1) Objective Lightning Probability Forecast: Phase I, (2) Severe Weather Forecast Decision Aid, (3) Hail Index, (4) Stable Low Cloud Evaluation, (5) Shuttle Ascent Camera Cloud Obstruction Forecast, (6) Range Standardization and Automation (RSA) and Legacy Wind Sensor Evaluation, (7) Advanced Regional Prediction System (ARPS) Optimization and Training Extension, and (8) User Control Interface for ARPS Data Analysis System (ADAS) Data Ingest

Description: Weather constraints to launching space vehicles are designed to prevent loss of the vehicle or mission due to weather hazards (See, e.g., Ref 1). Constraints include Lightning Launch Commit Criteria (LLCC) designed to avoid natural and triggered lightning. The LLCC currently in use at most American launch sites including the Eastern Range and Kennedy Space Center require the Launch Weather Officer to determine the height of cloud bases and tops, the location of cloud edges, and cloud transparency. The preferred method of making these determinations is visual observation, but when that isn't possible due to darkness or obscured vision, it is permissible to use radar. This note examines the relationship between visual and radar observations in three ways: A theoretical consideration of the relationship between radar reflectivity and optical transparency. An observational study relating radar reflectivity to cloud edge determined from in-situ measurements of cloud particle concentrations that determine the visible cloud edge. An observational study relating standard radar products to anvil cloud transparency. It is shown that these three approaches yield results consistent with each other and with the radar threshold specified in Reference 2 for LLCC evaluation.

Description: During August and September 2004, four hurricanes tested the mettle of Space Coast residents and the Kennedy Space Center (KSC) leadership and workforce. These threats underscored two important points: the very real vulnerability of KSC and its valuable space program assets to the devastating power of a hurricane, and the planning required to effectively deal with such threats. The damage was significant even though KSC did not experience sustained hurricane-force winds. To better understand and appreciate these points, this report provides an overview of the meteorological history of the Space Coast and what is involved in the planning, preparation, and recovery activities, as well as addressing the impacts of the 2004 hurricane season.

Description: The contour error map (CEM) algorithm and the software that implements the algorithm are means of quantifying correlations between sets of time-varying data that are binarized and registered on spatial grids. The present version of the software is intended for use in evaluating numerical weather forecasts against observational sea-breeze data. In cases in which observational data come from off-grid stations, it is necessary to preprocess the observational data to transform them into gridded data. First, the wind direction is gridded and binarized so that D(i,j;n) is the input to CEM based on forecast data and d(i,j;n) is the input to CEM based on gridded observational data. Here, i and j are spatial indices representing 1.25-km intervals along the west-to-east and south-to-north directions, respectively; and n is a time index representing 5-minute intervals. A binary value of D or d = 0 corresponds to an offshore wind, whereas a value of D or d = 1 corresponds to an onshore wind. CEM includes two notable subalgorithms: One identifies and verifies sea-breeze boundaries; the other, which can be invoked optionally, performs an image-erosion function for the purpose of attempting to eliminate river-breeze contributions in the wind fields.

Description: Presentation by Cristina Milesi, First Author, NASA Ames Research Center, Moffett Field, CA at the "Meeting the Challenge of Sea Level Rise in Santa Clara County" on June 19, 2005 Santa Clara County, bordering with the southern portion of the San Francisco Bay, is highly vulnerable to flooding and to sea level rise (SLR). In this presentation, the latest sea level rise projections for the San Francisco Bay will be discussed in the context of extreme water height frequency and extent of flooding vulnerability. I will also present preliminary estimations of levee requirements and possible mitigation through tidal restoration of existing salt ponds. The examples will draw mainly from the work done by the NASA Climate Adaptation Science Investigators at NASA Ames.

Description: This study analyzed daily satellite cloud observations and reanalysis dynamical parameters to determine how mid-tropospheric vertical velocity and advection over the sea surface temperature gradient control midlatitude North Pacific cloud properties. Optically thick clouds with high tops are generated by synoptic ascent, but two different cloud regimes occur under synoptic descent. When vertical motion is downward during summer, extensive stratocumulus cloudiness is associated with near surface northerly wind, while frequent cloudless pixels occur with southerly wind. Examinations of ship-reported cloud types indicates that midlatitude stratocumulus breaks up as the the boundary level decouples when it is advected equatorward over warmer water. Cumulus is prevalent under conditions of synoptic descent and cold advection during winter. Poleward advection of subtropical air over colder water causes stratification of the near-surface layer that inhibits upward mixing of moisture and suppresses cloudiness until a fog eventually forms. Averaging of cloud and radiation data into intervals of 500-hPa vertical velocity and advection over the SST gradient enables the cloud response to changes in temperature and the stratification of the lower troposphere to be investigated independent of the dynamics.

Description: This study investigated the spatial pattern of linear trends in surface-observed upper-level (combined mid-level and High-level) cloud cover, precipitation, and surface divergence over the tropical Indo-Pacific Ocean during 1952-1957. Cloud values were obtained from the Extended Edited Cloud Report Archive (EECRA), precipitation values were obtained from the Hulme/Climate Research Unit Data Set, and surface divergence was alternatively calculated from wind reported Comprehensive Ocean-Atmosphere Data Set and from Smith and Reynolds Extended Reconstructed sea level pressure data.

Description: The first paper was Multidecadal changes in near-global cloud cover and estimated cloud cover radiative forcing, by J. R. Norris (2005, J. Geophys. Res. - Atmos., 110, D08206, doi: lO.l029/2004JD005600). This study examined variability in zonal mean surface-observed upper-level (combined midlevel and high-level) and low-level cloud cover over land during 1971-1 996 and over ocean during 1952-1997. These data were averaged from individual synoptic reports in the Extended Edited Cloud Report Archive (EECRA). Although substantial interdecadal variability is present in the time series, long-term decreases in upper-level cloud cover occur over land and ocean at low and middle latitudes in both hemispheres. Near-global upper-level cloud cover declined by 1.5%-sky-cover over land between 1971 and 1996 and by 1.3%-sky-cover over ocean between 1952 and 1997. Consistency between EECRA upper-level cloud cover anomalies and those from the International Satellite Cloud Climatology Project (ISCCP) during 1984-1 997 suggests the surface-observed trends are real. The reduction in surface-observed upper-level cloud cover between the 1980s and 1990s is also consistent with the decadal increase in all-sky outgoing longwave radiation reported by the Earth Radiation Budget Satellite (EMS). Discrepancies occur between time series of EECRA and ISCCP low-level cloud cover due to identified and probable artifacts in satellite and surface cloud data. Radiative effects of surface-observed cloud cover anomalies, called "cloud cover radiative forcing (CCRF) anomalies," are estimated based on a linear relationship to climatological cloud radiative forcing per unit cloud cover. Zonal mean estimated longwave CCRF has decreased over most of the globe. Estimated shortwave CCRF has become slightly stronger over northern midlatitude oceans and slightly weaker over northern midlatitude land areas. A long-term decline in the magnitude of estimated shortwave CCRF occurs over low-latitude land and ocean, but comparison with EMS all-sky reflected shortwave radiation during 1985-1997 suggests this decrease may be underestimated.

Description: Numerical cloud models, which are based the non-hydrostatic equations of motion, have been extensively applied to cloud-scale and mesoscale processes during the past four decades. Because cloud-scale dynamics are treated explicitly, uncertainties stemming from convection that have to be parameterized in (hydrostatic) large-scale models are obviated, or at least mitigated, in cloud models. Global models will use the non-hydrostatic framework when their horizontal resolution becomes about 10 km, the theoretical limit for the hydrostatic approximation. This juncture will be reached one to two decades from now. Over the past generation, voluminous datasets on atmospheric convection have been accumulated from radar, instrumented aircraft, satellites, and rawinsonde measurements in field campaigns, enabling the detailed evaluation of models. Improved numerical methods have resulted in more accurate and efficient dynamical cores in models. Improvements have been made in the parameterizations of microphysical processes, radiation, boundary-layer effects, and turbulence; however, microphysical parameterizations remain a major source of uncertainty in all classes of atmospheric models. In recent years, exponentially increasing computer power has extended cloud-resolving-model integrations from hours to months, the number of computational grid points from less than a thousand to close to ten million. Three-dimensional models are now more prevalent. Much attention is devoted to precipitating cloud systems where the crucial 1-km scales are resolved in horizontal domains as large as 10,000 km in two-dimensions, and 1,000 x 1,000 square kilometers in three-dimensions. Cloud models now provide statistical information useful for developing more realistic physically-based parameterizations for climate models and numerical weather prediction models. A review of developments and applications of cloud models in the past, present and future will be presented in this talk. In particular, a new approach to using cloud-resolving models to study the interactions between clouds, precipitation and aerosols will be presented.

Description: Almost 20 years ago Olsen et al. [1987] reported Dynamics Explorer 1 Retarding Ion Mass Spectrometer observations of equatorially trapped, cold ions in the vicinity of the plasmapause. In that study the trapped population corresponded to a local minimum in density at the magnetic equator. During that time period observations were uncovered of a local maximum in plasma density at the equator. Until IMAGE there has been no good opportunity to experimentally revisit this topic, however until now no direct evidence of a relevant equatorial process has been recognized near the plasmapause during early recovery conditions. It appears that evidence has now been found in both the Extreme Ultraviolet Imager and Radio Plasma Imager observations. The observations, conditions, and properties of what appears to be an equatorially trapped and enhanced density near the magnetic equator will be presented and discussed.

Description: Accurate estimates of snow water equivalent and other properties play an important role in weather, natural hazard, and hydrological forecasting and climate modeling over a range of scales in space and time. Remote sensing-derived estimates have traditionally been of the "snapshot" type, but techniques involving models with assimilation are also being explored. In both cases, forward emission models are useful to understand the observed passive microwave signatures and developing retrieval algorithms. However, mismatches between passive microwave sensor resolutions and the scales of processes controlling subpixel heterogeneity can affect the accuracy of the estimates. Improving the spatial resolution of new passive microwave satellite sensors is a major desire in order to (literally) resolve such subpixel heterogeneity, but limited spacecraft and mission resources impose severe constraints and tradeoffs. In order to maximize science return while mitigating risk for a satellite concept, it is essential to understand the scaling behavior of snow in terms of what the sensor sees (brightness temperature) as well as in terms of the actual variability of snow. NASA's Cold Land Processes Experiment-1 (CLPX-1: Colorado, 2002 and 2003) was designed to provide data to measure these scaling behaviors for varying snow conditions in areas with forested, alpine, and meadow/pasture land cover. We will use observations from CLPX-1 ground, airborne, and satellite passive microwave sensors to examine and evaluate the scaling behavior of observed and modeled brightness temperatures and observed and retrieved snow parameters across scales from meters to 10's of kilometers. The conclusions will provide direct examples of the appropriate spatial sampling scales of new sensors for snow remote sensing. The analyses will also illustrate the effects and spatial scales of the underlying phenomena (e.g., land cover) that control subpixel heterogeneity.

Description: The NASA Columbia supercomputer was ranked second on the TOP500 List in November, 2004. Such a quantum jump in computing power provides unprecedented opportunities to conduct ultra-high resolution simulations with the finite-volume General Circulation Model (fvGCM). During 2004, the model was run in realtime experimentally at 0.25 degree resolution producing remarkable hurricane forecasts [Atlas et al., 2005]. In 2005, the horizontal resolution was further doubled, which makes the fvGCM comparable to the first mesoscale resolving General Circulation Model at the Earth Simulator Center [Ohfuchi et al., 2004]. Nine 5-day 0.125 degree simulations of three hurricanes in 2004 are presented first for model validation. Then it is shown how the model can simulate the formation of the Catalina eddies and Hawaiian lee vortices, which are generated by the interaction of the synoptic-scale flow with surface forcing, and have never been reproduced in a GCM before.)

Description: Considerable uncertainty surrounds the issue of whether precipitation over the tropical oceans (30deg N/S) systematically changes with interannual sea-surface temperature (SST) anomalies that accompany El Nino (warm) and La Nina (cold) events. Although it is well documented that El Nino-Southern Oscillation (ENSO) events with marked SST changes over the tropical oceans produce significant regional changes in precipitation, water vapor, and radiative fluxes in the tropics, we still cannot yet adequately quantify the associated net integrated changes to water and heat balance over the entire tropical oceanic or land sectors. Resolving this uncertainty is important since precipitation and latent heat release variations over land and ocean sectors are key components of the tropical heat balance in its most aggregated form. Rainfall estimates from the Version 5 Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR) averaged over the tropical oceans have not solved this issue and, in fact, show marked differences with estimates from two TRMM Microwave Imager (TMI) passive microwave algorithms. In this paper we will focus on findings that uncertainties in microphysical assumptions necessitated by the single-frequency PR measurement pose difficulties for detecting climate-related precipitation signals. Recent work has shown 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, yet the PR rainfall interannual variability (and attenuation derived predominantly from reflectivity) differs even in sign. We will explore these apparent inconsistencies and examine changes in new TRMM Version 6 retrievals. To place these results in a tropical water balance perspective we also examine interannual variations in evaporation over the tropical oceans made from TRMM and SSM/I (Special Sensor Microwave Imager) measurements of surface winds and humidity. Evaporation estimates from reanalysis and several global model experiments will also be compared to the TRMM findings and evaluated for consistency. The ability to detect regional shifts in freshwater flux over the oceans (equivalently, integrated moisture convergence) and moisture transport will be discussed.

Description: One notable aspect of Earth s climate is that although the planet appears to be very close to radiative balance at top-of-atmosphere (TOA), the atmosphere itself and underlying surface are not. Profound exchanges of energy between the atmosphere and oceans, land and cryosphere occur over a range of time scales. Recent evidence from broadband satellite measurements suggests that even these TOA fluxes contain some detectable variations. Our ability to measure and reconstruct radiative fluxes at the surface and at the top of atmosphere is improving rapidly. In this work we will evaluate two recently released estimates of radiative fluxes, focusing primarily on surface estimates. The International Satellite Cloud Climatology Project FD radiative flux profiles are available from mid-1 983 to near present and have been constructed by driving the radiative transfer physics from the Goddard Institute for Space Studies (GISS) global model with ISCCP clouds and TOVS (TIROS Operational Vertical Sounder)thermodynamic profiles. Fu!l and clear sky SW and LW fluxes are produced. A similar product from the NASA/GEWEX Surface Radiation Budget Project using different radiative flux codes and thermodynamics from the NAS/Goddard Earth Observing System (GEOS-1) assimilation model makes a similar calculation of surface fluxes. However this data set currently extends only through 1995. Significant differences in both interannual variability as well as trends are found between among these data sets. For radiative fluxes these differences are traced to TOVS thermodynamic soundings used to drive the ISCCP-FD calculations. Errors in near surface temperature and precipitable water cascade into ISCCP upward and downward IR flux components, demonstrably affecting interannual variability. Revised estimates of clear-sky fluxes over ocean are made using statistical algorithms and water vapor from the (SSM/I) Special Sensor Microwave Imager. These calculations show strong near-surface water vapor feedback over the tropical oceans in association with SST changes. However, it is also shown that ISCCP longwave cloud forcing, common to both the ISCCP-FD and GEWEX SRB retrievals, is the main driver of a long-term decrease in net LW flux to the surface during the near-20 year period covered by these revised estimates.

Description: The Goddard Chemistry Aerosol Radiation and Transport (GOCART) model is one of the fourteen global models participating in the International AEROCOM activity. The Oslo workshop is the 4th AEROCOM workshop. In this workshop, we are going to show updated model results of global aerosol extinction and absorption, vertical profiles, and direct and indirect effects. We will also discuss with other participants on model similarities and discrepancies, issues related to comparisons with in-situ and remote sensing data, and future directions.

Description: In the past few years, early solid precipitation detection and retrieval algorithms have been developed and shown to be applicable for snowing clouds and blizzards. NOAA has an operational snow versus rain classifier based on AMSU-B observations. Solid precipitation retrieval algorithms reported in the literature over the past two years include those that rely on neural nets, statistics, or physical relationships. All of the algorithms require the use of millimeter-wave radiometer observations. The millimeter-wave frequencies are especially sensitive to the scattering and emission properties of frozen particles due to the ice particle refractive index. Passive radiometric channels respond to both the integrated particle mass throughout the volume and field of view, and to the amount, location, and size distribution of the frozen (and liquid) particles with the sensitivity varying for different frequencies and hydrometeor types. This investigation probes the sensitivity of scattering and absorption coefficients, and hence computed brightness temperatures, resulting from variations in solid precipitation cloud profiles. The first study compares the single scattering, absorption, and asymmetry parameters associated with snow particles in clouds. Several methodologies are used to convert the physical characteristics (e.g., shape, size distributions, ice-air-water ratios) of ice particles to electromagnetic properties (e.g., absorption, scattering, and asymmetry factors). These methodologies include: conversion to solid ice particles, homogeneous dielectric mixing, or discrete dipole approximation. Changes in the conversion methodology can produce computed brightness temperature differences greater than 50 Kelvin.

Description: Recent GEWEX Cloud System Study (GCSS) model comparison projects have indicated that cloud-resolving models (CRMs) agree with observations better than traditional single-column models in simulating various types of clouds and cloud systems from different geographic locations. Current and future NASA satellite programs can provide cloud, precipitation, aerosol and other data at very fine spatial and temporal scales. It requires a coupled global circulation model (GCM) and cloud-scale model (termed a super-parameterization or multi-scale modeling framework, MMF) to use these satellite data to improve the understanding of the physical processes that are responsible for the variation in global and regional climate and hydrological systems. The use of a GCM will enable global coverage, and the use of a CRM will allow for better and ore sophisticated physical parameterization. NASA satellite and field campaign cloud related datasets can provide initial conditions as well as validation for both the MMF and CRMs. The Goddard MMF is based on the 2D Goddard Cumulus Ensemble (GCE) model and the Goddard finite volume general circulation model (fvGCM), and it has started production runs with two years results (1998 and 1999). Also, at Goddard, we have implemented several Goddard microphysical schemes (21CE, several 31CE), Goddard radiation (including explicity calculated cloud optical properties), and Goddard Land Information (LIS, that includes the CLM and NOAH land surface models) into a next generation regional scale model, WRF. In this talk, I will present: (1) A Brief review on GCE model and its applications on precipitation processes (microphysical and land processes), (2) The Goddard MMF and the major difference between two existing MMFs (CSU MMF and Goddard MMF), and preliminary results (the comparison with traditional GCMs), (3) A discussion on the Goddard WRF version (its developments and applications), and (4) The characteristics of the four-dimensional cloud data sets (or cloud library) stored at Goddard.

Description: Rainfall production is a fundamental process within the Earth;s hydrological cycle because it represents both a principal forcing term in surface water budgets, and its energetics corollary, latent heating, is the principal source of atmospheric diabatic heating. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The properties of the vertical distribution of latent heat release modulate large-scale meridional and zonal circulations with the Tropics - as well as modify the energetic efficiencies of mid-latitude weather systems. This paper highlights the retrieval of observatory, which was launched in November 1997 as a joint American-Japanese space endeavor. Since then, TRMM measurements have been providing an accurate four-dimensional amount of rainfall over the global Tropics and sub-tropics - information which can be used to estimate the spacetime structure of latent heating across the Earth's low latitudes. A set of algorithm methodologies has and continues to be developed to estimate latent heating based on rain rate profile retrievals obtained from TRMM measurements. These algorithms are briefly described followed by a discussion of the foremost latent heating products that can be generate from them. The investigation then provides an overview of how TRMM-derived latent heating information is currently being used in conjunction with global weather and climate models, concluding with remarks intended to stimulate further research on latent heating retrieval from satellites.

Description: Cloud-resolving case-study simulations of convective transport and lightning NO production have yielded results which are directly applicable to the design of lightning parameterizations for global chemical transport models. In this work we have used cloud-resolving models (the Goddard Cumulus Ensemble Model (GCE) and MMS) to drive an off-line cloud-scale chemical transport model (CSCTM). The CSCTM, in conjunction with aircraft measurements of NO(x) in thunderstorms and ground-based lightning observations, has been used to constrain the amount of NO produced per flash. Cloud and chemistry simulations for several case studies of storms in different environments will be presented. Observed lightning flash rates have been incorporated into the CSCTM, and several scenarios of NO production per intracloud (IC) and per cloud-to-ground (CG) flash have been tested for each storm. The resulting NO(x) mixing ratios are compared with aircraft measurements taken within the storm (typically the anvil region) to determine the most likely NO production scenario. The range of values of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced from observed NO spikes and from anvil flux calculations. Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. This result simplifies the lightning parameterization for global models (ie., an algorithm for estimating the ICKG ratio is not necessary). Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. Estimates of mean NO production per flash vary by a factor of three from one simulated storm to another. When combined with the global flash rate of 44 flashes per second from NASA's Optical Transient Detector (OTD) measurements, these estimates and the results from other techniques yield global NO production rates of 2-9 TgN/year. Simulations of the photochemistry over the 24 hours following a storm has been performed to determine the additional ozone production which can be attributed to lightning NO. Convective transport of HO(x), precursors leads to the generation of a HO(x), plume which substantially aids the downstream ozone production.

Description: The Lagrange multiplier theory and "pitch down method" developed in Part I of this study are applied to complete the calibration of a Citation aircraft that is instrumented with six field mill sensors. When side constraints related to average fields are used, the method performs well in computer simulations. For mill measurement errors of 1 V/m and a 5 V/m error in the mean fair weather field function, the 3-D storm electric field is retrieved to within an error of about 12%. A side constraint that involves estimating the detailed structure of the fair weather field was also tested using computer simulations. For mill measurement errors of 1 V/m, the method retrieves the 3-D storm field to within an error of about 8% if the fair weather field estimate is typically within 1 V/m of the true fair weather field. Using this side constraint and data from fair weather field maneuvers taken on 29 June 2001, the Citation aircraft was calibrated. The resulting calibration matrix was then used to retrieve storm electric fields during a Citation flight on 2 June 2001. The storm field results are encouraging and agree favorably with the results obtained from earlier calibration analyses that were based on iterative techniques.

Description: Effective management of energy resources is critical for the U.S. economy, the environment, and, more broadly, for sustainable development and alleviating poverty worldwide. The scope of energy management is broad, ranging from energy production and end use to emissions monitoring and mitigation and long-term planning. Given the extensive NASA Earth science research on energy and related weather and climate-related parameters, and rapidly advancing energy technologies and applications, there is great potential for increased application of NASA Earth science research to selected energy management issues and decision support tools. The NASA Energy Management Program Element is already involved in a number of projects applying NASA Earth science research to energy management issues, with a focus on solar and wind renewable energy and developing interests in energy modeling, short-term load forecasting, energy efficient building design, and biomass production.

Description: Using meteorology data, focusing on precipitable water (PW), obtained during the 2000-2003 thunderstorm seasons in Central Florida, this paper will, one, assess the skill and accuracy measurements of the current Mazany forecasting tool and, two, provide additional forecasting tools that can be used in predicting lightning. Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) are located in east Central Florida. KSC and CCAFS process and launch manned (NASA Space Shuttle) and unmanned (NASA and Air Force Expendable Launch Vehicles) space vehicles. One of the biggest cost impacts is unplanned launch scrubs due to inclement weather conditions such as thunderstorms. Each launch delay/scrub costs over a quarter million dollars, and the need to land the Shuttle at another landing site and return to KSC costs approximately $ 1M. Given the amount of time lost and costs incurred, the ability to accurately forecast (predict) when lightning will occur can result in significant cost and time savings. All lightning prediction models were developed using binary logistic regression. Lightning is the dependent variable and is binary. The independent variables are the Precipitable Water (PW) value for a given time of the day, the change in PW up to 12 hours, the electric field mill value, and the K-index value. In comparing the Mazany model results for the 1999 period B against actual observations for the 2000-2003 thunderstorm seasons, differences were found in the False Alarm Rate (FAR), Probability of Detection (POD) and Hit Rate (H). On average, the False Alarm Rate (FAR) increased by 58%, the Probability of Detection (POD) decreased by 31% and the Hit Rate decreased by 20%. In comparing the performance of the 6 hour forecast period to the performance of the 1.5 hour forecast period for the Mazany model, the FAR was lower by 15% and the Hit Rate was higher by 7%. However, the POD for the 6 hour forecast period was lower by 16% as compared to the POD of the 1.5 hour forecast period. Neither forecast period performed at the accuracy measures expected. A 2-Hr Forecasting Tool was developed to support a Phase I Lightning Advisory, which requires a 30-minute lead time for predicting lightning.

Description: The statistical distribution of the magnitude of the vector wind change over 0.25, 0.5, 1 and 2-h periods based on data from November 1999 through August 2001 is presented. The distributions of the 2-h u and v component wind changes are also presented for comparison. The wind changes at altitudes from 500 to 3000 m were measured using the Eastern Range network of five 915 MHz Doppler radar wind profilers. Quality controlled profiles were produced every 15 minutes for up to sixty gates, each representing 101 m in altitude over the range from 130 m to 6089 m. Five levels, each constituting three consecutive gates, were selected for analysis because of their significance to aerodynamic loads during the Space Shuttle ascent roll maneuver. The distribution of the magnitude of the vector wind change is found to be lognormal consistent with earlier work in the mid-troposphere. The parameters of the distribution vary with time lag, season and altitude. The component wind changes are symmetrically distributed with near-zero means, but the kurtosis coefficient is larger than that of a Gaussian distribution.

Description: The current understanding and modeling of the upper reaches of the atmosphere is incomplete. Upper atmospheric interactions with the lower atmosphere, effects of ionizing radiation, high altitude cloud phenomena, and the dynamical interaction with the magnetosphere require greater definition. The scientific objective of obtaining a greater understanding of the upper atmosphere can be achieved by designing, implementing, testing, and utilizing a facility that provides long period in-situ measurements of the mesosphere. Current direct sub-sonic measurements of the upper atmosphere are hampered by the approximately one minute sub-sonic observation window of a ballistic sounding rocket regardless of the launch angle. In-situ measurements at greater than transonic speeds impart energy into the molecular atmospheric system and distort the true atmospheric chemistry. A long duration, sub-sonic capability will significantly enhance our ability to observe and measure: (1) mesospheric lightning phenomena (sprites and blue jets) (2) composition, structure and stratification of noctilucent clouds (3) physics of seasonal radar echoes, gravity wave phenomena (4) chemistry of mesospheric gaseous ratio mixing (5) mesospheric interaction of ionizing radiation (6) dynamic electric and magnetic fields This new facility will also provide local field measurements which complement those that can be obtained through external measurements from satellite and ground-based platforms. The 400 foot (approximately 130 meter) diameter lightweight mega-mesospheric parachute system, deployed with a sounding rocket, is proposed herein as a method to increase sub-sonic mesospheric measurement time periods by more than an order of magnitude. The report outlines a multi-year evolving science instrumentation suite in parallel with the development of the mega meso-chute facility. The developmental issues surrounding the meso-chute are chiefly materials selection (thermal and structural) and deployment mechanism physics. Three mission cases were conceived and developed to include cost and schedules estimates. Each scenario has increasing scientific utility with paralleling launch weight, parachute hang-time, deployment altitude, and parachute size: (1) Case #1: $8.4M@24 months, 6kg payload, 20 min., 50km alt., 80 m. dia. (2) Case #2: $10.4M@24 months, 6kg payload, 20 min., 60km alt, 130m. dia. (3) Case #3: $13.6M@36 months, 30kg payload, 30 min., 90km alt., 200m. dia. The initial breakout cost for the parachute system is approximately $2M@24 months. This report identifies that although the challenges of the mega-meso-chute may be difficult, they can be surmounted and valuable results can be achieved.

Description: This paper is the second in a series in which kilometer-scale-resolving observations from the Atmospheric Radiation Measurement program and a cloud-resolving model (CRM) are used to evaluate the single-column model (SCM) version of the National Centers for Environmental Prediction Global Forecast System model. Part I demonstrated that kilometer-scale cirrus properties simulated by the SCM significantly differ from the cloud radar observations while the CRM simulation reproduced most of the cirrus properties as revealed by the observations. The present study describes an evaluation, through a comparison with the CRM, of the SCM's representation of detrainment from deep cumulus and ice-phase microphysics in an effort to better understand the findings of Part I. It is found that detrainment occurs too infrequently at a single level at a time in the SCM, although the detrainment rate averaged over the entire simulation period is somewhat comparable to that of the CRM simulation. Relatively too much detrained ice is sublimated when first detrained. Snow falls over too deep of a layer due to the assumption that snow source and sink terms exactly balance within one time step in the SCM. These characteristics in the SCM parameterizations may explain many of the differences in the cirrus properties between the SCM and the observations (or between the SCM and the CRM). A possible improvement for the SCM consists of the inclusion of multiple cumulus cloud types as in the original Arakawa-Schubert scheme, prognostically determining the stratiform cloud fraction and snow mixing ratio. This would allow better representation of the detrainment from deep convection, better coupling of the volume of detrained air with cloud fraction, and better representation of snow field.

Description: Citlalt6petl (Pico de Orizaba) is a dormant stratovolcano located at the eastern end of the trans-Mexican Volcanic Belt at approximately 19 degrees of latitude. It is one of the largest stratovolcanos in the world and at 5,630 meters above sea level, the highest mountain in Mexico and the third highest in North America. Situated on the summit cone and north face of the volcano is a permanent ice cap known as the Jamapa Glacier. Recent and historical studies of Citlaltepetl have been based primarily on volcanic risk assessment, in particular stability assessments of the summit cone. Relatively little work has been directed toward the glacial environment of the mountain, possibly due in part to its high altitude, steep slopes, and general inaccessibility. In addition to this glacier's potential to contribute to a better understanding of climate change, the Jamapa glacier and its environmental, cryologic and geologic setting could also serve as a valuable terrestrial analog to studies of Martian geology, hydrology, and subsurface ice.

Description: I have completed an update of global glacier volume change. All data of glacier annual mass balances, surface area over the period 1945/46 till 2004, outside the Greenland and Antarctic ice sheets were included in this update. As the result global glacier volume change have been calculated, also in terms of glacier contribution to sea level change. These results were sent to Working Group 1 and 2 of IPCC-4 as the basis for modeling of sea level towards the end of 2100. In this study I have concentrated on studying glacier systems of different scales, from primary (e.g. Devon ice cap) to regional (e.g. Canadian Arctic), continental scale (e,g., entire Arctic), and global (e.g., change in glacier volume and contribution to sea level rise).

Description: Until recently, the stable isotopic composition of chemically and datively important stratospheric species, such as ozone (O3), carbon dioxide (CO2), nitrous oxide (N2O), and methane (CH4), was largely unexplored, despite indications from the few measurements available and theoretical studies that global-scale isotopic variations will provide a unique tool for quantifying rates of global-scale mass transport into, within, and out of the stratosphere and for understanding the mechanisms of chemical reactions involved in ozone production. The number and geographical extent of observations are beginning to increase rapidly, however, as access to the stratosphere, both directly and by remote-sensing, has increased over the last 10 years and as new analytical techniques have been developed that make global-scale isotope measurements by whole-air sampling more feasible. The objective of this study, begun in April 1999, is to incorporate into the Livermore 2D model the likely photochemical fractionation processes that determine the isotopic compositions of stratospheric CO2, N2O, CH4, and O3, and to use the model results and new observations from NASA field campaigns in 1996 and 1997 to investigate stratospheric chemistry and mass transport. Additionally, since isotopic signatures from the stratosphere are transferred to the troposphere by downward transport at middle and high latitudes, the isotopic compositions may also serve as sensitive tracers of stratosphere-totroposphere transport. Comparisons of model results with stratospheric and upper tropospheric observations from these campaigns, as well as with ground-based observations from new NOAA and NSF-sponsored studies, will help determine whether the magnitudes of the stratospheric fractionation processes are large enough to use as global-scale tracers of transport into the troposphere and, if so, will be used to help constrain the degree of coupling between the troposphere and the stratosphere.

Description: A contrail tracking system was developed to help in the assessment of the effect of commercial jet contrails on the Earth's radiative budget. The tracking system was built by combining meteorological data from the Rapid Update Cycle (RUC) numerical weather prediction model with commercial air traffic flight track data and satellite imagery. A statistical contrail-forecasting model was created a combination of surface-based contrail observations and numerical weather analyses and forecasts. This model allows predictions of widespread contrail occurrences for contrail research on either a real-time basis or for long-term time scales. Satellite-derived cirrus cloud properties in polluted and unpolluted regions were compared to determine the impact of air traffic on cirrus.