ALBERT

Description: The vision document provides an overview of the Climate Change Science Program (CCSP) long-term strategic plan to enhance scientific understanding of global climate change.This document is a companion to the comprehensive Strategic Plan for the Climate Change Science Program. The report responds to the Presidents direction that climate change research activities be accelerated to provide the best possible scientific information to support public discussion and decisionmaking on climate-related issues.The plan also responds to Section 104 of the Global Change Research Act of 1990, which mandates the development and periodic updating of a long-term national global change research plan coordinated through the National Science and Technology Council.This is the first comprehensive update of a strategic plan for U.S. global change and climate change research since the origal plan for the U.S. Global Change Research Program was adopted at the inception of the program in 1989.

Description: The purpose of this Science and Implementation Plan is to describe GAPP science objectives and the activities required to meet these objectives, both specifically for the near-term and more generally for the longer-term. The GEWEX Americas Prediction Project (GAPP) is part of the Global Energy and Water Cycle Experiment (GEWEX) initiative that is aimed at observing, understanding and modeling the hydrological cycle and energy fluxes at various time and spatial scales. The mission of GAPP is to demonstrate skill in predicting changes in water resources over intraseasonal-to-interannual time scales, as an integral part of the climate system.

Description: Much of NASA's Arctic Research is run through its Cryospheric Sciences Program. Arctic research efforts to date have focused primarily on investigations of the mass balance of the largest Arctic land-ice masses and the mechanisms that control it, interactions among sea ice, polar oceans, and the polar atmosphere, atmospheric processes in the polar regions, energy exchanges in the Arctic. All of these efforts have been focused on characterizing, understanding, and predicting, changes in the Arctic. NASA's unique vantage from space provides an important perspective for the study of these large scale processes, while detailed process information is obtained through targeted in situ field and airborne campaigns and models. An overview of NASA investigations in the Arctic will be presented demonstrating how the synthesis of space-based technology, and these complementary components have advanced our understanding of physical processes in the Arctic.

Description: NASA conducted the Cirrus Regional Study of Tropical Anvils and Cirrus Layers (CRYSTAL) Florida Area Cirrus Experiment (FACE) during July 2002 for improved understanding of tropical cirrus. One of the goals was to improve the understanding of cirrus generation by convective updrafts. The reasons why some convective storms produce extensive cirrus anvils is only partially related to convective instability and the vertical transport ice mass by updrafts. Convective microphysics must also have an important role on cirrus generation, for example, there are hypotheses that homogeneous nucleation in convective updrafts is a major source of anvil ice particles. In this paper, we report on one intense CRYSTAL- FACE convective case on 16 July 2002 that produced extensive anvil.

Description: The capability for massively parallel programming (MPP) using a message passing interface (MPI) has been implemented into a three-dimensional version of the Goddard Cumulus Ensemble (GCE) model. The design for the MPP with MPI uses the concept of maintaining similar code structure between the whole domain as well as the portions after decomposition. Hence the model follows the same integration for single and multiple tasks (CPUs). Also, it provides for minimal changes to the original code, so it is easily modified and/or managed by the model developers and users who have little knowledge of MPP. The entire model domain could be sliced into one- or two-dimensional decomposition with a halo regime, which is overlaid on partial domains. The halo regime requires that no data be fetched across tasks during the computational stage, but it must be updated before the next computational stage through data exchange via MPI. For reproducible purposes, transposing data among tasks is required for spectral transform (Fast Fourier Transform, FFT), which is used in the anelastic version of the model for solving the pressure equation. The performance of the MPI-implemented codes (i.e., the compressible and anelastic versions) was tested on three different computing platforms. The major results are: 1) both versions have speedups of about 99% up to 256 tasks but not for 512 tasks; 2) the anelastic version has better speedup and efficiency because it requires more computations than that of the compressible version; 3) equal or approximately-equal numbers of slices between the x- and y- directions provide the fastest integration due to fewer data exchanges; and 4) one-dimensional slices in the x-direction result in the slowest integration due to the need for more memory relocation for computation.

Description: A detailed analysis of available in situ and remotely sensed N2O and CH4 data measured in the 1999/2000 winter Arctic vortex has been performed in order to quantify the temporal evolution of vortex descent. Differences in potential temperature (theta) among balloon and aircraft vertical profiles (an average of 19-23 K on a given N2O or CH4 isopleth) indicated significant vortex inhomogeneity in late fall as compared with late winter profiles. A composite fall vortex profile was constructed for 26 November 1999, whose error bars encompassed the observed variability. High-latitude extravortex profiles measured in different years and seasons revealed substantial variability in N2O and CH4 on theta surfaces, but all were clearly distinguishable from the first vortex profiles measured in late fall 1999. From these extravortex-vortex differences we inferred descent prior to 26 November: as much as 397 plus or minus 15 K (lsigma) at 30 ppbv N2O and 640 ppbv CH4, and falling to 28 plus or minus 13 K above 200 ppbv N2O and 1280 ppbv CH4. Changes in theta were determined on five N2O and CH4 isopleths from 26 November through 12 March, and descent rates were calculated on each N2O isopleth for several time intervals. The maximum descent rates were seen between 26 November and 27 January: 0.82 plus or minus 0.20 K/day averaged over 50- 250 ppbv N2O. By late winter (26 February to 12 March), the average rate had decreased to 0.10 plus or minus 0.25 K/day. Descent rates also decreased with increasing N2O; the winter average (26 November to 5 March) descent rate varied from 0.75 plus or minus 0.10 K/day at 50 ppbv to 0.40 plus or minus 0.11 K/day at 250 ppbv. Comparison of these results with observations and models of descent in prior years showed very good overall agreement. Two models of the 1999/2000 vortex descent, SLIMCAT and REPROBUS, despite theta offsets with respect to observed profiles of up to 20 K on most tracer isopleths, produced descent rates that agreed very favorably with the inferred rates from observation.

Description: The Goddard Cumulus Ensemble (GCE) model is used to examine the sensitivities of surface fluxes, explicit radiation, and ice microphysical processes on multi-day simulations of deep tropical convection over the Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE). The simulations incorporate large-scale advective temperature and moisture forcing, as well as large-scale momentum, that are updated every time step on a periodic lateral boundary grid. This study shows that when surface fluxes are eliminated, the mean atmosphere is much cooler and drier, convection and CAPE are much weaker, precipitation is less, and cloud coverage in stratiform regions much greater. Surface fluxes using the TOGA COARE flux algorithm are weaker than with the aerodynamic formulation, but closer to the observed fluxes. In addition, similar trends noted above for the case without surface fluxes are produced for the TOGA flux case, albeit to a much lesser extent. The elimination of explicit shortwave and longwave radiation is found to have only minimal effects on the mean thermodynamics, convection, and precipitation. However explicit radiation does have a significant impact on cloud temperatures and structure above 200 mb and on the overall mean vertical circulation. The removal of ice processes produces major changes in the structure of the cloud. Much of the liquid water is transported aloft and into anvils above the melting layer (600 mb), leaving narrow, but intense bands of rainfall in convective regions. The elimination of melting processes leads to greater hydrometeor mass below the melting layer, and produces a much warmer and moister boundary layer, leading to a greater mean CAPE. Finally, the elimination of the graupel species has only a small impact on mean total precipitation, thermodynamics, and dynamics of the simulation, but does produce much greater snow mass just above the melting layer. Some of these results differ from previous CRM studies of tropical systems, which is likely due to the type of simulated system, total time integration, and model setup.

Description: Observing system simulation experiments (OSSE's) provide an effective means to evaluate the potential impact of a proposed observing system, as well as to determine tradeoffs in their design, and to evaluate data assimilation methodology. Great care must be taken to ensure realism of the OSSE's, and in the interpretation of OSSE results. All of the OSSE's that have been conducted to date have demonstrated tremendous potential for space-based wind profile data to improve atmospheric analyses, forecasts, and research. This has been true for differing data assimilation systems, analysis methodology, and model resolutions. OSSE's clearly show much greater potential for observations of the complete wind profile than for single-level wind data or observations of the boundary layer alone.

Description: Spaceborne scatterometers are active microwave radar instruments designed to acquire near-simultaneous, spatially collocated measurements of the normalized radar backscattering cross section (sigma0) of the global surface from several azimuth and/or incidence angles. The primary objective of the scatterometer mission is to measure the near-surface wind speed and direction over the global ocean using sigma0 measurements together with a wind geophysical model function. However, since sigma0 measurements are collected globally all the time, sigma0 data can also be used for global land and ice applications. In this paper, we will first present the objectives of the QSCAT mission, the instrument design, and the unique features of the Ku-band scatterometer currently in operation, called SeaWinds on QuikSCAT (QSCAT). We will then present some emerging land and ocean applications of the QSCAT data, which include (1) global snow detection and monitoring, (2) melt region mapping on the Greenland ice sheet, (3) Monsoon flood detection and monitoring, (4) soil wetness application at large scale, and (5) hurricane monitoring and tracking.

Description: This paper studies the performance of a spaceborne precipitation radar in measuring vertical Doppler velocity of rainfall. As far as a downward pointing precipitation radar is concerned, one of the major problems affecting Doppler measurement at the nadir direction arises from the Non-Uniform Beam-Filling effect (NUBF). That is, when significant variation in rain rate is present within the radar IFOV (Instrument Field of View) in the along track direction. the Doppler shift caused by the radial component of the horizontal speed of the satellite is weighted differently among the portions of IFOV. The effects of this non-uniform weighting may dominate any other contribution. Under this condition, shape, average value and width of the Doppler spectrum may not be directly correlated with the vertical velocity of the precipitating particles. However, by using an inversion technique which over-samples the radar measurements in the along track direction, we show that the shift due to NUBF can be evaluated, and that the NUBF induced errors on average fall speed can be reduced.

Description: The conversion of cloud ice to snow by depositional growth, designated P(sub SFI), in the Goddard Cumulus Ensemble Model cloud physics parameterization is examined. The original formulation of P(sub SFI) is shown to produce excessive conversion of cloud ice to snow because of an implicit assumption that the relative humidity is 100% with respect to water even though the air may actually be quite less humid. Two possible corrections to this problem are proposed, the first involving application of a relative humidity dependent correction factor to the original formulation of P(sub SFI), and the second involving a new formulation of P(sub SFI) based on the equation for depositional growth of cloud ice.

Description: This project, supported by the NASA New Investigator Program, has primarily been funding a graduate student, Darren McKague. Since August 1999 Darren has been working part time at Raytheon, while continuing his PhD research. Darren is planning to finish his thesis work in May 2001, thus some of the work described here is ongoing. The proposed research was to use GOES visible and infrared imager data and SSM/I microwave data to obtain joint distributions of cirrus cloud ice mass and precipitation for a study region in the Eastern Tropical Pacific. These joint distributions of cirrus cloud and rainfall were to be compared to those from the CSU general circulation model to evaluate the cloud microphysical amd cumulus parameterizations in the GCM. Existing algorithms were to be used for the retrieval of cloud ice water path from GOES (Minnis) and rainfall from SSM/I (Wilheit). A theoretical study using radiative transfer models and realistic variations in cloud and precipitation profiles was to be used to estimate the retrieval errors. Due to the unavailability of the GOES satellite cloud retrieval algorithm from Dr. Minnis (a co-PI), there was a change in the approach and emphasis of the project. The new approach was to develop a completely new type of remote sensing algorithm - one to directly retrieve joint probability density functions (pdf's) of cloud properties from multi-dimensional histograms of satellite radiances. The usual approach is to retrieve individual pixels of variables (i.e. cloud optical depth), and then aggregate the information. Only statistical information is actually needed, however, and so a more direct method is desirable. We developed forward radiative transfer models for the SSM/I and GOES channels, originally for testing the retrieval algorithms. The visible and near infrared ice scattering information is obtained from geometric ray tracing of fractal ice crystals (Andreas Macke), while the mid-infrared and microwave scattering is computed with Mie scattering. The radiative transfer is performed with the Spherical Harmonic Discrete Ordinate Method (developed by the PI), and infrared molecular absorption is included with the correlated k-distribution method. The SHDOM radiances have been validated by comparison to version 2 of DISORT (the community "standard" discrete-ordinates radiative transfer model), however we use SHDOM since it is computationally more efficient.

Description: The microphysical parameterization of clouds and rain-cells plays a central role in atmospheric forward radiative transfer models used in calculating passive microwave brightness temperatures. The absorption and scattering properties of a hydrometeor-laden atmosphere are governed by particle phase, size distribution, aggregate density., shape, and dielectric constant. This study identifies the sensitivity of brightness temperatures with respect to the microphysical cloud parameterization. Cloud parameterizations for wideband (6-410 GHz observations of baseline brightness temperatures were studied for four evolutionary stages of an oceanic convective storm using a five-phase hydrometeor model in a planar-stratified scattering-based radiative transfer model. Five other microphysical cloud parameterizations were compared to the baseline calculations to evaluate brightness temperature sensitivity to gross changes in the hydrometeor size distributions and the ice-air-water ratios in the frozen or partly frozen phase. The comparison shows that, enlarging the rain drop size or adding water to the partly Frozen hydrometeor mix warms brightness temperatures by up to .55 K at 6 GHz. The cooling signature caused by ice scattering intensifies with increasing ice concentrations and at higher frequencies. An additional comparison to measured Convection and Moisture LA Experiment (CAMEX 3) brightness temperatures shows that in general all but, two parameterizations produce calculated T(sub B)'s that fall within the observed clear-air minima and maxima. The exceptions are for parameterizations that, enhance the scattering characteristics of frozen hydrometeors.

Description: A persistent, mesoscale region of intense eyewall convection contained within Hurricane Bonnie on 23 August 1998 is examined from multiple observations synthesized from the NASA ER-2 and DC-8 aircraft. The intense convection occurred late in the day as Bonnie was attaining its minimum central pressure and during a stage when the inner core featured a markedly asymmetric structure. The main purpose of this paper is to describe the internal structure of a convective burst and its relationship to the warm core using a synthesis of high-resolution satellite, aircraft radar, and in situ data. An exceptionally vigorous eyewall tower penetrating to nearly 18 km is described. A second intense eyewall tower, adjacent to the eye, is shown to be associated with a mesoscale subsiding current of air that extends horizontally nearly 20 km into the eye interior. The subsidence occupies a substantial depth within the eye and appears to be a much larger scale feature than the convectively-induced, symmetric overturning which commonly occurs on the upper-level flanks of convective towers in other tropical environments.

Description: Ensemble sets of simulation experiments were conducted with a single column model (SCM) using the Goddard GEOS II GCM physics containing a recent version of the Cumulus Scheme (McRAS) and a biosphere based land-fluxes scheme (SSiB). The study used the 18 July to 5 August 1995 ARM-CART (Atmospheric Radiation Measurement-Cloud Atmospheric Radiation Test-bed) data, which was collected at the ARM-CART site in the mid-western United States and analyzed for single column modeling (SCM) studies. The new findings affirm the earlier findings that the vegetation, which increases the solar energy absorption at the surface together with soil and soil-moisture dependent processes, which modulate the surface, fluxes (particularly evapotranspiration) together help to increase the local rainfall. In addition, the results also show that for the particular study period roughly 50% of the increased evaporation over the ARM-CART site would be converted into rainfall with the Column, while the remainder would be advected out to the large-scale. Notwithstanding the limitations of only one-way interaction (i.e., the large-scale influencing the regional physics and not vice versa), the current SCM simulations show a very robust relationship. The evaporation-precipitation relationship turns out to be independent of the soil types, and soil moisture; however, it is weakly dependent on the vegetation cover because of its surface-albedo effect. Clearly, these inferences are prone to weaknesses of the SCM physics, the assumptions of the large-scale being unaffected by gridscale (SCM-scale) changes in moist processes, and other limitations of the evaluation procedures.

Description: A response is provided to Diem and colleagues discussion on the whether urbanization-enhanced precipitation should be maximized in the South-Southeast of Atlanta, Georgia as documented in Shepherd et al. (2003). Diem and colleagues have offered a critical response to Shepherd et al. (2002). The response herein offers both general and specific responses to issues raised by Diem and colleagues.

Description: Measurements from the Halogen Occultation Experiment (HALOE) revealed the infrared signature of polar mesospheric clouds (PMCs), for the first time, HALOE PMC observations at eight wavelengths (2.45 - 10 microns) show remarkable agreement with model PMC spectra based on ice particle extinction, and thus provide the first confirmation that water ice is the primary component of PMCs. Because PMCs respond to changes in temperature and water vapor, they are considered an indicator of global climate change. We propose to further the understanding of PMCs using a decade of infrared measurements form HALOE. This effort will characterize PMC spectral properties, extinction profiles, and size distributions. Using this information, HALOE measurements will be used to make simultaneous retrievals of H2O3, and temperature, in the presence of PMCs. The simultaneous retrievals of particle properties, H2O3, and temperature will be used with HALOE NO data to provide a significant step forward in the knowledge of PMC characteristics and formation conditions. We will challenge fundamental theories of PMC formation, and investigate changes in PMC properties and related conditions over the length of the HALOE measurement record. HALOE has been operating without flaw since it was launched on October 11, 1991. Consequently, ten southern and ten northern PMC seasons have been observed thus far, providing a wealth of data for the study of PMC and related parameters.

Description: Heavy rainfall occurred over the western side of Taiwan's complex terrain from August 10 to 13, 1994 after Typhoon Doug moved northward from the East China Sea into Taiwan and on towards the Yellow Sea. On August 10, most of the rainfall fell over sloped areas. The heaviest daily rainfall totals were in excess of 200 mm over southwestern as well as central Taiwan. However, not much rainfall occurred over northern Taiwan. The lack of rainfall over northern Taiwan also occurred on August 11, 12 and 13. The larger rainfall amounts shifted westward from the sloped areas on August 10 toward lower terrain on August 11. On August 12 and 13, most of the higher rainfall amounts were found over the coastal area in southwestern Taiwan. Notably, about 300 to 400 mm per day fell over the coastal area in southwest Taiwan on August 12 and 13. The distribution of rainfall amount was different on August 10 and 11 (termed as Case 1) compared to August 12 and 13 (termed as Case 2). The environmental situation and precipitation characteristics are analyzed using EC/TOGA data, ground-based radar data, surface rainfall patterns, surface wind data, and upper air soundings. Chen at al. (2001) also categorized the precipitation pattern into two types, propagating and quasi-stationary. For the propagating type of precipitation, rainrates increased or remained the same as systems went from the plains to mountainous regions. With the quasi-stationary type of precipitation, however, rainrates decreased as precipitation propagated across the plains and into the mountains. The focus of this study is to understand what causes the h1aher amounts of rainfall over Taiwan, and what factors influence where the higher amounts of rainfall will occur, over sloped areas or over coastal areas.

Description: An observing system comprised of two lidars in geosynchronous orbit would enable the synoptic and meso-scale measurement of atmospheric winds and moisture, both of which are key first-order variables of the Earth's weather equation. Simultaneous measurement of these parameters at fast revisit rates promises large advancements in our weather prediction skills. Such capabilities would be unprecedented and yield greatly improved and finer resolution initial conditions for models, make existing costly and cumbersome measurement approaches obsolete, and obviate the use of numerical techniques needed to correct data obtained using present observing systems. Additionally, simultaneous synoptic wind and moisture observations would lead to improvements in model parameterizations, and in our knowledge of small-scale weather processes. Technology and science data product assessments are ongoing. Results will be presented during the conference.

Description: The development of a satellite infrared technique for estimating convective and stratiform rainfall and its application in studying the diurnal variability of rainfall in Amazonia are presented. The Convective-Stratiform. Technique, calibrated by coincident, physically retrieved rain rates from the Tropical Rain Measuring Mission (TRMM) Microwave Imager (TMI), is applied during January to April 1999 over northern South America. The diurnal cycle of rainfall, as well as the division between convective and stratiform rainfall is presented. Results compare well (a one-hour lag) with the diurnal cycle derived from Tropical Ocean-Global Atmosphere (TOGA) radar-estimated rainfall in Rondonia. The satellite estimates reveal that the convective rain constitutes, in the mean, 24% of the rain area while accounting for 67% of the rain volume. The effects of geography (rivers, lakes, coasts) and topography on the diurnal cycle of convection are examined. In particular, the Amazon River, downstream of Manaus, is shown to both enhance early morning rainfall and inhibit afternoon convection. Monthly estimates from this technique, dubbed CST/TMI, are verified over a dense rain gage network in the state of Ceara, in northeast Brazil. The CST/TMI showed a high bias equal to +33% of the gage mean, indicating that possibly the TMI estimates alone are also high. The root mean square difference (after removal of the bias) equaled 36.6% of the gage mean. The correlation coefficient was 0.77 based on 72 station-months.

Description: Global rainfall is the primary redistributor of earth's energy by the process of latent heat release. This forms the main driving force for the tropical circulation, which in turn impacts the global circulation .through transient events-such as El Nino. Hence, more precise and long-term time series of the rainfall and its variability is crucial to the understanding and prediction of the global climate and climate change. The Precipitation Radar (PR) abroad the US/Japan Tropical Rainfall Measuring Mission (TRMM) is the first radar ever launched into space that measures detailed vertical profiles of rain intensity over the tropics. One of the challenges in estimating rainfall from spaceborne radars is the presence of attenuation at frequencies, such as 14 GHz of the TRMM PR and future planned systems at this and higher frequencies. A common approach in current rainfall retrieval algorithms is to employed the path integrated attenuation (PIA) as a constraint to the retrieval, and hence overcome errors in the radar calibration or in the assumed rainfall parameters. PIA can either be derived from a radiometer or from the surface reference technique, in which a clear air radar measurement is compared with the measurement in the raining area. The current TRMM 2A21 PIA data product makes use of both a temporal and spatial clear air database for comparison to rainy measurements. In this paper we present results from analysis of TRMM surface backscatter cross-section (sigmaO) measurements from Nov 97-Feb99, and a comparison with sigmaO measurements obtained by the NASA Scatterometer (NSCAT) between Sept96-June97. Measurements for a given month from both instruments are compiled on a 1 deg. (lat.) x 1 deg. (lon.) x 1 hr. grid. This enables TRMM--NSCAT comparison and the investigation of seasonal and diurnal trends in both data sets. From preliminary analysis of TRMM sigmaO's we have decided not to treat the ocean as a single homogeneous region but to select a number of ocean sub-regions and individually analyze their trends. Likewise, and in a similar approach to previous studies of Seasat over-land data, we have selected a number of over-land regions for study. From said sigmaO maps and regional trend analysis we investigate possible sources of trends and variability. In addition, we study the effects of TRMM PR sensitivity through the PR "possible rain" class. Given NSCAT's inability to flag rain contaminated measurements we are able to gauge the impact of rain contamination on NSCAT monthly sigmaO maps, using TRMM measurements. The research described in this paper was carried out by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration (NASA), U.S.A.

Description: The NASA ER-2 and DC-8 aircraft collected remote sensing and in situ data sets from Hurricane Bonnie (23, 24, and 26 August 1998) during the Convection And Moisture Experimental-3 (CAMEX-3). Bonnie was an exceptional case where NASA and NOAA had five aircraft sampling both upper levels and lower altitudes. The ER-2 was instrumented with the ER-2 Doppler XBand radar (EDOP) and several radiometers ranging from visible to lower frequency microwaves. EDOP is a fixed dual-beam radar (nadir and forward-looking beams) which allows computation of both vertical and alongtrack horizontal winds. The hurricane secondary circulation is typically difficult to measure at upper levels due to aircraft altitude limitations and sensitivity of the lower altitude airborne radars. EDOP is in principle, well suited to measure these components of the wind. When ER-2 flies across the approximate center of the hurricane circulation, the along-track winds derived from EDOP, are approximately equal to the hurricane radial flow comprising the secondary circulation. Assuming that the hydrometeor fallspeeds can be approximated, the radial and vertical wind components of the secondary circulation can be measured. Since the hydrometeor motions can be estimated with more confidence in the higher altitude ice regions (i.e., graupel and mixed phase are complicated at lower altitudes), the derived radial and vertical winds have higher accuracy at upper levels. On the other hand, the reflectivities are extremely low at higher altitudes, resulting in fewer Doppler velocity estimates.

Description: The onset of the southeast Asian monsoon during 1997 and 1998 was simulated by coupling a mesoscale atmospheric model (MM5) and a detailed, land surface model, PLACE (the Parameterization for Land-Atmosphere-Cloud Exchange). The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The control simulation with the PLACE land surface model and variable sea surface temperature captured the basic signatures of the monsoon onset processes and associated rainfall statistics. Sensitivity tests indicated that simulations were significantly improved by including the PLACE land surface model. The mechanism by which the land surface processes affect the moisture transport and the convection during the onset of the southeast Asian monsoon were analyzed. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation: the southwest low-level flow over the Indo-China peninsula and the northern, cold frontal intrusion from southern China. The surface sensible and latent heat fluxes modified the low-level temperature distribution and merit, and therefore the low-level wind due to the thermal wind effect. The more realistic forcing of the sensible and latent heat fluxes from the detailed, land surface model improved the low-level wind simulation and associated moisture transport and convection.

Description: A cloud-resolving model is used to study sensitivities of two different microphysical schemes, one is the bulk type, and the other is an explicit bin scheme, in simulating a mid-latitude squall line case (PRE-STORM, June 10-11, 1985). Simulations using different microphysical schemes are compared with each other and also with the observations. Both the bulk and bin models reproduce the general features during the developing and mature stage of the system. The leading convective zone, the trailing stratiform region, the horizontal wind flow patterns, pressure perturbation associated with the storm dynamics, and the cool pool in front of the system all agree well with the observations. Both the observations and the bulk scheme simulation serve as validations for the newly incorporated bin scheme. However, it is also shown that, the bulk and bin simulations have distinct differences, most notably in the stratiform region. Weak convective cells exist in the stratiform region in the bulk simulation, but not in the bin simulation. These weak convective cells in the stratiform region are remnants of the previous stronger convections at the leading edge of the system. The bin simulation, on the other hand, has a horizontally homogeneous stratiform cloud structure, which agrees better with the observations. Preliminary examinations of the downdraft core strength, the potential temperature perturbation, and the evaporative cooling rate show that the differences between the bulk and bin models are due mainly to the stronger low-level evaporative cooling in convective zone simulated in the bulk model. Further quantitative analysis and sensitivity tests for this case using both the bulk and bin models will be presented in a companion paper.

Description: Despite the obvious notion that the presence of hail or graupel is a good indication of convection, the model results show this does not provide an objective benchmark partly due to the unrealistic presence of small amounts of hail or graupel throughout the anvil in the model but mainly because of the significant amounts of hail or graupel, especially in the tropical TOGA COARE simulation, in the transition zone. Without use of a "transition" category, it is open to debate as how this region should best be defined, as stratiform or as convective. So, the presence of significant hail or graupel contents in this zone significantly degrades its use an objective benchmark for convection. The separation algorithm comparison was done in the context of a cloud-resolving model. These models are widely used and serve a variety of purposes especially with regard to retrieving information that cannot be directly measured by providing synthetic data sets that are consistent and complete. Separation algorithms are regularly applied in these models. However, as with any modeling system, these types 'of models are constantly being improved to overcome any known deficiencies and make them more accurate representations of observed systems. The presence of hail and graupel in the anvil and the bias towards heavy rainfall rates are two such examples of areas that need improvement. Since, both of these can effect the perceived performance of the separation algorithms, the Lang et al. (2003) study did not want to overstate the relative performance of any specific algorithms.

Description: In this study, we investigate the seasonal variations in surface rainfall and associated large-scale processes in the tropical eastern Atlantic and West African region. The 5-yr (1998-2002) high-quality TRMM rainfall, sea surface temperature (SST), water vapor and cloud liquid water observations are applied along with the NCEP/NCAR reanalysis wind components and a 3-yr (2000-2002) Quickscat satellite-observed surface wind product. Major mean rainfall over West Africa tends to be concentrated in two regions and is observed in two different seasons, manifesting an abrupt shift of the mean rainfall zone during June-July. (i) Near the Gulf of Guinea (about 5 degN), intense convection and rainfall are seen during April-June and roughly follow the seasonality of SST in the tropical eastern Atlantic. (ii) Along the latitudes of about 10 deg. N over the interior West African continent, a second intense rain belt begins to develop from July and remains there during the later summer season. This belt co-exists with a northwardmoved African Easterly Jet (AEJ) and its accompanying horizonal and vertical shear zones, the appearance and intensification of an upper tropospheric Tropical Easterly Jet (TEJ), and a strong low-level westerly flow. Westward-propagating wave signals [ i e . , African easterly waves (AEWs)] dominate the synoptic-scale variability during July-September, in contrast to the evident eastward-propagating wave signals during May- June. The abrupt shift of mean rainfall zone thus turns out to be a combination of two different physical processes: (i) Evident seasonal cycles in the tropical eastern Atlantic ocean which modulate convection and rainfall in the Gulf of Guinea by means of SST thermal forcing and SST-related meridional gradient; (ii) The interaction among the AEJ, TEJ, low-level westerly flow, moist convection and AEWs during July-September which modulates rainfall variability in the interior West Africa, primarily within the ITCZ rain band. Evident seasonality in synoptic-scale wave signals is shown to be a good evidence for this seasonal evolution.

Description: We seek funding from NASA for the third year (2005) of the four-year period January 1, 2003 - December 31, 2006 for continued support of the MIT contributions to the multi-national global atmospheric trace species measurement program entitled Advanced Global Atmospheric Gases Experiment (AGAGE). The case for real-time high-frequency measurement networks like AGAGE is very strong and the observations and their interpretation are widely recognized for their importance to ozone depletion and climate change studies and to verification issues arising from the Montreal Protocol (ozone) and Kyoto Protocol (climate). The proposed AGAGE program is distinguished by its capability to measure over the globe at high frequency almost all of the important species in the Montreal Protocol and almost all of the significant non-CO2 gases in the Kyoto Protocol.

Description: Computational phenomena (i.e., spurious supersaturation and negative mixing ratio of cloud water) usually exist in cloud-resolving models when the time step for explicit integration is larger than a microphysical timescale in clouds. In this paper, the microphysical timescales in clouds are studied, showing that the timescale of water vapor condensation (or cloud water evaporation) is smaller than 10 s - the order of a typical time step for cloud-resolving models. To avoid spurious computational phenomena in cloud-resolving modeling, it is suggested that moist entropy be used as a prognostic thermodynamic variable, and temperature be diagnosed from that and other prognostic variables. A simple numerical model with moist entropy as a prognostic variable, for example, is presented to show that spurious computational phenomena are removed when moist entropy is used as a prognostic variable.

Description: Recently Grabowski and Khairoutdinov and Randall have proposed and demonstrated, the feasibility of using 2D CRMs (these did not include detailed land processes and used periodic lateral boundary conditions) as 'super-parameterizations or a multi-scale modeling framework MMF' for cloud processes within atmospheric generation circulation models (GCMs). In the MMF, a fine-resolution 2D CRM takes the place of the single-column parameterization that is used in conventional GCMs. An overview of this approach is given by Randall et al.. Since a CRM can explicitly simulate cloud processes at the natural space and time scales of clouddynamical processes, it's possible to compute statistics of cloud fields, including radiances and radar backscatters, directly simulated/produced by the CRM. A seed fund is available at NASA Goddard to build a MMF based on a 2D GCE model and the Goddard fvGCM. A prototype MMF will be developed at end of 2004. The major objectives are to test and evaluate its performance using the GCE model as the super parameterization and compare its results with other super parameterizations, conventional GCMs and observations. Recently, a detailed spectral-bin microphysical scheme was implemented into the Goddard Cumulus Ensemble (WE) model. The GCE model is a CRM and it has been improved and applied to study the precipitation processes in the past two decades. Atmospheric aerosols are also described using number density size-distribution functions. The model is tested by studying the evolution of deep tropical clouds in the west Pacific warm pool region and summertime convection over a mid-latitude continent with different concentrations of CCN: a low 'clean' concentration and a high 'dirty' concentration. The impact of atmospheric aerosol concentration on cloud and precipitation will be investigated. In addition, the similarities and differences between bulk- and spectral-bin microphysics will be presented,

Description: Fast amplitude variations due to atmospheric scintillation are the main concerns for the Deep Space Network (DSN) Ka-band downlink under clear weather conditions. A theoretical study of the amplitude scintillation variances for a finite aperture antenna is presented. Amplitude variances for weak scattering scenarios are examined using turbulence theory to describe atmospheric irregularities. We first apply the Kolmogorov turbulent spectrum to a point receiver for three different turbulent profile models, especially for an exponential model varying with altitude. These analytic solutions then are extended to a receiver with a finite aperture antenna for the three profile models. Smoothing effects of antenna aperture are expressed by gain factors. A group of scaling factor relations is derived to show the dependences of amplitude variances on signal wavelength, antenna size, and elevation angle. Finally, we use these analytic solutions to estimate the scintillation intensity for a DSN Goldstone 34-m receiving station. We find that the (rms) amplitude fluctuation is 0.13 dB at 20-deg elevation angle for an exponential model, while the fluctuation is 0.05 dB at 90 deg. These results will aid us in telecommunication system design and signal-fading prediction. They also provide a theoretical basis for further comparison with other measurements at Ka-band.

Description: A key aspect of this project is the establishment of a causal link from circulation anomalies in the lower mesosphere and stratopause region downward through the stratosphere to the troposphere. The observational link for stratospheric sudden warmings and surface climate is fairly clear. However, our understanding of the dynamics is incomplete. We have been making significant progress in the area of dynamical mechanisms by which circulation anomalies in the stratosphere affect the troposphere. We are trying to understand the details and sequence of events that occur when a middle atmosphere (wind) anomaly propagates downward to near the tropopause. The wind anomaly could be caused by a warming or solar variations in the low-latitude stratopause region, or could have other causes. The observations show a picture that is consistent with a circulation anomaly that descends to the tropopause region, and can be detected as low as the mid-troposphere. Processes near the stratopause in the tropics appear to be important precursors to the wintertime development of the northern polar vortex. This may affect significantly our understanding of the process by which low-latitude wind anomalies in the low mesosphere and upper stratosphere evolve through the winter and affect the polar vortex.

Description: The Pennsylvania State University-National Center for Atmospheric Research mesoscale model MM5 is used to simulate Hurricane Bonnie at high resolution (2-km spacing) in order to examine budgets of water vapor, cloud condensate, and precipitation. Virtually all budget terms are derived directly from the model (except for the effects of storm motion). The water vapor budget reveals that a majority of the condensation in the eyewall occurs in convective hot towers, while outside of the eyewall most of the condensation occurs in weak updrafts, indicative of a dominance of stratiform precipitation processes. The ocean source of water vapor in the eyewall region is only a very small fraction of that transported inward in the boundary layer inflow or that condensed in the updrafts. In contrast, in the outer regions, the ocean vapor source is larger owing to the larger area, counters the drying effect of low-level subsidence, and enhances the moisture transported in toward the eyewall. In this mature storm, cloud condensate is consumed as rapidly as it is produced. Cloud water peaks at the top of the boundary layer and within the melting layer, where. cooling from melting enhances condensation. Unlike squall lines, in the hurricane, very little condensate produced in the eyewall convection is transported outward into the surrounding precipitation area. Most of the mass ejected outward is likely in the form of small snow particles that seed the outer regions and enhance stratiform precipitation development through additional growth by vapor deposition and aggregation. In addition, artificial source terms for cloud and precipitation mass associated with setting negative mixing ratios, produced as a result of numerical advection errors, to zero are described. Although small at any given point and time, the cumulative effect of these terms contributes an amount of mass equivalent to 13% of the total condensation. Thus, these terms must be accounted for to balance the model budgets and the results suggest the need for improved model numerics.

Description: The notion that the continental-scale land-sea contrast is the main reason that monsoon circulation exists has been a long-held belief. The purpose of this paper is to point out that this notion should be substantially modified. The central idea of this notion states that in summer, radiative heating of the continent, say Asia, gives rise to a continental-scale thermal low and surrounding the thermal low in its southeast direction the low level wind flows in from south-west. This low-level inflow creates a convergence of moisture, which maintains the cumulus convection. And in winter, radiative cooling of continent gives rise to a thermal high and to its southeast the low-level wind is from northeast. The mechanism in this interpretation does undoubtedly exist. However, this mechanism, though believed to be the main driving force of monsoon, has not been tested in numerical experiments. There has been an increasing recognition in the recent years that monsoon is inextricably tied to the heating in the intertropical convergence zone (ITCZ). We propose that the main cause of monsoon is ITCZ's being substantially away from the equator. A brief qualitative explanation of why the ITCZ can be a source of monsoon circulation can be offered based on the circulation field forced by the ITCZ heating. The existence of the ITCZ's does not always have to rely on land-sea contrast on the continental scale. This is hinted in the fact that in February the ITCZ close to Australia (and its associated monsoon circulation) covers a longitudinal range several times as long as that of Australia and thus cannot possibly be caused mainly by the land-sea contrast associated with Australia. Yet, this cannot be used as a proof that the ITCZ in the Asian summer monsoon is not mainly due to land-sea contrast. One of the purposes of this work is to provide a convincing proof. In this work the role of land-sea contrast in the origin of monsoon is examined through numerical simulation with the Goddard general circulation model. The Asian and Australian monsoon circulations are obtained in a four-year integration and then the integration is repeated with Asia, the maritime continent, and Australia replaced by ocean. The sea surface temperature (SST) at each affected grid is specified as the SST at the first grid to the east that is an ocean grid in the first experiment. The latter integration shows that the monsoon circulation pattern over where south Asia and Australia were and the surrounding region has largely remained. The results discount land-sea contrast as the main cause of Asian monsoon. A third experiment is the same as the first except that the topography of Asia, the maritime continent, and Australia is reduced to zero. This experiment reveals that the difference between the first two experiments is due more to the removal of topography than to the removal of land-sea contrast. August precipitation is shown averaged over the last three years of each of the three experiments. They show that the Asian monsoon rainy region is largely intact in the second experiment and the difference between the second and the third experiment is mainly in the longitudinal location of the maximum precipitation. Additionally, in Asian and Australian winter monsoons land-sea contrast also plays only a modifying role. Although land-sea contrast plays only a modifying role in Asian and Australian (and Central American including Mexican) monsoons, it is the main reason that ITCZ (and thus monsoon) exists in Africa and South America. Thus, monsoons can be classified into two groups depending on whether land-sea contrast plays a major role.

Description: In dynamic meteorology, singular vectors (SVs) are the structures that maximize a given norm of a forecast perturbation given a tangent linear model and a quadratic constraint on the initial perturbation. That constraint is a prescription of the value of either the same or a different norm applied to the initial perturbations. In the sense that SVs maximize the forecast perturbations according to a specified measure, they may be considered as optimal perturbations. SVs are used to characterize predictability, to identify and correct possible initial errors given forecast errors, to create a set of significant perturbations for ensemble forecasting, or to determine locations for observation targeting. They are a specific application of generally defined singular vectors in mathematics.

Description: Land-atmosphere feedback, by which precipitation-induced soil moisture anomalies affect subsequent precipitation, may be an important element of Earth's climate system, but its very existence has never been demonstrated conclusively at regional to continental scales. Evidence for the feedback is sought in a 50-year observational precipitation dataset covering the United States. The precipitation variance and autocorrelation fields are characterized by features that agree (in structure, though not in magnitude) with those produced by an atmospheric general circulation model (AGCM). Because the model-generated features are known to result from land-atmosphere feedback alone, the observed features are highly suggestive of the existence of feedback in nature.

Description: The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China. Multiple observation platforms (e.g., upper-air soundings, Doppler radar, ships, wind profilers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convective storms and air pattern changes associated with monsoons over the South China Sea region. SCSMEX also provided rainfall estimates which allows for comparisons with those obtained from the Tropical Rainfall Measuring Mission (TRMM), a low earth orbit satellite designed to measure rainfall from space. The Goddard Cumulus Ensemble (GCE) model (with 1-km grid size) is used to understand and quantify the precipitation processes associated with the summer monsoon over the South China Sea. This is the first (loud-resolving model used to simulate precipitation processes in this particular region. The GCE-model results captured many of the observed precipitation characteristics because it used a fine grid size. For example, the temporal variation of the simulated rainfall compares quite well to the sounding-estimated rainfall variation. The time and domain-averaged temperature (heating/cooling) and water vapor (drying/ moistening) budgets are in good agreement with observations. The GCE-model-simulated rainfall amount also agrees well with TRMM rainfall data. The results show there is more evaporation from the ocean surface prior to the onset of the monsoon than after the on-~et of monsoon when rainfall increases. Forcing due to net radiation (solar heating minus longwave cooling) is responsible for about 25% of the precipitation in SCSMEX The transfer of heat from the ocean into the atmosphere does not contribute significantly to the rainfall in SCSMEX. Model sensitivity tests indicated that total rain production is reduced 17-18% in runs neglecting the ice phase. The SCSMEX results are compared to other GCE-model-simulated weather systems that developed during other field campaigns (i.e., west Pacific warm pool region, eastern Atlantic region and central USA). Large-scale forcing vie temperature and water vapor tendency, is the major energy source for net condensation in the tropical cases. The effects of large-scale cooling exceed that of large-scale moistening in the west pacific warm pool region and eastern Atlantic region. For SCSMEX, however, the effects of large-scale moistening predominate. Net radiation and sensible and latent hc,it fluxes play a much more important role in the central USA.

Description: CloudSat is a satellite experiment designed to measure the vertical structure of clouds from space. The expected launch of CloudSat is planned for 2004 and, once launched, CloudSat will orbit in formation as part of a constellation of satellites including NASA's Aqua and Aura satellites, a NASA-CNES lidar satellite (P-C) and a CNES satellite carrying a polarimeter (PARASOL). A unique feature that CloudSat brings to this constellation is the ability to fly a precise orbit enabling the fields of view of the CloudSat radar to be overlapped with the P-C lidar footprint and the other measurements of the EOS constellation. The precision of this overlap creates a unique multi-satellite observing system for studying the atmospheric processes essential to the hydrological cycle. The vertical profile of cloud properties provided by CloudSat fills a critical gap in the investigation of feedback mechanisms linking clouds to climate. Measuring the vertical profile of cloud properties requires a combination of active and passive instruments, and this will be achieved by combining the radar data of CloudSat with active and passive data from other sensors of the constellation. This paper describes the underpinning science, and gives an overview of the mission, and provides some idea of the expected products and anticipated application of these products. Notably, the CloudSat mission is expected to provide new knowledge about global cloudiness, stimulating new areas of research on clouds including data assimilation and cloud parameterization. The mission also provides an important opportunity to demonstrate active sensor technology for future scientific and tactical applications. The CloudSat mission is a partnership between NASA/JPL, the Canadian Space Agency, Colorado State University, the US Air Force, and the US Department of Energy.

Description: Several heavy precipitation episodes occurred over Taiwan from August 10 to 13, 1994. Precipitation patterns and characteristics are quite different between the precipitation events that occurred from August 10 and I I and from August 12 and 13. In Part I (Chen et al. 2001), the environmental situation and precipitation characteristics are analyzed using the EC/TOGA data, ground-based radar data, surface rainfall patterns, surface wind data, and upper air soundings. In this study (Part II), the Penn State/NCAR Mesoscale Model (MM5) is used to study the precipitation characteristics of these heavy precipitation events. Various physical processes (schemes) developed at NASA Goddard Space Flight Center (i.e., cloud microphysics scheme, radiative transfer model, and land-soil-vegetation surface model) have recently implemented into the MM5. These physical packages are described in the paper, Two way interactive nested grids are used with horizontal resolutions of 45, 15 and 5 km. The model results indicated that Cloud physics, land surface and radiation processes generally do not change the location (horizontal distribution) of heavy precipitation. The Goddard 3-class ice scheme produced more rainfall than the 2-class scheme. The Goddard multi-broad-band radiative transfer model reduced precipitation compared to a one-broad band (emissivity) radiation model. The Goddard land-soil-vegetation surface model also reduce the rainfall compared to a simple surface model in which the surface temperature is computed from a Surface energy budget following the "force-re store" method. However, model runs including all Goddard physical processes enhanced precipitation significantly for both cases. The results from these runs are in better agreement with observations. Despite improved simulations using different physical schemes, there are still some deficiencies in the model simulations. Some potential problems are discussed. Sensitivity tests (removing either terrain or radiative processes) are performed to identify the physical processes that determine the precipitation patterns and characteristics for heavy rainfall events. These sensitivity tests indicated that terrain can play a major role in determining the exact location for both precipitation events. The terrain can also play a major role in determining the intensity of precipitation for both events. However, it has a large impact on one event but a smaller one on the other. The radiative processes are also important for determining, the precipitation patterns for one case but. not the other. The radiative processes can also effect the total rainfall for both cases to different extents.

Description: Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. The simulations are made using two different numerical advection schemes: 4th order and positive definite advection. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain. The midlatitude PRESTORM simulation was found to be nearly invariant with respect to advection type for most quantities while for TOGA COARE fourth order advection produced numerous shallow convective cores and positive definite advection fewer cells that were both broader and deeper penetrating above the freezing level.

Description: This study used a two-dimensional coupled land/atmosphere (cloud-resolving) model to investigate the influence of land cover on the water budgets of squall lines in the Sahel. Study simulations used the same initial sounding and one of three different land covers, a sparsely vegetated semi-desert, a grassy savanna, and a dense evergreen broadleaf forest. All simulations began at midnight and ran for 24 hours to capture a full diurnal cycle. In the morning, the latent heat flux, boundary layer mixing ratio, and moist static energy in the boundary layer exhibited notable variations among the three land covers. The broadleaf forest had the highest latent heat flux, the shallowest, moistest, slowest growing boundary layer, and significantly more moist static energy per unit area than the savanna and semi-desert. Although all simulations produced squall lines by early afternoon, the broadleaf forest had the most intense, longest-lived squall lines with 29% more rainfall than the savanna and 37% more than the semi-desert. The sensitivity of the results to vegetation density, initial sounding humidity, and grid resolution was also assessed. There were greater differences in rainfall among land cover types than among simulations of the same land cover with varying amounts of vegetation. Small changes in humidity were equivalent in effect to large changes in land cover, producing large changes in the condensate and rainfall. Decreasing the humidity had a greater effect on rainfall volume than increasing the humidity. Reducing the grid resolution from 1.5 km to 0.5 km decreased the temperature and humidity of the cold pools and increased the rain volume.

Description: An extensive body of research this quarter is documented. Further methodical analysis of temperature residuals in Cryogenic Limb Array Etalon Spectrometer (CLAES) Version 8 level 3AT data show signatures during December 1992 at middle and high northern latitudes that, when compared to Naval Research Laboratory/Mountain Wave Forecast Model (NRL)/(MWFM) mountain wave hindcasts, reveal evidence of long mountain waves in these data over Eurasia, Greenland, Scandinavia and North America. The explicit detection of gravity waves in limb-scanned Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere (CRISTA) temperatures is modeled at length, to derive visibility functions. These insights are used to convert CRISTA gravity wave temperature residuals into data that more closely resemble gravity wave fluctuations detected in data from other satellite instruments, such as Microwave Limb Sounder (MLS), Limb Infrared Monitor of the Stratosphere (LIMS) and Global Positioning System/Meteorology (GPS)/(MET). Finally, newly issued mesospheric temperatures from inversion of CRISTA 15gin emissions are analyzed using a new method that uses separate Kalman fits to the ascending and descending node data. This allows us to study global gravity wave amplitudes at two local times, 12 hours apart. In the equatorial mesosphere, where a large diurnal tidal temperature signal exists, we see modulations of gravity wave activity that are consistent with gravity wave-tidal interactions produced by tidal temperature variability.

Description: Since the advent of meteorological satellites in the 1960's, numerous experiments have been conducted in order to evaluate the impact of these and other data on atmospheric analysis and prediction. Such studies have included both OSE'S and OSSE's. The OSE's were conducted to evaluate the impact of specific observations or classes of observations on analyses and forecasts. Such experiments have been performed for selected types of conventional data and for various satellite data sets as they became available. (See for example the 1989 ECMWF/EUMETSAT workshop proceedings on "The use of satellite data in operational numerical weather prediction" and the references contained therein.) The ODYSSEY were conducted to evaluate the potential for future observing systems to improve Numerical Weather Prediction NWP and to plan for the Global Weather Experiment and more recently for EVANS (Atlas et al., 1985a; Arnold and Day, 1986; Hoffman et al., 1990). In addition, OSSE's have been run to evaluate trade-offs in the design of observing systems and observing networks (Atlas and Emmitt, 1991; Rohaly and Krishnamurti, 1993), and to test new methodology for data assimilation (Atlas and Bloom, 1989).

Description: Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain.

Description: Hurricane Earl developed from a tropical wave that moved into the Gulf of Mexico, which triggered abundant convection. On 1 Sept. 1998, the wave was upgraded directly to a tropical storm. Earl reached hurricane status the next morning. The system moved erratically as it interacted with an upper level short wave trough rotating around a long wave trough to the northeast. The storm made landfall near 0600 UTC on 3 September near Panama City, FL. During August and September 1998, NASA conducted the Third Convection and Moisture Experiment (CAMEX-3). It focused on studying the intensity, track, and impacts at landfall of hurricanes. On the afternoon of 2 September 1998, the NASA ER2 high-altitude aircraft flying at 65,000 feet in tandem with the NASA DC-8 flying at 35,000 feet flew over and through, respectively, the eastern rainbands of Earl near the Florida Panhandle as the storm neared landfall in the region. Two approaches to studying Earl are undertaken here: first, an examination of the source and height of the dry air region using GOES-8 water vapor data and, second, a look into the impact of the dry air entrainment on the system using aircraft remote sensing data.

Description: The penetration of the MARSIS radar signal into the polar ice mass is modeled to determine the capability of the instrument to locate sub-glacial aquifers. As a ground penetrating radar, the orbiting MARSIS transmits a signal greater than 1 W between 1-5 MHz. In this work we will investigate the effect of ice conductive losses on the radar-detection of subsurface aquifers. Based on wave propagation analysis, it is found that for a bulk ice conductivity below 10-5 S/m, conductive losses in the medium are not significant. However, if the bulk ice conductivity is relatively large (greater than 10-5 S/m), the reflected signal from any deep aquifer will be absorbed as it propagates in the lossy ice medium limiting the probing depth.

Description: Summer tropical synoptic-scale waves over West Africa are quantified by the 850 mb meridional wind component from the NCEP/NCAR reanalysis project. Their relationships with surface precipitation patterns are further explored by applying the data from the Tropical Rainfall Measuring Mission (TRMM) satellite in combination with other satellite observations during 1998-2002. Evident wavelet spectral power peaks are seen within a period of 2.5 - 6 days in both meridional wind and precipitation. The most intense wave signals in meridional wind are concentrated along 15 deg N- 25 deg N. Wave signals in precipitation and corresponding wavelet cross-spectral signals between these two variables, however, are primarily located at 5 deg N- 15 deg N, the latitudes of major summer rain events. There is a tendency for the perturbations in meridional wind component to lag (lead) precipitation signals south (north) of 15 deg N. In some cases, either an in-phase or out-of-phase relationship can even be found between these two variables, suggesting a latitude-dependent horizontal structure for these waves and probably implying two distinct wave-convective coupling mechanisms. Moreover, the lagging relationship (and/or the out-of-phase tendency) is only observed south of 15 deg N during July-September, indicating a strong seasonal preference. This phase relationship is generally consistent with the horizontal wave structures from a composite analysis.

Description: The development of a ground based direct detection Doppler lidar based on the recently described aerosol double edge technique is reported. A pulsed, injection seeded Nd:YAG laser operating at 1064 nm is used to make range resolved measurements of atmospheric winds in the free troposphere. The wind measurements are determined by measuring the Doppler shift of the laser signal backscattered from atmospheric aerosols. The lidar instrument and double edge method are described and initial tropospheric wind profile measurements are presented. Wind profiles are reported for both day and night operation. The measurements extend to altitudes as high as 14 km and are compared to rawinsonde wind profile data from Dulles airport in Virginia. Vertical resolution of the lidar measurements is 330 m and the rms precision of the measurements is a low as 0.6 m/s.

Description: The SeaWinds on QuikSCAT scatterometer was developed by NASA JPL to measure the speed and direction of ocean surface winds. Simulations performed to estimate the performance of the instrument prior to its launch have indicated that the mid-swath accuracy is worse than that of the rest of the swath. This behavior is a general characteristic of scanning pencil beam scatterometers. For SeaWinds, the accuracy of the rest of the swath, and the size of the swath are such that the instrument meets its science requirements despite mid-swath shortcomings. However, by understanding the problem at mid-swath, we can improve the performance there as well. We discuss the underlying causes of the problem in detail and propose a new wind retrieval algorithm which improves mid-swath performance. The directional discrimination ability of the instrument varies with cross track distance wind speed, and direction. By estimating the range of likely wind directions for each measurement cell, one can optimally apply information from neighboring cells where necessary in order to reduce random wind direction errors without significantly degrading the resolution of the resultant wind field. In this manner we are able to achieve mid-swath RMS wind direction errors as low as 15 degrees for low winds and 10 degrees for moderate to high winds, while at the same time preserving high resolution structures such as cyclones and fronts.

Description: The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data fiom the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo- China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the lowlevel temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation.

Description: This quarter was largely devoted to a detailed study of temperature data acquired by the Cryogenic Limb Array Etalon Spectrometer (CLAES) on UARS. Our analysis used the same sequence of methods that have been developed, tested and refined on a more limited subset of temperature data acquired by the CRISTA instrument. We focused on a limited subset of our reasoning that geographical and vertical trends in the small-scale temperature variability could be compared with similar trends observed in November 1994 by the CRISTA-SPAS satellite. Results, backed up with hindcasts from the Mountain Wave Forecast Model (MWFM), reveal strong evidence of mountain waves, most persuasively in the Himalayas on 16-17 November, 1992. These CLAES results are coherent over the 30-50 km range and compare well with MWFM hindcasts for the same period. This constitutes, we believe, the first clear evidence that CLAES explicitly resolved long wavelength gravity waves in its CO2 temperature channel. A series of other tasks, related to mesoscale modeling of mountain waves in CRISTA data and fitting of ground-based and HRDI data on global scales, were seen through to publication stage in peer-reviewed journals.

Description: Ideal cloud-resolving models contain little-accumulative errors. When their domain is so large that synoptic large-scale circulations are accommodated, they can be used for the simulation of the interaction between convective clouds and the large-scale circulations. This paper sets up a framework for the models, using moist entropy as a prognostic variable and employing conservative numerical schemes. The models possess no accumulative errors of thermodynamic variables when they comply with a discrete constraint on entropy conservation and sound waves. Alternatively speaking, the discrete constraint is related to the correct representation of the large-scale convergence and advection of moist entropy. Since air density is involved in entropy conservation and sound waves, the challenge is how to compute sound waves efficiently under the constraint. To address the challenge, a compensation method is introduced on the basis of a reference isothermal atmosphere whose governing equations are solved analytically. Stability analysis and numerical experiments show that the method allows the models to integrate efficiently with a large time step.

Description: Reasor et al. (2004) examined the tilting of a vortex by shear and found that a damping mechanism intrinsic to the dry adiabatic dynamics suppresses departures from an upright state. This realignment occurs through projection of the tilt asymmetry onto two types of vortex Rossby waves: sheared vortex Rossby waves, in which the radial shear of the swirling flow axisymmetrizes tilt asymmetries, and a quasi-mode, or discrete, vortex Rossby wave that, in the absence of damping, causes precession of the upper vortex until it realigns with the lower one (thereafter undergoing repeated cycles of tilting and realignment). With damping, the vortex achieves a downshear-left equilibrium tilt. This paper examines the role of a damped quasi-mode in producing asymmetric vertical motions in a high-resolution simulation of Hurricane Bonnie (1998).

Description: Land-atmosphere feedback, by which (for example) precipitation-induced moisture anomalies at the land surface affect the overlying atmosphere and thereby the subsequent generation of precipitation, has been examined and quantified with many atmospheric general circulation models (AGCMs). Generally missing from such studies, however, is an indication of the extent to which the simulated feedback strength is model dependent. Four modeling groups have recently performed a highly controlled numerical experiment that allows an objective inter-model comparison of land-atmosphere feedback strength. The experiment essentially consists of an ensemble of simulations in which each member simulation artificially maintains the same time series of surface prognostic variables. Differences in atmospheric behavior between the ensemble members then indicates the degree to which the state of the land surface controls atmospheric processes in that model. A comparison of the four sets of experimental results shows that feedback strength does indeed vary significantly between the AGCMs.

Description: The One-Degree Daily (1DD) technique is described for producing globally complete daily estimates of precipitation on a 1 deg x 1 deg lat/long grid from currently available observational data. Where possible (40 deg N-40 deg S), the Threshold-Matched Precipitation Index (TMPI) provides precipitation estimates in which the 3-hourly infrared brightness temperatures (IR T(sub b)) are thresholded and all "cold" pixels are given a single precipitation rate. This approach is an adaptation of the Geostationary Operational Environmental Satellite (GOES) Precipitation Index (GPI), but for the TMPI the IR Tb threshold and conditional rain rate are set locally by month from Special Sensor Microwave/Imager (SSM/I)-based precipitation frequency and the Global Precipitation Climatology Project (GPCP) satellite-gauge (SG) combined monthly precipitation estimate, respectively. At higher latitudes the 1DD features a rescaled daily Television Infrared Observation Satellite (TIROS) Operational Vertical Sounder (TOVS) precipitation. The frequency of rain days in the TOVS is scaled down to match that in the TMPI at the data boundaries, and the resulting non-zero TOVS values are scaled locally to sum to the SG (which is a globally complete monthly product). The time series of the daily 1DD global images shows good continuity in time and across the data boundaries. Various examples are shown to illustrate uses. Validation for individual grid -box values shows a very high root-mean-square error but, it improves quickly when users perform time/space averaging according to their own requirements.

Description: A two-dimensional version of the Goddard Cumulus Ensemble (GCE) Model is used to simulate convective systems that developed in various geographic locations. Observed large-scale advective tendencies for potential temperature, water vapor mixing ratio, and horizontal momentum derived from field campaigns are used as the main forcing. By examining the surface energy budgets, the model results show that the two largest terms are net condensation (heating/drying) and imposed large-scale forcing (cooling/moistening) for tropical oceanic cases. These two terms arc opposite in sign, however. The contributions by net radiation and latent heat flux to the net condensation vary in these tropical cases, however. For cloud systems that developed over the South China Sea and eastern Atlantic, net radiation (cooling) accounts for about 20% or more of the net condensation. However, short-wave heating and long-wave cooling are in balance with each other for cloud systems over the West Pacific region such that the net radiation is very small. This is due to the thick anvil clouds simulated in the cloud systems over the Pacific region. Large-scale cooling exceeds large-scale moistening in the Pacific and Atlantic cases. For cloud systems over the South China Sea, however, there is more large-scale moistening than cooling even though the cloud systems developed in a very moist environment. though For three cloud systems that developed over a mid-latitude continent, the net radiation and sensible and latent heat fluxes play a much more important role. This means the accurate measurement of surface fluxes and radiation is crucial for simulating these mid-latitude cases.

Description: We compute upper tropospheric relative humidity profiles using water vapor profiles measured by an airborne DIAL and a ground-based Raman lidar. LASE water vapor and MTP temperature profiles acquired from the NASA DC-8 aircraft during the recent Pacific Exploratory Mission Tropics B (PEM Tropics B) field mission in the tropical Pacific and the SAGE-III Ozone Loss and Validation Experiment (SOLVE) in the Arctic as well as water vapor profiles derived from the ground-based DOE ARM Southern Great Plains (SGP) CART Raman lidar are used. Comparisons of the lidar water vapor measurements with available in situ measurements show reasonable agreement for water vapor mixing ratios above 0.05 g/kg. Relative humidity frequency distributions computed using LASE data indicate that ice supersaturation occurred about 5-11% of the time when temperatures were below -35 C. While a higher frequency of ice supersaturation was observed during SOLVE, higher peak values of relative humidity were observed during PEM Tropics B. The relative humidity fields associated with cirrus clouds are also examined.

Description: An observing system comprised of two lidars in geosychronous orbit would enable the synoptic and meso-scale measurement of atmospheric winds and moisture both of which are key first-order variables of the Earth's weather equation. Simultaneous measurement of these parameters at fast revisit rates promises large advancements in our weather prediction skills. Such capabilities would be unprecedented and a) yield greatly improved and finer resolution initial conditions for models, b) make existing costly and cumbersome measurement approaches obsolete, and c) obviate the use of numerical techniques needed to correct data obtained using present observing systems. Additionally, simultaneous synoptic wind and moisture observations would lead to improvements in model parameterizations, and in our knowledge of small-scale weather processes. Technology and science data product assessments are ongoing. Results will be presented during the conference.

Description: The goals of this research are: 1) To validate TIMED Doppler Interferometer (TIDI) winds using ground-based MF and meteor winds; and 2) To examine short-term (i. e., day-to-day and week-to-week) variability of the diurnal tide. This objective was to have originally been met using comparisons of short-term diurnal tidal determinations from ground-based (GB) winds with planetary-scale diurnal nonmigrating tidal definitions from TIDI winds.

Description: The first storm-scale, total lightning observations from space during tornadogenesis are presented. During the overpass of an Oklahoma supercell, just minutes prior to tornado touchdown on 17 April 1995, the NASA (National Aeronautics and Space Administration) OTD (Optical Transient Detector) detected a total of 143 flashes during approximately 3 minutes of observation time. The estimated total flash rate ranges from 45 (raw counts) to 78 (corrected for detection efficiency) flashes min-1. This total flash rate was at least 17 times greater than the cloud-to-ground lightning rate detected by the National Lightning Detection Network (NLDN), indicating most of the lightning was intracloud. Cloud-to-ground lightning at this time was also dominated by positive polarity flashes. In addition, total lightning rates were decreasing rapidly prior to touchdown. These OTD observations are consistent with the limited results from recent ground based measurements of total lightning activity in tornadic storms and corroborate that such storms have unusually high total flash rates, are dominated by intracloud lightning, and that the total flash rates are observed to decrease rapidly in the minutes prior to touchdown.

Description: The research facilities on the Svalbard archepelego provide a unique opportunity for observing the temporal and spatial characteristics of the interaction of the solar wind with the magnetosphere-ionosphere system. The first sounding rockets from the SvdRak range at Ny-Alesund have opened a new perspective for separating temporal and spatial effects. Finding a correlation between observations from the rockets and the Wind satellite located nearly 200 R(sub E) upstream in the solar wind that had a lag time less than the advection time forced consideration of tilted phase planes of the interp1anetary electric field. From this it was deduced that the interaction process has to involve high latitude merging, the interplanetary magnetic field (IMF) B(sub Y) bifurcates the cusp relative to the high latitude source regions in the Northern and Southern Hemispheres, IMF B(sub X) controls the interaction time in each hemisphere, and the small convection cell is driven by opposite hemisphere merging. Studies with 4 satellites in the solar wind have shown that the tilt of the phase plane varies on a minute-by-minute basis and change significantly on time scales of tens of minutes. Images of the cusp at 557.7 nm provide a nearly instantaneous picture of the temporal and spatial characteristics of merging at the magnetopause and show that the process can occur at multiple locations. The rate and location vary with time.

Description: This report summarizes the progress on our three-year program of research to refine the measurement capability for satellite-based instruments that monitor ozone and other trace species in the Earth's stratosphere and troposphere, to retrieve global distributions of these and other constituents h m the GOME and SCIAMACHY satellite instruments, and to conduct scientific studies for the ILAS instruments. This continues our involvements as a U.S. participant in GOME and SCIAMACHY since their inception, and as a member of the ILAS-II Science Team. These programs have led to the launch of the first satellite instrument specifically designed to measure height-resolved ozone, including the tropospheric component (GOME), and the development of the first satellite instrument that will measure tropospheric ozone simultaneously with NO2, CO, HCHO, N2O, H2O, and CH4 (SCIAMACHY). The GOME program now includes the GOME-2 instruments, to be launched on the Eumetsat Metop satellites, providing long-term continuity in European measurements of global ozone that complement the measurements of the TOMS, SBUV, OMI, OMPS instruments. The research primarily focuses on two areas: Data analysis, including algorithm development and validation studies that will improve the quality of retrieved data products, in support for future field campaigns (to complement in situ and airborne campaigns with satellite measurements), and scientific analyses to be interfaced to atmospheric modeling studies.

Description: The onset of the southeast Asian monsoon during 1997 and 1998 was simulated by coupling a mesoscale atmospheric model (MM5) and a detailed, land surface model, PLACE (the Parameterization for Land-Atmosphere-Cloud Exchange). The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The control simulation with the PLACE land surface model and variable sea surface temperature captured the basic signatures of the monsoon onset processes and associated rainfall statistics. Sensitivity tests indicated that simulations were sigmficantly improved by including the PLACE land surface model. The mechanism by which the land surface processes affect the moisture transport and the convection during the onset of the southeast Asian monsoon were analyzed. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation: the southwest low-level flow over the Indo-china peninsula and the northern, cold frontal intrusion from southern China. The surface sensible and latent heat fluxes modified the low-level temperature distribution and gradient, and therefore the low-level wind due to the thermal wind effect. The more realistic forcing of the sensible and latent heat fluxes from the detailed, land surface model improved the low-level wind simulation apd associated moisture transport and convection.

Description: On the windward side of southwestern Taiwan, about a quarter to a half of all rainfall during mid-July through August from 1994 to 2000 came from convective systems embedded in the southwesterly monsoon flow. k this study, the causes of two heavy rainfall events (daily rainfall exceeding 100 mm day over at least three rainfall stations) observed over the slopes and/or lowlands of southwestern Taiwan were examined. Data from European Center for Medium-Range Weather Forecasts /Tropical Ocean- Global Atmosphere (EC/TOGA) analyses, the rainfall stations of the Automatic Rainfall and Meteorological Telemetry System (ARMTS) and the conventional surface stations over Taiwan, and the simulation results from a regional-scale numerical model were used to accomplish the objectives. In one event (393 mm day on 9 August 1999), heavy rainfall was observed over the windward slopes of southern Taiwan in a potentially unstable environment with very humid air around 850 hPa. The extreme accumulation was simulated and attributed to orographic lifting effects. No preexisting convection drifted in from the Taiwan Strait into western Taiwan.

Description: Rainfall production is the fundamental variable within the Earth's hydrological cycle because it is both the principal forcing term in surface water budgets and its energetics corollary, latent heating, is the principal source of atmospheric diabatic heating. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water. The properties of the vertical distribution of latent heat release modulate large-scale meridional and zonal circulations within the tropics - as well as modifying the energetic efficiencies of midlatitude weather systems. This paper focuses on the retrieval of latent heat release from satellite measurements generated by the Tropical Rainfall Measuring Mission (TRMM) satellite observatory, which was launched in November 1997 as a joint American-Japanese space endeavor. Since then, TRMM measurements have been providing an accurate four-dimensional account of rainfall over the global tropics and sub-tropics, information which can be used to estimate the space-time structure of latent heating across the Earth's low latitudes. The paper examines how the observed TRMM distribution of rainfall has advanced an understanding of the global water and energy cycle and its consequent relationship to the atmospheric general circulation and climate via latent heat release. A set of algorithm methodologies that are being used to estimate latent heating based on rain rate retrievals from the TRMM observations are described. The characteristics of these algorithms and the latent heating products that can be generated from them are also described, along with validation analyses of the heating products themselves. Finally, the investigation provides an overview of how TRMM-derived latent heating information is currently being used in conjunction with global weather and climate models, concluding with remarks intended to stimulate further research on latent heating retrieval from satellites.

Description: Rainfall is the fundamental variable within the Earth's hydrological cycle because it is both the main forcing term leading to variations in continental and oceanic surface water budgets. The vertical distribution of latent heat release, which is accompanied with rain, modulates large-scale meridional and zonal circulations within the tropics as well as modifying the energetic efficiency of mid-latitude weather systems. Latent heat release itself is a consequence of phase changes between the vapor, liquid, and frozen states of water.This paper focuses on the retrieval of latent heat release from satellite measurements generated by the Tropical Rainfall Measuring Mission 0. The TRMM observatory, whose development was a joint US-Japan space endeavor, was launched in November 1997. TRMM measurements provide an accurate account of rainfall over the global tropics, information which can be .used to estimate the four-dimensional structure of latent heating across the entire tropical and sub-tropical regions. Various algorithm methodologies for estimating latent heating based on rain rate measurements from TRMM observations are described. The strengths and weaknesses of these algorithms, as well as the latent heating products generated by these algorithms, are also discussed along with validation analyses of the products. The investigation paper provides an overview of how TRMM-derived latent heating information is currently being used in conjunction with global weather and climate models, and concludes with remarks designed to stimulate further research on latent heating retrieval

Description: A high-resolution finite-volume General Circulation Model (fvGCM), resulting from a development effort of more than ten years, is now being run operationally at the NASA Goddard Space Flight Center and Ames Research Center. The model is based on a finite-volume dynamical core with terrain-following Lagrangian control-volume discretization and performs efficiently on massive parallel architectures. The computational efficiency allows simulations at a resolution of a quarter of a degree, which is double the resolution currently adopted by most global models in operational weather centers. Such fine global resolution brings us closer to overcoming a fundamental barrier in global atmospheric modeling for both weather and climate, because tropical cyclones and even tropical convective clusters can be more realistically represented. In this work, preliminary results of the fvGCM are shown. Fifteen simulations of four Atlantic tropical cyclones in 2002 and 2004 are chosen because of strong and varied difficulties presented to numerical weather forecasting. It is shown that the fvGCM, run at the resolution of a quarter of a degree, can produce very good forecasts of these tropical systems, adequately resolving problems like erratic track, abrupt recurvature, intense extratropical transition, multiple landfall and reintensification, and interaction among vortices.

Description: Rainfall is a key link in the hydrologic cycle and is a primary heat source for the atmosphere. The vertical distribution of latent-heat release, which is accompanied by rainfall, modulates the large-scale circulations of the tropics and in turn can impact midlatitude weather. This latent heat release is a consequence of phase changes between vapor, liquid, and solid water. The Tropical Rainfall Measuring Mission (TRMM), a joint U.S./Japan space project, was launched in November 1997. It provides an accurate measurement of rainfall over the global tropics which can be used to estimate the four-dimensional structure of latent heating over the global tropics. The distributions of rainfall and inferred heating can be used to advance our understanding of the global energy and water cycle. This paper describes several different algorithms for estimating latent heating using TRMM observations. The strengths and weaknesses of each algorithm as well as the heating products are also discussed. The validation of heating products will be exhibited. Finally, the application of this heating information to global circulation and climate models is presented.

Description: A short history of the development of the TRMM mission will be described focusing on the development of the science requirements, the close collaboration between the U.S. and Japan, and the technological and scientific success of the mission. The initial TRMM science workshop was held in 1986 and with the science requirements derived from that start a set of instruments were developed to measure tropical precipitation for the first mission with that parameter as the principal focus. Implementation of the mission concept was carried out by through close cooperation between U.S. and Japanese scientists, engineers and managers. The Tropical Rainfall Measuring Mission (TRMM) has now completed three years in orbit. A short summary of research highlights will be presented and these successes will be traced back to the scientific and programmatic approaches used to develop and manage the mission. Plans for a Global Precipitation Mission (GPM) based on the successful approach used for TRMM will also be outlined.

Description: The South China Sea Monsoon Experiment (SCSMEX) was conducted in May-June 1998. One of its major objectives is to better understand the key physical processes for the onset and evolution of the summer monsoon over Southeast Asia and southern China (Lau et al. 2000). Multiple observation platforms (e.g., soundings, Doppler radar, ships, wind seafarers, radiometers, etc.) during SCSMEX provided a first attempt at investigating the detailed characteristics of convection and circulation changes, associated with monsoons over the South China Sea region. SCSMEX also provided precipitation derived from atmospheric budgets (Johnson and Ciesielski 2002) and comparison to those obtained from the Tropical Rainfall Measuring Mission (TRMM). In this paper, a regional climate model and a cloud-resolving model are used to perform multi-day integrations to understand the precipitation processes associated with the summer monsoon over Southeast Asia and southern China. The regional climate model is used to understand the soil - precipitation interaction and feedback associated with a flood event that occurred in and around China's Atlantic River during SCSMEX. Sensitivity tests on various land surface models, cumulus parameterization schemes (CASE), sea surface temperature (SST) variations and midlatitude influences are also performed to understand the processes associated with the onset of the monsoon over the S. China Sea during SCSMEX. Cloud-resolving models (CRMs) use more sophisticated and physically realistic parameterizations of cloud microphysical processes with very fine spatial and temporal resolution. One of the major characteristics of CRMs is an explicit interaction between clouds, radiation and the land/ocean surface. It is for this reason that GEWEX (Global Energy and Water Cycle Experiment) has formed the GCSS (GEWEX Cloud System Study) expressly for the purpose of improving the representation of the moist processes in large-scale models using CRMs. The Goddard Cumulus Ensemble (GCE) model is a CRM and is used to simulate convective systems associated with the onset of the South China Sea monsoon in 1998. The BRUCE model includes the same land surface model, cloud physics, and radiation scheme used in the regional climate model. A comparison between the results from the GCE model and regional climate model is performed.

Description: We discuss in this paper some of the problems that exist today in the fall utilization of satellite data to improve weather forecasts and we propose specific recommendations to solve them. This discussion can be viewed as an aspect of the general debate on how best to organize the transition from research to operational satellites and how to evaluate the impact of a research instrument on numerical weather predictions. A method for providing this transition is offered by the National Polar-Orbiting Operational Environmental Satellite System (NPOESS) Preparatory Project (NPP). This mission will bridge the time between the present NOAA and Department of Defense (DOD) polar orbiting missions and the initiation of the converged NPOESS series and will evaluate some of the Earth Observing System (EOS) instruments as appropriate for operational missions. Thus, this mission can be viewed as an effort to meet the operational requirements of NOAA and DOD and the research requirements of NASA. More generally, however, it can be said that the process of going from the conception of new, more advanced instruments to their operational implementation and full utilization by the weather forecast communities is not optimal. Instruments developed for research purposes may have insufficient funding to explore their potential operational capabilities. Furthermore, instrument development programs designed for operational satellites typically have insufficient funding for assimilation algorithms needed to transform the satellite observations into data that can be used by sophisticated global weather forecast models. As a result, years often go by before satellite data are efficiently used for operational forecasts. NASA and NOAA each have unique expertise in the design of satellite instruments, their use for basic and applied research and their utilization in weather and climate research. At a time of limited resources, the two agencies must combine their efforts to work toward common goals of full utilization of satellite data. This is a challenge that requires the assimilation of myriad new data into increasingly sophisticated numerical forecast models that run on increasingly sophisticated computer systems. In section II, we briefly outline the impact of satellite data on the quality of the National Centers for Environmental Prediction (NCEP) forecasts. In section III, we describe the present status of the utilization of satellite data in NCEP models and the challenges that lie ahead. In section IV, we propose solutions whose goals are summarized in section V.

Description: The Intergovernmental Panel on Climate Change estimated a peak increase in ozone ranging from 7-12 ppbv (zonal and annual average, and relative to a baseline with no aircraft), due to the subsonic aircraft in the year 2015, corresponding to aircraft emissions of 1.3 TgN/year. This range of values presumably reflects differences in model input (e.g., chemical mechanism, ground emission fluxes, and meteorological fields), and algorithms. The model implemented by the Global Modeling Initiative allows testing the impact of individual model components on the assessment calculations. We present results of the impact of doubling the 1995 aircraft emissions of NOx, corresponding to an extra 0.56 TgN/year, utilizing meteorological data from NASA's Data Assimilation Office (DAO), the Goddard Institute for Space Studies (GISS), and the Middle Atmosphere Community Climate Model, version 3 (MACCM3). Comparison of results to observations can be used to assess the model performance. Peak ozone perturbations ranging from 1.7 to 2.2 ppbv of ozone are calculated using the different fields. These correspond to increases in total tropospheric ozone ranging from 3.3 to 4.1 Tg/Os. These perturbations are consistent with the IPCC results, due to the difference in aircraft emissions. However, the range of values calculated is much smaller than in IPCC.

Description: It is generally thought that, in order to compensate for lower solar flux and maintain liquid oceans on the early Earth, methane must have been an important greenhouse gas before approximately 2.2 billion years (Gyr) ago. This is based upon a simple thermodynamic calculation that relates the absence of siderite (FeCO3) in some pre-2.2-Gyr palaeosols to atmospheric CO2 concentrations that would have been too low to have provided the necessary greenhouse effect. Using multi-dimensional thermodynamic analyses and geological evidence, we show here that the absence of siderite in palaeosols does not constrain atmospheric CO2 concentrations. Siderite is absent in many palaeosols (both pre- and post-2.2-Gyr in age) because the O2 concentrations and pH conditions in well-aerated soils have favoured the formation of ferric (Fe3+)-rich minerals, such as goethite, rather than siderite. Siderite, however, has formed throughout geological history in subsurface environments, such as euxinic seas, where anaerobic organisms created H2-rich conditions. The abundance of large, massive siderite-rich beds in pre-1.8-Gyr sedimentary sequences and their carbon isotope ratios indicate that the atmospheric CO2 concentration was more than 100 times greater than today, causing the rain and ocean waters to be more acidic than today. We therefore conclude that CO2 alone (without a significant contribution from methane) could have provided the necessary greenhouse effect to maintain liquid oceans on the early Earth.

Description: The goal of this paper is to present an overview of the Convection and Moisture Experiment including CAMEX-3 and CAMEX-4 field campaigns including field operations, aircraft platforms and instrumentation, aircraft missions, and data acquired during 1998 and 2001 field phases. A total of eight tropical storms and hurricanes were investigated during the CAMEX field campaigns including Bonnie, Danielle, Earl, and Georges during 1998 and Chantal, Erin, Gabrielle, and Humberto during 2001. Most of these storms were sampled with aircraft over the open ocean, but Bonnie (1998), Georges (1998), and Gabrielle (2001) also provided opportunities to monitor landfalling impacts. A few of the storms were sampled on multiple occasions during a course of several days. Most notably of these was Hurricane Humberto, which was sampled on three consecutive days during a cycle of both increasing and decreasing intensity change. Information collected for each of the eight CAMEX tropical storms as well the TRMM validation activities have been archived and are readily available for distribution at the CAMEX web site.

Description: The lack of an adequate ancient analogue for future climates means that we ultimately must use and trust climate models, evaluated against modern observation and our best geologic records of warm and cold climates of the past. Armed with an elevated confidence in the models, we will then be able to make reliable predictions of the Earth's response to our risky experiment with the climate system.

Description: Over the past decade, the Goddard Mesoscale Modeling and Dynamics Group has used a popular regional scale model, MM5, to study precipitation processes. Our group is making contributions to the MM5 by incorporating the following physical and numerical packages: improved Goddard cloud processes, a land processes model (Parameterization for Land-Atmosphere-Cloud Exchange - PLACE), efficient but sophisticated radiative processes, conservation of hydrometeor mass (water budget), four-dimensional data assimilation for rainfall, and better computational methods for trace gas transport. At NASA Goddard, the MM5 has been used to study: (1) the impact of initial conditions, assimilation of satellite-derived rainfall, and cumulus parameterizations on rapidly intensifying oceanic cyclones, hurricanes and typhoons, (2) the dynamic and thermodynamic processes associated with the development of narrow cold frontal rainbands, (3) regional climate and water cycles, (4) the impact of vertical transport by clouds and lightning on trace gas distributiodproduction associated with South and North American mesoscale convective systems, (5) the development of a westerly wind burst (WWB) that occurred during the TOGA COARE and the diurnal variation of precipitation in the tropics, (6) a Florida sea breeze convective event and a Mid-US flood event using a sophisticated land surface model, (7) the influence of soil heterogeneity on land surface energy balance in the southwest GCIP region, (8) explicit simulations (with 1.33 to 4 km horizontal resolution) of hurricanes Bob (1991) and Bonnie (1998), (9) a heavy precipitation event over Taiwan, and (10) to make real time forecasts for a major NASA field program. In this paper, the modifications and simulated cases will be described and discussed.

Description: The NASA Short-term Prediction Research and Transition (SPoRT) Center in Huntsville, Alabama was created to accelerate the infusion of NASA earth science observations, data assimilation and modeling research into NWS forecast operations and decision-making. The principal focus of experimental products is on the regional scale with an emphasis on forecast improvements on a time scale of 0-24 hours. The SPoRT Center research is aligned with the regional prediction objectives of the US Weather Research Program dealing with 0-1 day forecast issues ranging from convective initiation to 24-hr quantitative precipitation forecasting. The SPoRT Center, together with its other interagency partners, universities, and the NASA/NOAA Joint Center for Satellite Data Assimilation, provides a means and a process to effectively transition NASA Earth Science Enterprise observations and technology to National Weather Service operations and decision makers at both the global/national and regional scales. This paper describes the process for the transition of experimental products into forecast operations, current products undergoing assessment by forecasters, and plans for the future.

Description: Building upon its successful collaboration with NASA, Cybernet introduced GestureStorm(TradeMark), a weather map management system that utilizes both body tracking and gesture recognition technology for televised weather reports. The GestureStorm software interface enables meteorologists to control the computerized visual effects on weather maps using hand gestures and body movements. Prior to this, meteorologists prepared scripted on-air reports to coincide with the timing of graphics on the computerized maps. They knew when and where weather and Doppler radar maps would appear on the monitor, but could not react in real time. Cybernet s new product frees forecasters from a script, enabling them to directly interact with the weather picture in real time, significantly reducing the preparation time for the broadcast. With this technology, meteorologists control the pace of the visuals and can incorporate spontaneous close-ups with simple hand movements. In cases of severe weather, meteorologists strive to report new information as quickly as possible. GestureStorm gives forecasters an edge as they track late-breaking storms. The reduced preparation and real-time interaction shave critical minutes from the time required to broadcast severe weather warnings.

Description: Since its founding in 1992, Global Science & Technology, Inc. (GST), of Greenbelt, Maryland, has been developing technologies and providing services in support of NASA scientific research. GST specialties include scientific analysis, science data and information systems, data visualization, communications, networking and Web technologies, computer science, and software system engineering. As a longtime contractor to Goddard Space Flight Center s Earth Science Directorate, GST scientific, engineering, and information technology staff have extensive qualifications with the synthesis of satellite, in situ, and Earth science data for weather- and climate-related projects. GST s experience in this arena is end-to-end, from building satellite ground receiving systems and science data systems, to product generation and research and analysis.

Description: The objective of this study was to conduct measurements of chemical kinetics parameters for reactions of importance in the stratosphere and upper troposphere, and to study the interaction of trace gases with ice surfaces in order to elucidate the mechanism of heterogeneous chlorine activation processes, using both a theoretical and an experimental approach. The measurements were carried out under temperature and pressure conditions covering those applicable to the stratosphere and upper troposphere. The main experimental technique employed was turbulent flow-chemical ionization mass spectrometry, which is particularly well suited for investigations of radical-radical reactions.

Description: The objective of this study is to conduct measurements of chemical kinetics parameters for heterogeneous reactions of importance in the stratosphere and the troposphere. It involves the elucidation of the mechanism of the interaction of HC1 vapor with ice surfaces, which is the first step in the heterogeneous chlorine activation processes, as well as the investigation of the atmospheric oxidation mechanism of soot particles emitted by biomass and fossil fuels. The techniques being employed include turbulent flow- chemical ionization mass spectrometry and optical ellipsometry, among others. The next section summarizes our research activities during the first year of the project, and the section that follows consists of the statement of work for the second year.

Description: QuikSCAT backscatter and DMSP SSM/I radiance data are used to derive sea ice motion for both the Arctic and Antarctic region using wavelet analysis method. This technique provides improved spatial coverage over the existing array of Arctic Ocean buoys and better temporal resolution over techniques utilizing satellite data from Synthetic Aperture Radar (SAR). Sea ice motion of the Arctic for the period from October 1999 to March 2000 derived from QuikSCAT and SSM/I data agrees well with that derived from ocean buoys quantitatively. Thus the ice tracking results from QuikSCAT and SSM/I are complement to each other, Then, three sea-ice drift daily results from QuikSCAT, SSM/I, and buoy data can be merged to generate composite maps with more complete coverage of sea ice motion than those from single data source. A series of composite sea ice motion maps for December 1999 show that the major circulation patterns of sea ice motion are changing and shifting significantly within every four days and they are dominated by wind forcing. Sea-ice drift in the summer can not be derived from NSCAT and SSM/I data. In later summer of 1999 (in September), however, QuikSCAT data can provide good sea ice motion information in the Arctic. QuiksCAT can also provide at least partial sea ice motion information until June 15 in early summer 1999. For the Antarctic, case study shows that sea ice motion derived from QuikSCAT data is predominantly forced by and is consistent with wind field derived from QuikSCAT around the polar region. These calibrated/validated results indicate that QuikSCAT, SSM/I, and buoy merged daily ice motion are suitably accurate to identify and closely locate sea ice processes, and to improve our current knowledge of sea ice drift and related processes through the data assimilation of ocean-ice numerical model.

Description: Understanding the climatically varying properties of water vapor transports from a robust observational perspective is an essential step in calibrating climate models. This is tantamount to measuring year-to-year changes of monthly- or seasonally-averaged, divergent water vapor transport distributions. This cannot be done effectively with conventional radiosonde data over ocean regions where sounding data are generally sparse. This talk describes how a methodology designed to derive atmospheric water vapor transports over the world oceans from satellite-retrieved precipitation (P) and evaporation (E) datasets circumvents the problem of inadequate sampling. Ultimately, the method is intended to take advantage of the relatively complete and consistent coverage, as well as continuity in sampling, associated with E and P datasets obtained from satellite measurements. Independent P and E retrievals from Special Sensor Microwave Imager (SSM/I) measurements, along with P retrievals from Tropical Rainfall Measuring Mission (TRMM) measurements, are used to obtain transports by solving a potential function for the divergence of water vapor transport as balanced by large scale E - P conditions.

Description: The 1993 US Midwest Flood produced record levels of flooding, in the Mississippi River Basin. This flooding resulted from repeated frontal passages and mesoscale convective complexes (MCCs) during the months of June and July. A better understanding of processes that influenced MCC development during the 1993 Flood may lead to improved forecasts of heavy precipitation and flooding. Here, we consider the impact of soil moisture on MCC development during a two-day period (June 23-24) of the 1993 US Midwest Flood. The purpose of this study is to assess the importance of soil moisture distribution on the timing, intensity, and location of heavy precipitation. In this study, the MM5-PLACE Atmosphere/Land-Surface Model is utilized. The atmospheric component consists of the Penn State/NCAR MM5 mesoscale model, and the land-surface component consists of the Goddard Parameterization for Land Atmosphere-Cloud Exchange (PLACE). Initial soil moisture is provided from two sources: 1) NCEP reanalysis, and 2) Antecedent Precipitation Index (API) using NOAA rain gauge measurements as a proxy for soil moisture. NCEP reanalysis provides coarse resolution initial soil moisture (2.5 degree), while API provides high resolution initial soil moisture (10-200 km depending on NOAA rain gauge spacing). Initial results indicate that the distribution of soil moisture has a significant impact on the timing and location of heavy precipitation during this two-day flood event. Precipitation in simulations with high resolution initial soil moisture agrees more closely with observed precipitation. These results suggest that high resolution soil moisture observations are necessary to accurately predict severe storm development, heavy precipitation, and subsequent flooding.

Description: Ship tracks represent a natural laboratory to study the effects of aerosols on clouds. A number of observations and simulations have shown that increased droplet concentrations in ship tracks increase their total cross-sectional area, thereby enhancing cloud albedo and providing a negative radiative forcing at the surface and the top of the atmosphere. In some cases, cloud water has been found to be enhanced in ship tracks, which has been attributed to suppression of drizzle and implies an enhanced susceptibility of cloud albedo to droplet concentrations. However, more recently compiled observations indicate that cloud water is instead reduced in daytime ship tracks on average. Such a response is consistent with cloud-burning due to solar absorption by soot (the semi-direct radiative forcing of aerosols), recently suggested to be suppressing trade cumulus cloud coverage over the Indian Ocean. We will summarize observational evidence and present large-eddy simulations that consider these competing mechanisms in the effects of aerosols on cloud albedo.

Description: Simultaneous measurement of the radiation and cloud fields on a global basis is a key component in the effort to understand and model the interaction between clouds and radiation at the top of the atmosphere, at the surface, and within the atmosphere. The NASA Clouds and Earth s Radiant Energy System (CERES) Project, begun in 1998, is meeting this need. Broadband shortwave (SW) and longwave radiance measurements taken by the CERES scanners at resolutions between 10 and 20 km on the Tropical Rainfall Measuring Mission (TRMM), Terra, and Aqua satellites are matched to simultaneous retrievals of cloud height, phase, particle size, water path, and optical depth OD from the TRMM Visible Infrared Scanner (VIRS) and the Moderate Resolution Imaging Spectroradiometer (MODIS) on Terra and Aqua. Besides aiding the interpretation of the broadband radiances, the CERES cloud properties are valuable for understanding cloud variations at a variety of scales. In this paper, the resulting CERES cloud data taken to date are averaged at several temporal scales to examine the temporal and spatial variability of the cloud properties on a global scale at a 1 resolution.

Description: The Goddard Earth Observing System Data Assimilation System (GEOS-DAS) timeseries is a globally gridded atmospheric data set for use in climate research. This near real-time data set is produced by the Data Assimilation Office (DAO) at the NASA Goddard Space Flight Center in direct support of the operational EOS instrument product generation from the Terra (12/1999 launch), Aqua (05/2002 launch) and Aura (01/2004 launch) spacecrafts. The data is archived in the EOS Core System (ECS) at the Goddard Earth Sciences Data and Information Services Center/Distributed Active Archive Center (GES DISC DAAC). The data is only a selection of the products available from the GEOS-DAS. The data is organized chronologically in timeseries format to facilitate the computation of statistics. GEOS-DAS data will be available for the time period January 1, 2000, through present.

Description: The application of hyper spectral radiometric data to climate research requires very high absolute radiometric accuracy and stability. We use cloud-free tropical ocean data from the Atmospheric InfraRed Sounder (AIR) Calibration Data Subset (ADCS) to show that the radiometric precision and stability required climate applications has been achieved. The sea surface skin temperatures derived from the AIRS 2616cm-1 super window channel are stable relative to the RTG.SST at the better than 8 mK/year level, and the spectral calibration is stable at the 1 ppm/year level. The excellent stability and accuracy are the result of the implementation of AIRS as a grating array spectrometer, which is cooled and stabilized within 10 mK at 155 K. Analysis of daily measurements of the temperature gradient between the surface and 7 km altitude show that the AIRS Calibration Data Subset has applications which extend its original intent for calibration support to climate research. The Atmospheric Infrared Sounder (AIRS) on the EOS Aqua satellite was launched into polar orbit in May 2002. AIRS covers the spectral region from 640 to 2700 cm-1 with 2378 independent channels and represents the first of a new generation of hyper spectral resolution sounders in support of global sounding data for weather forecasting and climate research.

Description: Microwave brightness temperatures at 18.7,36.5, and 89 GHz collected at the Local-Scale Observation Site (LSOS) of the NASA Cold-Land Processes Field Experiment in February, 2003 (third Intensive Observation Period) were simulated using a Dense Media Radiative Transfer model (DMRT), based on the Quasi Crystalline Approximation with Coherent Potential (QCA-CP). Inputs to the model were averaged from LSOS snow pit measurements, although different averages were used for the lower frequencies vs. the highest one, due to the different penetration depths and to the stratigraphy of the snowpack. Mean snow particle radius was computed as a best-fit parameter. Results show that the model was able to reproduce satisfactorily brightness temperatures measured by the University of Tokyo s Ground Based Microwave Radiometer system (CBMR-7). The values of the best-fit snow particle radii were found to fall within the range of values obtained by averaging the field-measured mean particle sizes for the three classes of Small, Medium and Large grain sizes measured at the LSOS site.

Description: The project focused on ice-ocean model development and in particular on the assimilation of ice motion data and ice concentration data into both regional and global models. Many of the resulting publications below deal with improvements made in the physics treated by the model and the procedures for assimilating data. Several papers examine how the ability of the model to simulate the past behavior of the ice cover, especially to represent the ice thickness and ice deformation, is improved by data assimilation. A second aspect of the work involved interpretation of modeled behavior. Resulting papers treat the decline of arctic ice thickness over the last thirty years, and how that decline was caused by a slight warming of the near-surface atmosphere, and also how large variation in ice thickness are due to changes in wind patterns associated with a well- known oscillation of the atmospheric circulation. The research resulted in over 20 published papers on these topics.

Description: This report summarizes our efforts on the funded project 'Quality Control and Analysis of Microphysical Data Collected in TRMM Airborne Validation Experiments', NASA NAG5-9663, Andrew Heymsfield, P. I. We begin this report by summarizing our activities in FY2000-FY2004. We then present some highlights of our work. The last part of the report lists the publications that have resulted from our funding through this grant.

Description: The Moderate Resolution Imaging Spectroradiometer (MODIS) on the NASA Earth Observing System (EOS) Terra and Aqua missions has shown considerable capability for mapping snowcover. The typical approach that has used, along with other criteria, the Normalized Snow Difference Index (NDSI) that takes the difference between 500 meter observations at 1.64 micrometers (MODIS band 6) and 0.555 micrometers (MODIS band 4) over the sum of these observations to determine whether MODIS pixels are snowcovered or not in mapping the extent of snowcover. For many hydrological and climate studies using remote sensing of snowcover, it is desirable to assess if the MODIS snowcover observations could not be enhanced by providing the fraction of snowcover in each MODIS observation (pixel). Pursuant to this objective studies have been conducted to assess whether there is sufficient "signal%o in the NDSI parameter to provide useful estimates of fractional snowcover in each MODIS 500 meter pixel. To accomplish this objective high spatial resolution (30 meter) Landsat snowcover observations were used and co-registered with MODIS 500 meter pixels. The NDSI approach was used to assess whether a Landsat pixel was or was not snowcovered. Then the number of snowcovered Landsat pixels within a MODIS pixel was used to determine the fraction of snowcover within each MODIS pixel. The e results were then used to develop statistical relationships between the NDSI value for each 500 meter MODIS pixel and the fraction of snowcover in the MODIS pixel. Such studies were conducted for three widely different areas covered by Landsat scenes in Alaska, Russia, and the Quebec Province in Canada. The statistical relationships indicate that a 10 percent accuracy can be attained. The variability in the statistical relationship for the three areas was found to be remarkably similar (-0.02 for mean error and less than 0.01 for mean absolute error and standard deviation). Independent tests of the relationships were accomplished by taking the relationship of fractional snow-cover to NDSI from one area (e.g., Alaska) and testing it on the other two areas (e.g. Russia and Quebec). Again the results showed that fractional snow-cover can be estimated to 10 percent. The results have been shown to have advantages over other published fractional snowcover algorithms applied to MODIS data. Most recently the fractional snow-cover algorithm has been applied using 500-meter observations over the state of Colorado for a period spanning 25 days. The results exhibit good behavior in mapping the spatial and temporal variability in snowcover over that 25-day period. Overall these studies indicate that robust estimates of fractional snow-cover can be attained using the NDSI parameter over areas extending in size from watersheds relatively large compared to MODIS pixels to global land cover. Other refinements to this approach as well as different approaches are being examined for mapping fractional snow-cover using MODIS observations.

Description: A taxonomy of tropical convective vertical structures is constructed through cluster analysis of three years of Tropical Rainfall Measuring Mission [TRMM] Precipitation Radar [PR] vertical profiles, their surface rainfall and associated radar-based classifiers (convective/stratiform and bright band existence). archetypal profile types are identified. These include nine convective types, divided into warm, "just cold", midlevel, deep and deep/wet-growth categories, seven stratiform types, divided into warm, "just cold", midlevel and deep categories, three "mixed" types (deep profiles with low reflectivity aloft), and six fragment types (non-precipitating anvils and sheared deep convective profiles). The taxonomy allows for description of any storm or local Convective spectrum by the nine primary convective and stratiform types, a significant reduction over full three-dimensional radar data which nonetheless retains vertical structure information. The analysis provides a quasi-independent corroboration of the TRMM 2A23 convective/stratiform classification. The global frequency of occurrence and contribution to rainfall for the profile types is presented, demonstrating primary rainfall contribution by midlevel glaciated convection and similar depth decaying/stratiform stages. Close correspondence is found between deep convective profile frequency and annualized lightning production. Passive microwave and lightning properties associated with the profiles are reported, and cases presented illustrating known nonuniqueness problems with 85 and 37 GHz brightness temperature pairs (the same pairs corresponding to both convective and stratiform profiles), and how supplementary lightning information might be used to mitigate these problems.

Description: NASA's Aura satellite is scheduled to launch in the second quarter of 2004 into a polar orbit. The Aura mission is designed to collect data to address three high priority environmental science questions: (1) Is the ozone layer recovering as expected? (2) What are the sources and processes that control tropospheric pollutants? And (3) what is the quantitative impact of constituents on climate change? Aura will answer these questions by globally measuring a comprehensive set of trace gases and aerosols in the troposphere and stratosphere. Aura data will also have applications for monitoring and predicting climate and air quality parameters. Aura s observations will continue the TOMS ozone trend record and provide an assessment as to whether the Montreal Protocol is achieving its objective. Aura will measure gases and aerosols in the upper troposphere and lower stratosphere that contribute to climate forcing. These data will be of sufficient coverage, vertical resolution, and accuracy to help constrain climate models. In addition, Aura observations of tropospheric ozone and its precursors will have regional as well as intercontinental coverage, which could improve emission inventories. Near real time data will tested for local air quality forecasts in collaboration with the US's Environmental Protection UV-B forecasts from Aura ozone and cloud cover data. An overview of Aura s instruments, data products, validation, and examples of data applications will be presented.

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

Description: The Global Precipitation Measurement (GPM) Mission will consist of a constellation of rain-measuring satellites, the main member of which (the core satellite) will serve as the measurement reference to the other members of the constellation. The core satellite is being developed jointly by the National Aeronautics and Space Administration (NASA) and the newly-named Japan Aerospace Exploration Agency (JAXA -- previously NASDA) along with its government partner, the Communications Research Laboratory (CRL). The GPM mission was proposed as a follow-up mission to the Tropical Rainfall Measuring Mission (TRMM) by both NASA and NASDA based on the unparalleled scientific success of TRMM, and has recently been joined by the European Space Agency (ESA) via its formulation of the European GPM mission (i.e., EGPM). GPM is an ambitious mission designed to produce accurate and frequent global observations of precipitation (both rain and snow) made possible by replacing the TRMM satellite with the new core satellite carrying an advanced radar-radiometer system, and serving as the centerpiece for the constellation of some eight (8) additional satellites being provided through international cooperation. The core satellite is to be flown up to high latitudes (inclined some 65-70 degrees), and will carry a Ku/Ka-band, nadir-scanning, dual-frequency precipitation radar (DPR) that is being developed by JAXA and CRL, along with a large aperture, extended frequency-range, conically-scanning passive microwave radiometer being developed by NASA and its industrial partners. Each constellation satellite will also carry some type of multi-channel passive microwave radiometer (as well as a multi-beam Ka-band radar in the case of EGPM) whose rain estimates will be calibrated and referenced to those made by the core satellite, producing for the first time fully-global, continuous, and bias-free precipitation datasets. GPM data will be delivered in near-realtime, taking a major step toward the operational use of precipitation information for model initialization and data assimilation in a number of application areas such as hazardous weather forecasting, flood warning, fresh water resource assessment, and crop growth prediction. In addition, GPM data will complement the now-existing global temperature record, allowing for improved assessments of climate change, particularly those processes in which the global water cycle both forces and responds to climatic drifts in global temperature conditions. A foremost element of this international constellation mission is a parallel international ground validation (GV) network. This GV network is needed to determine uncertainties in the rain retrievals, critical for application of the retrieval information in weather and hydrometeorological modeling and climate diagnostics, as well as assurances that the satellite retrievals of surface rainfall are consistent with those actually measured at the surface. The key aspects of this network is that it must be worldwide and created through the GPM partnership process. Therefore the network will consist of a confederation of government agencies, academic organizations, private institutions, and individual scientists from a collection of nations who have initiated the process by gathering in Abingdon to develop the fundamentals of the international GPM GV research programme. Therefore in keeping with our responsibilities as the front-runners of the programme, the main objectives of this workshop are: (1) to present and share opinions on interests, perspectives, and concerns about GPM GV research; (2) to examine the conceptual and/or planned GPM GV site templates from NASA, NASDA, ESA, and other partners; (3) to define the main scientific objectives of the international GPM GV research programme; (4) to formulate a preliminary set of international GPM GV science and measurement requirements; and (5) to convene a Steering Committee to aid the organization of the GPM GV program, to document its science implementation plans, and to aid planning for follow-up GPM GV meetings.

Description: The Earth Observing System (EOS) is a space-based observing system comprised of a series of satellite sensors by which scientists can monitor the Earth, a Data and Information System (EOSDIS) enabling researchers worldwide to access the satellite data, and an interdisciplinary science research program to interpret the satellite data. The Moderate Resolution Imaging Spectroradiometer (MODIS), developed as part of the Earth Observing System (EOS) and launched on Terra in December 1999 and Aqua in May 2002, is designed to meet the scientific needs for satellite remote sensing of clouds, aerosols, water vapor, and land and ocean surface properties. This sensor and multi-platform observing system is especially well suited to observing detailed interdisciplinary components of the Earth s surface and atmosphere in and around urban environments, including aerosol optical properties, cloud optical and microphysical properties of both liquid water and ice clouds, land surface reflectance, fire occurrence, and many other properties that influence the urban environment and are influenced by them. In this presentation I will summarize the current capabilities of MODIS and other EOS sensors currently in orbit to study human-induced climate variations.

Description: Over the past twenty years, rainfall retrieval algorithms have been developed to retrieve rainfall and vertical hydrometeor structures from passive microwave observations by making use of the fact that weighting functions for various frequencies peak at different levels within a rainy atmosphere. GPROF is one of two TMI rainfall algorithms. It is physically based retrieval that finds the vertical hydrometeor profile that best fits the brightness temperatures in the available passive radiometer channels. Matching is achieved using a library of hydrometeor profiles generated by cloud-resolving models (CRMs). The hydrometeor profiles have a corresponding surface precipitation rate. The algorithm retrieves the hydrometeor profiles and associated surface rainfall using a Bayesian approach that gives the estimated expected values. The ability of CRMs to produce cloud structures that are reliable and representative of observed storms is crucial for the success of GPROF. The cloud mycrophysics are one of the keys to achieving this. In addition, CRMs have been a very useful tool for GPM-algorithm developers through Cloud-Radiation Simulations (CRS), one of the nine GPM disciplinary research themes. This paper will discuss how to generate consistent and comprehensive 4D cloud datasets from an improved (i.e., in regard to bulk and multi-moment microphysics) CRM for TRMM and GPM rainfall retrieval algorithm developers. These cloud datasets include CRM-simulated clouds and cloud systems from different geographic locations in the tropics and midlatitudes. By linking the CRM with a passive microwave radiative-transfer model and using satellite and airborne data, the performance of the "cloud physics" can be assessed and in turn modified and improved. This paper will also address how to assess and improve the performance of various latent and diabatic heating algorithms and develop an algorithm to retrieve the vertical structure of apparent moistening (Q2). Considering that the GPM will produce high (temporal and spatial) resolution heating and rainfall data, these algorithms will be used to obtain the temporal and spatial distributions of surface rainfall and the associated vertical heating and moistening profiles throughout the subtropical and midlatitudes.

Description: The GPM mission is currently planned for start in the late 2007 - early 2008 time frame. Its main scientific goal is to help answer pressing scientific problems arising within the context of global and regional water cycling. These problems cut across a hierarchy of scales and include climate-water cycle interactions, techniques for improving weather and climate predictions, and better methods for combining observed precipitation with hydrometeorological prediction models for applications to hazardous flood-producing storms, seasonal flood draught conditions, and fresh water resource assessments. The GPM mission will expand the scope of precipitation measurement through the use of a constellation of some 9 satellites, one of which will be an advanced TRMM-like core satellite carrying a dual-frequency Ku-Ka band precipitation radar and an advanced, multifrequency passive microwave radiometer with vertical-horizontal polarization discrimination. The other constellation members will include new dedicated satellites and co-existing operational/research satellites carrying similar (but not identical) passive microwave radiometers. The goal of the constellation is to achieve approximately 3-hour sampling at any spot on the globe -- continuously. The constellation's orbit architecture will consist of a mix of sun-synchronous and non-sun-synchronous satellites with the core satellite providing measurements of cloud-precipitation microphysical processes plus calibration-quality rainrate retrievals to be used with the other retrieval information to ensure bias-free constellation coverage. A major requirement before the retrieved rainfall information generated by the GPM mission can be used effectively by prognostic models to improve weather forecasts, hydrometeorological forecasts, and climate model reanalysis simulations is a capability to quantify the error characteristics of the retrievals. A solution for this problem has been upheld in past precipitation missions because of the lack of suitable error modeling systems incorporated into the validation programs and data distribution systems. An overview of how NASA intends to overcome this problem for the GPM mission using a physically-based error modeling approach within a multi-faceted validation program is described. The solution is to first identify specific user requirements and then determine the most stringent of these requirements that embodies all essential error characterization information needed by the entire user community. In the context of NASA s scientific agenda for the GPM mission, the most stringent user requirement is found within the data assimilation community. The fundamental theory of data assimilation vis-a-vis ingesting satellite precipitation information into the pre-forecast initializations is based on quantifying the conditional bias and precision errors of individual rain retrievals, and the space-time structure of the precision error (i.e., the spatial-temporal error covariance). By generating the hardware and software capability to produce this information in a near real-time fashion, and to couple the derived quantitative error properties to the actual retrieved rainrates, all key validation users can be satisfied. The talk will describe the essential components of the hardware and software systems needed to generate such near real-time error properties, as well as the various paradigm shifts needed within the validation community to produce a validation program relevant to the precipitation user community.

Description: In situ measurements of water vapor concentration made during the CRYSTAL-FACE and Pre-AVE missions indicate higher than expected supersaturations in both clear and cloudy air near the cold tropical tropopause: (1) steady-state ice supersaturations of 20-30% were measured within cirrus at T 〈 200 K; (2) supersaturations exceeding 100% (near water saturation) were observed under cloud-free conditions near 187 K. The in-cloud measurements challenge the conventional belief that any water vapor in excess of ice saturation should be depleted by crystal growth given sufficient time. The high clear-sky supersaturations imply that thresholds for ice nucleation due to homogeneous freezing of aerosols (or any other mechanism) are much higher than those inferred from laboratory measurements. We will use simulations of Tropical Tropopause Layer (TTL) transport and cloud formation throughout the tropics to show that these effects have important implications for TTL cloud frequency and freeze-drying of air crossing the tropical tropopause cold trap.

Description: Asian dust typically originates in desert areas far from polluted urban regions. During transport, dust layers can interact with anthropogenic sulfate and soot aerosols from heavily polluted urban areas. Added to the complex effects of clouds and natural marine aerosols, dust particles reaching the marine environment can have drastically different properties than those from the source. Thus, understanding the unique temporal and spatial variations of Asian aerosols is of special importance in regional-to-global climate issues such as radiative forcing, the hydrological cycle, and primary biological productivity in the mid-Pacific Ocean. During ACE-Asia campaign, we have acquired ground- based (temporal) and satellite (spatial) measurements to infer aerosol physical/optical/radiative properties, column precipitable water amount, and surface reflectivity over this region. The inclusion of flux measurements permits the determination of aerosol radiative flux in addition to measurements of loading and optical depth. At the time of the Terra/MODIS, SeaWiFS, TOMS and other satellite overpasses, these ground-based observations can provide valuable data to compare with satellite retrievals over land. In this paper, we will demonstrate new capability of the Deep Blue algorithm to track the evolution of the Asian dust storm from sources to sinks. Although there are large areas often covered by clouds in the dust season in East Asia, this algorithm is able to distinguish heavy dust from clouds over the entire regions. Examination of the retrieved daily maps of dust plumes over East Asia clearly identifies the sources contributing to the dust loading in the atmosphe~~. We have compared the satellite retrieved aerosol optical thickness to the ground-based measurements and obtained a reasonable agreement between these two. Our results also indicate that there is a large difference in the retrieved value of spectral single scattering albedo of windblown dust between different sources in East Asia.