ALBERT

Description: The recent Department of Energy Atmospheric Radiation Measurement (ARM) Aerosol Intensive Operations Period (AIOP, 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(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 six different instruments: airborne Sun photometer (AATS-14), airborne nephelometer/absorption photometer, airborne cavity ring-down system, groundbased Raman lidar, and two ground-based elastic backscatter lidars. We find the in situ measured sigma(ep)(lambda) to be lower than the AATS-14 derived values. Bias differences are 0.002-0.004 Km!1 equivalent to 13-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(ep)(lambda) are higher: Bias differences are 0.004 Km(-1) (13%) and 0.007 Km(-1) (24%) for the two elastic backscatter lidars (MPLNET and MPLARM, lambda = 523 nm) and 0.029 Km(-1) (54%) for the Raman lidar (lambda = 355 nm). 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 six field campaigns conducted since 1996, airborne in situ measurements of sigma(ep)(lambda) tend to be biased slightly low (17% at visible wavelengths) when compared to airborne Sun photometer sigma(ep)(lambda). On the other hand, sigma(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: In May 2003, a Twin Otter airplane, equipped with instruments for making in situ measurements of aerosol optical properties, was deployed during the Atmospheric Radiation Measurements (ARM) Program s Aerosol Intensive Operational Period in Oklahoma. Several of the Twin Otter flights were flown in formation with an instrumented light aircraft (Cessna 172XP) that makes routine in situ aerosol profile flights over the site. This paper presents comparisons of measured scattering coefficients at 467 nm, 530 nm, and 675 nm between identical commercial nephelometers aboard each aircraft. Overall, the agreement between the two nephelometers decreases with longer wavelength. During the majority of the flights, the Twin Otter flew with a diffuser inlet while the Cessna had a 1 mm impactor, allowing for an estimation of the fine mode fraction aloft. The fine mode fraction aloft was then compared to the results of a ground-based nephelometer. Comparisons are also provided in which both nephelometers operated with identical 1 mm impactors. These scattering coefficient comparisons are favorable at the longer wavelengths (i.e., 530 nm and 675 nm), yet differed by approximately 30% at 467 nm. Mie scattering calculations were performed using size distribution measurements, made during the level flight legs. Results are also presented from Cadenza, a new continuous wave cavity ring-down (CW-CRD) instrument, which compared favorably (i.e., agreed within 2%) with data from other instruments aboard the Twin Otter. With this paper, we highlight the significant implications of coarse mode (larger than 1 mm) aerosol aloft with respect to aerosol optical properties.

Description: Using the efficient discrete-ordinate method, we present an analytical solution for radiative transfer in the coupled atmosphere-ocean system with rough air-water interface. The theoretical formulations of the radiative transfer equation and solution are described. The effects of surface roughness on radiation field in the atmosphere and ocean are studied and compared with measurements. The results show that ocean surface roughness has significant effects on the upwelling radiation in the atmosphere and the downwelling radiation in the ocean. As wind speed increases, the angular domain of sunglint broadens, the surface albedo decreases, and the transmission to ocean increases. The downward radiance field in the upper ocean is highly anisotropic, but this anisotropy decreases rapidly as surface wind increases and as depth in ocean increases. The effects of surface roughness on radiation also depend greatly on both wavelength and angle of incidence (i.e., solar elevation); these effects are significantly smaller throughout the spectrum at high sun. The model-observation discrepancies may indicate that the Cox-Munk surface roughness model is not sufficient for high wind conditions.

Description: Raman lidar water vapor and aerosol extinction profiles acquired during the daytime over the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) site in northern Oklahoma (36.606 N, 97.50 W, 315 m) are evaluated using profiles measured by in situ and remote sensing instruments deployed during the May 2003 Aerosol Intensive Operations Period (IOP). The automated algorithms used to derive these profiles from the Raman lidar data were first modified to reduce the adverse effects associated with a general loss of sensitivity of the Raman lidar since early 2002. The Raman lidar water vapor measurements, which are calibrated to match precipitable water vapor (PWV) derived from coincident microwave radiometer (MWR) measurements were, on average, 5-10% (0.3-0.6 g/m(exp 3) higher than the other measurements. Some of this difference is due to out-of-date line parameters that were subsequently updated in the MWR PWV retrievals. The Raman lidar aerosol extinction measurements were, on average, about 0.03 km(exp -1) higher than aerosol measurements derived from airborne Sun photometer measurements of aerosol optical thickness and in situ measurements of aerosol scattering and absorption. This bias, which was about 50% of the mean aerosol extinction measured during this IOP, decreased to about 10% when aerosol extinction comparisons were restricted to aerosol extinction values larger than 0.15 km(exp -1). The lidar measurements of the aerosol extinction/backscatter ratio and airborne Sun photometer measurements of the aerosol optical thickness were used along with in situ measurements of the aerosol size distribution to retrieve estimates of the aerosol single scattering albedo (omega(sub o)) and the effective complex refractive index. Retrieved values of omega(sub o) ranged from (0.91-0.98) and were in generally good agreement with omega(sub o) derived from airborne in situ measurements of scattering and absorption. Elevated aerosol layers located between about 2.6 and 3.6 km were observed by the Raman lidar on May 25 and May 27. The airborne measurements and lidar retrievals indicated that these layers, which were likely smoke produced by Siberian forest fires, were primarily composed of relatively large particles (r(sub eff) approximately 0.23 micrometers), and that the layers were relatively nonabsorbing (omega(sub o) approximately 0.96-0.98). Preliminary results show that major modifications that were made to the Raman lidar system during 2004 have dramatically improved the sensitivity in the aerosol and water vapor channels and reduced random errors in the aerosol scattering ratio and water vapor retrievals by an order of magnitude.

Description: Many of the clouds important to the Earth's energy balance, from the tropics to the Arctic, are optically thin and contain liquid water. Longwave and shortwave radiative fluxes are very sensitive to small perturbations of the cloud liquid water path (LWP) when the liquid water path is small (i.e., 〈 g/sq m) and, thus, the radiative properties of these clouds must be well understood to capture them correctly in climate models. We review the importance of these thin clouds to the Earth's energy balance, and explain the difficulties in observing them. In particular, because these clouds are optically thin, potentially mixed-phase, and often (i.e., have large 3-D variability), it is challenging to retrieve their microphysical properties accurately. We describe a retrieval algorithm intercomparison that was conducted to evaluate the issues involved. The intercomparison included eighteen different algorithms to evaluate their retrieved LWP, optical depth, and effective radii. Surprisingly, evaluation of the simplest case, a single-layer overcast cloud, revealed that huge discrepancies exist among the various techniques, even among different algorithms that are in the same general classification. This suggests that, despite considerable advances that have occurred in the field, much more work must be done, and we discuss potential avenues for future work.

Description: The large horizontal extent, location in the cold upper troposphere, and ice composition make cirrus clouds important modulators of the earth's radiation budget and climate. Cirrus cloud microphysical properties are difficult to measure and model because they are inhomogeneous in nature and their ice crystal size distribution and habit are not well characterized. Accurate retrievals of cloud properties are crucial for improving the representation of cloud scale processes in large-scale models and for accurately predicting the earth's future climate. A number of passive and active remote sensing retrieval algorithms exist for estimating the microphysical properties of upper tropospheric clouds. We believe significant progress has been made in the evolution of these retrieval algorithms in the last decade, however, there is room for improvement. Members of the Atmospheric Radiation measurement program (ARM) Cloud properties Working Group are involved in an intercomparison of optical depth(tau), ice water path, and characteristic particle size in clouds retrieved using ground-based instruments. The goals of this intercomparison are to evaluate the accuracy of state-of-the-art algorithms, quantify the uncertainties, and make recommendations for improvement.

Description: The Nimbus 7 Limb Infrared Monitor of the Stratosphere (LIMS) radiance profile dataset of 1978/79 was reconditioned and reprocessed to Version 6 (V6) profiles of temperature and species that are improved significantly over those from Version 5 (V5). The LIMS V6 dataset was archived for public use in 2002. Improvements for its ozone include: (1) a more accurate accounting for instrument and spacecraft motion effects in the radiances, (2) the use of better spectroscopic line parameters for its ozone forward model, (3) retrievals of all its scans, (4) more accurate and compatible temperature versus pressure profiles (or T(p)), which are needed for the registration of the ozone radiances and for the removal of temperature effects from them, and (5) a better accounting for interfering species in the lower stratosphere. The retrieved V6 ozone profiles extend from near cloud top altitudes to about 80 km and from 64S to 84N latitude with better sampling along the orbit than for the V5 dataset. Calculated estimates of the single-profile precision and accuracy are provided for the V6 ozone from this study. Precision estimates based on the data themselves are of order 3% or better from 1 to 30 hPa. Estimates of total systematic error for a single profile are hard to generalize because the separate sources of error may not all be of the same sign and they depend somewhat on the atmospheric state. It is estimated that the V6 zonal mean ozone distributions are accurate to within 9% to 7% from 50 hPa to 3 hPa, respectively. Effects of a temperature bias can be significant and may be present at 1 to 2 hPa though. There may be ozone biases of order 10% at those levels due to possible biases of up to +2 K, but there is no indication of a similar problem elsewhere in the stratosphere. Simulation studies show that the LIMS retrievals are also underestimating slightly the small amplitudes of the atmospheric temperature tides, which affect its retrieved day/night ozone differences. There are small biases in the middle to lower stratosphere for the ascending versus descending node LIMS ozone, due principally to not accounting for the asymmetric weighting of its radiances across the tangent layer. The estimates of total accuracy were assessed by comparing the daily zonal mean LIMS ozone distributions against those from the Nimbus 7 SBUV Version 8 (V8) dataset for the same period. Generally, the LIMS V6 ozone agrees well with SBUV, except perhaps in the tropical lower stratosphere where the LIMS ozone is less. Still, the accuracy for LIMS V6 ozone in the lower stratosphere is improved over that found for LIMS V5, as indicated by several LIMS comparisons with ECC ozonesonde profiles. The LIMS V6 ozone is considered especially suitable for detailed studies of large-scale stratospheric processes above the 100-hPa level. Comparison of diurnal, photochemical model calculations with the monthly-averaged, upper stratospheric ozone observed from LIMS V6 indicates only a slight ozone deficit for the model at about 2 hPa. However, that deficit exhibits little to no seasonal variation and is in good agreement with similar model comparisons for a seasonal time series of ozone obtained with ground-based microwave instruments. Because the LIMS V6 ozone in the lower stratosphere has improved accuracy and sampling versus that of V5, it should now be possible to conduct quantitative studies of ozone transport and chemistry for the northern hemisphere, polar stratospheric winter of 1978/79 a time period when the catalytic loss of ozone due to reactive chlorine should not have been a major factor for the Arctic region.

Description: The relationships of upper tropospheric water vapor (UTWV), cloud ice and sea surface temperature (SST) are examined in the annual cycles of ECMWF analyses and simulations from 15 atmosphere-ocean coupled models which were contributed to the IPCC AR4. The results are compared with the observed relationships based on UTWV and cloud ice measurements from MLS on Aura. It is shown that the ECMWF analyses produce positive correlations between UTWV, cloud ice and SST, similar to the MLS data. The rate of the increase of cloud ice and UTWV with SST is about 30% larger than that for MLS. For the IPCC simulations, the relationships between UTWV, cloud ice and SST are qualitatively captured. However, the magnitudes of the simulated cloud ice show a considerable disagreement between models, by nearly a factor of 10. The amplitudes of the approximate linear relations between UTWV, cloud ice and SST vary by a factor up to 4.

Description: We provide an overview of the nadir measurements of carbon monoxide (CO) obtained thus far by the Tropospheric Emission Spectrometer (TES). The instrument is a high resolution array Fourier transform spectrometer designed to measure infrared spectral radiances from low Earth orbit. It is one of four instruments successfully launched onboard the Aura platform into a sun synchronous orbit at an altitude of 705 km on July 15, 2004 from Vandenberg Air Force Base, California. Nadir spectra are recorded at 0.06/cm spectral resolution with a nadir footprint of 5 x 8 km. We describe the TES retrieval approach for the analysis of the nadir measurements, report averaging kernels for typical tropical and polar ocean locations, characterize random and systematic errors for those locations, and describe instrument performance changes in the CO spectral region as a function of time. Sample maps of retrieved CO for the middle and upper troposphere from global surveys during December 2005 and April 2006 highlight the potential of the results for measurement and tracking of global pollution and determining air quality from space.

Description: This paper discusses the performance of AIRS and examines how it is meeting its operational and research objectives based on the experience of more than 2 yr with AIRS data. We describe the science background and the performance of AIRS in terms of the accuracy and stability of its observed spectral radiances. We examine the validation of the retrieved temperature and water vapor profiles against collocated operational radiosondes, and then we assess the impact thereof on numerical weather forecasting of the assimilation of the AIRS spectra and the retrieved temperature. We close the paper with a discussion on the retrieval of several minor tropospheric constituents from AIRS spectra.

Description: In this paper, an analytical treatment of the atmospheric remote sensing problem of determining the raindrop size distribution (DSD) with a spaceborne multifrequency microwave nadir-looking radar system is presented. It is typically assumed that with two radar measurements at different frequencies one ought to be able to calculate two state variables of the DSD: a bulk quantity, such as the rain rate, and a distribution shape parameter. To determine if this nonlinear problem can indeed be solved, the DSD is modeled as a Gamma distribution and quadratic approximations to the corresponding radar-rain relations are used to examine the invertibility of the resulting system of equations in the case of two as well as three radar frequencies. From the investigation, it is found that for regions of DSD state space multiple solutions exist for two or even three different frequency radar measurements. This should not be surprising given the nonlinear coupled nature of the problem.

Description: Polarization-sensitive lidars have proven to be highly effective in discriminating between spherical and non-spherical particles in the atmosphere. These lidars use a linearly polarized laser and are equipped with a receiver that can separately measure the components of the return signal polarized parallel and perpendicular to the outgoing beam. In this work we describe a technique for calibrating polarization-sensitive lidars that was originally developed at NASA s Langley Research Center (LaRC) and has been used continually over the past fifteen years. The procedure uses a rotatable half-wave plate inserted into the optical path of the lidar receiver to introduce controlled amounts of polarization cross-talk into a sequence of atmospheric backscatter measurements. Solving the resulting system of nonlinear equations generates the system calibration constants (gain ratio, G, and offset angle, theta) required for deriving calibrated measurements of depolarization ratio from the lidar signals. In addition, this procedure also determines the mean depolarization ratio within the region of the atmosphere that is analyzed. Simulations and error propagation studies show the method to be both reliable and well behaved. Operational details of the technique are illustrated using measurements obtained as part of Langley Research Center s participation in the First ISCCP Regional Experiment (FIRE).

Description: The semi-direct effects of dust aerosols are analyzed over eastern Asia using 2 years (June 2002 to June 2004) of data from the Clouds and the Earth s Radiant Energy System (CERES) scanning radiometer and MODerate Resolution Imaging Spectroradiometer (MODIS) on the Aqua satellite, and 18 years (1984 to 2001) of International Satellite Cloud Climatology Project (ISCCP) data. The results show that the water path of dust-contaminated clouds is considerably smaller than that of dust-free clouds. The mean ice water path (IWP) and liquid water path (LWP) of dusty clouds are less than their dust-free counterparts by 23.7% and 49.8%, respectively. The long-term statistical relationship derived from ISCCP also confirms that there is significant negative correlation between dust storm index and ISCCP cloud water path. These results suggest that dust aerosols warm clouds, increase the evaporation of cloud droplets and further reduce cloud water path, the so-called semi-direct effect. The semi-direct effect may play a role in cloud development over arid and semi-arid areas of East Asia and contribute to the reduction of precipitation.

Description: A previously unreported phenomenon, a 'frozen-in' anticyclone (FrIAC) after the 2005 Arctic spring vortex breakup, was discovered in Earth Observing System (EOS) Microwave Limb Sounder (MLS) long-lived trace gas data. A tongue of low-latitude (high-N2O, low-H2O) air was drawn into high latitudes and confined in a tight anticyclone, then advected intact in the summer easterlies through late August. A similar feature in O3 disappeared by early April as a result of chemical processes. The FrIAC was initially advected upright at nearly the same speed at all levels from approx.660 to 1300 K (approx.25-45 km); increasing vertical wind shear after early June tilted the FrIAC and weakened it at higher levels. The associated feature in PV disappeared by early June; transport calculations fail to reproduce the remarkable persistence of the FrIAC, suggesting deficiencies in summer high-latitude winds. The historical PV record suggests that this phenomenon may have occurred several times before. The lack of a persistent signature in O3 or PV, along with its small size and rapid motion, make it unlikely that a FrIAC could have been reliably identified without hemispheric daily longlived trace gas profiles such as those from EOS MLS.

Description: The atmospheric moisture and temperature profiles from the Atmospheric Infrared Sounder (AIRS)/Advanced Microwave Sounding Unit on the NASA Aqua mission, in combination with the precipitation from the Tropical Rainfall Measuring Mission (TRMM), are employed to study the vertical moist thermodynamic structure and spatial-temporal evolution of the Madden-Julian oscillation (MJO). The AIRS data indicate that, in the Indian Ocean and western Pacific, the temperature anomaly exhibits a trimodal vertical structure: a warm (cold) anomaly in the free troposphere (800-250 hPa) and a cold (warm) anomaly near the tropopause (above 250 hPa) and in the lower troposphere (below 800 hPa) associated with enhanced (suppressed) convection. The AIRS moisture anomaly also shows markedly different vertical structures as a function of longitude and the strength of convection anomaly. Most significantly, the AIRS data demonstrate that, over the Indian Ocean and western Pacific, the enhanced (suppressed) convection is generally preceded in both time and space by a low-level warm and moist (cold and dry) anomaly and followed by a low-level cold and dry (warm and moist) anomaly. The MJO vertical moist thermodynamic structure from the AIRS data is in general agreement, particularly in the free troposphere, with previous studies based on global reanalysis and limited radiosonde data. However, major differences in the lower-troposphere moisture and temperature structure between the AIRS observations and the NCEP reanalysis are found over the Indian and Pacific Oceans, where there are very few conventional data to constrain the reanalysis. Specifically, the anomalous lower-troposphere temperature structure is much less well defined in NCEP than in AIRS for the western Pacific, and even has the opposite sign anomalies compared to AIRS relative to the wet/dry phase of the MJO in the Indian Ocean. Moreover, there are well-defined eastward-tilting variations of moisture with height in AIRS over the central and eastern Pacific that are less well defined, and in some cases absent, in NCEP. In addition, the correlation between MJO-related mid-tropospheric water vapor anomalies and TRMM precipitation anomalies is considerably more robust in AIRS than in NCEP, especially over the Indian Ocean. Overall, the AIRS results are quite consistent with those predicted by the frictional Kelvin-Rossby wave/conditional instability of the second kind (CISK) theory for the MJO.

Description: We describe the clustering algorithm used by the Lightning Imaging Sensor (LIS) and the Optical Transient Detector (OTD) for combining the lightning pulse data into events, groups, flashes, and areas. Events are single pixels that exceed the LIS/OTD background level during a single frame (2 ms). Groups are clusters of events that occur within the same frame and in adjacent pixels. Flashes are clusters of groups that occur within 330 ms and either 5.5 km (for LIS) or 16.5 km (for OTD) of each other. Areas are clusters of flashes that occur within 16.5 km of each other. Many investigators are utilizing the LIS/OTD flash data; therefore, we test how variations in the algorithms for the event group and group-flash clustering affect the flash count for a subset of the LIS data. We divided the subset into areas with low (1-3), medium (4-15), high (16-63), and very high (64+) flashes to see how changes in the clustering parameters affect the flash rates in these different sizes of areas. We found that as long as the cluster parameters are within about a factor of two of the current values, the flash counts do not change by more than about 20%. Therefore, the flash clustering algorithm used by the LIS and OTD sensors create flash rates that are relatively insensitive to reasonable variations in the clustering algorithms.

Description: ICE: The Nature, the History, and the Uses of an Astonishing Substance is a recently published book by Mariana Gosnell about ice. It covers not just the ice that is readily seen, such as sea ice, lake ice, icebergs, glaciers, ice sheets, and ice cubes, but also ice in the ground, in the atmosphere, inside plants and animals, and in outer space, plus new ice forms being created in scientific laboratories. Gosnell treats the reader to a well-written, easy-going mixture of science, adventure, history, applications, science methods and controversies, and philosophy, all centered in one way or another on ice. The book is 563 pages long and is filled with fascinating anecdotes and details, such as beetles in the Canadian Rockies that can supercool to 60 C below freezing and a lake in Minnesota where each winter typically 65,000 fishing shanties are set up on the lake's ice, many with couches, beds, television sets, and bathrooms. Gosnell also includes many practical suggestions. Among them: When driving on lake ice, keep your windows open, in case your vehicle breaks through the ice and you need to make a rapid exit.

Description: From analyses of blended space-based and ground-based global rainfall data, we found increasing trends in the occurrence of extreme heavy and light rain events, coupled to a decreasing trend in moderate rain events in the tropics during 1979-2003. The trends are consistent with a shift in the large-scale circulation associated with a) a relatively uniform increase in warm rain over the tropical oceans, b) enhanced ice-phase rain over the near-equatorial oceans, and c) reduced mixed-phase rain over the tropical ocean and land regions. Due to the large compensation among different rain categories, the total tropical rainfall trend remained undetectable.

Description: Aerosols over the Northeastern Pacific Ocean enhance the cloud drop number concentration and reduce the drop size for marine stratocumulus and cumulus clouds. These microphysical effects result in brighter clouds, as evidenced by a combination of aircraft and satellite observations. In-situ measurements from the Cloud Indirect Forcing Experiment (CIFEX) indicate that the mean cloud drop number concentration in low clouds over the polluted marine boundary layer is greater by 53/cu cm compared to clean clouds, and the mean cloud drop effective radius is smaller by 4 microns. We link these in-situ measurements of cloud modification by aerosols, for the first time, with collocated satellite broadband radiative flux observations from the Clouds and the Earth's Radiant Energy System (CERES) to show that these microphysical effects of aerosols enhance the top-of-atmosphere cooling by -9.9+/-4.3 W/sq m for overcast conditions.

Description: This study provides explanations for some of the experimental findings of Chao (2000) and Chao and Chen (2001) concerning the mechanisms responsible for the ITCZ in an aqua-planet model. These explanations are then applied to explain the origin of some of the systematic errors in the GCM simulation of ITCZ precipitatin over oceans. The ITCZ systematic errors are highly sensitive to model physics and by extension model horizontal resolution. The findings in this study along with those of Chao (2000) and Chao and Chen (2001, 2004) contribute to building a theoretical foundation for ITCZ study. A few possible methods of alleviating the systematic errors in the GCM simulaiton of ITCZ are discussed. This study uses a recent version of the Goddard Modeling and Assimilation Office's Goddard Earth Observing System (GEOS-5) GCM.

Description: Solar proton fluxes have been measured by satellites for over forty years (1963-2005). Several satellites, including the NASA Interplanetary Monitoring Platforms (1963-1993) and the NOAA Geostationary Operational Environmental Satellites (1994-2005), have been used to compile this long-term dataset. Some solar eruptions lead to solar proton events (SPEs) at the Earth, which typically last a few days. High energy solar protons associated with SPEs precipitate on the Earth's atmosphere and cause increases in odd hydrogen (HOx) and odd nitrogen (NOy) in the polar cap regions (greater than 60 degrees geomagnetic). The enhanced HOx leads to short-lived ozone depletion (~days) due to the short lifetime of HOx constituents. The enhanced NOy leads to long-lived ozone changes because of the long lifetime of the NOy family in the stratosphere and lower mesosphere. Very large SPEs occurred in 1972, 1989, 2000, 2001, and 2003 and were predicted to cause maximum total ozone depletions of 1-3%, which lasted for several months to years past the events. These long-term ozone changes caused by SPES are discussed.

Description: The SBUV/2 (Solar Backscattered Ultraviolet, model 2) instrument is designed to monitor ozone stratospheric profile and total column ozone using measurements of the Earth's backscattered ultraviolet albedo. We have previously demonstrated that the normal radiance measurements from SBUV/2 instruments, which sample 12 discrete wavelengths between 252 and 340 nm during each scan, can be used to identify polar mesospheric clouds (PMCs). Some SBUV/2 instruments also periodically view the earth in continuous scan mode, covering the wavelength range 160-400 nm with 0.15 nm sampling. Analysis of these data show PMC occurrence rates similar to the normal discrete scan results, although the observation technique reduces the number of daily measurements by a factor of six. PMC observed by SBUV/2 instruments show a monotonic variation in the residual spectral albedo over the wavelength range 250 300 nm, with maximum enhancements of 10 15% at 250 nm. This result is consistent with microphysical model predictions from Jensen [1989. A numerical model of polar mesospheric cloud formation and evolution, Ph. D. Thesis, University of Colorado]. We find no evidence for a systematic localized increase in PMC residual albedo for wavelengths near 260 nm, in contrast to the recently reported results from the MSX UVISI instrument [Carbary J.F., et al., 2004. Evidence for bimodal particle distribution from the spectra of polar mesospheric clouds. Geophysics Research. Letters 31, L13108]. This result is observed for three different SBUV/2 instruments in both Northern and Southern Hemisphere data over a 13-year span. Our Mie scattering calculations show that the location and magnitude of the 260 nm hump feature is dependent upon the specific scattering angles appropriate to the MSX measurements. Although it explains the MSX spectrum, the bimodal size distribution proposed by Carbary et al. (2004), cannot explain the lack of scattering angle dependence of the SBUV/2 spectral shapes. The spectral signature of the SBUV/2 continuous scan PMC data is thus inconsistent with the bimodal particle size distribution suggested by Carbary et al. (2004).

Description: We describe the operational algorithm for the retrieval of stratospheric, tropospheric, and total column densities of nitrogen dioxide NO2 from earthshine radiances measured by the Ozone Monitoring Instrument (OMI), aboard the EOS-Aura satellite. The algorithm uses the DOAS method for the retrieval of slant column NO densities. Air mass factors (AMFs) calculated from a stratospheric NO2 profile are used to make initial estimates of the vertical column density. Using data collected over a 24-h period, a smooth estimate of the global stratospheric field is constructed. Where the initial vertical column densities exceed the estimated stratospheric field, we infer the presence of tropospheric NO2, and recalculate the vertical column density (VCD) using an AMF calculated from an assumed tropospheric NO2 profile. The parameters that control the operational algorithm were selected with the aid of a set of data assembled from stratospheric and tropospheric chemical transport models. We apply the optimized algorithm to OMI data and present global maps of NO2 VCDs for the first time.

Description: Two approximations to convective transport have been implemented in an offline chemistry transport model (CTM) to explore the impact on calculated atmospheric CO2 distributions. GlobalCO2 in the year 2000 is simulated using theCTM driven by assimilated meteorological fields from the NASA s Goddard Earth Observation System Data Assimilation System, Version 4 (GEOS-4). The model simulates atmospheric CO2 by adopting the same CO2 emission inventory and dynamical modules as described in Kawa et al. (convective transport scheme denoted as Conv1). Conv1 approximates the convective transport by using the bulk convective mass fluxes to redistribute trace gases. The alternate approximation, Conv2, partitions fluxes into updraft and downdraft, as well as into entrainment and detrainment, and has potential to yield a more realistic simulation of vertical redistribution through deep convection. Replacing Conv1 by Conv2 results in an overestimate of CO2 over biospheric sink regions. The largest discrepancies result in a CO2 difference of about 7.8 ppm in the July NH boreal forest, which is about 30% of the CO2 seasonality for that area. These differences are compared to those produced by emission scenario variations constrained by the framework of Intergovernmental Panel on Climate Change (IPCC) to account for possible land use change and residual terrestrial CO2 sink. It is shown that the overestimated CO2 driven by Conv2 can be offset by introducing these supplemental emissions.

Description: The Atmospheric InfraRed Sounder (AIRS), flying aboard NASA's Earth Observing System (EOS) Aqua satellite with the Advanced Microwave Sounding Unit-A (AMSU-A), has been providing data for use in numerical weather prediction (NWP) and data assimilation systems (DAS) for over three years. The full AIRS data set is currently not transmitted in near-real-time (NRT) to the NWP centers. Instead, data sets with reduced spatial and spectral information are produced and made available in NRT. In this paper, we evaluate the use of different channel selections and error specifications. We achieved significant positive impact from the Aqua AIRS/AMSU-A combination in both hemispheres during our experimental time period of January 2003. The best results were obtained using a set of 156 channels that did not include any in the 6.7micron water vapor band. The latter have a large influence on both temperature and humidity analyses. If observation and background errors are not properly specified, the partitioning of temperature and humidity information from these channels will not be correct, and this can lead to a degradation in forecast skill. We found that changing the specified channel errors had a significant effect on the amount of data that entered into the analysis as a result of quality control thresholds that are related to the errors. However, changing the channel errors within a relatively small window did not significantly impact forecast skill with the 155 channel set. We also examined the effects of different types of spatial data reduction on assimilated data sets and NWP forecast skill. Whether we picked the center or the warmest AIRS pixel in a 3x3 array affected the amount of data ingested by the analysis but had a negligible impact on the forecast skill.

Description: Bridging the gap between current research and the classroom is a major challenge to today s instructor, especially in the sciences where progress happens quickly. NASA Goddard Space Flight Center and the University of Maryland teamed up in designing a graduate class project intended to provide a hands-on introduction to the physical basis for the retrieval of aerosol properties from state-of-the-art MODIS observations. Students learned to recognize spectral signatures of atmospheric aerosols and to perform spectral inversions. They became acquainted with the operational MODIS aerosol retrieval algorithm over oceans, and methods for its evaluation, including comparisons with groundbased AERONET sun-photometer data.

Description: The fifth-generation Pennsylvania State University National Center for Atmospheric Research (PSU NCAR) Mesoscale Model (MM5) is used to simulate Hurricane Erin (2001) at high resolution (4-km spacing) from its early development as a tropical depression on 7 September 2001, through a period of rapid intensification into a strong hurricane (8 9 September), and finally into a stage during which it maintains its intensity on 10 September. These three stages of formation, intensification, and maintenance in the simulation are in good agreement with the observed evolution of Erin. The simulation shows that during the formation and early portions of the intensification stages, intensification is favored because the environmental wind shear is weak and the system moves over a warm tongue of water. As Erin intensifies, the wind shear gradually increases with the approach of an upper-level trough and strengthening of a low-level high pressure system. By 10 September, the wind shear peaks and begins to decrease, the storm moves over slightly cooler waters, and the intensification ends. Important structural changes occur at this time as the outer precipitation shifts from the northeastern and eastern sides to the western side of the eye. A secondary wind maximum and an outer eyewall begin to develop as precipitation begins to surround the entire eye. The simulation is used to investigate the role of vertical wind shear in the changes of the precipitation structure that took place between 9 and 10 September by examining the effects of changes in storm-relative flow and changes in the shear-induced tilt. Qualitative agreement is found between the divergence pattern and advection of vorticity by the relative flow with convergence (divergence) generally associated with asymmetric inflow (outflow) in the eyewall region. The shift in the outer precipitation is consistent with a shift in the low-level relative inflow from the northeastern to the northwestern side of the storm. The changes in the relative flow are associated with changes in the environmental winds as the hurricane moves relative to the upper trough and the low-level high pressure system. Examination of the shear-induced tilt of the vortex shows that the change in the tilt direction is greater than that of the shear direction as the tilt shifts from a northerly orientation to northwesterly. Consistent with theory for adiabatic vortices, the maximum low-level convergence and upper-level divergence (and the maximum upward motion) occurs in the direction of tilt. Consequently, both mechanisms may play roles in the changes in the precipitation pattern.

Description: Precipitation estimation from satellite passive microwave radiometer observations is a problem that does not have a unique solution that is insensitive to errors in the input data. Traditionally, to make this problem well posed, a priori information derived from physical models or independent, high-quality observations is incorporated into the solution. In the present study, a database of precipitation profiles and associated brightness temperatures is constructed to serve as a priori information in a passive microwave radiometer algorithm. The precipitation profiles are derived from a Tropical Rainfall Measuring Mission (TRMM) combined radar radiometer algorithm, and the brightness temperatures are TRMM Microwave Imager (TMI) observed. Because the observed brightness temperatures are consistent with those derived from a radiative transfer model embedded in the combined algorithm, the precipitation brightness temperature database is considered to be physically consistent. The database examined here is derived from the analysis of a month-long record of TRMM data that yields more than a million profiles of precipitation and associated brightness temperatures. These profiles are clustered into a tractable number of classes based on the local sea surface temperature, a radiometer-based estimate of the echo-top height (the height beyond which the reflectivity drops below 17 dBZ), and brightness temperature principal components. For each class, the mean precipitation profile, brightness temperature principal components, and probability of occurrence are determined. The precipitation brightness temperature database supports a radiometer-only algorithm that incorporates a Bayesian estimation methodology. In the Bayesian framework, precipitation estimates are weighted averages of the mean precipitation values corresponding to the classes in the database, with the weights being determined according to the similarity between the observed brightness temperature principal components and the brightness temperature principal components of the classes. Because the classes are stratified by the sea surface temperature and the echo-top-height estimator, the number of classes that are considered for retrieval is significantly smaller than the total number of classes, making the algorithm computationally efficient. The radiometer-only algorithm is applied to TMI observations, and precipitation estimates are compared with combined TRMM precipitation radar (PR) TMI reference estimates. The TMI-only algorithm, supported by the empirically derived database, produces estimates that are more consistent with the reference values than the precipitation estimates from the version-6 TRMM facility TMI algorithm. Cloud-resolving model simulations are used to assign a latent heating profile to each precipitation profile in the empirically derived database, making it possible to estimate latent heating using the radiometer-only algorithm. Although the evaluation of latent heating estimates in this study is preliminary, because realistic conditional probability distribution functions are attached to latent heating structures in the algorithm s database, a generally positive impact on latent heating estimation from passive microwave observations is expected.

Description: In this paper we present results of a numerical study using the NASA finite-volume GCM to elucidate a plausible mechanism for aerosol impact on the Asian summer monsoon involving interaction with physical processes over the Tibetan Plateau (TP). During the premonsoon season of March April, dusts from the deserts of western China, Afghanistan/Pakistan, and the Middle East are transported into and stacked up against the northern and southern slopes of the TP. The absorption of solar radiation by dust heats up the elevated surface air over the slopes. On the southern slopes, the atmospheric heating is reinforced by black carbon from local emission. The heated air rises via dry convection, creating a positive temperature anomaly in the mid-to-upper troposphere over the TP relative to the region to the south. In May through early June in a manner akin to an elevated heat pump , the rising hot air forced by the increasing heating in the upper troposphere, draws in warm and moist air over the Indian subcontinent, setting the stage for the onset of the South Asia summer monsoon. Our results suggest that increased dust loading coupled with black carbon emission from local sources in northern India during late spring may lead to an advance of the rainy periods and subsequently an intensification of the Indian summer monsoon. The enhanced rainfall over India is associated with the development of an aerosol-induced large-scale sea level pressure anomaly pattern, which causes the East Asia (Mei-yu) rain belt to shift northwestward, suppressing rainfall over East Asia and the adjacent oceanic regions.

Description: Rainfall rate estimates from spaceborne microwave radiometers are generally accepted as reliable by a majority of the atmospheric science community. One of the Tropical Rainfall Measuring Mission (TRMM) facility rain-rate algorithms is based upon passive microwave observations from the TRMM Microwave Imager (TMI). In Part I of this series, improvements of the TMI algorithm that are required to introduce latent heating as an additional algorithm product are described. Here, estimates of surface rain rate, convective proportion, and latent heating are evaluated using independent ground-based estimates and satellite products. Instantaneous, 0.5 deg. -resolution estimates of surface rain rate over ocean from the improved TMI algorithm are well correlated with independent radar estimates (r approx. 0.88 over the Tropics), but bias reduction is the most significant improvement over earlier algorithms. The bias reduction is attributed to the greater breadth of cloud-resolving model simulations that support the improved algorithm and the more consistent and specific convective/stratiform rain separation method utilized. The bias of monthly 2.5 -resolution estimates is similarly reduced, with comparable correlations to radar estimates. Although the amount of independent latent heating data is limited, TMI-estimated latent heating profiles compare favorably with instantaneous estimates based upon dual-Doppler radar observations, and time series of surface rain-rate and heating profiles are generally consistent with those derived from rawinsonde analyses. Still, some biases in profile shape are evident, and these may be resolved with (a) additional contextual information brought to the estimation problem and/or (b) physically consistent and representative databases supporting the algorithm. A model of the random error in instantaneous 0.5 deg. -resolution rain-rate estimates appears to be consistent with the levels of error determined from TMI comparisons with collocated radar. Error model modifications for nonraining situations will be required, however. Sampling error represents only a portion of the total error in monthly 2.5 -resolution TMI estimates; the remaining error is attributed to random and systematic algorithm errors arising from the physical inconsistency and/or nonrepresentativeness of cloud-resolving-model-simulated profiles that support the algorithm.

Description: The direct effects of aerosols on global and regional climate during boreal spring are investigated based on numerical simulations with the NASA Global Modeling and Assimilation Office finite-volume general circulation model (fvGCM) with Microphyics of Clouds with the Relaxed Arakawa Schubert Scheme (McRAS), using aerosol forcing functions derived from the Goddard Ozone Chemistry Aerosol Radiation and Transport model (GOCART). The authors find that anomalous atmospheric heat sources induced by absorbing aerosols (dust and black carbon) excite 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, that is, South Asia, East Asia, and northern and western Africa. Significant atmospheric heating is found in regions with large loading of dust (over northern Africa and the Middle East) and black carbon (over Southeast 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 concentrations are low. Analyses of circulation anomalies show that the dipole anomaly is a part of an atmospheric teleconnection pattern 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. The surface temperature signature associated with the aerosol-induced teleconnection bears striking resemblance to the spatial pattern of observed long-term trend in surface temperature over Eurasia. Additionally, the boreal spring wave train pattern is similar to that reported by Fukutomi et al. associated with the boreal summer precipitation seesaw between eastern and western Siberia. The results of this study raise the possibility that global aerosol forcing during boreal spring may play an important role in spawning atmospheric teleconnections that affect regional and global climates.

Description: Recent advances in space-borne observations and numerical weather prediction models provide new opportunities for improving hurricane forecasts. In this study, state-of-the-art satellite observations are used to document the evolution of one of the most devastating tropical cyclones ever to hit the United States: Hurricane Katrina. The ECMWF and NASA global high-resolution forecasts, the latter being run in experimental mode, are compared with satellite observations, with a focus on precipitation and cloud processes. Future directions on modeling and observations are briefly discussed.

Description: The Earth continues to have a third of the ice that it had at the peak of the last ice age, although that ice continues to decrease, as it has, overall, for the past 18,000 years. Over the last 100 years, the retreat signal has been especially strong in ice shelves of the Arctic and along the Antarctic Peninsula, with a more mixed signal elsewhere. For instance, since the early 1990s the massive Greenland and Antarctic ice sheets have thinned along the coasts but thickened in the interior, and since the late 1970s sea ice has decreased in the Arctic but increased (slightly) in the Antarctic. Major difficulties in the interpretations of the climate record come from the high interannual variability of most cryosphere components and the lack of consistent long-term global data records, the latter problem now being slowly remedied, in part, through satellite technology.

Description: AIRS was launched on EOS Aqua on May 4, 2002, together with AMSU A and HSB, to form a next generation polar orbiting infrared and microwave atmospheric sounding system. The primary products of AIRS/AMSU are twice daily global fields of atmospheric temperature-humidity profiles, ozone profiles, sea/land surface skin temperature, and cloud related parameters including OLR. The sounding goals of AIRS are to produce 1 km tropospheric layer mean temperatures with an rms error of 1 K, and layer precipitable water with an rms error of 20 percent, in cases with up to 80 percent effective cloud cover. The basic theory used to analyze AIRS/AMSU/HSB data in the presence of clouds, called the at-launch algorithm, and a post-launch algorithm which differed only in the minor details from the at-launch algorithm, have been described previously. The post-launch algorithm, referred to as AIRS Version 4.0, has been used by the Goddard DAAC to analyze and distribute AIRS retrieval products. In this paper we show progress made toward the AIRS Version 5.0 algorithm which will be used by the Goddard DAAC starting late in 2006. A new methodology has been developed to provide accurate case by case error estimates for retrieved geophysical parameters and for the channel by channel cloud cleared radiances used to derive the geophysical parameters from the AIRS/AMSU observations. These error estimates are in turn used for quality control of the derived geophysical parameters and clear column radiances. Improvements made to the retrieval algorithm since Version 4.0 are described as well as results comparing Version 5.0 retrieval accuracy and spatial coverage with those obtained using Version 4.0.

Description: Polar bears are dependent on sea ice for survival. Climate warming in the Arctic has caused significant declines in coverage and thickness of sea ice in the polar basin and progressively earlier breakup in some areas. In four populations of polar bears in the eastern Canadian Arctic (including Western Hudson Bay), Inuit hunters report more bears near settlements during the open water period in recent years. These observations have been interpreted as evidence of increasing population size, resulting in increases in hunting quotas. However, long-term data on the population size and condition of polar bears in Western Hudson Bay, and population and harvest data from Baffin Bay, make it clear that those two populations at least are declining, not increasing. While the details vary in different arctic regions, analysis of passive-microwave satellite imagery, beginning in the late 1970s, indicates that the sea ice is breaking up at progressively earlier dates, so that bears must fast for longer periods during the open water season. Thus, at least part of the explanation for the appearance of more bears in coastal communities is likely that they are searching for alternative food sources because their stored body fat depots are being exhausted. We hypothesize that, if the climate continues to warm as projected by the IPCC, then polar bears in all five populations discussed in this paper will be stressed and are likely to decline in numbers, probably significantly so. As these populations decline, there will likely also be continuing, possibly increasing, numbers of problem interactions between bears and humans as the bears seek alternate food sources. Taken together, the data reported in this paper suggest that a precautionary approach be taken to the harvesting of polar bears and that the potential effects of climate warming be incorporated into planning for the management and conservation of this species throughout the Arctic.

Description: This special issue of the Journal of' Atmospheric and Solar-Terrestrial Physics is a compilation of 23 papers presented at The 2004 Huntsville Modeling Workshop: Challenges to Modeling thc San-Earth System held in Huntsville, AB on October 18-22, 2004. The title of the workshop appropriately captures the theme of what was presented and discussed by the 120 participants. Currently, end-to-end modeling of the Sun-Earth system is a major goal of the National Space Weather and NASA living with a star (LWS) programs. While profound advances have been made in modeling isolated regions of the Sun-Earth system, minimal progress has been achieved in modeling the end-to-end system. The transfer of mass, energy and momentum through the coupled Sun-Earth system spans a wide range of scales inn time and space. The uncertainty in the underlying physics responsible for coupling contiguous regions of the Sun-Earth system is recognized as a significant barrier to progress

Description: A multi-scale modeling framework (MMF), which replaces the conventional cloud parameterizations with a cloud-resolving model (CRM) in each grid column of a GCM, constitutes a new and promising approach. The MMF can provide for global coverage and two-way interactions between the CRMs and their parent GCM. The GCM allows global coverage and the CRM allows explicit simulation of cloud processes and their interactions with radiation and surface processes. A new MMF has been developed that is based the Goddard finite volume GCM (fvGCM) and the Goddard Cumulus Ensemble (GCE) model. This Goddard MMF produces many features that are similar to another MMF that was developed at Colorado State University (CSU), such as an improved .surface precipitation pattern, better cloudiness, improved diurnal variability over both oceans and continents, and a stronger, propagating Madden-Julian oscillation (MJO) compared to their parent GCMs using conventional cloud parameterizations. Both MMFs also produce a precipitation bias in the western Pacific during Northern Hemisphere summer. However, there are also notable differences between two MMFs. For example, the CSU MMF simulates less rainfall over land than its parent GCM. This is why the CSU MMF simulated less overall global rainfall than its parent GCM. The Goddard MMF overestimates global rainfall because of its oceanic component. Some critical issues associated with the Goddard MMF are presented in this paper.

Description: Snow fall and snow accumulation are key climate parameters due to the snow's high albedo, its thermal insulation, and its importance to the global water cycle. Satellite passive microwave radiometers currently provide the only means for the retrieval of snow depth and/or snow water equivalent (SWE) over land as well as over sea ice from space. All algorithms make use of the frequency-dependent amount of scattering of snow over a high-emissivity surface. Specifically, the difference between 37- and 19-GHz brightness temperatures is used to determine the depth of the snow or the SWE. With the availability of the Advanced Microwave Scanning Radiometer (AMSR-E) on the National Aeronautics and Space Administration's Earth Observing System Aqua satellite (launched in May 2002), a wider range of frequencies can be utilized. In this study we investigate, using model simulations, how snow depth retrievals are affected by the evolution of the physical properties of the snow (mainly grain size growth and densification), how they are affected by variations in atmospheric conditions and, finally, how the additional channels may help to reduce errors in passive microwave snow retrievals. The sensitivity of snow depth retrievals to atmospheric water vapor is confirmed through the comparison with precipitable water retrievals from the National Oceanic and Atmospheric Administration's Advanced Microwave Sounding Unit (AMSU-B). The results suggest that a combination of the 10-, 19-, 37-, and 89-GHz channels may significantly improve retrieval accuracy. Additionally, the development of a multisensor algorithm utilizing AMSR-E and AMSU-B data may help to obtain weather-corrected snow retrievals.

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 3-D) 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-dimentional(3-D) 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 microns) and one with liquid water efficient absorption of solar radiation (2.1 microns). In contrast to the plane-parallel approximation, a conventional approach to all current operational retrievals, 3-D 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: Three-dimensional reconstructions of six rocket-triggered lightning channels are derived from stereo photographs. These reconstructed channels are used to infer the behavior of the current in return strokes above the ground from current waveforms measured at the channel base and electric-field-change waveforms measured at a range of 5.2 kilometers for 24 return strokes in these channels. Streak photographs of 14 of the same strokes are analyzed to determine the rise times, propagation speeds, and amplitudes of relative light intensity for comparison with the electrical inferences. Results include the following: 1) The fine structure of the field-change waveforms that were radiated by these subsequent return strokes can be explained, in large part, by channel geometry. 2) The average 10 - 90% rise time of the stroke current increased by about a factor of seven in our sample, from an observed 0.31 plus or minus 0.17 microseconds at the surface to an inferred 2.2 plus or minus 0.5 microcseconds at 1 kilometer path length above the surface. 3) The three-dimensional propagation speed of the current front averaged 1.80 plus or minus 0.24 X 10(exp 8) meters per second over channel lengths typically greater than 1 kilometer. 4) Assuming that the measured current was entirely due to the return stroke forced an unreasonably large and abrupt reduction in inferred current amplitude over the first few tens of meters above the surface, especially in cases when the leader was bright relative to its stroke. Therefore, a significant fraction of the current at the surface was probably due to the leader, at least in such cases. 5) Peak return-stroke currents decreased by approximately 37 plus or minus 12% from 100 meters to 1 kilometer of path length above the surface. Because of uncertainty about how to partition the measured current between leader and return stroke, we are unable to infer the variation of current amplitude near the ground.

Description: The Earth's radiative energy balance and hydrological cycle are fundamentally coupled with the distribution and properties of clouds. Therefore, the ability to remotely infer cloud properties and their variation in space and time is crucial for establishing climatologies as a reference for validation of present-day climate models and in assessing future climate change. Remote cloud observations also provide data sets useful for testing and improving cloud model physics, and for assimilation into numerical weather prediction models. The MODerate Resolution Imaging Spectroradiometer (MODIS) imagers on the Terra and Aqua Earth Observing System (EOS) platforms provide the capability for globally retrieving these properties using passive solar reflectance and infrared techniques. In addition to providing measurements similar to those offered on a suite of historical operational weather platforms such as the Advanced Very High Resolution Radiometer (AVHRR), the High-resolution Infrared Radiation Sounder (HIRS), and the Geostationary Operational Environmental Satellite (GOES), MODIS provides additional spectral and/or spatial resolution in key atmospheric bands, along with on-board calibration, to expand the capability for global cloud property retrievals. The core MODIS operational cloud products include cloud top pressure, thermodynamic phase, optical thickness, particle size, and water path, and are derived globally at spatial resolutions of either 1- or 5-km (referred to as Level-2 or pixel-level products). In addition, the MODIS atmosphere team (collectively providing cloud, aerosol, and clear sky products) produces a combined gridded product (referred to as Level-3) aggregated to a 1 equal-angle grid, available for daily, eight-day, and monthly time periods. The wealth of information available from these products provides critical information for climate studies as well as the continuation and improved understanding of existing satellite-based cloud climatologies derived from heritage instruments. This chapter provides an overview of the MODIS Level-2 and -3 operational cloud products.

Description: Tropical Storm Chantal during August 2001 was a storm that failed to intensify over the few days prior to making landfall on the Yucatan Peninsula. An observational study of Tropical Storm Chantal is presented using a diverse dataset including remote and in situ measurements from the NASA ER-2 and DC-8 and the NOAA WP-3D N42RF aircraft and satellite. The authors discuss the storm structure from the larger-scale environment down to the convective scale. Large vertical shear (850-200-hPa shear magnitude range 8-15 m/s) plays a very important role in preventing Chantal from intensifying. The storm had a poorly defined vortex that only extended up to 5-6-km altitude, and an adjacent intense convective region that comprised a mesoscale convective system (MCS). The entire low-level circulation center was in the rain-free western side of the storm, about 80 km to the west-southwest of the MCS. The MCS appears to have been primarily the result of intense convergence between large-scale, low-level easterly flow with embedded downdrafts, and the cyclonic vortex flow. The individual cells in the MCS such as cell 2 during the period of the observations were extremely intense, with reflectivity core diameters of 10 km and peak updrafts exceeding 20 m/s. Associated with this MCS were two broad subsidence (warm) regions, both of which had portions over the vortex. The first layer near 700 hPa was directly above the vortex and covered most of it. The second layer near 500 hPa was along the forward and right flanks of cell 2 and undercut the anvil divergence region above. There was not much resemblance of these subsidence layers to typical upper-level warm cores in hurricanes that are necessary to support strong surface winds and a low central pressure. The observations are compared to previous studies of weakly sheared storms and modeling studies of shear effects and intensification. The configuration of the convective updrafts, low-level circulation, and lack of vertical coherence between the upper- and lower-level warming regions likely inhibited intensification of Chantal. This configuration is consistent with modeled vortices in sheared environments, which suggest the strongest convection and rain in the downshear left quadrant of the storm, and subsidence in the upshear right quadrant. The vertical shear profile is, however, different from what was assumed in previous modeling in that the winds are strongest in the lowest levels and the deep tropospheric vertical shear is on the order of 10-12 m/s.

Description: Ozone measurements from the OMI and MLS instruments on board the Aura satellite are used for deriving global distributions of tropospheric column ozone (TCO). TCO is determined using the tropospheric ozone residual method which involves subtracting measurements of MLS stratospheric column ozone (SCO) from OMI total column ozone after adjusting for intercalibration differences of the two instruments using the convective-cloud differential method. The derived TCO field, which covers one complete year of mostly continuous daily measurements from late August 2004 through August 2005, is used for studying the regional and global pollution on a timescale of a few days to months. The seasonal and zonal characteristics of the observed TCO fields are also compared with TCO fields derived from the Global Modeling Initiative's Chemical Transport Model. The model and observations show interesting similarities with respect to zonal and seasonal variations. However, there are notable differences, particularly over the vast region of the Saharan desert.

Description: We have developed an algorithm to retrieve scattering cloud pressures and other cloud properties with the Aura Ozone Monitoring Instrument (OMI). The scattering cloud pressure is retrieved using the effects of rotational Raman scattering (RRS). It is defined as the pressure of a Lambertian surface that would produce the observed amount of RRS consistent with the derived reflectivity of that surface. The independent pixel approximation is used in conjunction with the Lambertian-equivalent reflectivity model to provide an effective radiative cloud fraction and scattering pressure in the presence of broken or thin cloud. The derived cloud pressures will enable accurate retrievals of trace gas mixing ratios, including ozone, in the troposphere within and above clouds. We describe details of the algorithm that will be used for the first release of these products. We compare our scattering cloud pressures with cloud-top pressures and other cloud properties from the Aqua Moderate-Resolution Imaging Spectroradiometer (MODIS) instrument. OMI and MODIS are part of the so-called A-train satellites flying in formation within 30 min of each other. Differences between OMI and MODIS are expected because the MODIS observations in the thermal infrared are more sensitive to the cloud top whereas the backscattered photons in the ultraviolet can penetrate deeper into clouds. Radiative transfer calculations are consistent with the observed differences. The OMI cloud pressures are shown to be correlated with the cirrus reflectance. This relationship indicates that OMI can probe through thin or moderately thick cirrus to lower lying water clouds.

Description: A four year record of MODIS spaceborne data provides a new measurement tool to assess the aerosol direct radiative effect at the top of the atmosphere. MODIS derives the aerosol optical thickness and microphysical properties from the scattered sunlight at 0.55-2.1 microns. The monthly MODIS data used here are accumulated measurements across a wide range of view and scattering angles and represent the aerosol s spectrally resolved angular properties. We use these data consistently to compute with estimated accuracy of +/-0.6W/sq m the reflected sunlight by the aerosol over global oceans in cloud free conditions. The MODIS high spatial resolution (0.5 km) allows observation of the aerosol impact between clouds that can be missed by other sensors with larger footprints. We found that over the clear-sky global ocean the aerosol reflected 5.3+/-0.6W/sq m with an average radiative efficiency of 49+/-2W/sq m per unit optical thickness. The seasonal and regional distribution of the aerosol radiative effects are discussed. The analysis adds a new measurement perspective to a climate change problem dominated so far by models.

Description: An objective assessment of the impact of a new cloud scheme, called Microphysics of Clouds with Relaxed Arakawa-Schubert Scheme (McRAS) (together with its radiation modules), on the finite volume general circulation model (fvGCM) was made with a set of ensemble forecasts that invoke performance evaluation over both weather and climate timescales. The performance of McRAS (and its radiation modules) was compared with that of the National Center for Atmospheric Research Community Climate Model (NCAR CCM3) cloud scheme (with its NCAR physics radiation). We specifically chose the boreal summer months of May and June 2003, which were characterized by an anomalously wet eastern half of the continental United States as well as northern regions of Amazonia. The evaluation employed an ensemble of 70 daily 10-day forecasts covering the 61 days of the study period. Each forecast was started from the analyzed initial state of the atmosphere and spun-up soil moisture from the first-day forecasts with the model. Monthly statistics of these forecasts with up to 10-day lead time provided a robust estimate of the behavior of the simulated monthly rainfall anomalies. Patterns of simulated versus observed rainfall, 500-hPa heights, and top-of-the-atmosphere net radiation were recast into regional anomaly correlations. The correlations were compared among the simulations with each of the schemes. The results show that fvGCM with McRAS and its radiation package performed discernibly better than the original fvGCM with CCM3 cloud physics plus its radiation package. The McRAS cloud scheme also showed a reasonably positive response to the observed sea surface temperature on mean monthly rainfall fields at different time leads. This analysis represents a method for helpful systematic evaluation prior to selection of a new scheme in a global model.

Description: There are several dozen papers that study the effects of cloud horizontal inhomogeneity on the retrievals of cloud optical thickness, but only a few of them deal with cloud droplet sizes. This paper is one of the first comprehensive attempts to fill this gap: It takes a close theoretical look at the radiative effects of cloud 3-D structure in retrievals of droplet effective radii. Under some general assumptions, it was found that ignoring subpixel (unresolved) variability produces a negative bias in the retrieved effective radius, while ignoring cloud inhomogeneity at scales larger than a pixel scale (resolved variability), on the contrary, leads to overestimation of the domain average droplet size. The theoretical results are illustrated with examples from Large Eddy Simulations (LES) of cumulus (Cu) and stratocumulus (Sc) cloud fields. The analysis of cloud drop size distributions retrieved from both LES fields confirms that ignoring shadowing in 1-D retrievals results in substantial overestimation of effective radii which is more pronounced for broken Cu than for Sc clouds. Collocated measurements of broken Cu clouds by Moderate Resolution Imaging Spectrometer (MODIS) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) are used to check simulations and theory with observations. The analysis of ASTER and MODIS data and associated derived products recommends against blindly using retrieved effective radii for broken cloud fields, especially if one wants to relate aerosol amounts to cloud droplet sizes.

Description: Two 20-day, continental midlatitude cases are simulated with a three-dimensional (3D) cloud-resolving model (CRM) and compared to Atmospheric Radiation Measurement (ARM) data. This evaluation of long-term cloud-resolving model simulations focuses on the evaluation of clouds and surface fluxes. All numerical experiments, as compared to observations, simulate surface precipitation well but over-predict clouds, especially in the upper troposphere. The sensitivity of cloud properties to dimensionality and other factors is studied to isolate the origins of the over prediction of clouds. Due to the difference in buoyancy damping between 2D and 3D models, surface precipitation fluctuates rapidly with time, and spurious dehumidification occurs near the tropopause in the 2D CRM. Surface fluxes from a land data assimilation system are compared with ARM observations. They are used in place of the ARM surface fluxes to test the sensitivity of simulated clouds to surface fluxes. Summertime simulations show that surface fluxes from the assimilation system bring about a better simulation of diurnal cloud variation in the lower troposphere.

Description: There is an apparent inconsistency between the total column ozone derived from the total ozone mapping spectrometer (TOMS) and aircraft observations within the eye region of tropical cyclones. The higher spectral resolution, coverage, and sampling of the ozone monitoring instrument (OMI) on NASA s Aura satellite as compared with TOMS allows for improved ozone retrievals by including estimates of cloud pressure derived simultaneously using the effects of rotational Raman scattering. The retrieved cloud pressures from OM1 are more appropriate than the climatological cloud-top pressures based on infrared measurements used in the TOMS and initial OM1 algorithms. We find that total ozone within the eye of hurricane Katrina is significantly overestimated when we use climatological cloud pressures. Using OMI-retrieved cloud pressures, total ozone in the eye is similar to that in the surrounding area. The corrected total ozone is in better agreement with aircraft measurements that imply relatively small or negligible amounts of stratospheric intrusion into the eye region of tropical cyclones.

Description: In mid-March 2005 the northern lower stratospheric polar vortex experienced a severe stretching episode, bringing a large polar filament far south of Alaska toward Hawaii. This meridional intrusion of rare extent, coinciding with the polar vortex final warming and breakdown, was followed by a zonal stretching in the wake of the easterly propagating subtropical main flow. This caused polar air to remain over Hawaii for several days before diluting into the subtropics. After being successfully forecasted to pass over Hawaii by the high-resolution potential vorticity advection model Modele Isentrope du transport Meso-echelle de l'Ozone Stratospherique par Advection (MIMOSA), the filament was observed on isentropic surfaces between 415 K and 455 K (17-20 km) by the Jet Propulsion Laboratory stratospheric ozone lidar measurements at Mauna Loa Observatory, Hawaii, between 16 and 19 March 2005. It was materialized as a thin layer of enhanced ozone peaking at 1.6 ppmv in a region where the climatological values usually average 1.0 ppmv. These values were compared to those obtained by the three dimensional Chemistry-Transport Model MIMOSA-CHIM. Agreement between lidar and model was excellent, particularly in the similar appearance of the ozone peak near 435 K (18.5 km) on 16 March, and the persistence of this layer at higher isentropic levels for the following three days. Passive ozone, also modeled by MIMOSA-CHIM, was at about 3-4 ppmv inside the filament while above Hawaii. A detailed history of the modeled chemistry inside the filament suggests that the air mass was still polar ozone- depleted when passing over Hawaii. The filament quickly separated from the main vortex after its Hawaiian overpass. It never reconnected and, in less than 10 days, dispersed entirely in the subtropics.

Description: Clouds exert an important influence on tropospheric photochemistry through modification of solar radiation that determines photolysis frequencies (J-values). We assess the radiative effect of clouds on photolysis frequencies and key oxidants in the troposphere with a global three-dimensional (3-D) chemical transport model (GEOS-CHEM) driven by assimilated meteorological observations from the Goddard Earth Observing System data assimilation system (GEOS DAS) at the NASA Global Modeling and Assimilation Office (GMAO). We focus on the year of 2001 with the GEOS-3 meteorological observations. Photolysis frequencies are calculated using the Fast-J radiative transfer algorithm. The GEOS-3 global cloud optical depth and cloud fraction are evaluated and generally consistent with the satellite retrieval products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the International Satellite Cloud Climatology Project (ISCCP). Results using the linear assumption, which assumes linear scaling of cloud optical depth with cloud fraction in a grid box, show global mean OH concentrations generally increase by less than 6% because of the radiative effect of clouds. The OH distribution shows much larger changes (with maximum decrease of approx.20% near the surface), reflecting the opposite effects of enhanced (weakened) photochemistry above (below) clouds. The global mean photolysis frequencies for J[O1D] and J[NO2] in the troposphere change by less than 5% because of clouds; global mean O3 concentrations in the troposphere increase by less than 5%. This study shows tropical upper tropospheric O3 to be less sensitive to the radiative effect of clouds than previously reported (approx.5% versus approx.20-30%). These results emphasize that the dominant effect of clouds is to influence the vertical redistribution of the intensity of photochemical activity while global average effects remain modest, again contrasting with previous studies. Differing vertical distributions of clouds may explain part, but not the majority, of these discrepancies between models. Using an approximate random overlap or a maximum-random overlap scheme to take account of the effect of cloud overlap in the vertical reduces the impact of clouds on photochemistry but does not significantly change our results with respect to the modest global average effect.

Description: A high-resolution numerical simulation of Hurricane Erin (2001) is used to examine the organization of vertical motion in the eyewall and how that organization responds to a large and rapid increase in the environmental vertical wind shear and subsequent decrease in shear. During the early intensification period, prior to the onset of significant shear, the upward motion in the eyewall was concentrated in small-scale convective updrafts that formed in association with regions of concentrated vorticity (herein termed mesovortices) with no preferred formation region in the eyewall. Asymmetric flow within the eye was weak. As the shear increased, an azimuthal wavenumber 1 asymmetry in storm structure developed with updrafts tending to form on the downshear to downshear-left side of the eyewall. Continued intensification of the shear led to increasing wavenumber 1 asymmetry, large vortex tilt, and a change in eyewall structure and vertical motion organization. During this time, the eyewall structure was dominated by a vortex couplet with a cyclonic (anticyclonic) vortex on the downtilt-left (downtilt-right) side of the eyewall and strong asymmetric flow across the eye that led to strong mixing of eyewall vorticity into the eye. Upward motion was concentrated over an azimuthally broader region on the downtilt side of the eyewall, upstream of the cyclonic vortex, where low-level environmental inflow converged with the asymmetric outflow from the eye. As the shear diminished, the vortex tilt and wavenumber 1 asymmetry decreased, while the organization of updrafts trended back toward that seen during the weak shear period.

Description: Snow depth on sea ice plays a critical role in the heat exchange between ocean and atmosphere because of its thermal insulation property. Furthermore, a heavy snow load on the relatively thin Southern Ocean sea-ice cover submerges the ice floes below sea level, causing snow-to-ice conversion. Snowfall is also an important freshwater source into the weakly stratified ocean. Snow depth on sea-ice information can be used as an indirect measure of solid precipitation. Satellite passive microwave data are used to investigate the interannual and regional variability of the snow cover on sea ice. In this study we make use of 12 years (1992-2003) of Special Sensor Microwave/Imager (SSM/I) radiances to calculate average monthly snow depth on the Antarctic sea-ice cover. The results show a slight increase in snow depth and a partial eastward propagation of maximum snow depths, which may be related to the Antarctic Circumpolar Wave.

Description: The terminology "convectively-coupled Kelvin waves" has been used frequently in the literature to refer to the 15 m/s eastward-moving planetary and large-scale waves in the tropics. This note points out that this terminology is not appropriate, since these waves contain Rossby waves and mixed Rossby-gravity waves also. The significance of pointing out this misnomer is that a better understanding of these waves may contribute to the search for their cause.

Description: The sensitivity of Arctic ozone loss to polar stratospheric cloud volume (V(sub PSC)) and chlorine and bromine loading is explored using chemistry and transport models (CTMs). A simulation using multi-decadal output from a general circulation model (GCM) in the Goddard Space Flight Center (GSFC) CTM complements one recycling a single year s GCM output in the Global Modeling Initiative (GMI) CTM. Winter polar ozone loss in the GSFC CTM depends on equivalent effective stratospheric chlorine (EESC) and polar vortex characteristics (temperatures, descent, isolation, polar stratospheric cloud amount). Polar ozone loss in the GMI CTM depends only on changes in EESC as the dynamics repeat annually. The GSFC CTM simulation reproduces a linear relationship between ozone loss and Vpsc derived from observations for 1992 - 2003 which holds for EESC within approx.85% of its maximum (approx.1990 - 2020). The GMI simulation shows that ozone loss varies linearly with EESC for constant, high V(sub PSC).

Description: Global microwave rainfall retrievals from a 5-satellite constellation, including TMI from TRMM, SSWI from DMSP F13, F14 and F15, and AMSR-E from EOS-AQUA, are assimilated into the NASA Goddard Earth Observing System (GEOS) Data Assimilation System (DAS) using a 1-D variational continuous assimilation (VCA) algorithm. The physical and dynamical impact of rainfall assimilation on GEOS analyses and forecasts is examined at various temporal and spatial scales. This study demonstrates that the 1-D VCA algorithm, which was originally developed and evaluated for rainfall assimilations over tropical oceans, can effectively assimilate satellite microwave rainfall retrievals and improve GEOS analyses over both the Tropics and the extratropics where the atmospheric processes are dominated by different large-scale dynamics and moist physics, and also over the land, where rainfall estimates from passive microwave radiometers are believed to be less accurate. Results show that rainfall assimilation renders the GEOS analysis physically and dynamically more consistent with the observed precipitation at the monthly-mean and 6-hour time scales. Over regions where the model precipitation tends to misbehave in distinctly different rainy regimes, the 1-D VCA algorithm, by compensating for errors in the model s moist time-tendency in a 6-h analysis window, is able to bring the rainfall analysis closer to the observed. The radiation and cloud fields also tend to be in better agreement with independent satellite observations in the rainfall-assimilation m especially over regions where rainfall analyses indicate large improvements. Assimilation experiments with and without rainfall data for a midlatitude frontal system clearly indicates that the GEOS analysis is improved through changes in the thermodynamic and dynamic fields that respond to the rainfall assimilation. The synoptic structures of temperature, moisture, winds, divergence, and vertical motion, as well as vorticity are more realistically captured across the front. Short-term forecasts using initial conditions assimilated with rainfall data also show slight improvements. 1

Description: The Antarctic ozone hole demonstrates large-scale, man-made affects on our atmosphere. Surface observations now show that human produced ozone depleting substances (ODSs) are declining. The ozone hole should soon start to diminish because of this decline. Herein we demonstrate an ozone hole parametric model. This model is based upon: 1) a new algorithm for estimating C1 and Br levels over Antarctica and 2) late-spring Antarctic stratospheric temperatures. This parametric model explains 95% of the ozone hole area s variance. We use future ODS levels to predict ozone hole recovery. Full recovery to 1980 levels will occur in approximately 2068. The ozone hole area will very slowly decline over the next 2 decades. Detection of a statistically significant decrease of area will not occur until approximately 2024. We further show that nominal Antarctic stratospheric greenhouse gas forced temperature change should have a small impact on the ozone hole.

Description: Global and large regional rainfall variations and possible long-term changes are examined using the 26-year (1979-2004) GPCP monthly dataset (Adler et al., 2003). Our emphasis is to discriminate among variations due to ENSO, volcanic events, and possible long-term climate changes in the tropics. Although the global linear change of precipitation in the data set is near zero during the time period, an increase in tropical rainfall is noted, with a weaker decrease over northern hemisphere middle latitudes. Focusing on the tropics (25degS-25degN), the data set indicates an upward trend (0.06 mm/day/decade) and a downward trend (-0.02 mm/day/decade) over tropical ocean and land, respectively. This corresponds to an about 4.9% increase (ocean) and 1.6% decrease (land) during the entire 26-year time period. Techniques are applied to isolate and quantify variations due to ENSO and two major volcanic eruptions (El Chichon, March 1982; Pinatubo, June 1991) in order to examine longer time-scale changes. The ENSO events generally do not impact the tropical total rainfall, but, of course, induce significant anomalies with opposite signs over tropical land and ocean. The impact of the two volcanic eruptions is estimated to be about a 5% reduction in tropical rainfall over both land and ocean. A modified data set (with ENSO and volcano effects removed) retains the same approximate linear change slopes, but with reduced variance, thereby increasing the confidence levels associated with the long-term rainfall changes in the tropics 2

Description: The NPOESS Preparatory Project (NPP) provides remotely-sensed land, ocean, atmospheric, ozone, and sounder data that will serve the meteorological and global climate change scientific communities while also providing risk reduction for the National Polar-orbiting Operational Environmental Satellite System (NPOESS), the U.S. Government s future low-Earth orbiting satellite system monitoring global weather and environmental conditions. NPOESS and NPP are a new era, not only because the sensors will provide unprecedented quality and volume of data but also because it is a joint mission of three federal agencies, NASA, NOAA, and DoD. NASA's primary science role in NPP is to independently assess the quality of the NPP science and environmental data records. Such assessment is critical for making NPOESS products the best that they can be for operational use and ultimately for climate studies. The Science Data Segment (SDS) supports science assessment by assuring the timely provision of NPP data to NASA s science teams organized by climate measurement themes. The SDS breaks down into nine major elements, an input element that receives data from the operational agencies and acts as a buffer, a calibration analysis element, five elements devoted to measurement based quality assessment, an element used to test algorithmic improvements, and an element that provides overall science direction. This paper will describe how the NPP SDS will leverage on NASA experience to provide a mission-reliable research capability for science assessment of NPP derived measurements.

Description: An intrusion of vortex edge air in D the interior of the Arctic polar vortex was observed on the January 31,2005 flight of the NASA DC-8 aircraft. This intrusion was identified as anomalously high values of ozone by the AROTAL and DIAL lidars. Our analysis shows that this intrusion formed when a blocking feature near Iceland collapsed, allowing edge air to sweep into the vortex interior. along the DC-8 flight track also shows the intrusion in both ozone and HNO3. Polar Stratospheric Clouds (PSCs) were observed by the DIAL lidar on the DC-8. The spatial variability of the PSCs can be explained using MLS HNO3 and H2O observations and meteorological analysis temperatures. We also estimate vortex denitrification using the relationship between N2O and HNO3. Reverse domain fill back trajectory calculations are used to focus on the features in the MLS data. The trajectory results improve the agreement between lidar measured ozone and MLS ozone and also improve the agreement between the HNO3 measurements PSC locations. The back trajectory calculations allow us to compute the local denitrification rate and reduction of HCl within the filament. We estimate a denitrification rate of about lO%/day after exposure to below PSC formation temperature. Analysis of Aura MLS observations made

Description: The TRMM Multi-satellite Precipitation Analysis (TMPA) provides a calibration-based sequential scheme for combining multiple precipitation estimates from satellites, as well as gauge analyses where feasible, at fine scales (0.25 degrees x 0.25 degrees and 3-hourly). It is available both after and in real time, based on calibration by the TRMM Combined Instrument and TRMM Microwave Imager precipitation products, respectively. Only the after-real-time product incorporates gauge data at the present. The data set covers the latitude band 50 degrees N-S for the period 1998 to the delayed present. Early validation results are as follows: The TMPA provides reasonable performance at monthly scales, although it is shown to have precipitation rate dependent low bias due to lack of sensitivity to low precipitation rates in one of the input products (based on AMSU-B). At finer scales the TMPA is successful at approximately reproducing the surface-observation-based histogram of precipitation, as well as reasonably detecting large daily events. The TMPA, however, has lower skill in correctly specifying moderate and light event amounts on short time intervals, in common with other fine-scale estimators. Examples are provided of a flood event and diurnal cycle determination.

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 (pathlength and width) and an estimation of the tornado intensity. This study explores the possibility of using near-real-time medium and high spatial resolution satellite imagery from the NASA Earth Observing System satellites to provide additional information for the surveys. Moderate Resolution Imaging Spectroradiometer (MODIS) and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) data were used to study the damage tracks from three tornadic storms: the La Plata, Maryland, storm of 28 April 2002 and the Ellsinore and Marquand, Missouri, storms of 24 April 2002. These storms varied in intensity and occurred over regions with significantly different land cover. It was found that, depending on the nature of the land cover, tornado damage tracks from intense storms (F1 or greater) and hail storms may be evident in ASTER, Landsat, and MODIS satellite imagery. In areas where the land cover is dominated by forests, the scar patterns can show up very clearly, while in areas of grassland and regions with few trees, scar patterns are not as obvious or cannot be seen at all in the satellite imagery. The detection of previously unidentified segments of a damage track caused by the 24 April 2002 Marquand, Missouri, tornado demonstrates the utility of satellite imagery for damage surveys. However, the capability to detect tornado tracks in satellite imagery depends on the ability to observe the ground without obstruction from space and appears to be as much dependent on the nature of the underlying surface and land cover as on the severity of the tornadic storm.

Description: When radiometers on a satellite are pointed towards the planet with the goal of understanding a phenomenon quantitatively, rather than just creating a pleasing image, the task at hand is often problematic. The signal at the detector can be affected by scattering, absorption, and emission; and these can be due to atmospheric constituents (gases, clouds, and aerosols), the earth's surface, and subsurface features. When targeting surface phenomena, the remote sensing algorithm needs to account for the radiation associated with the atmospheric constituents. Likewise, one needs to correct for the radiation leaving the surface, when atmospheric phenomena are of interest. Rigorous validation of such remote sensing products is a real challenge. In visible and near infrared wavelengths, the jumble of effects on atmospheric radiation are best accomplished over dark surfaces with fairly uniform reflective properties (spatial homogeneity) in the satellite instrument's field of view (FOV). The ocean's surface meets this criteria; land surfaces - which are brighter, more spatially inhomogeneous, and more changeable with time - generally do not. NASA's Clouds and the Earth's Radiant Energy System (CERES) project has used this backdrop to establish a radiation monitoring site in Virginia's coastal Atlantic Ocean. The project, called the CERES Ocean Validation Experiment (COVE), is located on a rigid ocean platform allowing the accurate measurement of radiation parameters that require precise leveling and pointing unavailable from ships or buoys. The COVE site is an optimal location for verifying radiative transfer models and remote sensing algorithms used in climate research; because of the platform's small size, there are no island wake effects; and suites of sensors can be simultaneously trained both on the sky and directly on ocean itself. This paper describes the site, the types of measurements made, multiple years of atmospheric and ocean surface radiation observations, and satellite validation results.

Description: We examine differences in total precipitable water vapor (PWV) from the Atmospheric Infrared Sounder (AIRS) and the Advanced Microwave Scanning Radiometer (AMSR-E) experiments sharing the Aqua spacecraft platform. Both systems provide estimates of PWV over water surfaces. We compare AIRS and AMSR-E PWV to constrain AIRS retrieval uncertainties as functions of AIRS retrieved infrared cloud fraction. PWV differences between the two instruments vary only weakly with infrared cloud fraction up to about 70%. Maps of AIRS-AMSR-E PWV differences vary with location and season. Observational biases, when both instruments observe identical scenes, are generally less than 5%. Exceptions are in cold air outbreaks where AIRS is biased moist by 10-20% or 10-60% (depending on retrieval processing) and at high latitudes in winter where AIRS is dry by 5-10%. Sampling biases, from different sampling characteristics of AIRS and AMSR-E, vary in sign and magnitude. AIRS sampling is dry by up to 30% in most high-latitude regions but moist by 5-15% in subtropical stratus cloud belts. Over the northwest Pacific, AIRS samples conditions more moist than AMSR-E by a much as 60%. We hypothesize that both wet and dry sampling biases are due to the effects of clouds on the AIRS retrieval methodology. The sign and magnitude of these biases depend upon the types of cloud present and on the relationship between clouds and PWV. These results for PWV imply that climatologies of height-resolved water vapor from AIRS must take into consideration local meteorological processes affecting AIRS sampling.

Description: Temperature versus pressure or T(p) time series from the Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research Satellite (UARS) have been extended and re-analyzed for the period of 1991-2005 and for the upper stratosphere and mesosphere in 10-degree wide latitude zones from 60S to 60N. Even though sampling from a solar occultation experiment is somewhat limited, it is shown to be quite adequate for developing both the seasonal and longer-term variations in T(p). Multiple linear regression (MLR) techniques were used in the re-analyses for the seasonal and the significant interannual, solar cycle (SC-like or decadal-scale), and linear trend terms. A simple SC-like term of 11-yr period was fitted to the time series residuals after accounting for the seasonal and interannual terms. Highly significant SC-like responses were found for both the upper mesosphere and the upper stratosphere. The phases of these SC-like terms were checked for their continuity with latitude and pressure-altitude, and in almost all cases they are directly in-phase with that of standard proxies for the solar flux variations. The analyzed, max minus min, responses at low latitudes are of order 1 K, while at middle latitudes they are as large as 3 K in the upper mesosphere. Highly significant, linear cooling trends were found at middle latitudes of the middle to upper mesosphere (about -2 K/decade), at tropical latitudes of the middle mesosphere (about -1 K/decade), and at 2 hPa (or order -1 K/decade).

Description: In the current version of the Numerical Spectral Model (NSM), the Quasi-biennial Oscillation (QBO) is generated primarily by small-scale gravity waves (GW) from Hines' Doppler Spread Parameterization (DSP). The model does not have topography, and the planetary waves are solely generated by instabilities. We discuss a 3D modeling study that describes the QBO extending from the stratosphere into the upper mesosphere, where the oscillation produces significant inter-annual variations in the diurnal tide. The numerical results are compared with temperature measurements from the SABER (TIMED) and MLS (UARS) instruments obtained by Huang et al. (2006). With a GW source that peaks at the Equator and is taken to be isotropic and independent of season, the NSM generates a QBO with variable periods around 26 months and zonal wind amplitudes of almost 25 m/s at 30 km. As reported earlier, the NSM reproduces the observed equinoctial maxima in the diurnal tide at altitudes around 95 km. The modeled QBO modulates the tide such that the seasonal amplitude maxima can vary from one year to another by as much as 30%. To shed light on the underlying mechanisms, the relative importance of the advection terms are discussed, and they are shown to be important in the stratosphere. At altitudes above 80 km, however, the QBO-related inter-annual variations of the tide are generated primarily by GW momentum deposition. In qualitative agreement with the SABER measurements, the model generates distinct zonal-mean QBO temperature variations in the stratosphere and mesosphere. In the stratosphere, the computed amplitudes are not much smaller than those observed, and the rate of downward propagation at the Equator is reproduced. The modeled temperature amplitudes in the mesosphere, however, are much smaller than those observed. The observed and computed temperature variations of the QBO peak at the Equator but extend with phase reversals to high latitudes, in contrast to the zonal winds that are confined to equatorial latitudes. Hemispherical asymmetries also appear in both the model results and the observations. The temperature amplitudes outside the equatorial region however tend to occur at lower latitudes in the model results. While there is qualitative agreement between the TIMED measurements and the model prediction, there are some areas of significant disagreement that require us to reexamine the present version of the NSM. The numerical results critically depend on the chosen parameters that determine the wave forcing, and there are a number of avenues to improve the performance of the model that had not been tuned to fit the observations. The GW spectrum and its latitude dependence in the troposphere are not well known, and numerical experiments are discussed that describe the related model response. While it appears that eastward propagating Kelvin waves and westward propagating Rossby gravity waves are not the primary source to generate the QBO, the GW forcing can seed the oscillation and act as a catalyst to enhance effectiveness of these planetary waves.

Description: One of the challenges facing atmospheric scientists is to interpret trends in multi-decadal data records. Although data records from satellite instruments are not as long as some ground-based records, global coverage and resolved vertical profiles provide unique information for identifying signatures of climate change. For example, the Halogen Occultation Experiment (HALOE) on the Upper Atmosphere Research Satellite provided profiles of O3, H2O, HC1, HF, CH4 from October 1991 until November 2005. There are also multi-annual ground based measurements of the column HCl. Middle latitude ground-based measurements show a seasonal cycle, and the HALOE profiles show that this is driven by the seasonal change in the composition and mass of the region between the tropopause and 380K surface (the lowermost stratosphere). Understanding the processes that produce the seasonal cycle makes it possible to interpret a future change in the seasonal cycle as a marker of a change in the stratospheric residual circulation. Satellite observations have also provided key information for improving the physical basis of models used to predict future composition and climate circulation. An example is the "tape recorder" signature in tropical stratospheric water vapor, i.e., the slow ascent of high and low water vapor anomalies roughly corresponding to the tropopause temperature at the time air entered the stratosphere. This signature has become a key diagnostic of performance for climate models.

Description: Seasonal snow and soil moisture (SM) are important elements of the Earth's water cycle, acting as storage elements and by modulating water and energy exchanges between the land surface and the atmosphere and space. Accurate estimates of snow water equivalent (SWE) and SM also play an important role in weather and climate models over a range of scales. There are many advantages for the retrieval of snow parameters using passive microwave sensors. However, the current scale mismatch between many high-resolution applications and available satellite snow sensor footprints (approx.25km) presents challenges. There is also a scale mismatch for planned L-band soil moisture satellite measurements with their large (50-100 km) passive microwave footprints. Combined radiometer-radar schemes are being explored as possible aids, but this is still an area of ongoing research. Retrieval accuracies are limited in part by complex, and potentially nonlinear relationships between SWE or SM and microwave emission signatures. We will use airborne and satellite observations from two field campaigns to examine the scaling behavior of brightness temperatures and SWE and SM across scales from 100 meters to 10's of kilometers. The results will provide examples of the appropriate spatial sampling scales for new sensors, and guidance for downscaling schemes. The two campaigns are NASA's Cold Land Processes Experiment- 1 (CLPX- 1 : Colorado, 2002 & 2003), and the Australian National Airborne Field Experiment (NAFE, Nov, 2005). The analyses will also illustrate the roles and spatial scales of the underlying phenomena (e.g., vegetation, topography) that control subpixel heterogeneity.

Description: More than three years of temperature observations from the SABER (TIMED) and MLS WARS) instruments are analyzed to study the annual and inter-annual variations extending from the stratosphere into the upper mesosphere. The SABER measurements provide data from a wide altitude range (15 to 95 km) for the years 2002 to 2004, while the MLS data were taken in the 16 to 55 km altitude range a decade earlier. Because of the sampling properties of SABER and MLS, the variations with local solar time must be accounted for when estimating the zonal mean variations. An algorithm is thus applied that delineates with Fourier analysis the year-long variations of the migrating tides and zonal mean component. The amplitude of the diurnal tide near the equator shows a strong semiannual periodicity with maxima near equinox, which vary from year to year to indicate the influence from the Quasi-biennial Oscillation (QBO) in the zonal circulation. The zonal mean QBO temperature variations are analyzed over a range of latitudes and altitudes, and the results are presented for latitudes from 48"s to 48"N. New results are obtained for the QBO, especially in the upper stratosphere and mesosphere, and at mid-latitudes. At Equatorial latitudes, the QBO amplitudes show local peaks, albeit small, that occur at different altitudes. From about 20 to 40 km, and within about 15" of the Equator, the amplitudes can approach 3S K for the stratospheric QBO or SQBO. For the mesospheric QBO or MQBO, we find peaks near 70 km, with temperature amplitudes reaching 3.5"K, and near 85 km, the amplitudes approach 2.5OK. Morphologically, the amplitude and phase variations derived from the SABER and MLS measurements are in qualitative agreement. The QBO amplitudes tend to peak at the Equator but then increase again pole-ward of about 15" to 20'. The phase progression with altitude varies more gradually at the Equator than at mid-latitudes. A comparison of the observations with results from the Numerical Spectral Model (NSM) reveals that there is qualitative agreement. The NSM generates the QBO extending from the stratosphere into the upper mesosphere, with temperature variations extending to mid latitudes, but the predicted amplitudes are smaller than those observed.

Description: The observed reduction in cloud-to-ground lightning in the new-equational zone is examined from the perspective of the width of the main negative charge region. Thermodynamic observations of cloud base height also show a climatological minimum value in the near-equatorial region. The association of low cloud base with both narrow updrafts and narrow changing zones may impede the bridging of the large air gap to ground, and thereby suppress cloud-to-ground lightning activity. This width dependence may be more important than the approx. 10% variation in height of the freezing level in the encouraging flashes to ground.

Description: The additional heating of the air over the city is the result of the replacement of naturally vegetated surfaces with those composed of asphalt, concrete, rooftops and other manmade materials. The temperatures of these artificial surfaces can be 20 to 40 0 C higher than vegetated surfaces. Materials such as asphalt store much of the sun s energy and remains hot long after sunset. This produces a dome of elevated air temperatures 5 to 8 C greater over the city, compared to the air temperatures over adjacent rural areas. This effect is called the "urban heat island". Urban landscapes are a complex mixture of vegetated and nonvegetated surfaces. It is difficult to take enough temperature measurements over a large city area to characterize the complexity of urban radiant surface temperature variability. However, the use of remotely sensed thermal data from airborne scanners are ideal for the task. The NASA Airborne Thermal and Land Applications Sensor (ATLAS) operates in the visual and IR bands was used in February 2004 to collect data from San Juan, Puerto Rico with the main objective of investigating the Urban Heat Island (UHI) in tropical cities. In this presentation we will examine the techniques of analyzing remotely sensed data for measuring the effect of various urban surfaces on their contribution to the urban heat island effect. Results from data collected from other US cities of Sacramento, Salt Lake City and Baton Rouge will be used to compare the "urban fabric" among the cities.

Description: During early July, 1999, the DROPPS (Distribution and Role of Particles in the Polar Summer Mesosphere) campaign launched two rocket payloads whose purpose was to study the polar summer MLT (mesosphere and lower thermosphere), particularly PMSEs (polar mesospheric summer echoes) and PMCs (polar mesospheric clouds). The rockets were launched from the Anderya Rocket Range in Norway the nights of the 5th and 14th of July. Both payloads included a front-mounted PID (Particle Impact Detector) consisting of charge and mass telescopes to measure aerosol and dust mass distributions. Ice particles of nanometer size are believed to be responsible for PMSEs through the process of electron scavenging. Evidence for this process is suggested, for example, by the presence of an electron "biteout" simultaneously measured by several instruments at an altitude of approx. 82 - 87km during the first DROPPS launch. This presentation will characterize similarities and differences between both flights as seen by the charge and mass telescopes, starting at launch until the loss of data on the downleg of each flight. Various stages of the flights will be considered in detail, such as the PMSE layer and the apogee at 117 km, as well as the calibration of the data before launch.

Description: Status and plans for GPCP are presented along with scientific findings from the current data set. Global and large regional rainfall variations and possible long-term changes are examined using the 26-year (1979-2004) monthly dataset from the Global Precipitation Climatology Project (GPCP). One emphasis is to discriminate among the variations due to ENSO, volcanic events and possible long-term changes. Although the global change of precipitation in the data set is near zero, the data set does indicate an upward trend (0.13 mm/day/25yr) and a downward trend (-0.06 mm/day/25yr) over tropical oceans and lands (25S-25N), respectively. This corresponds to a 4% increase (ocean) and 2% decrease (land) during this time period. Simple techniques are derived to attempt to eliminate variations due to ENSO and major volcanic eruptions in the Tropics. Using only annual values two "volcano years" are determined by examining ocean-land coupled variations in precipitation related to ENSO and other phenomena. The outlier years coincide with Pinatubo and El Chicon eruptions. The ENSO signal is reduced by deriving mean ocean and land values for El Nino, La Nina and neutral conditions based on Nino 3.4 SST and normalizing the annual ocean and land precipitation to the neutral set of cases. The impact of the two major volcanic eruptions over the past 25 years is estimated to be about a 5% reduction in tropical rainfall. The modified data set (with ENSO and volcano effect at least partially removed) retains the same approximate linear change slopes over the data set period, but with reduced variance leading to significance tests with results in the 90-95% range. Intercomparisons between the GPCP, SSM/I (1988-2004), and TRMM (1998-2004) satellite rainfall products and alternate gauge analyses over land are made to attempt to increase or decrease confidence in the changes seen in the GPCP analysis.

Description: Persistent and strong dependence of rain rate on the day of the week has been found in Tropical Rainfall Measuring Mission (TRMM) satellite estimates of summer afternoon rainfall over the southeast U.S. and the nearby Atlantic from 1998 to 2005. Midweek (Tue-Thu) rain rates and rain area appear to increase over land, and this increase is accompanied by a corresponding diminution of rainfall over nearby waters. Reanalysis data from atmospheric models, suggest that there is a corresponding weekly variation in atmospheric winds consistent with the changes in rainfall. These variations are almost certainly caused by weekly variations in human activity. The most likely cause of the observed changes in rainfall is the well documented weekly variation in atmospheric pollution. Particulate pollution is highest in the middle of the week. Considerable observational and modeling evidence has accumulated concerning the effects of aerosols on precipitation. Most of this evidence relates to the suppression of precipitation by aerosols, but it has been argued that storms in highly unstable moist environments can be invigorated by aerosols, and some modeling studies seem to confirm this. The strong weekly cycle in rainfall observed over the southeast U.S. along with what appears to be dynamical suppression of rainfall over the nearby Atlantic, and the lack of an observable cycle over the southwest U.S., are consistent with this theory.

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 CFWs. 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 (1 998 and 1999). 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), and (3) A discussion on the Goddard WRF version (its developments and applications).

Description: Parameterization of lightning NO(x) in global chemical transport models requires a method to specify the temporal and geographic distribution of flash rates, an estimate of NO production per flash, and a method to specify the effective vertical distribution of NO production. Flash rate parameterizations using the upward convective cloud mass flux have been developed and tested for several GCM and assimilated meteorological data sets. This parameterization and others involving the cloud-top height and convective precipitation provide a generally acceptable climatological distribution of flashes after normalization with satellite data. However, considerable improvement in terms of geographic and temporal variability may be achieved in the future using cloud microphysics data available from next-generation GCMs. Cloud-resolving case-study simulations of convective transport and lightning NO production have yielded results which are directly applicable to the NO production per flash and the vertical NO distribution portions of the global model lightning parameterizations. In this work we have used cloud-resolving models (the Goddard Cumulus Ensemble Model (GCE) and MM5) 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. 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 in each of several case study storms. 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. Estimates of mean NO production per flash vary by a factor of three from one simulated storm to another. 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. These profiles show a larger percentage of the total lightning NO(x) mass in the free troposphere than those suggested by Pickering et al. (1998).

Description: Spectral land surface albedo is an important parameter for describing the radiative properties of the Earth. Accordingly it reflects the consequences of natural and human interactions, such as anthropogenic, meteorological, and phenological effects, on global and local climatological trends. Consequently, albedos are integral parts in a variety of research areas, such as general circulation models (GCMs), energy balance studies, modeling of land use and land use change, and biophysical, oceanographic, and meteorological studies. , Over five years of land surface anisotropy, diffuse bihemispherical (white-sky) albedo and direct beam directional hemispherical (black-sky) albedo from observations acquired by the MODIS instruments aboard NASA s Terra and Aqua satellite platforms have provided researchers with unprecedented spatial, spectral, and temporal information on the land surface s radiative characteristics. However, roughly 30% of the global land surface, on an annual equal-angle basis, is obscured due to persistent and transient cloud cover, while another 207% is obscured due to ephemeral and seasonal snow effects. This precludes the MOD43B3 albedo products from being directly used in some remote sensing and ground-based applications, climate models, and global change research projects. To provide researchers with the requisite spatially complete global snow-free land surface albedo dataset, an ecosystem-dependent temporal interpolation technique was developed to fill missing or lower quality data and snow covered values from the official MOD43B3 dataset with geophysically realistic values. The method imposes pixel-level and local regional ecosystem-dependent phenological behavior onto retrieved pixel temporal data in such a way as to maintain pixel-level spatial and spectral detail and integrity. The phenological curves are derived from statistics based on the MODIS MOD12Q1 IGBP land cover classification product geolocated with the MOD43B3 data.

Description: In 2010, the National Aeronautics and Space Agency (NASA) of the U.S.A. and the Japanese Exploration Agency (JAXA) of Japan in cooperation with other U.S. and international partners will launch the Global Precipitation Measurement (GPM) mission. The mission center-piece is a core U.S. provided satellite holding a scanning microwave imager provided by the U.S. A. and a dual-frequency precipitation radar provided by Japan. The core satellite is in a 65 deg inclination (the current Tropical Rainfall Measuring Mission -TRMM is in a 35 deg inclination). Joining the core satellite are a constellation of approximately 8 satellites containing scanning radiometers. The purpose of the constellation is to increase the global sampling capability of the mission. One constellation satellite will be provided by NASA. Some of the constellation radiometers will provided by other U.S. agencies with existing (e.g. SSMI/S) or planned (NPOESS) radiometer assets. International groups have also expressed interest in contributing to the GPM mission including providing radiometer data for the GPM constellation. The use of a heterogeneous group of scanning radiometers each with its own unique purpose, characteristics and calibration offers a significant challenge for combining brightness temperatures or rain retrievals to create meaningful combined global radiometer products. However, the availability of active dual precipitation radar on the GPM core in combination with a well-calibrated radiometer on the same platform offers the possibility of inter-calibrating the constellation radiometers using the core satellite as a calibrator. This paper describes a joint NASA/GSFC and Colorado State University prototype effort at inter-calibrating existing radiometers using such a core calibrator approach. In the prototype, existing radiometers (i.e. SSM/I and AMSR-E) are intercalibrated, as required, using the TRMM Precipitation Radar (PR) and TRMM Microwave Imager (TMI) as the calibration core. This prototype includes four basic steps: provide the data in a normalized-common logical radiometer format, convert all data to a granule that stops and starts at the southernmost latitude rather than at the equator, perform basic processing and finally inter-calibrate the radiometers with the TRMM core as appropriate. Basic processing includes quality control, cross-track bias correction as well as the core action of converting T(sub a), to T(sub b). The purpose of the inter-calibration step is to provide consistent brightness temperatures for the various SSMA sensors, TMI and AMSR-E. The C appended to the level of the product indicates common brightness temperature (T(sub c)). The paper provides a summary of each of these actions. The details of the processing as well as data format and tools can be obtained by connecting to the website at mrain.atmos.colostate.eduLEVEL1C .

Description: The availability of climatological chlorophyll-a concentration data products from the SeaWiFS mission spanning the eight-year mission period allowed the creation of a climatological anomaly analysis function in Giovanni, the GES DISC Interactive Online Visualization and ANalysis Infrastructure. This study utilizes the Giovanni anomaly analysis function to examine mesoscale anomalies in the North Atlantic Ocean during the springtime North Atlantic Bloom. This examination indicates that areas exhibiting positive anomalies and areas exhibiting negative anomalies are coherent over significant spatial scales, with relatively abrupt boundaries between areas with positive and negative anomalies. Year-to-year variability in anomaly "intensity" can be caused by either variability in the temporal occurrence of the bloom peak or by variability in the peak chlorophyll concentration in a particular area. The study will also discuss the feasibility of combining chlorophyll anomaly analysis with other data types.

Description: The GCE (Goddard Cumulus Ensemble) model, which has been developed and improved at NASA Goddard Space Flight Center over the past two decades, is considered as one of the finer and state-of-the-art CRMs (Cloud Resolving Models) in the research community. As the chosen CRM for a NASA Interdisciplinary Science (IDS) Project, GCE has recently been successfully upgraded into an MPI (Message Passing Interface) version with which great improvement has been achieved in computational efficiency, scalability, and portability. By basically using the large-scale temperature and moisture advective forcing, as well as the temperature, water vapor and wind fields obtained from TRMM (Tropical Rainfall Measuring Mission) field experiments such as SCSMEX (South China Sea Monsoon Experiment) and KWAJEX (Kwajalein Experiment), our recent 2-D and 3-D GCE simulations were able to capture detailed convective systems typical of the targeted (simulated) regions. The GEOS-3 [Goddard EOS (Earth Observing System) Version-3] reanalysis data have also been proposed and successfully implemented for usage in the proposed/performed GCE long-term simulations (i.e., aiming at producing massive simulated cloud data -- Cloud Library) in compensating the scarcity of real field experimental data in both time and space (location). Preliminary 2-D or 3-D pilot results using GEOS-3 data have generally showed good qualitative agreement (yet some quantitative difference) with the respective numerical results using the SCSMEX observations. The first objective of this paper is to ensure the GEOS-3 data quality by comparing the model results obtained from several pairs of simulations using the real observations and GEOS-3 reanalysis data. The different large-scale advective forcing obtained from these two kinds of resources (i.e., sounding observations and GEOS-3 reanalysis) has been considered as a major critical factor in producing various model results. The second objective of this paper is therefore to investigate and present such an impact of large-scale forcing on various modeled quantities (such as hydrometeors, rainfall, and etc.). A third objective is to validate the overall GCE 3-D model performance by comparing the numerical results with sounding observations, as well as available satellite retrievals.

Description: Long recognized as a major environmental hazard, aerosol is now known to have strong impacts on both regional and global climate. It has been estimated that aerosol may reduce by up to 10% of the seasonal mean solar radiation reaching the earth surface, producing a global cooling effect that opposes global warming (Climate Change 2001). This means that the potential perils that humans have committed to global warming may be far greater than what we can detect at the present. As a key component of the Earth climate system, the water cycle is profoundly affected by the presence of aerosols in the atmosphere. Through the so-called direct effect , aerosol scatters and/or absorbs solar radiation, thus cooling the earth surface and changing the horizontal and vertical radiational heating contrast in the atmosphere. The heating contrast drives anomalous atmospheric circulation, resulting in changes in convection, clouds, and rainfall. Another way aerosol can affect the water cycle is through the so-called indirect effects, whereby aerosol increases the number of cloud condensation nuclei, prolongs life time of clouds, and inhibits the growth of cloud drops to raindrops. This leads to more clouds, and increased reflection of solar radiation, and further cooling at the earth surface. In monsoon regions, the response of the water cycle to aerosol forcing is especially complex, not only because of presence of diverse mix of aerosol species with vastly different radiative properties, but also because the monsoon is strongly influenced by ocean and land surface processes, land use, land change, as well as regional and global greenhouse warming effects. Thus, sorting out the impacts of aerosol forcing, and interaction with the monsoon water cycle is a very challenging problem. In this talk, I will offer some insights into how aerosols may impact the Asian monsoon based on preliminary results from satellite observations and climate model experiments. Specifically, I will discuss the elevated heat pump hypothesis, involving atmospheric heating by absorbing aerosols (dust and black carbon) over the southern slopes of the Himalayas, and feedback with the deep convection, in modifying monsoon water cycle over South .and East Asia. The role of aerosol forcing relative to those due to sea surface temperature and land surface processes, as well as observation requirements to verify such a hypothesis will also be discussed.

Description: Ozone is an atmospheric trace gas with multiple impacts on the environment. Global ozone fields are needed for air quality predictions, estimation of the ultraviolet radiation reaching the surface, climate-radiation studies, and may also have an impact on longer-term weather predictions. We estimate global ozone fields in the stratosphere and troposphere by combining the data from EOS Aura satellite with an ozone model using data assimilation. Ozone exhibits a large temporal variability in the lower stratosphere. Our previous work showed that assimilation of satellite data from limb-sounding geometry helps constrain ozone profiles in that region. We assimilated ozone data from the Aura Microwave Limb Sounder (MLS) and the Ozone Monitoring Instrument (OMI) into the ozone system at NASA's Global Modeling and Assimilation Office (GMAO). Ozone is transported within a general circulation model (GCM) which includes parameterizations for stratospheric photochemistry, tropospheric chemistry, and a simple scheme for heterogeneous ozone loss. The focus of this study is on the representation of ozone in the lower stratosphere and tropospheric ozone columns. We plan to extend studies of tropospheric ozone distribution through assimilation of ozone data from the Tropospheric Emission Spectrometer (TES). Comparisons with ozone sondes and occultation data show that assimilation of Aura data reproduces ozone gradients and variability in the lower stratosphere well. We proceed by separating the contributions to temporal changes in the ozone field into those that are due to the model and those that are due to the assimilation of Aura data. The impacts of Aura data are illustrated and their role in the representation of ozone variability in the lower stratosphere and troposphere is shown.

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. 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 explicitly 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), and (3) A discussion on the Goddard WRF version (its developments and applications).

Description: In writing the integral equations for the median mass diameter and particle concentration, or comparable parameters of the raindrop size distribution, it is apparent that when attenuation effects are included, the forms of the equations for polarimetric and dual wavelength radars are identical. In both sets of equations, differences in the backscattering and extinction cross sections appear: in the polarimetric equations, the differences are taken with respect polarization at a fixed frequency while for the dual wavelength equations, the differences are taken with respect to wavelength at a fixed polarization. Because the forms of the equations are the same, the ways in which they can be solved are similar as well. To avoid instabilities in the forward recursion procedure, the equations can be expressed in the form of a final-value. Solving the equations in this way traditionally has required estimates of the path attenuations to the final gate: either the attenuations at horizontal and vertical polarizations at the same frequency or attenuations at two frequencies with the same polarization. This has been done for dual-frequency (air/spaceborne case) and polarimetric radars by the respective use of the surface reference technique and the differential phase shift. An alternative to solving the constrained version of the equations is an iterative procedure recently proposed in which independent estimates of path attenuation are not required. Although the procedure has limitations, it appears to be quite useful. Simulations of the retrievals help clarify the relationship between the constrained and unconstrained approaches and their application to the polarimetric and dual-wavelength equations.

Description: In the presence of rain, 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. Conversely, the difference in radar reflectivity factors (in dB) between the upper and lower frequencies is independent of water vapor absorption and can be used to estimate the median mass diameter of the hydrometeors. For a down-looking radar, path-integrated estimates of water vapor absorption may be possible under rain-free as well as raining conditions by using the surface returns at the three frequencies. Cross-talk or interference between the precipitation and water vapor estimates depends on the frequency separation of the channels as well as on the phase state and the median mass diameter of the hydrometeors. Simulations of the retrieval of water vapor absorption show that the largest source of variability arises from the variance in the measured radar return powers while the largest biases occur in the mixed-phase region. Use of high pulse repetition frequencies and signal whitening methods may be needed to obtain the large number of independent samples required. Measurements over a fractional bandwidth, defined as the ratio of the difference between the upper and lower frequencies to the center frequency, up to about 0.2 should be passible in a differential frequency mode, where a single transceiver and antenna are used. Difficulties in frequency allocation may require alternative choices of frequency where the water vapor absorptions at the low and high frequencies are unequal. We consider the degradation in the retrieval accuracy when the frequencies are not optimum.

Description: Real clouds and cloud systems are inherently three-dimensional (3D). Because of the limitations in computer resources, however, most cloud-resolving models (CRMs) today are still two-dimensional (2D). A few 3D CRMs have been used to study the response of clouds to large-scale forcing. In these 3D simulations, the model domain was small, and the integration time was 6 hours. Only recently have 3D experiments been performed for multi-day periods for tropical cloud systems with large horizontal domains at the National Center for Atmospheric Research (NCAR), NOAA GFDL, the U.K. Met. Office, Colorado State University and NASA Goddard Space Flight Center. An improved 3D Goddard Cumulus Ensemble (GCE) model was recently used to simulate periods during TOGA COARE (December 19-27, 1992), GATE (september 1-7, 1974), SCSMEX (May 18-26, June 2-11, 1998) and KWAJEX (August 7-13, August 18-21, and August 29-September 12, 1999) using a 512 by 512 km domain and 41 vertical layers. The major objectives of this paper are: (1) to identify the differences and similarities in the simulated precipitation processes and their associated surface and water energy budgets in TOGA COARE, GATE, KWAJEX, and SCSMEX, and (2) to asses the impact of microphysics, radiation budget and surface fluxes on the organization of convection in tropics.

Description: Cloud physics is inevitably affected by the smoke particle (CCN, cloud condensation nuclei) size distribution below the clouds. Therefore, the size distributions parameterized as spectral bin microphysics are needed to explicitly study the effects 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 Goddard Cumulus Ensemble (GCE) model. The formulation for the explicit spectral-bin 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: The approach to create a comprehensive emission inventory for the time period 1980 to 2000 is described in this paper. We have recently compiled an emission database, which we will use for a 21 year simulation of tropospheric aerosols with the GOCART model. Particular attention was paid to the time-dependent SO2, black carbon and organic carbon aerosol emissions. For the emission of SO2 from sporadically erupting volcanoes, we assembled emission data from the Global Volcanism Program of the Smithsonian Institution, using the VEI to derive the volcanic cloud height and the SO2 amount, and amended this dataset by the SO2 emission data from the TOMS instrument when available. 3-dimensional aircraft emission data was obtained for a number of years from the AEAP project, converted from burned fuel to SO2 and interpolated to each year, taking the sparsity of the flight patterns into account. Other anthropogenic SO2 emissions are based on gridded emissions from the EDGAR 2000 database (excluding sources from aircraft, biomass burning and international ship traffic), which were scaled to individual years with country/regional based emission inventories. Gridded SO2 emissions from international ship traffic for 2000 and the scaling factors for other years are from [Eyring et al., 2005]. We used gridded anthropogenic black and organic carbon emissions for 1996 [Bond et al., 2005], again excluding aircraft, biomass burning and ship sources. These emissions were scaled with regional based emission inventories from 1980 to 2000 to derive gridded emissions for each year. The biomass burning emissions are based on a climatology, which is scaled with regional scaling factors derived from the TOMS aerosol index and the AVHRR/ASTR fire counts to each year [Duncan et al., 2003]. Details on the integration of the information from the various sources will be provided and the distribution patterns and total emissions in the final product will be discussed.

Description: Develop climatologies of gridded CG lightning densities and frequencies of occurrence for the Melbourne, FL National Weather Service (NWS MLB) county warning area. These grids are used to create a first-guess field for the lightning threat index map that is available on the NWS MLB NASA KSCIKT website. Forecasters previously created this map from scratch. Having the climatologies as a background field will increase consistency between forecasters and decrease their workload. Delivered all files containing the lightning climatologies, the data, and the code used to create the climatologies to NWS MLB. Completed and distributed a final memorandum describing how the climatologies were created. All the files were installed on the NWS MLB computer system, and then the code was compiled and tested to ensure that it worked properly on their operating system. The climatologies and their descriptions are posted on the NWS MLB website. Forecasting Low-Level Convergent Bands Under Southeast Flow Provide guidance to operational personnel that will help improve their forecasts of cloud bands under large-scale southeast flow. When these bands occur, they can lead to cloud, rain, and thunderstorm occurrences that adversely affect launch, landing, and ground operations at Kennedy Space Center/Cape Canaveral Air Force Station (KSC/CCAFS). Completed the first draft of the final report. The conclusions from this task indicated low-level wind speed and direction, low-level high pressure ridge position, east coast sea breeze front activity and upper-level jet streak position have the greatest influence on convergent band formation and movement during southeasterly flow.

Description: No abstract available

Description: This report summarizes the Applied Meteorology Unit (AMU) activities for the first quarter of Fiscal Year 2006 (October - December 2005). A detailed project schedule is included in the Appendix.

Description: Accurate vertical wind profiles are essential to successful launch or landing. Wind changes can make it impossible to fly a desired trajectory or avoid dangerous vehicle loads, possibly resulting in loss of mission. Balloons take an hour to generate a profile up to 20 km, but major wind changes can occur in 20 minutes. Wind profilers have the temporal response to detect such last minute hazards. They also measure the winds directly overhead while balloons blow downwind. At the Eastern Range (ER), altitudes from 2 to 20 km are sampled by a 50-MHz profiler every 4 minutes. The surface to 3 km is sampled by five 915-MHz profilers every 15 minutes. The Range Safety office assesses the risk of potential toxic chemical dispersion. They use observational data and model output to estimate the spatial extent and concentration of substances dispersed within the boundary layer. The ER uses 915-MHz profilers as both a real time observation system and as input to dispersion models. The WR has similar plans. Wind profilers support engineering analyses for the Space Shuttle. The 50-IVl11z profiler was used recently to analyze changes in the low frequency wind and low vertical wavenumber content of wind profiles in the 3 to 15 km region of the atmosphere. The 915-MHz profiler network was used to study temporal wind change within the boundary layer.

Description: The purpose of the work was to determine minimum atmospheric electric fields required for lightning initiation from an airborne vehicle at various altitudes up to 10 km. The problem was reduced to the determination of a condition for initiation of a viable positive leader from a conductive object in an ambient electric field. It was shown that, depending on air density and shape and dimensions of the object, critical atmospheric fields are governed by the condition for leader viability or that for corona onset. To establish quantitative criteria for reduced air densities, available observations of spark discharges in long laboratory gaps were analyzed, the effect of air density on leader velocity was discussed and evolution in time of the properties of plasma in the leader channel was numerically simulated. The results obtained were used to evaluate the effect of pressure on the quantitative relationships between the potential difference near the leader tip, leader current and its velocity; based on these relationships, criteria for steady development of a leader were determined for various air pressures. Atmospheric electric fields required for lightning initiation from rods and ellipsoidal objects of various dimensions were calculated at different air densities. It was shown that there is no simple way to extend critical ambient fields obtained for some given objects and pressures to other objects and pressures.

Description: A coordinated aircraft - radar project that investigated the electric fields, cloud microphysics and radar reflectivity of thunderstorm anvils near Kennedy Space Center is described. Measurements from two cases illustrate the extensive nature of the microphysics and electric field observations. As the aircraft flew from the edges of anvils into the interior, electric fields very frequently increased abruptly from approx.1 to 〉10 kV/m even though the particle concentrations and radar reflectivity increased smoothly. The abrupt increase in field usually occurred when the aircraft entered regions with a reflectivity of 10 to 15 dBZ. It is suggested that the abrupt increase in electric field may be because the charge advection from the storm core did not occur across the entire breadth of the anvil and was not constant in time. Screening layers were not detected near the edges of the anvils. Some long-lived anvils showed subsequent enhancement of electric field and reflectivity and growth of particles, which if localized, might be a factor in explaining the abrupt change of field in some cases. Comparisons of electric field magnitude with particle concentration or reflectivity for a combined data set that included all anvil measurements showed a threshold behavior. When the average reflectivity, such as in a 3-km cube, was less than approximately 5 dBZ, the electric field magnitude was 〈3 kV/m. Based on these findings, the Volume Averaged Height Integrated Radar Reflectivity (VAHIRR) is now being used by NASA, the Air Force and Federal Aviation Administration in new Lightning Launch Commit Criteria as a diagnostic for high electric fields in anvils.

Description: Forecasters at the Space Meteorology Group (SMG) issue 30 to 90 minute forecasts for low cloud ceilings at the Space Shuttle Landing Facility (TTS) to support Space Shuttle landings. Mission verification statistics have shown ceilings to be the number one forecast challenge for SMG. More specifically, forecasters at SMG are concerned with any rapidly developing clouds/ceilings below 8000 ft in a stable, capped thermodynamic environment. Therefore, the Applied Meteorology Unit (AMU) was tasked to examine archived events of rapid stable cloud formation resulting in ceilings below 8000 ft, and document the atmospheric regimes favoring this type of cloud development. The AMU examined the cool season months of November to March during the years of 1993-2003 for days that had low-level inversions and rapid, stable low cloud formation that resulted in ceilings violating the Space Shuttle Flight Rules. The AMU wrote and modified existing code to identify inversions from the morning (-10 UTC) Cape Canaveral, FL rawinsonde (XMR) during the cool season and output pertinent sounding information. They parsed all days with cloud ceilings below 8000 ft at TTS, forming a database of possible rapidly-developing low ceiling events. Days with precipitation or noticeable fog burn-off situations were excluded from the database. In the first phase of this work, only the daytime hours were examined for possible ceiling development events since low clouds are easier to diagnose with visible satellite imagery. Phase II of this work includes expanding the database to include nighttime cases which is underway as this abstract is being written. For the nighttime cases, the AMU will analyze both the 00 UTC soundings and the 10 UTC soundings to examine those data for the presence of a low-level inversion. The 00 UTC soundings will probably not have a surface-based inversion, but the presence of inversions or "neutral" layers aloft and below 8,000 ft will most likely help define the stable regime, being a thermodynamically "capped" environment. Occurrences of elevated low-level inversions or stable layers will be highlighted in conjunction with nights that experienced a possible development or onset of cloud ceilings below 8,000 ft. Using these criteria to narrow down the database, the AMU will then use archived IR satellite imagery for these possible events. This presentation summarizes the composite meteorological conditions for 20 daytime event days with rapid low cloud ceiling formation and 48 non-events days consisting of advection or widespread low cloud ceilings and describes two sample cases of daytime rapidly-developing low cloud ceilings. The authors will also summarize the work from the nighttime cases and describe a representative sample case from this data set.

Description: The Applied Meteorology Unit (AMU) located at the Kennedy Space Center (KSC)/Cape Canaveral Air Force Station (CCAFS) implemented an operational configuration of the Advanced Regional Prediction System (ARPS) Data Analysis System (ADAS), as well as the ARPS numerical weather prediction (NWP) model. Operational, high-resolution ADAS analyses have been produced from this configuration at the National Weather Service in Melbourne, FL (NWS MLB) and the Spaceflight Meteorology Group (SMG) over the past several years. Since that time, ADAS fields have become an integral part of forecast operations at both NWS MLB and SMG. To continue providing additional utility, the AMU has been tasked to implement visualization products to assess the potential for supercell thunderstorms and significant tornadoes, and to improve assessments of short-term cloud-to-ground (CG) lightning potential. This paper and presentation focuses on the visualization products developed by the AMU for the operational high-resolution ADAS and AR.PS at the NWS MLB and SMG. The two severe weather threat graphics implemented within ADAS/ARPS are the Supercell Composite Parameter (SCP) and Significant Tornado Parameter (SIP). The SCP was designed to identify areas with supercell thunderstorm potential through a combination of several instability and shear parameters. The SIP was designed to identify areas that favor supercells producing significant tornadoes (F2 or greater intensity) versus non-tornadic supercells. Both indices were developed by the NOAAINWS Storm Prediction Center (SPC) and were normalized by key threshold values based on previous studies. The indices apply only to discrete storms, not other convective modes. In a post-analysis mode, the AMU calculated SCP and SIP for graphical output using an ADAS configuration similar to the operational set-ups at NWS MLB and SMG. Graphical images from ADAS were generated every 15 minutes for 13 August 2004, the day that Hurricane Charley approached and made landfall on the Florida peninsula. Several tornadoes struck the interior of the Florida peninsula in advance of Hurricane Charley's landfall during the daylight hours of 13 August. Since SPC had previously examined this case using SCP and SIP graphics generated from output of the Rapid Update Cycle (RUC) model, this day served as a good benchmark to compare and validate the high-resolution ADAS graphics against the smoother RUC analyses, which serves as background fields to the ADAS analyses. The ADAS-generated SCP and STP graphics have been integrated into the suite of products examined operationally by NWS MLB forecasters and are used to provide additional guidance for assessment of the near-storm environment during convective situations.

Description: This study compared the performance of mechanical and ultrasonic anemometers at the Eastern Range (ER; Kennedy Space Center and Cape Canaveral Air Force Station on Florida's Atlantic coast) and the Western Range (WR; Vandenberg Air Force Base on California's Pacific coast). Launch Weather Officers, forecasters, and Range Safety analysts need to understand the performance of wind sensors at the ER and WR for weather warnings, watches, advisories, special ground processing operations, launch pad exposure forecasts, user Launch Commit Criteria (LCC) forecasts and evaluations, and toxic dispersion support. The current ER and WR weather tower wind instruments are being changed from the current propeller-and-vane (ER) and cup-and-vane (WR) sensors to ultrasonic sensors through the Range Standardization and Automation (RSA) program. The differences between mechanical and ultrasonic techniques have been found to cause differences in the statistics of peak wind speed in previous studies. The 45th Weather Squadron (45 WS) and the 30th Weather Squadron (30 WS) requested the Applied Meteorology Unit (AMU) to compare data between RSA and current sensors to determine if there are significant differences. Approximately 3 weeks of Legacy and RSA wind data from each range were used in the study, archived during May and June 2005. The ER data spanned the full diurnal cycle, while the WR data was confined to 1000-1600 local time. The sample of 1-minute data from numerous levels on 5 different towers on each range totaled more than 500,000 minutes of data (482,979 minutes of data after quality control). The 10 towers were instrumented at several levels, ranging from 12 ft to 492 ft above ground level. The RSA sensors were collocated at the same vertical levels as the present sensors and typically within 15 ft horizontally of each another. Data from a total of 53 RSA ultrasonic sensors, collocated with present sensors were compared. The 1-minute average wind speed/direction and the 1-second peak wind speed/direction were compared.

Description: The atmospheric angular momentum is closely related to variations in the Earth rotation. The atmospheric excitation function (AEF), or namely atmospheric effective angular momentum function, is introduced in studying the atmospheric excitation of the Earth's variable rotation. It may be separated into two portions, i.e, the "wind" terms due to the atmospheric motion relative to the mantle and the "pressure" terms due to the variations of atmospheric mass distribution evident through surface pressure changes. The AEF wind terms during the period of 1948-2004 are re-processed from the NCEP/NCAR (National Centers for Environmental Prediction-National Center for Atmospheric Research) reanalysis 6-hourly wind and pressure fields. Some previous calculations were approximate, in that the wind terms were integrated from an isobaric lower boundary of 1000 hPa. To consider the surface topography effect, however, the AEF is computed by integration using the winds from the Earth's surface to 10 hPa, the top atmospheric model level, instead of from 1000 hPa. For these two cases, only a minor difference, equivalent to approx. 0.004 milliseconds in length-of-day variation, exists with respect to the axial wind term. However, considerable differences, equivalent to 5-6 milliarcseconds in polar motion, are found regarding equatorial wind terms. We further compare the total equatorial AEF (with and without the topographic effect) with the polar motion excitation function (PMEF) during the period of 1980-2003. The equatorial AEF gets generally closer to the PMEF, and improved coherences are found between them when the topography effect is included. Keywords: Atmospheric angular momentum, Atmospheric excitation function, Earth rotation, Topography, Wind, Pressure.

Description: The Intercontinental Chemical Transport Experiment-North America (INTEX-NA) is an international field campaign envisioned to investigate the transport and transformation of gases and aerosols on transcontinental/intercontinental scales and assess their impact on air quality and climate. Phase B (INTEX-B) of the mission was conducted during a 10- week period from March 1 to May 15, 2006 and focused initially on pollution outflow from the Mexico City Metropolitan Area, later addressing the transport of pollution from Asia to North America during springtime meteorological conditions. During the deployment, fast-response (1-s resolution) CO2 measurements were recorded aboard the NASA DC-8 providing valuable regional-scale information on carbon sources and sinks over sparsely sampled areas of North America and adjacent ocean basins. When coupled with the enormously sophisticated chemistry payload on the DC-8, these measurements collectively afford extremely powerful multi-tracer constraints for carbon source/sink attribution. Preliminary examination of the two data sets from the INTEX-B campaign, acquired one month apart, reveals not only the influence of the CO2 seasonal cycle, but also the preponderance of human population and industrial activity in the northern hemisphere. In this presentation, a synergy of the ensemble of airborne and surface observations, bottomup emission inventories, as well as transport history are invoked in a GIS framework to elucidate the source/sink processes reflected in the observations. The airborne CO2 data, along with simultaneous surface measurements (e.g. NOAA ESRL), are examined to establish the vertical distribution and variability of CO2 as a function of location. The role of localized sources, long-range transport, the biosphere, stratospheric exchange, and dynamical processes on the CO2 spatial variability observed throughout the tropospheric column will be discussed.

Description: The CERES Flight Models 1 through 4 instruments were launched aboard NASA's Earth Observing System (EOS) Terra and Aqua Spacecraft into 705 Km sun-synchronous orbits with 10:30 a.m. and 1:30 p.m. equatorial crossing times. These instruments supplement measurements made by the CERES Proto Flight Model (PFM) instrument launched aboard NASA's Tropical Rainfall Measuring Mission (TRMM) into a 350 Km, 38-degree mid-inclined orbit. CERES Climate Data Records consist of geolocated and calibrated instantaneous filtered and unfiltered radiances through temporally and spatially averaged TOA, Surface and Atmospheric fluxes. CERES filtered radiance measurements cover three spectral bands including shortwave (0.3 to 5 microns), total (0.3 to 100 microns) and an atmospheric window channel (8 to 12 microns). The CERES Earth Radiation Budget measurements represent a new era in radiation climate data, realizing a factor of 2 to 4 improvement in calibration accuracy and stability over the previous ERBE climate records, while striving for the next goal of 0.3-percent per decade absolute stability. The current improvement is derived from two sources: the incorporation of lessons learned from the ERBE mission in the design of the CERES instruments and the development of a rigorous and comprehensive radiometric validation protocol consisting of individual studies covering different spatial, spectral and temporal time scales on data collected both pre and post launch. Once this ensemble of individual perspectives is collected and organized, a cohesive and highly rigorous picture of the overall end-to-end performance of the CERES instrument's and data processing algorithms may be clearly established. This approach has resulted in unprecedented levels of accuracy for radiation budget instruments and data products with calibration stability of better than 0.2-percent and calibration traceability from ground to flight of 0.25-percent. The current work summarizes the development, philosophy and implementation of the protocol designed to rigorously quantify the quality of the data products as well as the level of agreement between the CERES TRMM, Terra and Aqua climate data records.