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 re-evaporation 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.

The climate research community uses atmospheric reanalysis data sets to understand a wide range of processes and variability in the atmosphere, yet different reanalyses may give very different results for the same diagnostics. The Stratosphere–troposphere Processes And their Role in Climate (SPARC) Reanalysis Intercomparison Project (S-RIP) is a coordinated activity to compare reanalysis data sets using a variety of key diagnostics. The objectives of this project are to identify differences among reanalyses and understand their underlying causes, to provide guidance on appropriate usage of various reanalysis products in scientific studies, particularly those of relevance to SPARC, and to contribute to future improvements in the reanalysis products by establishing collaborative links between reanalysis centres and data users. The project focuses predominantly on differences among reanalyses, although studies that include operational analyses and studies comparing reanalyses with observations are also included when appropriate. The emphasis is on diagnostics of the upper troposphere, stratosphere, and lower mesosphere. This paper summarizes the motivation and goals of the S-RIP activity and extensively reviews key technical aspects of the reanalysis data sets that are the focus of this activity. The special issue The SPARC Reanalysis Intercomparison Project (S-RIP) in this journal serves to collect research with relevance to the S-RIP in preparation for the publication of the planned two (interim and full) S-RIP reports.

Global modeling of atmospheric chemistry is a grand computational challenge because of the need to simulate large coupled systems of ∼100–1000 chemical species interacting with transport on all scales. Offline chemical transport models (CTMs), where the chemical continuity equations are solved using meteorological data as input, have usability advantages and are important vehicles for developing atmospheric chemistry knowledge that can then be transferred to Earth system models. However, they have generally not been designed to take advantage of massively parallel computing architectures. Here, we develop such a high-performance capability for GEOS-Chem (GCHP), a CTM driven by meteorological data from the NASA Goddard Earth Observation System (GEOS) and used by hundreds of research groups worldwide. GCHP is a grid-independent implementation of GEOS-Chem using the Earth System Modeling Framework (ESMF) that permits the same standard model to operate in a distributed-memory framework for massive parallelization. GCHP also allows GEOS-Chem to take advantage of the native GEOS cubed-sphere grid for greater accuracy and computational efficiency in simulating transport. GCHP enables GEOS-Chem simulations to be conducted with high computational scalability up to at least 500 cores, so that global simulations of stratosphere–troposphere oxidant–aerosol chemistry at C180 spatial resolution (∼0.5∘×0.625∘) or finer become routinely feasible.

The Global Modeling and Assimilation Office (GMAO) has recently released a new version of the Goddard Earth Observing System (GEOS) Subseasonal to Seasonal prediction (S2S) system, GEOS‐S2S‐2, that represents a substantial improvement in performance and infrastructure over the previous system. The system is described here in detail, and results are presented from forecasts, climate equillibrium simulations, and data assimilation experiments. The climate or equillibrium state of the atmosphere and ocean showed a substantial reduction in bias relative to GEOS‐S2S‐1. The GEOS‐S2S‐2 coupled reanalysis also showed substantial improvements, attributed to the assimilation of along‐track absolute dynamic topography. The forecast skill on subseasonal scales showed a much improved prediction of the Madden‐Julian Oscillation in GEOS‐S2S‐2, and on a seasonal scale the tropical Pacific forecasts show substantial improvement in the east and comparable skill to GEOS‐S2S‐1 in the central Pacific. GEOS‐S2S‐2 anomaly correlations of both land surface temperature and precipitation were comparable to GEOS‐S2S‐1 and showed substantially reduced root‐mean‐square error of surface temperature. The remaining issues described here are being addressed in the development of GEOS‐S2S Version 3, and with that system GMAO will continue its tradition of maintaining a state‐of‐the‐art seasonal prediction system for use in evaluating the impact on seasonal and decadal forecasts of assimilating newly available satellite observations, as well as evaluating additional sources of predictability in the Earth system through the expanded coupling of the Earth system model and assimilation components.

Global simulations of atmospheric chemistry are commonly conducted with off-line chemical transport models (CTMs) driven by archived meteorological data from general circulation models (GCMs). The off-line approach has the advantages of simplicity and expediency, but it incurs errors due to temporal averaging in the meteorological archive and the inability to reproduce the GCM transport algorithms exactly. The CTM simulation is also often conducted at coarser grid resolution than the parent GCM. Here we investigate this cascade of CTM errors by using ²²²Rn–²¹⁰Pb–⁷Be chemical tracer simulations off-line in the GEOS-Chem CTM at rectilinear 0.25^∘ × 0.3125^∘ (≈ 25 km) and 2^∘ × 2.5^∘ (≈ 200 km) resolutions and online in the parent GEOS-5 GCM at cubed-sphere c360 (≈ 25 km) and c48 (≈ 200 km) horizontal resolutions. The c360 GEOS-5 GCM meteorological archive, updated every 3 h and remapped to 0.25^∘ × 0.3125^∘, is the standard operational product generated by the NASA Global Modeling and Assimilation Office (GMAO) and used as input by GEOS-Chem. We find that the GEOS-Chem ²²²Rn simulation at native 0.25^∘ × 0.3125^∘ resolution is affected by vertical transport errors of up to 20 % relative to the GEOS-5 c360 online simulation, in part due to loss of transient organized vertical motions in the GCM (resolved convection) that are temporally averaged out in the 3 h meteorological archive. There is also significant error caused by operational remapping of the meteorological archive from a cubed-sphere to a rectilinear grid. Decreasing the GEOS-Chem resolution from 0.25^∘ × 0.3125^∘ to 2^∘ × 2.5^∘ induces further weakening of vertical transport as transient vertical motions are averaged out spatially and temporally. The resulting ²²²Rn concentrations simulated by the coarse-resolution GEOS-Chem are overestimated by up to 40 % in surface air relative to the online c360 simulations and underestimated by up to 40 % in the upper troposphere, while the tropospheric lifetimes of ²¹⁰Pb and ⁷Be against aerosol deposition are affected by 5–10 %. The lost vertical transport in the coarse-resolution GEOS-Chem simulation can be partly restored by recomputing the convective mass fluxes at the appropriate resolution to replace the archived convective mass fluxes and by correcting for bias in the spatial averaging of boundary layer mixing depths.

We examine differences among reanalysis high-cloud products in the tropics, assess the impacts of these differences on radiation budgets at the top of the atmosphere and within the tropical upper troposphere and lower stratosphere (UTLS), and discuss their possible origins in the context of the reanalysis models. We focus on the ERA5 (fifth-generation European Centre for Medium-range Weather Forecasts – ECMWF – reanalysis), ERA-Interim (ECMWF Interim Reanalysis), JRA-55 (Japanese 55-year Reanalysis), MERRA-2 (Modern-Era Retrospective Analysis for Research and Applications, Version 2), and CFSR/CFSv2 (Climate Forecast System Reanalysis/Climate Forecast System Version 2) reanalyses. As a general rule, JRA-55 produces the smallest tropical high-cloud fractions and cloud water contents among the reanalyses, while MERRA-2 produces the largest. Accordingly, long-wave cloud radiative effects are relatively weak in JRA-55 and relatively strong in MERRA-2. Only MERRA-2 and ERA5 among the reanalyses produce tropical-mean values of outgoing long-wave radiation (OLR) close to those observed, but ERA5 tends to underestimate cloud effects, while MERRA-2 tends to overestimate variability. ERA5 also produces distributions of long-wave, short-wave, and total cloud radiative effects at the top of the atmosphere that are very consistent with those observed. The other reanalyses all exhibit substantial biases in at least one of these metrics, although compensation between the long-wave and short-wave effects helps to constrain biases in the total cloud radiative effect for most reanalyses. The vertical distribution of cloud water content emerges as a key difference between ERA-Interim and other reanalyses. Whereas ERA-Interim shows a monotonic decrease of cloud water content with increasing height, the other reanalyses all produce distinct anvil layers. The latter is in better agreement with observations and yields very different profiles of radiative heating in the UTLS. For example, whereas the altitude of the level of zero net radiative heating tends to be lower in convective regions than in the rest of the tropics in ERA-Interim, the opposite is true for the other four reanalyses. Differences in cloud water content also help to explain systematic differences in radiative heating in the tropical lower stratosphere among the reanalyses. We discuss several ways in which aspects of the cloud and convection schemes impact the tropical environment. Discrepancies in the vertical profiles of temperature and specific humidity in convective regions are particularly noteworthy, as these variables are directly constrained by data assimilation, are widely used, and feed back to convective behaviour through their relationships with thermodynamic stability.