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.

The Subseasonal Experiment (SubX) is a multimodel subseasonal prediction experiment designed around operational requirements with the goal of improving subseasonal forecasts. Seven global models have produced 17 years of retrospective (re)forecasts and more than a year of weekly real-time forecasts. The reforecasts and forecasts are archived at the Data Library of the International Research Institute for Climate and Society, Columbia University, providing a comprehensive database for research on subseasonal to seasonal predictability and predictions. The SubX models show skill for temperature and precipitation 3 weeks ahead of time in specific regions. The SubX multimodel ensemble mean is more skillful than any individual model overall. Skill in simulating the Madden-Julian oscillation (MJO) and the North Atlantic Oscillation (NAO), two sources of subseasonal predictability, is also evaluated, with skillful predictions of the MJO 4 weeks in advance and of the NAO 2 weeks in advance. SubX is also able to make useful contributions to operational forecast guidance at the Climate Prediction Center. Additionally, SubX provides information on the potential for extreme precipitation associated with tropical cyclones, which can help emergency management and aid organizations to plan for disasters.

Atmosphere: NASA GEOS-5 (internal tag: Heracles h54p1c180o05_2) run at 1⁄2 degree horizontal resolution and 72 vertical layers Ocean: Modular Ocean Model Version 5 (MOM5) of Geophysical Fluid Dynamics Laboratory run at 1⁄2 degree horizontal resolution and 40 vertical layers Land: Catchment land surface model Sea ice: Los Alamos Sea Ice model (CICE) Initialization of the reforecasts: Atmosphere from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2); ocean and the sea ice from GMAO ocean analysis; and land from MERRA-2 land fields which are based on observations-corrected precipitation. The reforecasts are initialized every 5 days with a total of 73 start dates per year, with 7 ensemble members per start date for 19992015. Each reforecast is run for 45 days. Prediction and Predictability of the Madden Julian Oscillation in the NASA GEOS-5 Seasonal-to-Subseasonal System Deepthi Achuthavarier1,2, Randal Koster2, Jelena Marshak2, Siegfried Schubert3 and Andrea Molod2 1Universities Space Research Association, Columbia, MD, 2Global Modeling and Assimilation Office, NASA Goddard Space Flight Center, Greenbelt, MD, 3Science Systems and Applications, Lanham, MD.

Tropical-extratropical teleconnections are considered key to advancing subseasonal prediction. The Madden Julian oscillation (MJO), characterized by large scale convective envelopes propagating along the tropical Indo-Pacific sector, is known to modulate midlatitude circulation and associated weather patterns. Although there is a general consensus on the MJO's influence on the midlatitude circulation, which is thought to be due to modulations of the North Atlantic Oscillation (NAO) and the Pacific North America (PNA) pattern, relatively less is known about the predictability of these teleconnection signals in dynamical forecast models. The composite evolution of the midlatitude circulation anomalies and associated wave train structure as the delayed response to tropical heating are reported in many studies that have examined reanalyses and long model simulations. However, it is yet to be determined whether they lend any beneficial subseasonal forecast skill, especially to weekly mean surface temperature and precipitation over North America. Investigating useful predictable signals from the MJO teleconnections is also complicated by the fact that the MJO is a moving heat source with an approximate periodicity of 30-60 days, and that the structure and amplitude of the midlatitude response can be sensitive to the longitudinal positioning of the heating anomaly as well as the propagation speed of the MJO. The objective of this study is to investigate the impact of MJO teleconnections on forecast accuracy at 2-3 week lead over North America, with an emphasis on the above-mentioned lesser known aspects of these teleconnections. To this end, we utilize a suite of subseasonal reforecasts performed with the latest NASA GEOS-5 seasonal-to-subseasonal (S2S) system. These reforecasts were performed as part of the NOAA SubX project, wherein the NASA GEOS-5 atmosphere-ocean coupled model was run at degree horizontal resolution, initialized every 5 days for the period 1999-2016. The GEOS-5 model shows skillful predictions of the MJO, with the correlation coefficient based on the real-time multivariate MJO (RMM) index staying at or above 0.5 up to forecast lead 26-36 days. The system is thus a useful tool for investigating MJO teleconnection processes.

Atmosphere: NASA GEOS run at 1⁄2 degree horizontal resolution and 72 vertical layers Ocean: Modular Ocean Model Version 5 (MOM5) of Geophysical Fluid Dynamics Laboratory run at 1⁄2 degree horizontal resolution and 40 vertical layers Land: Catchment land surface model Sea ice: Los Alamos Sea Ice model (CICE) Initialization of the reforecasts: Atmosphere from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2); ocean and the sea ice from GMAO ocean analysis; and land from MERRA-2 land fields which are based on observations-corrected precipitation. The reforecasts are initialized every 5 days with a total of 73 start dates per year. Prediction skill of the MJO teleconnection signals in the NASA GEOS subseasonal reforecasts Deepthi Achuthavarier1,2, Randal Koster2, Jelena Marshak2, Siegfried Schubert3 and Andrea Molod2 1Universities Space Research Association, Columbia, MD, 2Global Modeling and Assimilation Office, NASA Goddard Space Flight Center, Greenbelt, MD, 3Science Systems and Applications, Lanham, MD.

In this presentation we present an overview of the GMAO Sub-Seasonal and Seasonal Prediction System with a focus on the computing time and resources and actual time it takes to complete a full set of hindcasts. The goal is to come up with some solutions to allow us to run more ensemble members for the next version of the system which will be higher resolution and take many more resources.

In this study, we examine the prediction skill and predictability of the Madden Julian Oscillation (MJO) in a recent version of the NASA GEOS-5 atmosphere-ocean coupled model run at at 1/2 degree horizontal resolution. The results are based on a suite of hindcasts produced as part of the NOAA SubX project, consisting of seven ensemble members initialized every 5 days for the period 1999-2015. The atmospheric initial conditions were taken from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), and the ocean and the sea ice were taken from a GMAO ocean analysis. The land states were initialized from the MERRA-2 land output, which is based on observation-corrected precipitation fields. We investigated the MJO prediction skill in terms of the bivariate correlation coefficient for the real-time multivariate MJO (RMM) indices. The correlation coefficient stays at or above 0.5 out to forecast lead times of 26-36 days, with a pronounced increase in skill for forecasts initialized from phase 3, when the MJO convective anomaly is located in the central tropical Indian Ocean. A corresponding estimate of the upper limit of the predictability is calculated by considering a single ensemble member as the truth and verifying the ensemble mean of the remaining members against that. The predictability estimates fall between 35-37 days (taken as forecast lead when the correlation reaches 0.5) and are rather insensitive to the initial MJO phase. The model shows slightly higher skill when the initial conditions contain strong MJO events compared to weak events, although the difference in skill is evident only from lead 1 to 20. Similar to other models, the RMM-index-based skill arises mostly from the circulation components of the index. The skill of the convective component of the index drops to 0.5 by day 20 as opposed to day 30 for circulation fields. The propagation of the MJO anomalies over the Maritime Continent does not appear problematic in the GEOS-5 hindcasts implying that the Maritime Continent predictability barrier may not be a major concern in this model. Finally, the MJO prediction skill in this version of GEOS-5 is superior to that of the current seasonal prediction system at the GMAO; this could be partly attributed to a slightly better representation of the MJO in the free running version of this model and partly to the improved atmospheric initialization from MERRA-2.

This study examines the mechanisms of the Pacific decadal oscillation (PDO) in the NASA GEOS-5 general circulation model (GCM). Similar to several other state-of-the-art GCMs, the El Niño-Southern Oscillation (ENSO) simulated by the GEOS-5 has a strong biennial periodicity. Since this is a model bias that precludes a strong role of ENSO, it provides a unique environment to assess the other leading mechanisms of North Pacific decadal variability. Despite the biennial ENSO periodicity, the model simulates a realistic PDO pattern in the North Pacific that is resolved as the first empirical orthogonal function (EOF) of winter mean sea surface temperature (SST). The spectrum of the PDO indicates no preferred periodicity. The SST anomalies associated with the PDO, particularly its basin wide structure, are primarily forced by the Aleutian low through Ekman transport. The slow geostrophic transport in association with the meridional adjustment of the subtropical gyre is limited to a narrow region in the Kuroshio–Oyashio extension, north of 40°N. The atmosphere’s response to the PDO, while weak, projects onto the North Pacific Oscillation (NPO), a meridional dipole in sea level pressure. Both the lack of preferred periodicity and the weak atmospheric response indicate an air-sea coupled oscillation is an unlikely mechanism in this model. In agreement with recent studies, the NPO is correlated with the North Pacific Gyre Oscillation, which is another leading EOF of North Pacific SST variability. The results emphasize the role of atmospheric variability in the North Pacific SST modes, thereby bringing into question the potential for their predictability.

The daily surface and atmospheric radiative fluxes from NASA Clouds and the Earth's Radiant Energy System (CERES) Synoptic 1 degree (SYN1deg) Ed3A are among the most widely used data to study cloud-radiative feedback. The CERES SYN1deg data are based on Fu-Liou radiative transfer computations that use specific humidity (Q) and air temperature (T) from NASA Global Modeling and Assimilation Office (GMAO) reanalyses as inputs and are therefore subject to the quality of those fields. This study uses in situ Q and T observations collected during the Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign to augment the input stream used in the NASA GMAO reanalysis and assess the impact on the CERES daily surface and atmospheric longwave estimates. The results show that the assimilation of DYNAMO observations considerably improves the vertical profiles of analyzed Q and T over and near DYNAMO stations by moistening and warming the lower troposphere and upper troposphere and drying and cooling the mid-upper troposphere. As a result of these changes in Q and T, the computed CERES daily surface downward longwave flux increases by about 5Wm(-2), due mainly to the warming and moistening in the lower troposphere; the computed daily top-of-atmosphere (TOA) outgoing longwave radiation increases by 2-3Wm(-2) during dry periods only. Correspondingly, the estimated local atmospheric longwave radiative cooling enhances by about 5Wm(-2) (7-8Wm(-2)) during wet (dry) periods. These changes reduce the bias in the CERES SYN1deg-like daily longwave estimates at both the TOA and surface and represent an improvement over the DYNAMO region.

The National Aeronautics and Space Administration (NASA) nature run (NR), released for use in observing system simulation experiments (OSSEs), is a 2-yr-long global nonhydrostatic free-running simulation at a horizontal resolution of 7 km, forced by observed sea surface temperatures (SSTs) and sea ice, and inclusive of interactive aerosols and trace gases. This article evaluates the NR with respect to tropical cyclone (TC) activity. It is emphasized that to serve as an NR, a long-term simulation must be able to produce realistic TCs, which arise out of realistic large-scale forcings. The presence in the NR of the relevant dynamical features over the African monsoon region and the tropical Atlantic is confirmed, along with realistic African easterly wave activity. The NR Atlantic TC seasons, produced with 2005 and 2006 SSTs, show interannual variability consistent with observations, with much stronger activity in 2005. An investigation of TC activity over all the other basins (eastern and western North Pacific Ocean, north and south Indian Ocean, and Australian region), together with important elements of the atmospheric circulation, such as the Somali jet and westerly bursts, reveals that the model captures the fundamental aspects of TC seasons in every basin, producing a realistic number of TCs with realistic tracks, life spans, and structures. This confirms that the NASA NR is a very suitable tool for OSSEs targeting TCs and represents an improvement with respect to previous long simulations that have served the global atmospheric OSSE community.

This study examines the impact of the Dynamics of the Madden-Julian Oscillation (DYNAMO) campaign in situ observations on NASA Goddard Earth Observing System version 5 (GEOS-5) reanalyses and the improvements gained thereby in the representation of the Madden-Julian Oscillation (MJO) initiation processes. To this end, we produced a global, high-resolution (1/4 degrees spatially) reanalysis that assimilates the level-4, quality-controlled DYNAMO upper air soundings from about 87 stations in the equatorial Indian Ocean region along with a companion data-denied control reanalysis. The DYNAMO reanalysis produces a more realistic vertical structure of the temperature and moisture in the central tropical Indian Ocean by correcting the model biases, namely, the cold and dry biases in the lower troposphere and warm bias in the upper troposphere. The reanalysis horizontal winds are substantially improved, in that, the westerly acceleration and vertical shear of the zonal wind are enhanced. The DYNAMO reanalysis shows enhanced low-level diabatic heating, moisture anomalies and vertical velocity during the MJO initiation. Due to the warmer lower troposphere, the deep convection is invigorated, which is evident in convective cloud fraction. The GEOS-5 atmospheric general circulation model (AGCM) employed in the reanalysis is overall successful in assimilating the additional DYNAMO observations, except for an erroneous model response for medium rain rates, between 700 and 600hPa, reminiscent of a bias in earlier versions of the AGCM. The moist heating profile shows a sharp decrease there due to the excessive convective rain re-evaporation, which is partly offset by the temperature increment produced by the analysis.

This report documents an evaluation by the Global Modeling and Assimilation Office (GMAO) of a two-year 7-km-resolution non-hydrostatic global mesoscale simulation produced with the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model. The simulation was produced as a Nature Run for conducting observing system simulation experiments (OSSEs). Generation of the GEOS-5 Nature Run (G5NR) was motivated in part by the desire of the OSSE community for an improved high-resolution sequel to an existing Nature Run produced by the European Centre for Medium-Range Weather Forecasts (ECMWF), which has served the community for several years. The intended use of the G5NR in this context is for generating simulated observations to test proposed observing system designs regarding new instruments and their deployments. Because NASA's interest in OSSEs extends beyond traditional weather forecasting applications, the G5NR includes, in addition to standard meteorological components, a suite of aerosol types and several trace gas concentrations, with emissions downscaled to 10 km using ancillary information such as power plant location, population density and night-light information. The evaluation exercise described here involved more than twenty-five GMAO scientists investigating various aspects of the G5NR performance, including time mean temperature and wind fields, energy spectra, precipitation and the hydrological cycle, the representation of waves, tropical cyclones and midlatitude storms, land and ocean surface characteristics, the representation and forcing effects of clouds and radiation, dynamics of the stratosphere and mesosphere, and the representation of aerosols and trace gases. Comparisons are made with observational data sets when possible, as well as with reanalyses and other long model simulations. The evaluation is broad in scope, as it is meant to assess the overall realism of basic aspects of the G5NR deemed relevant to the conduct of OSSEs. However, because of the relatively short record and other practical considerations, these comparisons cannot provide a definitive, statistically sound assessment of all model deficiencies, or guarantee the G5NR's suitability for all OSSE applications. Differences between the observed and simulated behavior also must be judged in the context of basic internal atmospheric variability which can introduce variations that are not necessarily controlled by the prescribed sea surface temperatures used in generating the G5NR. The results show that the G5NR performs well as measured by the majority of metrics applied in this evaluation. Particular benefits derived from the 7-km resolution of G5NR include realistic representations of extreme weather events in both the tropics and extratropics including tropical cyclones, Nor'easters and mesoscale convective complexes; improved representation of the diurnal cycle of precipitation over land; well-resolved surface-atmosphere interactions such as katabatic wind flows over Antarctica and Greenland; and resolution of orographically generated gravity waves that propagate into the upper atmosphere and influence the large scale circulation. Obvious deficiencies in the G5NR include a "splitting" of the inter-tropical convergence zone, which leads to a weaker-than-observed Hadley circulation and related deficiencies in the depiction of stationary wave patterns. Also, while the G5NR captures global cloud features and radiative effects well in general, close comparison with observations reveals higher-than-observed cloud brightness, likely due to an overabundance of cloud condensate; less distinct cloud minima in subtropical subsidence zones, consistent with a weak Hadley circualtion; and too few near-coastal marine stratocumulus clouds.

How tropical cyclone (TC) activity in the northwestern Pacific might change in a future climate is assessed using multidecadal Atmospheric Model Intercomparison Project (AMIP)-style and time-slice simulations with the ECMWF Integrated Forecast System (IFS) at 16-km and 125-km global resolution. Both models reproduce many aspects of the present-day TC climatology and variability well, although the 16-km IFS is far more skillful in simulating the full intensity distribution and genesis locations, including their changes in response to El Nino-Southern Oscillation. Both IFS models project a small change in TC frequency at the end of the twenty-first century related to distinct shifts in genesis locations. In the 16-km IFS, this shift is southward and is likely driven by the southeastward penetration of the monsoon trough/subtropical high circulation system and the southward shift in activity of the synoptic-scale tropical disturbances in response to the strengthening of deep convective activity over the central equatorial Pacific in a future climate. The 16-km IFS also projects about a 50% increase in the power dissipation index, mainly due to significant increases in the frequency of the more intense storms, which is comparable to the natural variability in the model. Based on composite analysis of large samples of supertyphoons, both the development rate and the peak intensities of these storms increase in a future climate, which is consistent with their tendency to develop more to the south, within an environment that is thermodynamically more favorable for faster development and higher intensities. Coherent changes in the vertical structure of supertyphoon composites show system-scale amplification of the primary and secondary circulations with signs of contraction, a deeper warm core, and an upward shift in the outflow layer and the frequency of the most intense updrafts. Considering the large differences in the projections of TC intensity change between the 16-km and 125-km IFS, this study further emphasizes the need for high-resolution modeling in assessing potential changes in TC activity.