Validation of the Aura Microwave Limb Sounder Temperature and Geopotential Height Measurements

Global satellite observations of temperature and geopotential height (GPH) from the Microwave Limb Sounder (MLS) on the EOS Aura spacecraft are discussed. The precision, resolution, and accuracy of the data produced by the MLS version 2.2 processing algorithms are quantified, and recommendations for data screening are made. Temperature precision is 1 K or better from 316 hPa to 3.16 hPa, degrading to ∼3 K at 0.001 hPa. The vertical resolution is 3 km at 31.6 hPa, degrading to 6 km at 316 hPa and to ∼13 km at 0.001 hPa. Comparisons with analyses (Goddard Earth Observing System version 5.0.1 (GEOS-5), European Centre for Medium-range Weather Forecasts (ECMWF), Met Office (MetO)) and other observations (CHAllenging Minisatellite Payload (CHAMP), Atmospheric Infrared Sounder/Advanced Microwave Sounder Unit (AIRS/AMSU), Sounding of the Atmosphere using Broadband Radiometry (SABER), Halogen Occultation Experiment (HALOE), Atmospheric Chemistry Experiment (ACE), radiosondes) indicate that MLS temperature has persistent, pressure-dependent biases which are between −2.5 K and +1 K between 316 hPa and 10 hPa. The 100-hPa MLS v2.2 GPH surface has a bias of ∼150 m relative to the GEOS-5 values. These biases are compared to modeled systematic uncertainties. GPH biases relative to correlative measurements generally increase with height owing to an overall cold bias in MLS temperature relative to correlative temperature measurements in the upper stratosphere and mesosphere

The Atmospheric Infrared Sounder Version 6 Cloud Products

The version 6 cloud products of the Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Sounding Unit (AMSU) instrument suite are described. The cloud top temperature, pressure, and height and effective cloud fraction are now reported at the AIRS field-of-view (FOV) resolution. Significant improvements in cloud height assignment over version 5 are shown with FOV-scale comparisons to cloud vertical structure observed by the CloudSat 94 GHz radar and the Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP). Cloud thermodynamic phase (ice, liquid, and unknown phase), ice cloud effective diameter D(sub e), and ice cloud optical thickness () are derived using an optimal estimation methodology for AIRS FOVs, and global distributions for 2007 are presented. The largest values of tau are found in the storm tracks and near convection in the tropics, while D(sub e) is largest on the equatorial side of the midlatitude storm tracks in both hemispheres, and lowest in tropical thin cirrus and the winter polar atmosphere. Over the Maritime Continent the diurnal variability of tau is significantly larger than for the total cloud fraction, ice cloud frequency, and D(sub e), and is anchored to the island archipelago morphology. Important differences are described between northern and southern hemispheric midlatitude cyclones using storm center composites. The infrared-based cloud retrievals of AIRS provide unique, decadal-scale and global observations of clouds over portions of the diurnal and annual cycles, and capture variability within the mesoscale and synoptic scales at all latitudes

EOS Data and Information System (EOSDIS)

In the past decade, science and technology have reached levels that permit assessments of global environmental change. Scientific success in understanding global environmental change depends on integration and management of numerous data sources. The Global Change Data and Information System (GCDIS) must provide for the management of data, information dissemination, and technology transfer. The Earth Observing System Data and Information System (EOSDIS) is NASA's portion of this global change information system

Validation of the Aura Microwave Limb Sounder Temperature and Geopotential Height Measurements

This paper describes the retrievals algorithm used to determine temperature and height from radiance measurements by the Microwave Limb Sounder on EOS Aura. MLS is a "limbscanning" instrument, meaning that it views the atmosphere along paths that do not intersect the surface - it actually looks forwards from the Aura satellite. This means that the temperature retrievals are for a "profile" of the atmosphere somewhat ahead of the satellite. Because of the need to view a finite sample of the atmosphere, the sample spans a box about 1.5km deep and several tens of kilometers in width; the optical characteristics of the atmosphere mean that the sample is representative of a tube about 200-300km long in the direction of view. The retrievals use temperature analyses from NASA's Goddard Earth Observing System, Version 5 (GEOS-5) data assimilation system as a priori states. The temperature retrievals are somewhat deperrde~zt on these a priori states, especially in the lower stratosphere. An important part of the validation of any new dataset involves comparison with other, independent datasets. A large part of this study is concerned with such comparisons, using a number of independent space-based measurements obtained using different techniques, and with meteorological analyses. The MLS temperature data are shown to have biases that vary with height, but also depend on the validation dataset. MLS data are apparently biased slightly cold relative to correlative data in the upper troposphere and slightly warm in the middle stratosphere. A warm MLS bias in the upper stratosphere may be due to a cold bias in GEOS-5 temperatures

Preliminary Evaluation of Influence of Aerosols on the Simulation of Brightness Temperature in the NASA's Goddard Earth Observing System Atmospheric Data Assimilation System

This document reports on preliminary results obtained when studying the impact of aerosols on the calculation of brightness temperature (BT) for satellite infrared (IR) instruments that are currently assimilated in a 3DVAR configuration of Goddard Earth Observing System (GEOS)-atmospheric data assimilation system (ADAS). A set of fifteen aerosol species simulated by the Goddard Chemistry Aerosol Radiation and Transport (GOCART) model is used to evaluate the influence of the aerosol fields on the Community Radiative Transfer Model (CRTM) calculations taking place in the observation operators of the Gridpoint Statistical Interpolation (GSI) analysis system of GEOSADAS. Results indicate that taking aerosols into account in the BT calculation improves the fit to observations over regions with significant amounts of dust. The cooling effect obtained with the aerosol-affected BT leads to a slight warming of the analyzed surface temperature (by about 0:5oK) in the tropical Atlantic ocean (off northwest Africa), whereas the effect on the air temperature aloft is negligible. In addition, this study identifies a few technical issues to be addressed in future work if aerosol-affected BT are to be implemented in reanalysis and operational settings. The computational cost of applying CRTM aerosol absorption and scattering options is too high to justify their use, given the size of the benefits obtained. Furthermore, the differentiation between clouds and aerosols in GSI cloud detection procedures needs satisfactory revision

Remote sensing of tropical tropopause layer radiation balance using A-train measurements

Determining the level of zero net radiative heating (LZH) is critical to understanding parcel trajectory in the Tropical Tropopause Layer (TTL) and associated stratospheric hydration processes. Previous studies of the TTL radiative balance have focused on using radiosonde data, but remote sensing measurements from polar-orbiting satellites may provide the relevant horizontal and vertical information for assessing TTL solar heating and infrared cooling rates, especially across the Pacific Ocean. CloudSat provides a considerable amount of vertical information about the distribution of cloud properties relevant to heating rate analysis. The ability of CloudSat measurements and ancillary information to constrain LZH is explored. We employ formal error propagation analysis for derived heating rate uncertainty given the CloudSat cloud property retrieval algorithms. Estimation of the LZH to within approximately 0.5 to 1 km is achievable with CloudSat, but it has a low-altitude bias because the radar is unable to detect thin cirrus. This can be remedied with the proper utilization of Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) lidar backscatter information. By utilizing an orbital simulation with the GISS data set, we explore the representativeness of non-cross-track scanning active sounders in terms of describing the LZH distribution. In order to supplement CloudSat, we explore the ability of Atmospheric Infrared Sounder (AIRS) and Advanced Microwave Scanning Radiometer-EOS (AMSR-E) to constrain LZH and find that these passive sounders are useful where the cloud top height does not exceed 7 km. The spatiotemporal distributions of LZH derived from CloudSat and CALIPSO measurements are presented which suggest that thin cirrus have a limited effect on LZH mean values but affect LZH variability

Estimation of the vertical profile of sulfur dioxide injection into the atmosphere by a volcanic eruption using satellite column measurements and inverse transport modeling

International audienceAn analytical inversion method has been developed to estimate the vertical profile of SO2 emissions from volcanic eruptions. The method uses satellite-observed total SO2 columns and an atmospheric transport model (FLEXPART) to exploit the fact that winds change with altitude ? thus, the position and shape of the volcanic plume bear information on its emission altitude. The method finds the vertical emission distribution which minimizes the total difference between simulated and observed SO2 columns while also considering a priori information. We have tested the method with the eruption of Jebel at Tair on 30 September 2007 for which a comprehensive observational data set from various satellite instruments (AIRS, OMI, SEVIRI, CALIPSO) is available. Using satellite data from the first 24 h after the eruption for the inversion, we found an emission maximum near 16 km above sea level (asl), and secondary maxima near 5, 9, 12 and 14 km a.s.l. 60% of the emission occurred above the tropopause. The emission profile obtained in the inversion was then used to simulate the transport of the plume over the following week. The modeled plume agrees very well with SO2 total columns observed by OMI, and its altitude and width agree mostly within 1?2 km with CALIPSO observations of stratospheric aerosol produced from the SO2. The inversion result is robust against various changes in both the a priori and the observations. Even when using only SEVIRI data from the first 15 h after the eruption, the emission profile was reasonably well estimated. The method is computationally very fast. It is therefore suitable for implementation within an operational environment, such as the Volcanic Ash Advisory Centers, to predict the threat posed by volcanic ash for air traffic. It could also be helpful for assessing the sulfur input into the stratosphere, be it in the context of volcanic processes or also for proposed geo-engineering techniques to counteract global warming

Contrail cirrus supporting areas in model and observations

Contrails form and persist dependent on the surrounding moisture, temperature and pressure fields and on fuel and aircraft specific variables. After formation, contrail persistence requires only supersaturation relative to ice. The fractional area in which contrails can form is called potential contrail coverage. We introduce a potential contrail cirrus coverage equivalent to the cloud free supersaturated area. This field, simulated by the ECHAM4 climate model, agrees fairly well with estimates of supersaturation frequency as inferred from aircraft and satellite measurements. In areas where the two potential coverages are different, especially at lower flight levels, potential contrail coverage is not a valid estimate of maximum attainable contrail cirrus coverage. We parameterize both potential coverages consistently with the ECHAM4 cloud cover parameterization. A comparison of the potential contrail coverage with an earlier estimate reveals substantial differences especially at upper height levels in the tropics

The GEOS-5 Data Assimilation System-Documentation of Versions 5.0.1, 5.1.0, and 5.2.0

This report documents the GEOS-5 global atmospheric model and data assimilation system (DAS), including the versions 5.0.1, 5.1.0, and 5.2.0, which have been implemented in products distributed for use by various NASA instrument team algorithms and ultimately for the Modem Era Retrospective analysis for Research and Applications (MERRA). The DAS is the integration of the GEOS-5 atmospheric model with the Gridpoint Statistical Interpolation (GSI) Analysis, a joint analysis system developed by the NOAA/National Centers for Environmental Prediction and the NASA/Global Modeling and Assimilation Office. The primary performance drivers for the GEOS DAS are temperature and moisture fields suitable for the EOS instrument teams, wind fields for the transport studies of the stratospheric and tropospheric chemistry communities, and climate-quality analyses to support studies of the hydrological cycle through MERRA. The GEOS-5 atmospheric model has been approved for open source release and is available from: http://opensource.gsfc.nasa.gov/projects/GEOS-5/GEOS-5.php

Observed and Modeled Mountain Waves from the Surface to the Mesosphere Near the Drake Passage

Four state-of-the-science numerical weather prediction (NWP) models were used to perform mountain wave- (MW) resolving hind-casts over the Drake Passage of a 10-day period in 2010 with numerous observed MW cases. The Integrated Forecast System (IFS) and the Icosahedral Nonhydrostatic (ICON) model were run at Δx ≈ 9 and 13 km globally. TheWeather Research and Forecasting (WRF) model and the Met Office Unified Model (UM) were both configured with a Δx = 3 km regional domain. All domains had tops near 1 Pa (z ≈ 80 km). These deep domains allowed quantitative validation against Atmospheric InfraRed Sounder (AIRS) observations, accounting for observation time, viewing geometry, and radiative transfer. All models reproduced observed middle-atmosphere MWs with remarkable skill. Increased horizontal resolution improved validations. Still, all models underrepresented observed MW amplitudes, even after accounting for model effective resolution and instrument noise, suggesting even at Δx ≈ 3 km resolution, small-scale MWs are under-resolved and/or over-diffused. MWdrag parameterizations are still necessary in NWP models at current operational resolutions of Δx ≈ 10 km. Upper GW sponge layers in the operationally configured models significantly, artificially reduced MW amplitudes in the upper stratosphere and mesosphere. In the IFS, parameterized GW drags partly compensated this deficiency, but still, total drags were ≈ 6 time smaller than that resolved at Δx ≈ 3 km. Meridionally propagating MWs significantly enhance zonal drag over the Drake Passage. Interestingly, drag associated with meridional fluxes of zonal momentum (i.e. u'v') were important; not accounting for these terms results in a drag in the wrong direction at and below the polar night jet