The Modern-Era Retrospective Analysis for Research and Applications (MERRA) was undertaken by NASA’s Global Modeling and Assimilation Office with two primary objectives: to place observations from NASA’s Earth Observing System satellites into a climate context and to improve upon the hydrologic cycle represented in earlier generations of reanalyses. Focusing on the satellite era, from 1979 to the present, MERRA has achieved its goals with significant improvements in precipitation and water vapor climatology. Here, a brief overview of the system and some aspects of its performance, including quality assessment diagnostics from innovation and residual statistics, is given. By comparing MERRA with other updated reanalyses [the interim version of the next ECMWF Re-Analysis (ERA-Interim) and the Climate Forecast System Reanalysis (CFSR)], advances made in this new generation of reanalyses, as well as remaining deficiencies, are identified. Although there is little difference between the new reanalyses in many aspects of climate variability, substantial differences remain in poorly constrained quantities such as precipitation and surface fluxes. These differences, due to variations both in the models and in the analysis techniques, are an important measure of the uncertainty in reanalysis products. It is also found that all reanalyses are still quite sensitive to observing system changes. Dealing with this sensitivity remains the most pressing challenge for the next generation of reanalyses. Production has now caught up to the current period and MERRA is being continued as a near-real-time climate analysis. The output is available online through the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC).

An overview of the Community Earth System Model Version 2 (CESM2) is provided, including a discussion of the challenges encountered during its development and how they were addressed. In addition, an evaluation of a pair of CESM2 long preindustrial control and historical ensemble simulations is presented. These simulations were performed using the nominal 1° horizontal resolution configuration of the coupled model with both the “low‐top” (40 km, with limited chemistry) and “high‐top” (130 km, with comprehensive chemistry) versions of the atmospheric component. CESM2 contains many substantial science and infrastructure improvements and new capabilities since its previous major release, CESM1, resulting in improved historical simulations in comparison to CESM1 and available observations. These include major reductions in low‐latitude precipitation and shortwave cloud forcing biases; better representation of the Madden‐Julian Oscillation; better El Niño‐Southern Oscillation‐related teleconnections; and a global land carbon accumulation trend that agrees well with observationally based estimates. Most tropospheric and surface features of the low‐ and high‐top simulations are very similar to each other, so these improvements are present in both configurations. CESM2 has an equilibrium climate sensitivity of 5.1–5.3 °C, larger than in CESM1, primarily due to a combination of relatively small changes to cloud microphysics and boundary layer parameters. In contrast, CESM2's transient climate response of 1.9–2.0 °C is comparable to that of CESM1. The model outputs from these and many other simulations are available to the research community, and they represent CESM2's contributions to the Coupled Model Intercomparison Project Phase 6.

Abstract. The Modern-Era Retrospective Analysis for Research and Applications-2 (MERRA2) version of the Goddard Earth Observing System-5 (GEOS-5) atmospheric general circulation model (AGCM) is currently in use in the NASA Global Modeling and Assimilation Office (GMAO) at a wide range of resolutions for a variety of applications. Details of the changes in parameterizations subsequent to the version in the original MERRA reanalysis are presented here. Results of a series of atmosphere-only sensitivity studies are shown to demonstrate changes in simulated climate associated with specific changes in physical parameterizations, and the impact of the newly implemented resolution-aware behavior on simulations at different resolutions is demonstrated. The GEOS-5 AGCM presented here is the model used as part of the GMAO MERRA2 reanalysis, global mesoscale simulations at 10 km resolution through 1.5 km resolution, the real-time numerical weather prediction system, and for atmosphere-only, coupled ocean-atmosphere and coupled atmosphere-chemistry simulations.The seasonal mean climate of the MERRA2 version of the GEOS-5 AGCM represents a substantial improvement over the simulated climate of the MERRA version at all resolutions and for all applications. Fundamental improvements in simulated climate are associated with the increased reevaporation of frozen precipitation and cloud condensate, resulting in a wetter atmosphere. Improvements in simulated climate are also shown to be attributable to changes in the background gravity wave drag, and to upgrades in the relationship between the ocean surface stress and the ocean roughness. The series of resolution-aware parameters related to the moist physics was shown to result in improvements at higher resolutions and result in AGCM simulations that exhibit seamless behavior across different resolutions and applications.

[1] Present-day shortcomings in the representation of upper tropospheric ice clouds in general circulation models (GCMs) lead to errors in weather and climate forecasts as well as account for a source of uncertainty in climate change projections. An ongoing challenge in rectifying these shortcomings has been the availability of adequate, high-quality, global observations targeting ice clouds and related precipitating hydrometeors. In addition, the inadequacy of the modeled physics and the often disjointed nature between model representation and the characteristics of the retrieved/observed values have hampered GCM development and validation efforts from making effective use of the measurements that have been available. Thus, even though parameterizations in GCMs accounting for cloud ice processes have, in some cases, become more sophisticated in recent years, this development has largely occurred independently of the global-scale measurements. With the relatively recent addition of satellite-derived products from Aura/Microwave Limb Sounder (MLS) and CloudSat, there are now considerably more resources with new and unique capabilities to evaluate GCMs. In this article, we illustrate the shortcomings evident in model representations of cloud ice through a comparison of the simulations assessed in the Intergovernmental Panel on Climate Change Fourth Assessment Report, briefly discuss the range of global observational resources that are available, and describe the essential components of the model parameterizations that characterize their ''cloud'' ice and related fields. Using this information as background, we (1) discuss some of the main considerations and cautions that must be taken into account in making model-data comparisons related to cloud ice, (2) illustrate present progress and uncertainties in applying satellite cloud ice (namely from MLS and CloudSat) to model diagnosis, (3) show some indications of model improvements, and finally (4) discuss a number of remaining questions and suggestions for pathways forward.

[1] The Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere (CRISTA) instrument measured stratospheric temperatures and trace species concentrations with high precision and spatial resolution during two missions. The measuring technique is infrared limb-sounding of optically thin emissions. In a general approach, we investigate the applicability of the technique to measure gravity waves (GWs) in the retrieved temperature data. It is shown that GWs with wavelengths of the order of 100-200 km horizontally can be detected. The results are applicable to any instrument using the same technique. We discuss additional constraints inherent to the CRISTA instrument. The vertical field of view and the influence of the sampling and retrieval imply that waves with vertical wavelengths $3-5 km or larger can be retrieved. Global distributions of GW fluctuations were extracted from temperature data measured by CRISTA using Maximum Entropy Method (MEM) and Harmonic Analysis (HA), yielding height profiles of vertical wavelength and peak amplitude for fluctuations in each scanned profile. The method is discussed and compared to Fourier transform analyses and standard deviations. Analysis of data from the first mission reveals large GW amplitudes in the stratosphere over southernmost South America. These waves obey the dispersion relation for linear two-dimensional mountain waves (MWs). The horizontal structure on 6 November 1994 is compared to temperature fields calculated by the Pennsylvania State University (PSU)/National Center for Atmospheric Research (NCAR) mesoscale model (MM5). It is demonstrated that precise knowledge of the instrument's sensitivity is essential. Particularly good agreement is found at the southern tip of South America where the MM5 accurately reproduces the amplitudes and phases of a large-scale wave with 400 km horizontal wavelength. Targeted ray-tracing simulations allow us to interpret some of the observed wave features. A companion paper will discuss MWs on a global scale and estimates the fraction that MWs contribute to the total GW energy (Preusse et al., in preparation, 2002).

International audienceFor the first time, a formal comparison is made between gravity wave momentum fluxes in models and those derived from observations. Although gravity waves occur over a wide range of spatial and temporal scales, the focus of this paper is on scales that are being parameterized in present climate models, sub-1000-km scales. Only observational methods that permit derivation of gravity wave momentum fluxes over large geographical areas are discussed, and these are from satellite temperature measurements, constant-density long-duration balloons, and high-vertical-resolution radiosonde data. The models discussed include two high-resolution models in which gravity waves are explicitly modeled, Kanto and the Community Atmosphere Model, version 5 (CAM5), and three climate models containing gravity wave parameterizations, MAECHAM5, Hadley Centre Global Environmental Model 3 (HadGEM3), and the Goddard Institute for Space Studies (GISS) model. Measurements generally show similar flux magnitudes as in models, except that the fluxes derived from satellite measurements fall off more rapidly with height. This is likely due to limitations on the observable range of wavelengths, although other factors may contribute. When one accounts for this more rapid fall off, the geographical distribution of the fluxes from observations and models compare reasonably well, except for certain features that depend on the specification of the nonorographic gravity wave source functions in the climate models. For instance, both the observed fluxes and those in the high-resolution models are very small at summer high latitudes, but this is not the case for some of the climate models. This comparison between gravity wave fluxes from climate models, high-resolution models, and fluxes derived from observations indicates that such efforts offer a promising path toward improving specifications of gravity wave sources in climate models

During the 1930s, the United States experienced one of the most devastating droughts of the past century. The drought affected almost two-thirds of the country and parts of Mexico and Canada and was infamous for the numerous dust storms that occurred in the southern Great Plains. In this study, we present model results that indicate that the drought was caused by anomalous tropical sea surface temperatures during that decade and that interactions between the atmosphere and the land surface increased its severity. We also contrast the 1930s drought with other North American droughts of the 20th century.

[1] To assess the current status of climate models in simulating clouds, basic cloud climatologies from ten atmospheric general circulation models are compared with satellite measurements from the International Satellite Cloud Climatology Project (ISCCP) and the Clouds and Earth's Radiant Energy System (CERES) program. An ISCCP simulator is employed in all models to facilitate the comparison. Models simulated a four-fold difference in high-top clouds. There are also, however, large uncertainties in satellite high thin clouds to effectively constrain the models. The majority of models only simulated 30-40% of middle-top clouds in the ISCCP and CERES data sets. Half of the models underestimated low clouds, while none overestimated them at a statistically significant level. When stratified in the optical thickness ranges, the majority of the models simulated optically thick clouds more than twice the satellite observations. Most models, however, underestimated optically intermediate and thin clouds. Compensations of these clouds biases are used to explain the simulated longwave and shortwave cloud radiative forcing at the top of the atmosphere. Seasonal sensitivities of clouds are also analyzed to compare with observations. Models are shown to simulate seasonal variations better for high clouds than for low clouds. Latitudinal distribution of the seasonal variations correlate with satellite measurements at >0.9, 0.6-0.9, and À0.2-0.7 levels for high, middle, and low clouds, respectively. The seasonal sensitivities of cloud types are found to strongly depend on the basic cloud climatology in the models. Models that systematically underestimate middle clouds also underestimate seasonal variations, while those that overestimate optically thick clouds also overestimate their seasonal sensitivities. Possible causes of the systematic cloud biases in the models are discussed.