ALBERT

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

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

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

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

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

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

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

Description: This paper presents a detailed study of a single thunderstorm anvil cirrus cloud measured on 21 July 2002 near southern Florida during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers--Florida Area Cirrus Experiment (CRYSTAL-FACE). NASA WB-57F and University of North Dakota Citation aircraft tracked the microphysical and radiative development of the anvil for 3 h. Measurements showed that the cloud mass that was advected downwind from the thunderstorm was separated vertically into two layers: a cirrus anvil with cloud-top temperatures of -45?C lay below a second, thin tropopause cirrus (TTC) layer with the same horizontal dimensions as the anvil and temperatures near -70?C. In both cloud layers, ice crystals smaller than 50 ?m across dominated the size distributions and cloud radiative properties. In the anvil, ice crystals larger than 50 ?m aggregated and precipitated while small ice crystals increasingly dominated the size distributions; as a consequence, measured ice water contents and ice crystal effective radii decreased with time. Meanwhile, the anvil thinned vertically and maintained a stratification similar to its environment. Because effective radii were small, radiative heating and cooling were concentrated in layers approximately 100 m thick at the anvil top and base. A simple analysis suggests that the anvil cirrus spread laterally because mixing in these radiatively driven layers created horizontal pressure gradients between the cloud and its stratified environment. The TTC layer also spread but, unlike the anvil, did not dissipate--perhaps because the anvil shielded the TTC from terrestrial infrared heating. Calculations of top-of-troposphere radiative forcing above the anvil and TTC showed strong cooling that tapered as the anvil evolved.

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