ALBERT

Description: During the Rumphius Biohistorical Expedition (4.xi-17.xii.1990) and some additional field trips to Ambon, a total of 26 species of asteroids were collected. Seven species were new to Ambon, and four species new to Indonesia. The asteroid fauna of Ambon is now represented by 44 species. Anthenea difficilis is transferred to the genus Gymnanthenea. Taxonomical notes are given on the collected specimens for some species.

Description: A new sponge species Callyspongia (Euplacella) biru spec. nov. is described from various locations within Indonesia (NE Bali, SW & NE Sulawesi and NE Kalimantan). The species is assigned to the subgenus Euplacella of the genus Callyspongia because it has brushes of spicules on the nodes of its double-meshed ectosomal skeleton. The species is discussed and compared with relevant Callyspongia records from other Indonesian and neighbouring Indo-Pacific areas.

Description: This manuscript describes a method to class@ cirrus cloud ice particle shape using lidar depolarization measurements as a basis for segregating different particle shape regimes. Measurements from the ER-2 Cloud Physics Lidar (CPL) system during CRYSTAL-FACE provide the basis for this work. While the CPL onboard the ER-2 aircraft was providing remote sensing measurements of cirrus clouds, the Cloud Particle Imager (CPI) onboard the WB-57 aircraft was flying inside those same clouds to sample particle sizes. The results of classifying particle shapes using the CPL data are compared to the in situ measurements made using the CPI , and there is found to be good agreement between the particle shape inferred from the CPL data and that actually measured by the CPI. If proven practical, application of this technique to spaceborne observations could lead to large-scale classification of cirrus cloud particle shapes.

Description: During the ACE Asia intensive field campaign conducted in the spring of 2001 aerosol properties were measured onboard the R/V Ronald H. Brown to study the effects of the Asian aerosol on atmospheric chemistry and climate in downwind regions. Aerosol properties measured in the marine boundary layer included chemical composition; number size distribution; and light scattering, hemispheric backscattering, and absorption coefficients. In addition, optical depth and vertical profiles of aerosol 180 deg backscatter were measured. Aerosol within the ACE Asia study region was found to be a complex mixture resulting from marine, pollution, volcanic, and dust sources. Presented here as a function of air mass source region are the mass fractions of the dominant aerosol chemical components, the fraction of the scattering measured at the surface due to each component, mass scattering efficiencies of the individual components, aerosol scattering and absorption coefficients, single scattering albedo, Angstrom exponents, optical depth, and vertical profiles of aerosol extinction. All results except aerosol optical depth and the vertical profiles of aerosol extinction are reported at a relative humidity of 55 +/- 5%. An over-determined data set was collected so that measured and calculated aerosol properties could be compared, internal consistency in the data set could be assessed, and sources of uncertainty could be identified. By taking into account non-sphericity of the dust aerosol, calculated and measured aerosol mass and scattering coefficients agreed within overall experimental uncertainties. Differences between measured and calculated aerosol absorption coefficients were not within reasonable uncertainty limits, however, and may indicate the inability of Mie theory and the assumption of internally mixed homogeneous spheres to predict absorption by the ACE Asia aerosol. Mass scattering efficiencies of non-sea salt sulfate aerosol, sea salt, submicron particulate organic matter, and dust found for the ACE Asia aerosol are comparable to values estimated for ACE 1, Aerosols99, and INDOEX. Unique to the ACE Asia aerosol was the large mass fractions of dust, the dominance of dust in controlling the aerosol optical properties, and the interaction of dust with soot aerosol.

Description: The terminology East Asian summer monsoon is used to refer to the heavy rainfall in southeast China including the Yangtze River Valley starting in May and ending in August (e.g., Chen and Chang 1980, Tao and Chen 1987, Ding 1992, Chang et al. 2000a.) This rainfall region is associated with the Mei-Yu front, which extends to Japan and its neighborhood and is called Baiu there. The Mei-Yu front becomes prominent in May and has a slow northward movement. From May to July the elongated rain belt moves from the southeast coast of China to the Yangtze River Valley. The rain belt extends north-east-ward to south of Japan in May and later covers Korea also. The purpose of this note is to point out that the terminology of East Asian summer monsoon is a misnomer to refer to the portion of this rainbelt residing over East Asia, in the sense that it is not a monsoon.

Description: During the period December 23,1997 and December August 31,2004, we accomplished the development of 2 CFD codes for DNS/LES/RANS simulation of turbine cascade flows, namely LESTool and UNCLE. LESTool is a structured code making use of 5th order upwind differencing scheme and UNCLE is a second-order-accuracy unstructured code. LESTool has both Dynamic SGS and Spalart's DES models and UNCLE makes use of URANS and DES models. The current report provides a description of methodologies used in the codes.

Description: NASA's 4th Convection and Moisture Experiment (CAMEX-4) focused on Atlantic hurricanes during the 2001 hurricane season and it involved both NASA and NOAA participation. The NASA ER-2 and DC-8 aircraft were instrumented with unique remote sensing instruments to help increase the overall understanding of hurricanes. This paper is concerned about one of the storms studied, Tropical Storm Chantal, that was a weak storm which failed to intense into a hurricane. One of the practical questions of high importance is why some tropical sto~ins intensify into hurricanes, and others remain weak or die altogether. The magnitude of the difference between the horizontal winds at lower levels and upper altitudes in a tropical storm, i.e., the wind shear, is one important quantity that can affect the intensification of a tropical storm. Strong shear as was present during Tropical Storm Chantal s lifetime and it was detrimental to its intensification. The paper presents an analysis of unique aircraft observations collected from Chantal including an on-board radar, radiometers, dropsondes, and flight level measurements. These measurements have enabled us to examine the internal structure of the winds and thermal structure of Chantal. Most of the previous studies have involved intense hurricanes that overcame the effects of shear and this work has provided new insights into what prevents a weaker storm from intensifying. The storm had extremely intense thunderstorms and rainfall, yet its main circulation was confined to low levels of the atmosphere. Chantal's thermal structure was not configured properly for the storm to intensify. It is most typical that huricanes have a warm core structure where warm temperatures in upper levels of a storm s circulation help intensify surface winds and lower its central pressure. Chantal had two weaker warm layers instead of a well-defined warm core. These layers have been related to the horizontal and vertical winds and precipitation structure and have helped us learn more about why this storm didn't develop.

Description: Understanding the exchange of gases between the stratosphere and the troposphere is important for determining how pollutants enter the stratosphere and how they leave. This study does a global analysis of that the exchange of mass between the stratosphere and the troposphere. While the exchange of mass is not the same as the exchange of constituents, you can t get the constituent exchange right if you have the mass exchange wrong. Thus this kind of calculation is an important test for models which also compute trace gas transport. In this study I computed the mass exchange for two assimilated data sets and a GCM. The models all agree that amount of mass descending from the stratosphere to the troposphere in the Northern Hemisphere extra tropics is approx. 10(exp 10) kg/s averaged over a year. The value for the Southern Hemisphere by about a factor of two. ( 10(exp 10) kg of air is the amount of air in 100 km x 100 km area with a depth of 100 m - roughly the size of the D.C. metro area to a depth of 300 feet.) Most people have the idea that most of the mass enters the stratosphere through the tropics. But this study shows that almost 5 times more mass enters the stratosphere through the extra-tropics. This mass, however, is quickly recycled out again. Thus the lower most stratosphere is a mixture of upper stratospheric air and tropospheric air. This is an important result for understanding the chemistry of the lower stratosphere.

Description: The knowledge of the atmospherics phenomenon is an important part in the communication system. The principal factor that contributes to the attenuation in a Ka band communication system is the rain attenuation. We have four years of tropical region observations. The data in the tropical region was taken in Humacao, Puerto Rico. Previous data had been collected at various climate regions such as desserts, template area and sub-tropical regions. Figure 1 shows the ITU-R rain zone map for North America. Rain rates are important to the rain attenuation prediction models. The models that predict attenuation generally are of two different kinds. The first one is the regression models. By using a data set these models provide an idea of the observed attenuation and rain rates distribution in the present, past and future. The second kinds of models are physical models which use the probability density functions (PDF).

Description: The Hubble Space Telescope (HST) was launched April 24, 1990, and was deployed April 25 into low Earth orbit (LEO). It was soon discovered that the metal poles holding the solar arrays were expanding and contracting as the telescope orbited the Earth passing between the sunlight and the Earth s shadow. The expansion and contraction, although very small, was enough to cause the telescope to shake because of thermal-induced jitters, a detrimental effect when trying to take pictures millions of miles away. Therefore, the European Space Agency (ESA, the provider of the solar arrays) built new solar arrays (SA-11) that contained bi-stem thermal shields which insulated the solar array metal poles. These thermal shields were made of 2 mil thick aluminized-Teflon fluorinated ethylene propylene (FEP) rings fused together into a circular bellows shape. The new solar arrays were put on the HST during an extravehicular activity (EVA), also called an astronaut space walk, during the first servicing mission (SM1) in December 1993. An on-orbit photograph of the HST with the SA-11, and a close up of the bellows-like structure of the thermal shields is provided in Figure 1.

Description: In the polar region of the upper mesosphere, horizontal wind oscillations have been observed with periods around 10 hours. Waves with such a period are generated in our Numerical Spectral Model (NSM), and they are identified as planetary-scale inertio gravity waves (IGW). These IGWs have periods between 9 and 11 hours and appear above 60 km in the zonal mean (m = 0), as well as in zonal wavenumbers m = 1 to 4. The waves can propagate eastward and westward and have vertical wavelengths around 25 km. The amplitudes in the wind field are typically between 10 and 20 m/s and can reach 30 m/s in the westward propagating component for m = 1 at the poles. In the temperature perturbations, the wave amplitudes above 100 km are typically 5 K and as large as 10 K for m = 0 at the poles. The IGWs are intermittent but reveal systematic seasonal variations, with the largest amplitudes occurring generally in late winter and spring. In the NSM, the IGW are generated like the planetary waves (PW). They are produced apparently by the instabilities that arise in the zonal mean circulation. Relative to the PWs, however, the IGWs propagate zonally with much larger velocities, such that they are not affected much by interactions with the background zonal winds. Since the IGWs can propagate through the mesosphere without much interaction, except for viscous dissipation, one should then expect that they reach the thermosphere with significant and measurable amplitudes.

Description: Tropical deep convection and its dynamical effect on the tropopause and stratosphere are investigated using a suite of data from the Upper Atmospheric Research Satellite (UARS) Microwave Limb Sounder (MLS), including upper tropospheric humidity, cloud radiance, and gravity wave measurements.

Description: The nature of observed variations in temperature-salinity (T-S) relationship between El Nino and non-El Nino years in the pycnocline of the eastern equatorial Pacific Ocean (NINO3 region, 5(deg)S-5(deg)N, 150(deg)W-90(deg)W) is investigated using an ocean general circulation model. The origin of the subject water mass is identified using the adjoint of a simulated passive tracer. The higher salinity during El Nino is attributed to larger convergence of saltier water from the Southern Hemisphere and smaller convergence of fresher water from the Northern Hemisphere.

Description: Comparisons are performed between spatially averaged sea surface temperatures (ASST2) as derived from the second Along-Track Scanning Radiometer (ATSR-2) on board the second European Remote Sensing Satellite (ERS-2) and the NOAA-NASA Advanced Very High Resolution Radiometer (AVHRR) Oceans Pathfinder dataset (MPFSST). Difference maps, MPFSST 2 ASST2, along with the application of a simple statistical regression model to aerosol and cloud data from the Total Ozone Mapping Spectrometer ( TOMS), are used to examine the impact of possible aerosol and cloud contamination. Differences varied regionally, but the largest biases were seen off western Africa. Nighttime and daytime differences off western Africa were reduced from -0.5degrees to -0.2degreesC and from -0.1degrees to 0degreesC, respectively. Significant cloud flagging, based on the model, occurred in the Indian Ocean, the equatorial Pacific, and in the vicinity of the Gulf Stream. Comparisons of the MPFSST and the ASST2 with in situ data from the 2002 version of the World Oceanic Database (WOD02) off western Africa show larger mean differences for the MPFSST. The smallest mean differences occurred for nighttime ASST2 - WOD02 with a value of 0.0degrees +/- 0.4degreesC.

Description: A computational heat transfer design methodology was developed to study the dual-engine linear aerospike plume-induced base-heating environment during one power-pack out, in ascent flight. It includes a three-dimensional, finite volume, viscous, chemically reacting, and pressure-based computational fluid dynamics formulation, a special base-bleed boundary condition, and a three-dimensional, finite volume, and spectral-line-based weighted-sum-of-gray-gases absorption computational radiation heat transfer formulation. A separate radiation model was used for diagnostic purposes. The computational methodology was systematically benchmarked. In this study, near-base radiative heat fluxes were computed, and they compared well with those measured during static linear aerospike engine tests. The base-heating environment of 18 trajectory points selected from three power-pack out scenarios was computed. The computed asymmetric base-heating physics were analyzed. The power-pack out condition has the most impact on convective base heating when it happens early in flight. The source of its impact comes from the asymmetric and reduced base bleed.

Description: The authors use satellite data to examine the relationship between lightning and upper-level radar reflectivity. They find correlations between average flash rates and upper-level reflectivities over both land and ocean, although both flash rates and reflectivities are much lower over ocean than land. Analysis of the data using Empirical Orthogonal Functions (EOFs) shows similar EOFs for averaged lightning and reflectivity. In contrast, the EOFs of the anomalies of lightning and reflectivity have different spatial patterns; however, both have principal component time series that are correlated with the Southern Oscillation Index and, hence, El Nino. Differences in behavior of the lightning and reflectivity anomaly EOFs and principal components suggest that El Nino plays a smaller role in lightning anomaly than precipitation anomaly.

Description: This paper details a novel three-dimensional (3D) hydro-focusing micro cell sorter for micro flow cytometry applications. The unit was microfabricated by means of SU-8 3D lithography. The 3D microstructure for coaxial sheathing was designed, microfabricated, and tested. Three-dimensional hydro-focusing capability was demonstrated with an experiment to sort labeled tanned sheep erythrocytes (red blood cells). This polymer hydro-focusing microstructure is easily microfabricated and integrated with other polymer microfluidic structures.

Description: In this study, a framework is discussed to apply air/space-borne dual-wavelength radar for the estimation of characteristic parameters of hydrometeors. The focus of our study is on the Global Precipitation Measurements (GPM) precipitation radar, a dual-wavelength radar that operates at Ku (13.8 GHz) and Ka (35 GHz) bands. As the droplet size distributions (DSD) of rain are expressed as the Gamma function, a procedure is described to derive the median volume diameter (D(sub 0)) and particle number concentration (N(sub T)) of rain. The correspondences of an important quantity of dual-wavelength radar, defined as deferential frequency ratio (DFR), to the D(sub 0) in the melting region are given as a function of the distance from the 0 C isotherm. A self-consistent iterative algorithm that shows a promising to account for rain attenuation of radar and infer the DSD without use of surface reference technique (SRT) is examined by applying it to the apparent radar reflectivity profiles simulated from the DSD model and then comparing the estimates with the model (true) results. For light to moderate rain the self-consistent rain profiling approach converges to unique and correct solutions only if the same shape factors of Gamma functions are used both to generate and retrieve the rain profiles, but does not converges to the true solutions if the DSD form is not chosen correctly. To further examine the dual-wavelength techniques, the self-consistent algorithm, along with forward and backward rain profiling algorithms, is then applied to the measurements taken from the 2nd generation Precipitation Radar (PR-2) built by Jet Propulsion Laboratory. It is found that rain profiles estimated from the forward and backward approaches are not sensitive to shape factor of DSD Gamma distribution, but the self-consistent method is.

Description: A revised Bayesian algorithm for estimating surface rain rate, convective rain proportion, and latent heating/drying profiles from satellite-borne passive microwave radiometer observations over ocean backgrounds is described. The algorithm searches a large database of cloud-radiative model simulations to find cloud profiles that are radiatively consistent with a given set of microwave radiance measurements. The properties of these radiatively consistent profiles are then composited to obtain best estimates of the observed properties. The revised algorithm is supported by an expanded and more physically consistent database of cloud-radiative model simulations. The algorithm also features a better quantification of the convective and non-convective contributions to total rainfall, a new geographic database, and an improved representation of background radiances in rain-free regions. Bias and random error estimates are derived from applications of the algorithm to synthetic radiance data, based upon a subset of cloud resolving model simulations, and from the Bayesian formulation itself. Synthetic rain rate and latent heating estimates exhibit a trend of high (low) bias for low (high) retrieved values. The Bayesian estimates of random error are propagated to represent errors at coarser time and space resolutions, based upon applications of the algorithm to TRMM Microwave Imager (TMI) data. Errors in instantaneous rain rate estimates at 0.5 deg resolution range from approximately 50% at 1 mm/h to 20% at 14 mm/h. These errors represent about 70-90% of the mean random deviation between collocated passive microwave and spaceborne radar rain rate estimates. The cumulative algorithm error in TMI estimates at monthly, 2.5 deg resolution is relatively small (less than 6% at 5 mm/day) compared to the random error due to infrequent satellite temporal sampling (8-35% at the same rain rate).

Description: The water cycle is the key circuit moving water through the Earth's system. This large system, powered by energy from the sun, is a continuous exchange of moisture between the oceans, the atmosphere, and the land. Precipitation (including rain, snow, sleet, freezing rain, and hail), is the primary mechanism for transporting water from the atmosphere back to the Earth's surface and is the key physical process that links aspects of climate, weather, and the global water cycle. Global precipitation and associate cloud processes are critical for understanding the water cycle balance on a global scale and interactions with the Earth's climate system. However, unlike measurement of less dynamic and more homogenous meteorological fields such as pressure or even temperature, accurate assessment of global precipitation is particularly challenging due to its highly stochastic and rapidly changing nature. It is not uncommon to observe a broad spectrum of precipitation rates and distributions over very localized time scales. Furthermore, precipitating systems generally exhibit nonhomogeneous spatial distributions of rain rates over local to global domains.

Description: Precise repeat airborne laser surveys were conducted over the major ice caps in the Canadian Arctic Archipelago in the spring of 1995 and 2000 in order to measure elevation changes in the region. Our measurements reveal thinning at lower elevations (below 1600 m) on most of the ice caps and glaciers, but either very little change or thickening at higher elevations in the ice cap accumulation zones. Recent increases in precipitation in the area can account for the slight thickening where it was observed, but not for the thinning at lower elevations. For the northern ice caps on the Queen Elizabeth Islands, thinning was generally less than 0.5 m/yr , which is consistent with what would be expected from the warm temperature anomalies in the region for the 5-year period between surveys and appears to be a continuation of a trend that began in the mid 1980s. Further south, however, on the Barnes and Penny ice caps on Baffin Island, this thinning was much more pronounced at over 1 m/yr in the lower elevations. Here temperature anomalies were very small, and the thinning at low elevations far exceeds any associated enhanced ablation. The observations on Barnes, and perhaps Penny are consistent with the idea that the observed thinning is part of a much longer term deglaciation, as has been previously suggested for Barnes Ice Cap. Based on the regional relationships between elevation and elevation-change in our data, the 1995-2000 mass balance for the region is estimated to be 25 cu km/yr of ice, which corresponds to a sea level increase of 0.064 mm/ yr . This places it among the more significant sources of eustatic sea level rise, though not as substantial as Greenland ice sheet, Alaskan glaciers, or the Patagonian ice fields.

Description: A unique set of measurements of wind, water vapor mixing ratio and boundary layer height variability was observed during the first MOP dryline mission of 22 May 2002. Water vapor mixing ratio from the Scanning Raman Lidar (SRL), high-resolution profiles of aerosol backscatter from the HARLIE and wind profiles from the GLOW are combined with the vertical velocity derived from the NCAR/ISS/MAPR and the high-resolution FMCW radar to reveal the convective variability of the cumulus cloud-topped boundary layer. A combined analysis of the in-situ and remote sensing data from aircraft, radiosonde, lidars, and radars reveals moisture variability within boundary layer updraft and downdraft regions as well as characterizes the boundary layer height variability in the dry and moist sides of the dryline. The profiler site measurements will be tied to aircraft data to reveal the relative intensity and location of these updrafts to the dry line. This study provides unprecedented high temporal and spatial resolution measurements of wind, moisture and backscatter within a dryline and the associated convective boundary layer.

Description: Although dry/moist potential vorticity is a useful physical quantity for meteorological analysis, it cannot be applied to the analysis of 2D simulations. A convective vorticity vector (CVV) is introduced in this study to analyze 2D cloud-resolving simulation data associated with 2D tropical convection. The cloud model is forced by the vertical velocity, zonal wind, horizontal advection, and sea surface temperature obtained from the TOGA COARE, and is integrated for a selected 10-day period. The CVV has zonal and vertical components in the 2D x-z frame. Analysis of zonally-averaged and mass-integrated quantities shows that the correlation coefficient between the vertical component of the CVV and the sum of the cloud hydrometeor mixing ratios is 0.81, whereas the correlation coefficient between the zonal component and the sum of the mixing ratios is only 0.18. This indicates that the vertical component of the CVV is closely associated with tropical convection. The tendency equation for the vertical component of the CVV is derived and the zonally-averaged and mass-integrated tendency budgets are analyzed. The tendency of the vertical component of the CVV is determined by the interaction between the vorticity and the zonal gradient of cloud heating. The results demonstrate that the vertical component of the CVV is a cloud-linked parameter and can be used to study tropical convection.

Description: Previous estimates of land-atmosphere interaction (the impact of soil moisture on precipitation) have been limited by a severe paucity of relevant observational data and by the model-dependence of the various computational estimates. To counter this limitation, a dozen climate modeling groups have recently performed the same highly-controlled numerical experiment as part of a coordinated intercomparison project. This allows, for the first time ever, a superior multi-model approach to the estimation of the regions on the globe where precipitation is affected by soil moisture anomalies during Northern Hemisphere summer. Such estimation has many potential benefits; it can contribute, for example, to seasonal rainfall prediction efforts.

Description: Over the tropical land regions observations of the 85 GHz brightness temperature (T(sub 85v)) made by the TRMM Microwave Imager (TMI) radiometer when analyzed with the help of rain rate (R(sub pR)) deduced from the TRMM Precipitation Radar (PR) indicate that there are two maxima in rain rate. One strong maximum occurs when T(sub 85) has a value of about 220 K and the other weaker one when T(sub 85v) is much colder approx. 150 K. Together with the help of earlier studies based on airborne Doppler Radar observations and radiative transfer theoretical simulations, we infer the maximum near 220 K is a result of relatively weak scattering due to super cooled rain drops and water coated ice hydrometeors associated with a developing thunderstorm (Cb) that has a strong updraft. The other maximum is associated with strong scattering due to ice particles that are formed when the updraft collapses and the rain from the Cb is transit2oning from convective type to stratiform type. Incorporating these ideas and with a view to improve the estimation of rain rate from existing operational method applicable to the tropical land areas, we have developed a rain retrieval model. This model utilizes two parameters, that have a horizontal scale of approx. 20km, deduced from the TMI measurements at 19, 21 and 37 GHz (T(sub 19v), T(sub 21v), T(sub 37v). The third parameter in the model, namely the horizontal gradient of brightness temperature within the 20 km scale, is deduced from TMI measurements at 85 GHz. Utilizing these parameters our retrieval model is formulated to yield instantaneous rain rate on a scale of 20 km and seasonal average on a mesoscale that agree well with that of the PR.

Description: Linear regression trends for the years 1979-2003 were computed using the new Version 8 merged Solar Backscatter Ultraviolet (SBUV) data set of ozone profiles. These trends were compared to trends computed using ozone profiles from the Goddard Space Flight Center (GSFC) zonally averaged coupled model. Observed and modeled annual trends between 50 N and 50 S were a maximum in the higher latitudes of the upper stratosphere, with southern hemisphere (SH) trends greater than northern hemisphere (NH) trends. The observed upper stratospheric maximum annual trend is -5.5 +/- 0.9 % per decade (1 sigma) at 47.5 S and -3.8 +/- 0.5 % per decade at 47.5 N, to be compared with the modeled trends of -4.5 +/- 0.3 % per decade in the SH and -4.0 +/- 0.2% per decade in the NH. Both observed and modeled trends are most negative in winter and least negative in summer, although the modeled seasonal difference is less than observed. Model trends are shown to be greatest in winter due to a repartitioning of chlorine species and the increasing abundance of chlorine with time. The model results show that trend differences can occur depending on whether ozone profiles are in mixing ratio or number density coordinates, and on whether they are recorded on pressure or altitude levels.

Description: The data from twenty years of the NCEP numerical weather model have been used to calculate the IMF hydrostatic mapping function for several sites distributed in latitude from -66 degrees to +78 degrees. Comparison of heights estimated with the NMF hydrostatic mapping function demonstrates that using NMFh results in height errors at annual and semi-annual periods with amplitudes as large as approximately 8 mm and 4 mm, respectively, when data down to 5 degrees are included. The errors are smallest at the equator and increase towards the poles.

Description: Spatially detailed satellite data of mean color, sea ice concentration, surface temperature, clouds, and wind have been analyzed to quantify and study the large scale regional and temporal variability of phytoplankton blooms in the Arctic and peripheral seas from 1998 to 2002. In the Arctic basin, phytoplankton chlorophyll displays a large symmetry with the Eastern Arctic having about fivefold higher concentrations than those of the Western Arctic. Large monthly and yearly variability is also observed in the peripheral seas with the largest blooms occurring in the Bering Sea, Sea of Okhotsk, and the Barents Sea during spring. There is large interannual and seasonal variability in biomass with average chlorophyll concentrations in 2002 and 2001 being higher than earlier years in spring and summer. The seasonality in the latitudinal distribution of blooms is also very different such that the North Atlantic is usually most expansive in spring while the North Pacific is more extensive in autumn. Environmental factors that influence phytoplankton growth were examined, and results show relatively high negative correlation with sea ice retreat and strong positive correlation with temperature in early spring. Plankton growth, as indicated by biomass accumulation, in the Arctic and subarctic increases up to a threshold surface temperature of about 276-277 degree K (3-4 degree C) beyond which the concentrations start to decrease suggesting an optimal temperature or nutrient depletion. The correlation with clouds is significant in some areas but negligible in other areas, while the correlations with wind speed and its components are generally weak. The effects of clouds and winds are less predictable with weekly climatologies because of unknown effects of averaging variable and intermittent physical forcing (e.g. over storm event scales with mixing and upwelling of nutrients) and the time scales of acclimation by the phytoplankton.

Description: The NASA/GSFC Scanning Raman Lidar (SRL) was stationed on Andros Island, Bahamas for the third Convection and Moisture Experiment (CAMEX 3) held in August - September, 1998 and acquired an extensive set of water vapor and cirrus cloud measurements (Whiteman et al., 2001). The cirrus data studied here have been segmented by generating mechanism. Distinct differences in the optical properties of the clouds are found when the cirrus are hurricane-induced versus thunderstom-induced. Relationships of cirrus cloud optical depth, mean cloud temperature, and layer mean extinction-to-backscatter ratio (S) are presented and compared with mid-latitude and tropical results. Hurricane-induced cirrus clouds are found to generally possess lower values of S than thunderstorm induced clouds. Comparison of these measurements of S are made with other studies revealing at times large differences in the measurements. Given that S is a required parameter for spacebased retrievals of cloud optical depth using backscatter lidar, these large diffaences in S measurements present difficulties for space-based retrievals of cirrus cloud extinction and optical depth.

Description: Scattering properties of the Ocean surface have been widely used as a calibration reference for airborne and spaceborne microwave sensors. However, at millimeter-wave frequencies, the ocean surface backscattering mechanism is still not well understood, in part, due to the lack of experimental measurements. During the Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE), measurements of ocean surface backscattering were made using a 94-GHz (W-band) cloud radar onboard a NASA ER-2 high-altitude aircraft. The measurement set includes the normalized Ocean surface cross section over a range of the incidence angles under a variety of wind conditions. Analysis of the radar measurements shows good agreement with a quasi-specular scattering model. This unprecedented dataset enhances our knowledge about the Ocean surface scattering mechanism at 94 GHz. The results of this work support the proposition of using the Ocean surface as a calibration reference for airborne millimeter-wave cloud radars and for the ongoing NASA CloudSat mission, which will use a 94-GHz spaceborne cloud radar for global cloud measurements.

Description: Urbanization is one of the extreme cases of land use change. Most of world s population has moved to urban areas. Although currently only 1.2% of the land is considered urban, the spatial coverage and density of cities are expected to rapidly increase in the near future. It is estimated that by the year 2025, 60% of the world s population will live in cities. Human activity in urban environments also alters atmospheric composition; impacts components of the water cycle; and modifies the carbon cycle and ecosystems. However, our understanding of urbanization on the total Earth-climate system is incomplete. Better understanding of how the Earth s atmosphere-ocean-land-biosphere components interact as a coupled system and the influence of the urban environment on this climate system is critical. The goal of the 2003 AGU Union session Human-induced climate variations on urban areas: From observations to modeling was to bring together scientists from interdisciplinary backgrounds to discuss the data, scientific approaches and recent results on observing and modeling components of the urban environment with the intent of sampling our current stand and discussing future direction on this topic. Herein, a summary and discussion of the observations component of the session are presented.

Description: This paper describes ozone profiles from sonde data during the period of NASA s TRACE-A and the more recent SHADOZ (Southern Hemisphere Additional Ozonesondes) period. The data were taken by the South African Weather Service at the Irene (25 deg.54 min S; 28 deg. 13 min. E) station near Pretoria, South Africa, an area that is a unique mixture of local industry, heavy biofuels use and importation of biomass burning ozone from neighboring countries to the north. The main findings are: (1) With its geographical position at the edge of the subtropical transition zone, mid- latitude dynamical influences are evident at Irene, predominantly in winter when upper tropospheric ozone is enhanced as a result of stratospheric-tropospheric exchange. (2) There has been an increase in the near-surface ozone amount between the early 1990s and a decade later, presumably due to an influx of rural population toward the Johannesburg-Pretoria area, as well as with industrial growth and development. (3) Most significant for developing approaches for satellite ozone profile climatologies, cluster analysis has enabled the delineation of a background and "most polluted" profile. Enhancements of at least 30% occur throughout the troposphere in spring and in certain layers increases of 100 % are observed.

Description: Results provided by two different assimilation methodologies involving data from passive and active space-borne microwave instruments are presented. The impact of the precipitation estimates produced by the TRMM Microwave Imager (TMI) and Special Sensor Microwave/Imager (SSM/I) in a previously developed 1D variational continuous assimilation algorithm for assimilating tropical rainfall is shown on two hurricane cases. Results on the impact of the SeaWinds scatterometer on the intensity and track forecast of a mid-Atlantic hurricane are also presented. This work is the outcome of a collaborative effort between NASA and NOAA and indicates the substantial improvement in tropical cyclone forecasting that can result from the assimilation of space-based data in global atmospheric models.

Description: The Earth Observing System (EOS) Aura satellite is scheduled to launch in the second quarter of 2004. The Aura mission is designed to attack three science questions: (1) Is the ozone layer recovering as expected? (2) What are the sources and processes that control tropospheric pollutants? (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 at high vertical and horizontal resolution. Fig. 1 shows the Aura spacecraft and its four instruments.

Description: Howard (1833a) made the first documented observation of a temperature difference between an urban area and its rural environment. Manley (1958) termed this contrast the "urban heat island (UHI)". The UHI has now become a widely acknowledged, observed, and researched phenomenon because of its broad implications. It is estimated that by the year 2025, 60% of the world's population will live in cities (UNFP, 1999). In the United States, the current urban growth rate is approximately 12.5%, with 80% currently living in urban areas. As cities continue to grow, urban sprawl creates unique problems related to land use, transportation, agriculture, housing, pollution, and development for policymakers. Urban expansion and its associated urban heat islands also have measurable impacts on weather and climate processes.

Description: The July 2004 CRYSTAL-FACE field program, based from Key West, Florida, showed that long-lived thin tropopause cirrus (TTC) layers were common above thunderstorm anvils. This paper investigates the origins of these cloud using airborne measurements. The horizontal dimensions of the TTC were nearly identical to the convectively formed anvil beneath. However, the TTC did not appear to have originated from convective detrainment. Rather it appears to have formed in stably stratified air derived from high altitudes near the tropopause. The TTC was separated from the anvil by approx. 1 km, it lacked precipitation particles, and it was strongly depleted in HDO. Nonetheless, compared to surrounding clear air near the tropopause, the TTC was enriched in moisture and trace gases in a manner consistent with it having mixed with the same convective airmass that produced the anvil. Unlike surrounding air, the TTC had embedded a monochromatic gravity wave with a wavelength of 2 km and an amplitude of several hundred meters. Combined, this evidence, supported by a photograph from CRYSTAL-FACE, leads to the conjecture that the TTC originated as a pileus cloud layer, which formed near the tropopause ahead of vigorous convective uplift. We hypothesize that the pileus was penetrated by the convection, moistened through mixing, and once the convection subsided, it was sustained by radiative cooling due to the presence of the anvil layer beneath.

Description: The 6-year (1998-2003) rainfall products from the Tropical Rainfall Measuring Mission (TRMM) are used to quantify the Intertropical Convergence Zone (ITCZ) in the eastern Pacific (defined by longitudinal averages over 90 degrees W-130 degrees W) during boreal spring (March-April). The double ITCZ phenomenon, represented by the occurrence of two maxima with respect to latitude in monthly mean rainfall, is observed in most but not all of the years studied. The relative spatial locations of maxima in sea surface temperature (SST), rainfall, and surface pressure are examined. Interannual and weekly variability are characterized in SST, rainfall, surface convergence, total column water vapor, and cloud water. There appears to be a competition for rainfall between the two hemispheres during this season. When one of the two rainfall maxima is particularly strong, the other tends to be weak, with the total rainfall integrated over the two varying less than does the difference between the rainfall integrated over each separately. There is some evidence for a similar competition between the SST maxima in the two hemispheres, but this is more ambiguous, and there is evidence that some variations in the relative strengths of the two rainfall maxima may be independent of SST. Using a 25-year (1979-2003) monthly rainfall dataset from the Global Precipitation Climatology Project (GPCP), four distinct ITCZ types during March-April are defined, based on the relative strengths of rainfall peaks north and south of, and right over the equator. Composite meridional profiles and spatial distributions of rainfall and SST are documented for each type. Consistent with previous studies, an equatorial cold tongue is essential to the existence of the double ITCZs. However, too strong a cold tongue may dampen either the southern or northern rainfall maximum, depending on the magnitude of SST north of the equator.

Description: Remote sensing of snow depth has been used to infer snow depth for many years. Passive microwave remote sensing of snow depth is compared with the snow gauge data.

Description: The Tropical Rainfall Measuring Mission (TRMM) provides monthly rainfall estimates using data collected by the TRMM satellite. These estimates cover a substantial fraction of the earth's surface. The physical validation of TRMM estimates involves corroborating the accuracy of spaceborne estimates of areal rainfall by inferring errors and biases from ground-based rain estimates. The TRMM error budget consists of two major sources of error: retrieval and sampling. Sampling errors are intrinsic to the process of estimating monthly rainfall and occur because the satellite extrapolates monthly rainfall from a small subset of measurements collected only during satellite overpasses. Retrieval errors, on the other hand, are related to the process of collecting measurements while the satellite is overhead. One of the big challenges confronting the TRMM validation effort is how to best estimate these two main components of the TRMM error budget, which are not easily decoupled. This four-year study computed bulk sampling and retrieval errors for the TRMM microwave imager (TMI) and the precipitation radar (PR) by applying a technique that sub-samples gauge data at TRMM overpass times. Gridded monthly rain estimates are then computed from the monthly bulk statistics of the collected samples, providing a sensor-dependent gauge rain estimate that is assumed to include a TRMM equivalent sampling error. The sub-sampled gauge rain estimates are then used in conjunction with the monthly satellite and gauge (without sub- sampling) estimates to decouple retrieval and sampling errors. The computed mean sampling errors for the TMI and PR were 5.9% and 7.796, respectively, in good agreement with theoretical predictions. The PR year-to-year retrieval biases exceeded corresponding TMI biases, but it was found that these differences were partially due to negative TMI biases during cold months and positive TMI biases during warm months.

Description: Large errors in numerical weather prediction are often associated with explosive cyclogenesis. Most studes focus on the under-forecasting error, i.e. cases of rapidly developing cyclones which are poorly predicted in numerical models. However, the over-forecasting error (i.e., to predict an explosively developing cyclone which does not occur in reality) is a very common error that severely impacts the forecasting skill of all models and may also present economic costs if associated with operational forecasting. Unnecessary precautions taken by marine activities can result in severe economic loss. Moreover, frequent occurrence of over-forecasting can undermine the reliance on operational weather forecasting. Therefore, it is important to understand and reduce the prdctions of extreme weather associated with explosive cyclones which do not actually develop. In this study we choose a very prominent case of over-forecasting error in the northwestern Pacific. A 960 hPa cyclone develops in less than 24 hour in the 5-day forecast, with a deepening rate of about 30 hPa in one day. The cyclone is not versed in the analyses and is thus a case of severe over-forecasting. By assimilating AIRS data, the error is largely eliminated. By following the propagation of the anomaly that generates the spurious cyclone, it is found that a small mid-tropospheric geopotential height negative anomaly over the northern part of the Indian subcontinent in the initial conditions, propagates westward, is amplified by orography, and generates a very intense jet streak in the subtropical jet stream, with consequent explosive cyclogenesis over the Pacific. The AIRS assimilation eliminates this anomaly that may have been caused by erroneous upper-air data, and represents the jet stream more correctly. The energy associated with the jet is distributed over a much broader area and as a consequence a multiple, but much more moderate cyclogenesis is observed.

Description: The study employs a 108-year precipitation data record to identify statistically significant anomalies in rainfall downwind of the Phoenix urban region. The analysis reveals that during the monsoon season locations northeastern suburbs and exurbs of the Phoenix metropolitan area have experienced statistically significant increases in mean precipitation of 12 to 14 percent from a pre-urban (1895-1949) to post-urban (1950-2003) period. Mean and median post-urban precipitation totals in the anomaly region are significantly greater, in the statistical sense, than regions west of the city and in nearby mountainous regions of similar or greater topography. Further analysis of satellite-based rainfall totals for the summer of 2003 also reveal the existence of the anomaly region during a severe drought period. The anomaly can not simply be attributed to maximum topographic relief and is hypothesize to be related to urban-topographic interactions.

Description: The possible impact of an increase in global temperatures of about 2 C, as may be caused by a doubling of atmospheric CO2, is studied using historical satellite records of surface temperatures and sea ice from late 1970s to 2003. Updated satellite data indicate that the perennial ice continued to decline at an even faster rate of 9.2 % per decade than previously reported while concurrently, the surface temperatures have steadily been going up in most places except for some parts of northern Russia. Surface temperature is shown to be highly correlated with sea ice concentration in the seasonal sea ice regions. Results of regression analysis indicates that for every 1 C increase in temperature, the perennial ice area decreases by about 1.48 x 10(exp 6) square kilometers with the correlation coefficient being significant but only -0.57. Arctic warming is estimated to be about 0.46 C per decade on average in the Arctic but is shown to be off center with respect to the North Pole, and is prominent mainly in the Western Arctic and North America. The length of melt has been increasing by 13 days per decade over sea ice covered areas suggesting a thinning in the ice cover. The length of melt also increased by 5 days per decade over Greenland, 7 days per decade over the permafrost areas of North America but practically no change in Eurasia. Statistically derived projections indicate that the perennial sea ice cover would decline considerably in 2025, 2035, and 2060 when temperatures are predicted by models to reach the 2 C global increase.

Description: Ozone measurements from ozonesondes, AROTAL, DIAL, and POAM III instruments during the SOLVE-2/VINTERSOL period are composited in a time-varying, flow-following quasi-conservative (PV-6) coordinate space; the resulting composites from each instrument are mapped onto the other instruments locations and times. The mapped data are then used to intercompare data from the different instruments. Overall, the four data sets are found to be in good agreement. AROTAL shows somewhat lower values below 16 km, and DIAL has a positive bias at the upper limits of its altitude range. These intercomparisons are consistent with those obtained from more conventional near-coincident profiles, where available. Although the PV-theta mapping technique entails larger uncertainties of individual profile differences compared to direct near-coincident comparisons, the ability to include much larger numbers of comparisons can make this technique advantageous.

Description: Conventional wisdom is that lidar pulses do not significantly penetrate clouds having optical thickness exceeding about tau = 2, and that no returns are detectable from more than a shallow skin depth. Yet optically thicker clouds of tau much greater than 2 reflect a larger fraction of visible photons, and account for much of Earth s global average albedo. As cloud layer thickness grows, an increasing fraction of reflected photons are scattered multiple times within the cloud, and return from a diffuse concentric halo that grows around the incident pulse, increasing in horizontal area with layer physical thickness. The reflected halo is largely undetected by narrow field-of-view (FoV) receivers commonly used in lidar applications. THOR - Thickness from Off-beam Returns - is an airborne wide-angle detection system with multiple FoVs, capable of observing the diffuse halo, detecting wide-angle signal from which physical thickness of optically thick clouds can be retrieved. In this paper we describe the THOR system, demonstrate that the halo signal is stronger for thicker clouds, and validate physical thickness retrievals for clouds having z 〉 20, from NASA P-3B flights over the Department of Energy/Atmospheric Radiation Measurement/Southern Great Plains site, using the lidar, radar and other ancillary ground-based data.

Description: In this chapter, we describe the suite of Earth Observing System (EOS) Moderate-Resolution Imaging Spectroradiometer (MODIS) Terra and Aqua snow and sea ice products. Global, daily products, developed at Goddard Space Flight Center, are archived and distributed through the National Snow and Ice Data Center at various resolutions and on different grids useful for different communities Snow products include binary snow cover, snow albedo, and in the near future, fraction of snow in a 5OO-m pixel. Sea ice products include ice extent determined with two different algorithms, and sea ice surface temperature. The algorithms used to develop these products are described. Both the snow and sea ice products, available since February 24,2000, are useful for modelers. Validation of the products is also discussed.

Description: We report results from a study with the Numerical Spectral Model (NSM), which produces in the mesosphere significant inter-annual variations in the diurnal tide. Applying Hines Doppler Spread Parameterization (DPS), small-scale gravity waves (GW) drive the Quasi-biennial Oscillation (QBO) and Semi-annual Oscillation (SAO). With a GW source that peaks at the equator and is taken to be isotropic and independent of season, the NSM generates near the equator a QBO with variable periods around 27 months and zonal wind amplitudes close to 20 m / s at 30 Ism. As reported earlier, the NSM reproduces the observed equinoctial maxima in the diurnal tide at altitudes around 95 km. In the present paper it is shown that the QBO modulates the tide such that the seasonal amplitude maxima can vary from one year to another by as much as 30%. Since the period of the QBO is variable, its phase relative to the seasonal cycle changes. The magnitude of the QBO modulation of the tide thus varies considerably as our long-term model simulation shows. To shed light on the underlying mechanism, the relative importance of the linearized advection terms are discussed that involve the meridional and vertical winds of the diurnal tide.

Description: Sensitivity experiments with an atmospheric general circulation model (AGCM) show that parameterized rain re-evaporation has a large impact on simulated precipitation patterns in the tropical Pacific, especially on the configuration of the model s intertropical convergence zone (ITCZ). Weak re-evaporation leads t o the formation of a "double ITCZ" during the northern warm season. The double ITCZ is accompanied by strong coupling between precipitation and high-frequency vertical motion in the planetary boundary layer (PBL). Strong reevaporation leads to a better overall agreement of simulated precipitation with observations. The model s double ITCZ bias is reduced. At the same time, correlation between high-frequency vertical motion in the PBL and precipitation is reduced. Experiments with modified physics suggest that evaporative cooling by rain near the PBL top weakens the coupling between precipitation and vertical motion. This may reduce the model s tendency to form double ITCZs. The strength of high-frequency vertical motions in the PBL was also reduced directly through the introduction of a diffusive cumulus momentum transport (DCMT) parameterization. The DCMT had a visible impact on simulated precipitation in the tropics, but did not reduce the model s double bias in all cases.

Description: The 1930s was characterized by a decade of rainfall deficits and high temperatures that desiccated much of the United States Great Plains. Numerous dust storms created one of the most severe environmental catastrophes in U.S. history and led to the popular characterization of much of the southern Great Plains as the Dust Bowl . In this study, we show that the origin of the drought was in the anomalous tropical sea surface temperatures that occurred during that decade. We further show that interactions between the atmosphere and the land surface were essential to the development of the severe drought conditions. The results are based on simulations with the NASA Seasonal-to-Interannual Prediction Project general circulation model forced with observed and idealized sea surface temperatures. We contrast the 1930s drought with other major droughts of the 20th century, and speculate on the possibility of another Dust Bowl developing in the foreseeable future.

Description: The Sensors and Electronics Technology Branch of the NASA Glenn Research Center is developing thin-film-based sensors for surface measurement in propulsion system research. Thin-film sensors do not require special machining of the components on which they are mounted, and they are considerably thinner than wire- or foil-based sensors. One type of sensor being advanced is the thin-film thermocouple, specifically for applications in high-temperature combustion environments. Ceramics are being demonstrated as having the potential to meet the demands of thin-film thermocouples in advanced aerospace environments. The maximum-use temperature of noble metal thin-film thermocouples, 1500 C (2700 F), may not be adequate for components used in the increasingly harsh conditions of advanced aircraft and next-generation launch vehicles. Ceramic-based thermocouples are known for their high stability and robustness at temperatures exceeding 1500 C, but are typically in the form of bulky rods or probes. As part of ASTP, Glenn's Sensors and Electronics Technology Branch is leading an in-house effort to apply ceramics as thin-film thermocouples for extremely high-temperature applications as part of ASTP. Since the purity of the ceramics is crucial for the stability of the thermocouples, Glenn's Ceramics Branch and Case Western Reserve University are developing high-purity ceramic sputtering targets for fabricating high-temperature sensors. Glenn's Microsystems Fabrication Laboratory, supported by the Akima Corporation, is using these targets to fabricate thermocouple samples for testing. The first of the materials used were chromium silicide (CrSi) and tantalum carbide (TaC). These refractory materials are expected to survive temperatures in excess of 1500 C. Preliminary results indicate that the thermoelectric voltage output of a thin-film CrSi versus TaC thermocouple is 15 times that of the standard type R (platinum-rhodium versus platinum) thermocouple, producing 20 mV with a 200 C temperature gradient. The photograph on the left shows the CrSi-TaC thermocouple in a test fixture at Glenn, and the resulting output signal is shown on the right. The temperature differential across the sample, from the center of the sample inside the oven to the sample mount outside the oven, is measured using a type R thermocouple on the sample.

Description: The unsteady, incompressible Navier-Stokes equations are used for the direct numerical simulation (DNS) of spatially evolving disturbances in a three-dimensional (3-D) attachment-line boundary layer. Two-dimensional (2-D) disturbances are introduced either by forcing at the in ow or by harmonic-source generators at the wall; 3-D disturbances are introduced by harmonic-source generators at the wall. The DNS results are in good agreement with both 2-D non-parallel theory (for small-amplitude disturbances) and weakly nonlinear theory (for finite-amplitude disturbances), which validates the two theories. The 2-D DNS results indicate that nonlinear disturbance growth occurs near branch II of the neutral stability curve; however, steady suction can be used to stabilize this disturbance growth. For 3-D instabilities that are generated o the attachment line, spreading both toward and away from the attachment line causes energy transfer to the attachment-line and downstream instabilities; suction stabilizes these instabilities. Furthermore, 3-D instabilities are more stable than 2-D or quasi-2-D instabilities.

Description: A fast multigrid solver for the steady incompressible Euler equations is presented. Unlike time-marching schemes, this approach uses relaxation of the steady equations. Application of this method results in a discretization that correctly distinguishes between the advection and elliptic parts of the operator, allowing efficient smoothers to be constructed. Solvers for both unstructured triangular grids and structured quadrilateral grids have been written. Computations for channel flow and flow over a nonlifting airfoil have computed. Using Gauss-Seidel relaxation ordered in the flow direction, textbook multigrid convergence rates of nearly one order-of-magnitude residual reduction per multigrid cycle are achieved, independent of the grid spacing. This approach also may be applied to the compressible Euler equations and the incompressible Navier-Stokes equations.

Description: The structural analysis results for a graphite/epoxy quasi-isotropic circular plate subjected to a forced rotation at the boundary and pressure is presented. The analysis is to support a specialized material characterization test for composite cryogenic tanks. Finite element models were used to ensure panel integrity and determine the pressure necessary to achieve a predetermined equal biaxial strain value. The displacement results due to the forced rotation at the boundary led to a detailed study of the bending stiffness matrix [D]. The variation of the bending stiffness terms as a function of angular position is presented graphically, as well as, an illustrative technique of considering the laminate as an I-beam.

Description: Data from a single Weather Surveillance Radar-1988 Doppler (WSR-88D) and the National Lightning Detection Network are used to examine the characteristics of the convective storms that produced a severe tornado outbreak, including three tornadoes that reached F3 intensity, within Tropical Storm Beryl s remnants on 16 August 1994. Comparison of the radar data with reports of tornadoes suggests that only 13 cells produced the 29 tornadoes that were documented in Georgia and the Carolinas on that date. Six of these cells spawned multiple tornadoes, and the radar data confirm the presence of miniature supercells. One of the cells was identifiable on radar for 11 h. spawning tornadoes over a time period spanning approximately 6.5 h. Several other tornadic cells also exhibited great longevity, with cell lifetimes longer than ever previously documented in a landfalling tropical cyclone (TC) tornado event. This event is easily the most intense TC tornado outbreak yet documented with WSR-88Ds. Time-height analyses of the three strongest tornadic supercells are presented in order to document storm kinematic structure and to show how these storms appear at different ranges from a WSR-88D. In addition, cloud-to-ground (CG) lightning data are examined in Beryl s remnants. Although the tornadic cells were responsible for most of Beryl's CG lightning, their flash rates were only weak to moderate, and in all the tornadic storms the lightning flashes were almost entirely negative in polarity. A few of the single-tornado storms produced no detectable CG lightning at all. There is evidence that CG lightning rates decreased during the tornadoes, compared to 30-min periods before the tornadoes. A number of the storms spawned tornadoes just after producing their final CG lightning flashes. Contrary to the findings for flash rates, both peak currents and positive flash percentages were larger in Beryl's nontornadic storms than in the tornadic ones.

Description: Various artificial compressibility methods for calculating the three-dimensional incompressible Navier-Stokes equations are compared. Each method is described and numerical solutions to test problems are conducted. A comparison based on convergence behavior, accuracy, and robustness is given.

Description: This study provides systematic method for reducing power consumption in reduced gravity systems by adopting minimum velocity required to provide adequate CHF and preclude detrimental effects of reduced gravity . This study proves it is possible to use existing 1 ge flow boiling and CHF correlations and models to design reduced gravity systems provided minimum velocity criteria are met

Description: Recent workshops to define strategic research on critical issues in microgravity fluids and transport phenomena in support of mission orientated needs of NASA and many technical conferences over the years in support of fundamental research targeting NASA's long range missions goal have identified several phase change processes needed to design advanced space and planetary based systems for long duration operations Recommendation noted that phase change processes are profoundly affected by gravitational environment.

Description: We present experimental data on flow pattern transitions, pressure drop and flow characteristics for cocurrent gas-liquid flow through packed columns in microgravity. The flow pattern transition data indicates that the pulse flow regime exists over a wider range of gas and liquid flow rates under microgravity conditions compared to 1-g and the widely used Talmor map in 1-g is not applicable for predicting the transition boundaries. A new transition criterion between bubble and pulse flow in microgravity is proposed and tested using the data. Since there is no static head in microgravity, the pressure drop measured is the true frictional pressure drop. The pressure drop data, which has much smaller scatter than most reported 1-g data clearly shows that capillary effects can enhance the pressure drop (especially in the bubble flow regime) as much as 200% compared to that predicted by the single phase Ergun equation. The pressure drop data are correlated in terms of a two-phase friction factor and its dependence on the gas and liquid Reynolds numbers and the Suratman number. The influence of gravity on the pulse amplitude and frequency is also discussed and compared to that under normal gravity conditions. Experimental work is planned to determine the gas-liquid and liquid-solid mass transfer coefficients. Because of enhanced interfacial effects, we expect the gas-liquid transfer coefficients kLa and kGa (where a is the gas-liquid interfacial area) to be higher in microgravity than in normal gravity at the same flow conditions. This will be verified by gas absorption experiments, with and without reaction in the liquid phase, using oxygen, carbon dioxide, water and dilute aqueous amine solutions. The liquid-solid mass transfer coefficient will also be determined in the bubble as well as the pulse flow regimes using solid benzoic acid particles in the packing and measuring their rate of dissolution. The mass transfer coefficients in microgravity will be compared to those in normal gravity cocurrent flow to determine the mass transfer enhancement and propose new mass transfer correlations for two-phase gas-liquid flows through packed beds in microgravity.

Description: One particular characteristic observed in unsteady shear layers is the phase shift relative to the main flow. In attached boundary layers this will have an effect both on the instantaneous skin friction and heat transfer. In separation bubbles the contribution to the drag is dominated by the pressure distribution. However, the most significant effect appears to be the phase shift on the transition process. Unsteady transition behaviour may determine the bursting of the bubble resulting in an un-recoverable full separation. An early analysis of the phase shift was performed by Stokes for the incompressible boundary layer of an oscillating wall and an oscillating main flow. An amplitude overshoot within the shear layer as well as a phase shift were observed that can be attributed to the relatively slow diffusion of viscous stresses compared to the fast change of pressure. Experiments in a low speed facility with the boundary layer of a flat plate were evaluated in respect to phase shift. A pressure distribution similar to that on the suction surface of a turbomachinery aerofoil was superimposed generating a typical transitional separation bubble. A periodically unsteady main flow in the suction type wind tunnel was introduced via a rotating flap downstream of the test section. The experiments covered a range of the three similarity parameters of momentum-loss-thickness Reynolds-number of 92 to 226 and Strouhal-number (reduced frequency) of 0.0001 to 0.0004 at the separation point, and an amplitude range up to 19 %. The free stream turbulence level was less than 1% .Upstream of the separation point the phase shift in the laminar boundary layer does not appear to be affected significantly bay either of the three parameters. The trend perpendicular to the wall is similar to the Stokes analysis. The problem scales well with the wave velocity introduced by Stokes, however, the lag of the main flow near the wall is less than indicated analytically. The separation point comes closest to the Stokes analysis but the phase is still 20 degrees lower at the wall.

Description: Among the numerous causes for unsteadiness in turbo machinery flows are turbulence and flow environment, wakes from stationary and rotating vanes, boundary layer separation, boundary layer/shear layer instabilities, presence of shock waves and deliberate unsteadiness for flow control purposes. These unsteady phenomena may lead to flow-structure interactions such as flutter and forced vibration as well as system instabilities such as stall and surge. A major issue of unsteadiness relates to the fact that a fundamental understanding of unsteady flow physics is lacking and requires continued attention. Accurate simulations and sufficient high fidelity experimental data are not available. The Glenn Research Center plan for Engine Component Flow Physics Modeling is part of the NASA 21st Century Aircraft Program. The main components of the plan include Low Pressure Turbine National Combustor Code. The goals, technical output and benefits/impacts of each element are described in the presentation. The specific areas selected for discussion in this presentation are blade wake interactions, flow control, and combustor exit turbulence and modeling.

Description: The Effects of elevated free-stream turbulence (FST) on the natural and periodically excited separation bubbles were studied experimentally, due to the relevance of this flow to low-pressure turbine blades at low Reynolds numbers. A bubble was formed at the leading edge of a flat plate and the FST level was altered by placing a grid across the flow at different locations upstream of the plate. The mixing across the separated shear-layer, forming the free boundary of the bubble, increased due to the elevated FST and due to nominally two-dimensional periodic excitation, both flattening and shortening the bubble. Periodic excitation at frequencies that were at least an order of magnitude lower than those associated with the initial shear-layer instability, were very effective at low FST, because the amplitudes of the excitation frequency and its harmonic were amplified over the bubble. High frequency excitation (F+ 3, based on the length of the baseline low FST bubble) had a major effect close to the separation location, while farther downstream the excited fluctuations rapidly decayed in the reattachment region. Low frequency excitation, that generated waves comparable to the length of the unperturbed bubble (F+ 1) were less effective and their magnitude decayed at a slower rate downstream of reattachment. An increase in the level of the FST reduced the net effect of the periodic excitation on the mixing enhancement and subsequent reattachment process, probably due to a destructive interference between the nominally 2D excitation and the random (in space and time) FST, reducing the spanwise coherence and therefore the effectiveness of the current control strategy. However, even at the reduced effectiveness of 2D periodic excitation at elevated FST, it accelerated the reattachment process and the recovery rate of the reattached boundary layer, enhancing the boundary layer resistance to repeat separation and reducing its momentum loss further downstream.

Description: In this brief, we investigate the turbulent condensation of a population of droplets by means of a direct numerical simulation. To that end, a coupled Navier-Stokes/Lagrangian solver is used where each particle is tracked and its growth by water vapor condensation is monitored exactly. The main goals of the study are to find out whether turbulence broadens the droplet size distribution, as observed in in situ measurements. The second issue is to understand if and for how long a correlation between the droplet radius and the local supersaturation exists for the purpose of modeling sub-grid scale microphysics in cloud-resolving codes. This brief is organized as follows. In Section 2 the governing equations are presented, including the droplet condensation model. The implementation of the forcing procedure is described in Section 3. The simulation results are presented in Section 4 together with a sketch of a simple stochastic model for turbulent condensation. Conclusions and the main outcomes of the study are given in Section 5.

Description: Accurate, global profiling of wind velocity is highly desired by NASA, NOAA, the DOD/DOC/NASA Integrated Program Office (IPO)/NPOESS, DOD, and others for many applications such as validation and improvement of climate models, and improved weather prediction. The most promising technology to deliver this measurement from space is Doppler Wind Lidar (DWL). The NASA/NOAA Global Tropospheric Wind Sounder (GTWS) program is currently in the process of generating the science requirements for a space-based sensor. In order to optimize the process of defining science requirements, it is important for the scientific and user community to understand the nature of the wind measurements that DWL can make. These measurements are very different from those made by passive imaging sensors or by active radar sensors. The purpose of this paper is to convey the sampling characteristics and data product trade-offs of an orbiting DWL.

Description: Large-eddy simulation (LES) is a widely used technique in armospheric modeling research. In LES, large, unsteady, three dimensional structures are resolved and small structures that are not resolved on the computational grid are modeled. A filtering operation is applied to distinguish between resolved and unresolved scales. We present two near-surface models that have found use in atmospheric modeling. We also suggest a simpler eddy viscosity model that adopts Prandtl's mixing length model (Prandtl 1925) in the vicinity of the surface and blends with the dynamic Smagotinsky model (Germano et al, 1991) away from the surface. We evaluate the performance of these surface models by simulating a neutraly stratified atmospheric boundary layer.

Description: Three-dimensional computational techniques, in particular the uncoupled CFD-DSMC of the present study, are available to be applied to problems such as jet interactions with variable density regions ranging from a continuum jet to a rarefied free stream. When the value of the jet to free stream momentum flux ratio approximately greater than 2000 for a sharp leading edge flat plate forward separation vortices induced by the jet interaction are present near the surface. Also as the free stream number density n (infinity) decreases, the extent and magnitude of normalized pressure increases and moves upstream of the nozzle exit. Thus for the flat plate model the effect of decreasing n (infinity) is to change the sign of the moment caused by the jet interaction on the flat plate surface.

Description: NASA has developed an Earth Observing System (EOS) consisting of a series of satellites designed to study global change from space. The EOS flagship is the EOS TERRA satellite, launched in December 1999, equipped with five unique sensors to monitor and study the Earth s heat budget and many of the key controlling variables governing the Earth's climate system. CLAMS, the Chesapeake Lighthouse and Aircraft Measurements for Satellites field campaign was conducted from NASA Wallops Flight Facility and successfully executed over the middle Atlantic eastern seaboard from July 10 August 2, 2001. CLAMS is primarily a shortwave closure experiment designed to validate and improve EOS TERRA satellite data products being derived from three sensors: CERES (Clouds and Earth's Radiant Energy System), MISR (Multi-angle Imaging Spectro-Radiometer) and MODIS (MODerate Resolution Imaging Spectroradiometer). CLAMS is jointly sponsored by the CERES, MISR and MODIS instrument teams and the NASA GEWEX Global Aerosol Climatology Project (GACP). CLAMS primary objectives are to validate satellite-based retrievals of aerosol properties and vertical profiles of radiative flux, temperature and water vapor. Central to CLAMS measurement strategy is the Chesapeake Lighthouse, a stable sea platform located in the Atlantic Ocean, 13 miles east of Virginia Beach near the mouth of the Chesapeake Bay and the site of an ongoing CERES Ocean Validation Experiment (COVE). Six research aircraft were deployed to make detailed measurements of the atmosphere and ocean surface in the vicinity of COVE, over the surrounding ocean, over nearby NOAA buoys and over a few land sites. The measurements are used to validate and provide ground truth for simultaneous products being derived from TERRA data, a key step toward an improved understanding and ability to predict changes in the Earth's climate. One of the two CERES instruments on-board TERRA was programmed for Rotating Azimuth Plane Scans (RAPS) during CLAMS, increasing the CERES coverage over COVE by a factor of 10. Nine coordinated aircraft missions and numerous additional sorties were flown under a variety of atmospheric conditions and aerosol loadings. On one golden day, July 17, all six aircraft flew coordinated patterns, vertically stacked between 100 ft and 65,000 ft over the COVE site as the TERRA satellite orbited overhead. A summary of CLAMS measurement campaign and a description of the platforms and measurements is given.

Description: InSPACE is a microgravity fluid physics experiment that was operated on the International Space Station (ISS) in the Microgravity Science Glovebox from late March 2003 through early July 2003. (InSPACE is an acronym for Investigating the Structure of Paramagnetic Aggregates From Colloidal Emulsions.) The purpose of the experiment is to obtain fundamental data of the complex properties of an exciting class of smart materials termed magnetorheological (MR) fluids. MR fluids are suspensions, or colloids, comprised of small (micrometer-sized) superparamagnetic particles in a nonmagnetic medium. Colloids are suspensions of very small particles suspended in a liquid. (Examples of other colloids are blood, milk, and paint.) These controllable fluids can quickly transition into a nearly solid state when exposed to a magnetic field and return to their original liquid state when the magnetic field is removed. Controlling the strength of the magnetic field can control the relative stiffness of these fluids. MR fluids can be used to improve or develop new seat suspensions, robotics, clutches, airplane landing gear, and vibration damping systems. The principal investigator for InSPACE is Professor Alice P. Gast of the Massachusetts Institute of Technology (MIT). The InSPACE hardware was developed at the NASA Glenn Research Center. The InSPACE samples were delivered to the ISS in November 2002, on the Space Shuttle Endeavour, on Space Station Utilization Flight UF-2/STS113. Operations began on March 31, 2003, with the processing of three different particle size samples at multiple test parameters. This investigation focused on determining the structural organization of MR colloidal aggregates when exposed to a pulsing magnetic field. On Earth, the aggregates take the shape of footballs with spiky tips. This characteristic shape may be influenced by the pull of gravity, which causes most particles initially suspended in the fluid to sediment, (i.e., settle and collect at the bottom of the cell). In the absence of sedimentation effects on the ISS, the behavior and shape of these MR aggregate structures are dominated exclusively by magnetic and surface tension forces. The microscopic detail of these structures was imaged under two orthogonal camera views. The video was downlinked to the InSPACE team at Glenn's Telescience Support Center and to MIT and also recorded onboard the ISS on videotapes that will be brought back to the ground by the space shuttles. The study examined the effect on the structure formation by varying the magnetic field strength and pulse frequency, and particle size. Fundamental data that characterized the structure formation were obtained. InSPACE completed its last planned test run on July 2, 2003. Operations occurred on 21 days over approximately a 3-month period. Forty-one test points were completed during 26 test runs. During the initial testing, the procedures followed by the crew were modified to maximize the observation of some unexpected and interesting aggregate behavior. As a result Dr. Gast has reported on the formation of aggregate shapes that are more extended and diverse than those observed on the ground. Sheets of magnetic material folded over in a labyrinth pattern and large columnar aggregates with complex interfaces with the surrounding fluid are examples of the interesting structures that have been observed on the ISS. In light of these early findings, the understanding of the fundamental properties of MR fluids on the basis of ground-based observations may need to be reconsidered.The experiments on the ISS have provided a vast amount of video data for analysis. While this analysis is ongoing, plans are being made for additional experimental runs. For this purpose, additional hardware and cells containing samples of different magnetic particles and sizes are being fabricated for a future launch to the ISS. The InSPACE hardware will remain on orbit until this testing is completed.

Description: Leakage and wear are two fundamental problems in all traditional turbine seals that contribute to an engine's inefficiency. The solutions to seal leakage and wear conflict in the conventional design space. Reducing the clearance between the seal and rotating shaft reduces leakage but increases wear because of increased contact incidents. Increasing the clearance to reduce the contact between parts reduces wear but increases parasitic leakage. The goal of this effort is to develop a seal that restricts leakage flow using acoustic pressure while operating in a noncontacting manner, thereby increasing life. In 1996, Dr. Timothy Lucas announced his discovery of a method to produce shock-free high-amplitude pressure waves. For the first time, the formation of large acoustic pressures was possible using dissonant resonators. A pre-prototype acoustic seal developed at the NASA Glenn Research Center exploits this fundamental acoustic discovery: a specially shaped cavity oscillated at the contained fluid's resonant frequency produces high-amplitude acoustic pressure waves of a magnitude approaching those required of today's seals. While the original researchers are continuing their development of acoustic pumps, refrigeration compressors, and electronic thermal management systems using this technology, the goal of researchers at Glenn is to apply these acoustic principles to a revolutionary sealing device. When the acoustic resonator shape is optimized for the sealing device, the flow from a high-pressure cavity to a low-pressure cavity will be restricted by a series of high-amplitude standing pressure waves of higher pressure than the pressure to be sealed. Since the sealing resonator cavity will not touch the adjacent sealing structures, seal wear will be eliminated, improving system life. Under a cooperative agreement between Glenn and the Ohio Aerospace Institute (OAI), an acoustic-based pre-prototype seal was demonstrated for the first time. A pressurized cavity was attached to one end of the resonator while the other end remained open to ambient pressure. Measurements were taken at several values of applied pressure with the assembly stationary, oscillated at an off-resonance frequency, and then oscillated on-resonance. The three cases show that the flow through the conical resonator can be reduced by oscillating the resonator at the resonance frequency of the air contained within the cavity. The results are currently being compared with results obtained from a commercial computational fluid dynamics code. The objective is to improve the design through numerical simulation before fabricating a next-generation prototype sealing device. Future work is aimed at implementing acoustic seal design improvements to further reduce the leakage flow rate through the device and at reducing the device's overall size.

Description: A new technology for reducing turbulent skin friction, called the Microblowing Technique (MBT), is presented. Results from proof-of-concept experiments show that this technology could potentially reduce turbulent skin friction by more than 50% of the skin friction of a solid flat plate for subsonic and supersonic flow conditions. The primary purpose of this review paper is to provide readers with information on the turbulent skin friction reduction obtained from many experiments using the MBT. Although the MBT has a penalty for obtaining the microblowing air associated with it, some combinations of the MBT with suction boundary layer control methods are an attractive alternative for a real application. Several computational simulations to understand the flow physics of the MBT are also included. More experiments and computational fluid dynamics (CFD) computations are needed for the understanding of the unsteady flow nature of the MBT and the optimization of this new technology.

Description: Solve the time-dependent inviscid flow equations for this geometry subject to the specified inflow/outflow mean conditions and the fluctuating inflow velocity distortion. (1) Compute the unsteady solution until periodicity in pressure is achieved by showing that at least two successive periods are identical. Periodicity must be achieved on both the airfoil surface and the inflow/outflow boundaries. (2) Once periodicity is achieved, compute the pressure frequency spectra on the reference airfoil on both the upper and lower surfaces at x=(-0.25c,0.00, +0.25c), on the inflow boundary at (x,y)={1.5c,-0.3c), (-1.5c,0.0),(-1.5c,0.3c)} and on the outflow boundary at (x,y)= {(1.5c,-0.3c),(1.5c,0.0), (1.5c,0.3c)}. Express the spectral results in dB using the standard definition 20 log(P(sub(r.m.s)/P(sub ref), where p(sub ref) == 20 microPa. (3) Extract the harmonic pressure distributions on the inflow and outflow boundaries (i.e., on x= -/+ 1.5c lines) at the fundamental frequency omega and apply a Fourier transform in y direction to identify the spatial (i.e., mode order) structure of the pressure perturbations. Express the result in dB for each mode order. Repeat the process for the frequencies 2omega and 3omega.

Description: Three benchmark problems from the current and previous CAA workshops involving tone noise generated in viscous flows are investigated using the CE/SE finite volume method. The CE/SE method is first briefly reviewed. Then, the benchmark problems, namely, flow past a single cylinder (CAA Workshop II problem), flow past twin cylinders (from the current CAA Workshop IV, Category 5, Problem 1) and flow past a deep cavity with overhang (CAA Workshop III problem) are investigated. Generally good results are obtained in comparison with the experimental data.

Description: The effects of gravity on the bubble formation in an annular jet were studied. The experiments were conducted in the 2.2-second drop tower at the NASA Glenn Research Center. Terrestrial gravity experiments were conducted at the Fluid Dynamics Research Laboratory at the University of Oklahoma. Stainless steel tubing with inner diameters of 1/8" (gas inner annulus) and 5/16" (liquid outer annulus) served as the injector. A rectangular test section, 6" x 6" x 14" tall, made out of half-inch thick Lexan was used. Images of the annular jet were acquired using a high-speed camera. The effects of gravity and varying liquid and gas flow rates on bubble size, wavelength, and breakup length were documented. In general, the bubble diameter was found to be larger in terrestrial gravity than in microgravity for varying Weber numbers (0.05 - 0.16 and 5 - 11) and liquid flow rates (1.5 ft/s - 3.0 ft/s). The wavelength was found to be larger in terrestrial gravity than in microgravity, but remained constant for varying Weber numbers. For low Weber numbers (0.05 - 0.16), the breakup length in microgravity was significantly higher than in terrestrial gravity. Comparison with linear stability analysis showed estimated bubble sizes within 9% of experimental bubble sizes. Bubble size compared to other terrestrial gravity experiments with same flow conditions showed distinct differences in bubble size, which displayed the importance of injector geometry on bubble formation.

Description: The objective of the present study is to develop a two-fluid model formulation with interfacial area transport equation applicable for microgravity conditions. The new model is expected to make a leapfrog improvement by furnishing the constitutive relations for the interfacial interaction terms with the interfacial area transport equation, which can dynamically model the changes of the interfacial structures. In the first year of this three-year project supported by the U.S. NASA, Office of Biological and Physics Research, the primary focus is to design and construct a ground-based, microgravity two-phase flow simulation facility, in which two immiscible fluids with close density will be used. In predicting the two-phase flow behaviors in any two-phase flow system, the interfacial transfer terms are among the most essential factors in the modeling. These interfacial transfer terms in a two-fluid model specify the rate of phase change, momentum exchange, and energy transfer at the interface between the two phases. For the two-phase flow under the microgravity condition, the stability of the fluid particle interface and the interfacial structures are quite different from those under normal gravity condition. The flow structure may not reach an equilibrium condition and the two fluids may be loosely coupled such that the inertia terms of each fluid should be considered separately by use of the two-fluid model. Previous studies indicated that, unless phase-interaction terms are accurately modeled in the two-fluid model, the complex modeling does not necessarily warrant an accurate solution.

Description: Flow regime and pressure drop data was obtained and analyzed. Pulse flow exists at lower liquid flow rates in 0-g compared to 1-g. 1-g flow regime maps do not apply in microgravity. Pressure drop is higher in microgravity (enhanced interfacial effects).

Description: This report describes preliminary results of work done by JHU/APL under contract to the NASA Glenn Research Center to support flight testing of the Universal Access Transceiver (UAT) data link as a medium for weather data exchange. It presents a high level architectural description of the use of UAT to meet the program objectives with an identification of issues associated with the use of this data link, including a high level definition of the changes required to UAT avionics and ground-based receivers and supporting ground infrastructure to support implementation of the recommended architecture with focus on the issues associated with these changes.

Description: Automated cloud detection and tracking is an important step in assessing global climate change via remote sensing. Cloud masks, which indicate whether individual pixels depict clouds, are included in many of the data products that are based on data acquired on- board earth satellites. Many cloud-mask algorithms have the form of decision trees, which employ sequential tests that scientists designed based on empirical astrophysics studies and astrophysics simulations. Limitations of existing cloud masks restrict our ability to accurately track changes in cloud patterns over time. In this study we explored the potential benefits of automatically-learned decision trees for detecting clouds from images acquired using the Advanced Very High Resolution Radiometer (AVHRR) instrument on board the NOAA-14 weather satellite of the National Oceanic and Atmospheric Administration. We constructed three decision trees for a sample of 8km-daily AVHRR data from 2000 using a decision-tree learning procedure provided within MATLAB(R), and compared the accuracy of the decision trees to the accuracy of the cloud mask. We used ground observations collected by the National Aeronautics and Space Administration Clouds and the Earth s Radiant Energy Systems S COOL project as the gold standard. For the sample data, the accuracy of automatically learned decision trees was greater than the accuracy of the cloud masks included in the AVHRR data product.

Description: In fiscal year 2003, test cell 23 of the Research Combustion Laboratory (RCL 23) at the NASA Glenn Research Center was upgraded with the addition of gaseous hydrogen as a working propellant and the addition of a 450-psig air-supply system. Test flexibility was further enhanced by upgrades to the facility control systems. RCL 23 can now test with gaseous hydrogen flow rates up to 0.05 lbm/sec and jet fuel flow rates up to 0.62 lbm/sec. Research airflow rates up to 3 lbm/sec are possible with the 450-psig supply system over a range of inlet temperatures. Nonvitiated, heated air is supplied from a shell and tube heat exchanger. The maximum nonvitiated facility air temperature is 1100 F at 1.5 lbm/sec. Research-section exhaust temperatures are limited to 3200 F because of material and cooling capacity limits. A variety of support systems are available depending on the research hardware configuration. Test section ignition can be provided via either a hydrogen air torch system or an electronic spark system. Emissions measurements are obtained with either pneumatically or electromechanically actuated gas sample probes, and the electromechanical system allows for radial measurements at a user-specified axial location for measurement of emissions profiles. Gas analysis data can be obtained for a variety of species, including carbon monoxide (CO), carbon dioxide (CO2), nitrogen oxides (NO and NOx), oxygen (O2), unburnt hydrocarbons, and unburnt hydrogen. Facility control is accomplished with a programmable logic control system. Facility operations have been upgraded to a system based on graphical user interface control screens. A data system is available for real-time acquisition and monitoring of both measurements in engineering units and performance calculations. The upgrades have made RCL 23 a highly flexible facility for research into low emissions gas turbine combustor concepts, and the flame tube configuration inherently allows for a variety of fuel nozzle configurations to be tested in a cost-effective manner. RCL 23 is poised to be a leading facility for developing modern low-emission fuel nozzles for use with jet fuel and alternative fuels.

Description: The operation of a packed bed reactor (PBR) involves gas and liquid flowing simultaneously through a fixed-bed of solid particles. Depending on the application, the particles can be various shapes and sizes but are generally designed to force the two fluid phases through a tortuous route of narrow channels connecting the interstitial space. The PBR is the most common type of reactor in industry because it provides for intimate contact and high rates of transport between the phases needed to sustain chemical or biological reactions. The packing may also serve as either a catalyst or as a support for growing biological material. Furthermore, this type of reactor is relatively compact and requires minimal power to operate. This makes it an excellent candidate for unit operations in support of long-duration human space activities.

Description: Recent research efforts within the Microgravity Fluid Physics Branch of the NASA Glenn Research Center have necessitated the development of a microscope capable of high-resolution, three-dimensional imaging of intracellular structure and tissue morphology. Standard optical microscopy works well for thin samples, but it does not allow the imaging of thick samples because of severe degradation caused by out-of-focus object structure. Confocal microscopy, which is a laser-based scanning microscopy, provides improved three-dimensional imaging and true optical sectioning by excluding the out-of-focus light. However, in confocal microscopy, out-of-focus object structure is still illuminated by the incoming beam, which can lead to substantial photo-bleaching. In addition, confocal microscopy is plagued by limited penetration depth, signal loss due to the presence of a confocal pinhole, and the possibility of live-cell damage. Two-photon microscopy is a novel form of laser-based scanning microscopy that allows three-dimensional imaging without many of the problems inherent in confocal microscopy. Unlike one-photon microscopy, it utilizes the nonlinear absorption of two near-infrared photons. However, the efficiency of two-photon absorption is much lower than that of one-photon absorption because of the nonlinear (i.e., quadratic) electric field dependence, so an ultrafast pulsed laser source must typically be employed. On the other hand, this stringent energy density requirement effectively localizes fluorophore excitation to the focal volume. Consequently, two-photon microscopy provides optical sectioning and confocal performance without the need for a signal-limiting pinhole. In addition, there is a reduction in photo-damage because of the longer excitation wavelength, a reduction in background fluorescence, and a 4 increase in penetration depth over confocal methods because of the reduction in Rayleigh scattering.

Description: Many valuable advances in combustion science have come from observations of microgravity flames. This research is contributing to the improved efficiency and reduced emissions of practical combustors and is benefiting terrestrial and spacecraft fire safety. Unfortunately, difficulties associated with microgravity have prevented many types of measurements in microgravity flames. In particular, temperature measurements in flames are extremely important but have been limited in microgravity. A novel method of measuring temperatures in microgravity flames is being developed in-house at the National Center for Microgravity Research and the NASA Glenn Research Center and is described here. Called thin-filament pyrometry, it involves using a camera to determine the local gas temperature from the intensity of inserted fibers glowing in a flame. It is demonstrated here to provide accurate measurements of gas temperatures in a flame simultaneously at many locations. The experiment is shown. The flame is a laminar gas jet diffusion flame fueled by methane (CH4) flowing from a 14-mm round burner at a pressure of 1 atm. A coflowing stream of air is used to prevent flame flicker. Nine glowing fibers are visible. These fibers are made of silicon carbide (SiC) and have a diameter of 15 m (for comparison, the average human hair is 75 m in diameter). Because the fibers are so thin, they do little to disturb the flame and their temperature remains close to that of the local gas. The flame and glowing filaments were imaged with a digital black-and-white video camera. This camera has an imaging area of 1000 by 1000 pixels and a wide dynamic range of 12 bits. The resolution of the camera and optics was 0.1 mm. Optical filters were placed in front of the camera to limit incoming light to 750, 850, 950, and 1050 nm. Temperatures were measured in the same flame in the absence of fibers using 50-m Btype thermocouples. These thermocouples provide very accurate temperatures, but they generally are not useful in microgravity tests because they measure temperature at only one location at a time. Thermocouple measurements at a height of 11 mm above the burner were used to calibrate the thin-filament pyrometry system at all four wavelengths. This calibration was used to perform thin-filament pyrometry at other heights above the burner. One such profile is shown in this graph; this is for a height of 21 mm. The agreement between the pyrometry measurements and thermocouple results at this height is excellent in the range of 1000 to 2000 K, with an estimated uncertainty of 50 K and an estimated upper limit of 2500 K. Neither the thermocouple nor the thin-filament pyrometry temperatures have been corrected for radiation, but the correction is expected to be nearly the same for both methods. We anticipate that thin-filament pyrometry similar to that developed here will become an important diagnostic for studies of microgravity flames owing to its accuracy and its ability to simultaneously measure finely spaced temperatures.

Description: When discharged into an International Space Station (ISS) payload rack, a carbon dioxide (CO2) portable fire extinguisher (PFE) must extinguish a fire by decreasing the oxygen in the rack by 50 percent within 60 sec. The length of time needed for this oxygen reduction throughout the rack and the length of time that the CO2 concentration remains high enough to prevent the fire from reigniting is important when determining the effectiveness of the response and postfire procedures. Furthermore, in the absence of gravity, the local flow velocity can make the difference between a fire that spreads rapidly and one that self-extinguishes after ignition. A numerical simulation of the discharge of CO2 from PFE into the Combustion Integrated Rack (CIR) in microgravity was performed to obtain the local velocity and CO2 concentration. The complicated flow field around the PFE nozzle exits was modeled by sources of equivalent mass and momentum flux at a location downstream of the nozzle. The time for the concentration of CO2 to reach a level that would extinguish a fire anywhere in the rack was determined using the Fire Dynamics Simulator (FDS), a computational fluid dynamics code developed by the National Institute of Standards and Technology specifically to evaluate the development of a fire and smoke transport. The simulation shows that CO2, as well as any smoke and combustion gases produced by a fire, would be discharged into the ISS cabin through the resource utility panel at the bottom of the rack. These simulations will be validated by comparing the results with velocity and CO2 concentration measurements obtained during the fire suppression system verification tests conducted on the CIR in March 2003. Once these numerical simulations are validated, portions of the ISS labs and living areas will be modeled to determine the local flow conditions before, during, and after a fire event. These simulations can yield specific information about how long it takes for smoke and combustion gases produced by a fire to reach a detector location, how large the fire would be when the detector alarms, and the behavior of the fire until it has been extinguished. This new capability could then be used to optimize the location of fire detectors and fire-suppression ports as well as to evaluate the effectiveness of fire suppressants and response strategies. Numerical data collected from these simulations could also be used to develop a virtual reality fire event for crew training and fire safety awareness. This work is funded by NASA's Bioastronautics Initiative, which has the objective of ensuring and enhancing the health, safety, and performance of humans in space. As part of this initiative, the Microgravity Combustion Science Branch at the NASA Glenn Research Center is conducting spacecraft fire safety research to significantly improve fire safety on inhabited spacecraft.

Description: A fast-running unsteady aerodynamics code, LINFLUX, was previously developed for predicting turbomachinery flutter. This linearized code, based on a frequency domain method, models the effects of steady blade loading through a nonlinear steady flow field. The LINFLUX code, which is 6 to 7 times faster than the corresponding nonlinear time domain code, is suitable for use in the initial design phase. Earlier, this code was verified through application to a research fan, and it was shown that the predictions of work per cycle and flutter compared well with those from a nonlinear time-marching aeroelastic code, TURBO-AE. Now, the LINFLUX code has been applied to real configurations: fans developed under the Energy Efficient Engine (E-cubed) Program and the Quiet Aircraft Technology (QAT) project. The LINFLUX code starts with a steady nonlinear aerodynamic flow field and solves the unsteady linearized Euler equations to calculate the unsteady aerodynamic forces on the turbomachinery blades. First, a steady aerodynamic solution is computed for given operating conditions using the nonlinear unsteady aerodynamic code TURBO-AE. A blade vibration analysis is done to determine the frequencies and mode shapes of the vibrating blades, and an interface code is used to convert the steady aerodynamic solution to a form required by LINFLUX. A preprocessor is used to interpolate the mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh. Then, LINFLUX is used to calculate the unsteady aerodynamic pressure distribution for a given vibration mode, frequency, and interblade phase angle. Finally, a post-processor uses the unsteady pressures to calculate the generalized aerodynamic forces, eigenvalues, an esponse amplitudes. The eigenvalues determine the flutter frequency and damping. Results of flutter calculations from the LINFLUX code are presented for (1) the E-cubed fan developed under the E-cubed program and (2) the Quiet High Speed Fan (QHSF) developed under the Quiet Aircraft Technology project. The results are compared with those obtained from the TURBO-AE code. A graph of the work done per vibration cycle for the first vibration mode of the E-cubed fan is shown. It can be seen that the LINFLUX results show a very good comparison with TURBO-AE results over the entire range of interblade phase angle. The work done per vibration cycle for the first vibration mode of the QHSF fan is shown. Once again, the LINFLUX results compare very well with the results from the TURBOAE code.

Description: A new in-house test capability has been developed at the NASA Glenn Research Center, where a critical component of the Stirling Radioisotope Generator (SRG) is undergoing extensive testing to aid the development of analytical life prediction methodology and to experimentally aid in verification of the flight-design component's life. The new facility includes two test rigs that are performing creep testing of the SRG heater head pressure vessel test articles at design temperature and with wall stresses ranging from operating level to seven times that (see the following photograph).

Description: The TopMaker technique was developed in an effort to reduce the time required for grid generation in complex numerical studies. Topology generation accounts for much of the man-hours required for structured multiblock grids. With regard to structured multiblock grids, topology refers to how the blocks are arranged and connected. A two-dimensional multiblock topology generation technique has been developed at the NASA Glenn Research Center. Very general configurations can be addressed by the technique. A configuration is defined by a collection of non-intersecting closed curves, which will be referred to as loops. More than a single loop implies that holes exist in the domain, which poses no problem. This technique requires only the medial vertices and the touch points that define each vertex. From the information about the medial vertices, the connectivity between medial vertices is generated. The physical shape of the medial edge is not required. By applying a few simple rules to each medial edge, a multiblock topology can be generated without user intervention. The resulting topologies contain only the level of complexity dictated by the configurations. Grid lines remain attached to the boundary except at sharp concave turns, where a change in index family is introduced as would be desired. Keeping grid lines attached to the boundary is especially important in computational fluid dynamics, where highly clustered grids are used near no-slip boundaries. This technique is simple and robust and can easily be incorporated into the overall grid-generation process.

Description: Over the past thirty years, numerical methods and simulation tools for incompressible flows have been advanced as a subset of the computational fluid dynamics (CFD) discipline. Although incompressible flows are encountered in many areas of engineering, simulation of compressible flow has been the major driver for developing computational algorithms and tools. This is probably due to the rather stringent requirements for predicting aerodynamic performance characteristics of flight vehicles, while flow devices involving low-speed or incompressible flow could be reasonably well designed without resorting to accurate numerical simulations. As flow devices are required to be more sophisticated and highly efficient CFD took become increasingly important in fluid engineering for incompressible and low-speed flow. This paper reviews some of the successes made possible by advances in computational technologies during the same period, and discusses some of the current challenges faced in computing incompressible flows.

Description: Transport of solid-liquid slurries in pipeline transport over short and medium distances is very important in many industries, including mining related processes. The particle image velocimetry technique was successfully utilized to investigate the velocities and kinetic energy fluctuations of slurry particles at the tongue region of an optically-clear centrifugal pump. The experiments were conducted using 500 micron glass beads at volumetric Concentrations of 2.5% and 5% and at pump speeds of 725 rpm and 1000 rpm. The fluctuation kinetic energy increased approximately 200% to 500% as the pump speed was increased from 725 rpm to IO00 rpm. The directional impingement mechanism is more significant at the pressure side of the blade, tongue and the casing. This mechanism becomes more important as the speed increases. This suggests that the impeller; tongue and the casing of the slurry pump can wear out quickly, especially with an increase in speed. In this paper the emphasis is on the tongue region. The random impingement mechanism caused by the fluctuation kinetic energy of the solids can play an important role on the erosion of the tongue area.

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

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

Description: Based on an empirical analysis of measurements with the High Resolution Doppler Imager (HRDI) on the UARS spacecraft in the upper mesosphere (95 km), persistent and regular intra-seasonal oscillations (ISO) with periods of about 2 to 4 months have recently been reported in the zonal-mean meridional winds. Similar oscillations have also been discussed independently in a modeling study, and they were attributed to wave-mean-flow interactions. The observed and modeled meridional wind ISOs were largely confined to low latitudes. We report here an analysis of concurrent temperature measurements on UARS, which produces oscillations similar to those seen in the meridional winds. Although the temperature oscillations are observed at lower altitudes (55 km), their phase variations with latitude are qualitatively consistent with the inferred properties seen in the meridional winds and thus provide independent evidence for the existence of ISOs in the mesosphere.

Description: Formation and evolution of cirrostratus in response to weak, uniform and constant synoptic forcing is simulated using a one-dimensional numerical model with explicit microphysics, in which the particle size distribution in each grid box is fully resolved. A series of tests of the model response to nucleation modes (homogeneous-freezing-only/heterogeneous nucleation) and heterogeneous nucleation parameters are performed. In the case studied here, nucleation is first activated in the prescribed moist layer. A continuous cloud-top nucleation zone with a depth depending on the vertical humidity gradient and one of the nucleation parameters is developed afterward. For the heterogeneous nucleation cases, intermittent nucleation zones in the mid-upper portion of the cloud form where the relative humidity is on the rise, because existent ice crystals do not uptake excess water vapor efficiently, and ice nuclei (IN) are available. Vertical resolution as fine as 1 m is required for realistic simulation of the homogeneous-freezing-only scenario, while the model resolution requirement is more relaxed in the cases where heterogeneous nucleation dominates. Bulk microphysical and optical properties are evaluated and compared. Ice particle number flux divergence, which is due to the vertical gradient of the gravity-induced particle sedimentation, is constantly and rapidly changing the local ice number concentration, even in the nucleation zone. When the depth of the nucleation zone is shallow, particle number concentration decreases rapidly as ice particles grow and sediment away from the nucleation zone. When the depth of the nucleation zone is large, a region of high ice number concentration can be sustained. The depth of nucleation zone is an important parameter to be considered in parametric treatments of ice cloud generation.

Description: The Goddard Lidar Observatory for Winds (GLOW) is a mobile direct detection Doppler lidar system which uses the double edge technique to measure the Doppler shift of the molecular backscattered laser signal at a wavelength of 355 nm. In the spring of 2002 GLOW was deployed to the western Oklahoma profiling site (36 deg 33.500 min. N, 100 deg. 36.371 min. W) to participate in the International H2O Project (IHOP). During the IHOP campaign over 240 hours of wind profiles were obtained with the GLOW lidar in support of a variety of scientific investigations.

Description: This study presents a systematic and quantitative analysis of the effect of inhomogeneous surface albedo on shortwave cloud absorption estimates. We use 3D radiative transfer modeling with gradually complex clouds over a simplified surface to calculate cloud absorption. We find that averaging surface albedo always underestimates cloud absorption, and thus accounting for surface heterogeneity always enhances cloud absorption. However, the impact on cloud absorption estimates is not enough to explain the discrepancy between measured and model calculated shortwave cloud absorptions.

Description: High-resolution mesoscale model simulations of the 6-7 May 2000 Missouri flash flood event were performed to test the impact of model initialization and land surface treatment on timing, intensity, and location of extreme precipitation. In this flash flood event, a mesoscale convective system (MCS) produced over 340 mm of rain in roughly 9 hours in some locations. Two different types of model initialization were employed: 1) NCEP global reanalysis with 2.5-degree grid spacing and 12-hour temporal resolution, and 2) Eta reanalysis with 40- km grid spacing and $hour temporal resolution. In addition, two different land surface treatments were considered. A simple land scheme. (SLAB) keeps soil moisture fixed at initial values throughout the simulation, while a more sophisticated land model (PLACE) allows for r interactive feedback. Simulations with high-resolution Eta model initialization show considerable improvement in the intensity of precipitation due to the presence in the initialization of a residual mesoscale convective vortex (hlCV) from a previous MCS. Simulations with the PLACE land model show improved location of heavy precipitation. Since soil moisture can vary over time in the PLACE model, surface energy fluxes exhibit strong spatial gradients. These surface energy flux gradients help produce a strong low-level jet (LLJ) in the correct location. The LLJ then interacts with the cold outflow boundary of the MCS to produce new convective cells. The simulation with both high-resolution model initialization and time-varying soil moisture test reproduces the intensity and location of observed rainfall.

Description: The influence of hurricane-ocean coupling on intensity and track of tropical cyclones (TCs) is investigated through idealized numerical experiments using a coupled hurricane-ocean model. The focus is placed on how air-sea interaction affects TC tracks and intensity. It is found that the symmetric sea surface temperature (SST) cooling is primarily responsible for the TC weakening in the coupled experiments because the induced asymmetric circulation associated with the asymmetric SST anomalies is weak and shallow. The track difference between the coupled and fixed SST experiments is generally small because of the competing processes. One is associated with the modified TC asymmetries. The asymmetric SST anomalies - weaken the surface fluxes in the rear and enhance the fluxes in the front. As a result, the enhanced diabatic heating is located on the southern side for a westward-moving TC, tending to shift the TC southward. The symmetric SST anomalies weakens the TC intensity and thus the dymmetrization process, leading to more prominent TC asymmetries. The other is associated with the weakening of the beta drift resulting from the weakening of the TC outer strength. In the coupled experiment, the weakening of the beta drift leads to a more northward shift. By adjusting the vortex outer strength of the initial vortices, the beta drift can vary while the effect of air-sea interaction changes little. Two types of track differences simulated in the previous numerical studies are obtained.

Description: The association between tall precipitation and tropical cyclone intensification may have implications for the difficult task of forecasting the destructive potential of these storms. We propose a novel way to use radar-observed rain height to help predict tropical cyclone intensity. Then, we adapt this technique for use on the much more plentiful data from infrared and microwave instruments.

Description: Abstract: Since 1998 the Southern Hemisphere ADditional OZonesondes (SHADOZ) project has collected more than 2000 ozone profiles from a dozen tropical and subtropical sites using balloon-borne electrochemical concentration cell (ECC) ozonesondes. The data (with accompanying pressure-temperature-humidity soundings) are archived. Analysis of ozonesonde imprecision within the SHADOZ dataset revealed that variations in ozonesonde technique could lead to station-to-station biases in the measurements. In this paper imprecisions and accuracy in the SHADOZ dataset are examined in light of new data. When SHADOZ total ozone column amounts are compared to version 8 TOMS (2004 release), discrepancies between sonde and satellite datasets decline 1-2 percentage points on average, compared to version 7 TOMS. Variability among stations is evaluated using total ozone normalized to TOMS and results of laboratory tests on ozonesondes (JOSE-2O00, Julich Ozonesonde Intercomparison Experiment). Ozone deviations from a standard instrument in the JOSE flight simulation chamber resemble those of SHADOZ station data relative to a SHADOZ-defined climatological reference. Certain systematic variations in SHADOZ ozone profiles are accounted for by differences in solution composition, data processing and instrument (manufacturer). Instrument bias leads to a greater ozone measurement above 25 km over Nairobi and to lower total column ozone at three Pacific sites compared to other SHADOZ stations at 0-20 deg.S.

Description: Observational records in the past 50 years show an upward trend of boreal-summer precipitation over central eastern China and a downward trend over northern China. During boreal spring, the trend is upward over southeastern China and downward over central eastern China. This study explores the forcing mechanism of these trends in association with the global sea-surface temperature (SST) variations on the interannual and inter-decadal timescales. Results based on Singular Value Decomposition analyses (SVD) show that the interannual variability of China precipitation in boreal spring and summer can be well defined by two centers of actions for each season, which are co-varying with two interannual modes of SSTs. The first SVD modes of precipitation in spring and summer, which are centered in southeastern China and northern China, respectively, are linked to an ENSO-like mode of SSTs. The second SVD modes of precipitation in both seasons are confined to central eastern China, and are primarily linked to SST variations over the warm pool and Indian Ocean. Features of the anomalous 850-hPa winds and 700-Wa geopotential height corresponding to these modes support a physical mechanism that explains the causal links between the modal variations of precipitation and SSTs. On the decadal and longer timescale, similar causal links are found between the same modes of precipitation and SSTs, except for the case of springtime precipitation over central eastern China. For this case, while the interannual mode of precipitation is positively correlated with the interannual variations of SSTs over the warm pool and Indian Ocean; the inter-decadal mode is negatively correlated with a different SST mode, the North Pacific mode. The later is responsible for the observed downward trend of springtime precipitation over central eastern China. For all other cases, both the interannual and inter-decadal variations of precipitation can be explained by the same mode of SSTs. The upward trend of springtime precipitation over southeastern China and downward trend of summertime precipitation over northern China are attributable to the warming trend of the ENSO-like mode. The recent frequent summertime floods over central eastern China are linked to the warming trend of SSTs over the warm pool and Indian Ocean.

Description: A short course on off-beam cloud lidar is given. Specific topics addressed include: motivation and goal of off-beam cloud lidar; diffusion physics; numeric amalysis; and validity of the diffusion approximation. A demo of the process is included.

Description: This Symposium is intended to bring together the often distinct cultures of the Stability and Control (S&C) community and the Computational Fluid Dynamics (CFD) community. The COMSAC program is itself a new effort by NASA Langley to accelerate the application of high end CFD methodologies to the demanding job of predicting stability and control characteristics of aircraft. This talk is intended to set the stage for needing a program like COMSAC. It is not intended to give details of the program itself. The topics include: 1) S&C Challenges; 2) Aero prediction methodology; 3) CFD applications; 4) NASA COMSAC planning; 5) Objectives of symposium; and 6) Closing remarks.

Description: The motivation behind the inclusion of unsteady aerodynamics and aeroelastic effects in the computation of stability and control (S&C) derivatives will be discussed as they pertain to aeroelastic and aeroservoelastic analysis. This topic will be addressed in the context of two applications, the first being the estimation of S&C derivatives for a cable-mounted aeroservoelastic wind tunnel model tested in the NASA Langley Research Center (LaRC) Transonic Dynamics Tunnel (TDT). The second application will be the prediction of the nonlinear aeroservoelastic phenomenon known as Residual Pitch Oscillation (RPO) on the B-2 Bomber. Techniques and strategies used in these applications to compute S&C derivatives and perform flight simulations will be reviewed, and computational results will be presented.