Aimed at reducing deficiencies in representing the Madden-Julian oscillation (MJO) in general circulation models (GCMs), a global model evaluation project on vertical structure and physical processes of the MJO was coordinated. In this paper, results from the climate simulation component of this project are reported. It is shown that the MJO remains a great challenge in these latest generation GCMs. The systematic eastward propagation of the MJO is only well simulated in about one fourth of the total participating models. The observed vertical westward tilt with altitude of the MJO is well simulated in good MJO models but not in the poor ones. Damped Kelvin wave responses to the east of convection in the lower troposphere could be responsible for the missing MJO preconditioning process in these poor MJO models. Several process-oriented diagnostics were conducted to discriminate key processes for realistic MJO simulations. While large-scale rainfall partition and low-level mean zonal winds over the Indo-Pacific in a model are not found to be closely associated with its MJO skill, two metrics, including the low-level relative humidity difference between high-and low-rain events and seasonal mean gross moist stability, exhibit statistically significant correlations with the MJO performance. It is further indicated that increased cloud-radiative feedback tends to be associated with reduced amplitude of intraseasonal variability, which is incompatible with the radiative instability theory previously proposed for the MJO. Results in this study confirm that inclusion of air-sea interaction can lead to significant improvement in simulating the MJO.

The current literature provides compelling evidence suggesting that an eddy-resolving (as opposed to eddy-permitting or eddy-parameterized) ocean component model will significantly impact the simulation of the large-scale climate, although this has not been fully tested to date in multi-decadal global coupled climate simulations. The purpose of this paper is to document how increased ocean model resolution impacts the simulation of large-scale climate variability. The model used for this study is the NCAR Community Climate System Model version 3.5 (CCSM3.5) -the forerunner to CCSM4. Two experiments are reported here. The first experiment (i.e., control) is a 155-year present-day climate simulation using a 0.5º atmosphere component (zonal resolution 0.625º meridional resolution 0.5º) coupled to ocean and sea-ice components with zonal resolution of 1.2º and meridional resolution varying from 0.27º at the equator to 0.54º in the mid-latitudes. The second simulation uses the same atmospheric model coupled to 0.1º ocean and sea-ice component models. The simulations are compared in terms of how the representation of smaller scale features in the time mean ocean circulation and ocean eddies impact the mean and variable climate. In terms of the global mean surface temperature, the enhanced ocean resolution leads to a ubiquitous surface warming of about 0.2 o C. The warming is largest in the Arctic and regions of strong ocean fronts and ocean eddy activity (i.e., Southern Ocean, western boundary currents). The Arctic warming is associated with significant losses of sea-ice in the high-resolution simulation. The SST gradients in the North Atlantic, in particular, are better resolved in the high-resolution model leading to significantly sharper temperature gradients and associated large-scale shifts in the rainfall. In the extra-tropics, the interannual sea surface temperature anomaly (SSTA) variability is increased with the resolved eddies, but decreases in the deep tropics (i.e., the variance of El Niño and the Southern Oscillation is reduced). Changes in global SSTA teleconnections and local air-sea feedback are also documented and show large changes in ocean-atmosphere coupling.3

The sensitivity to the horizontal resolution of the climate, anthropogenic climate change, and seasonal predictive skill of the ECMWF model has been studied as part of Project Athena—an international collaboration formed to test the hypothesis that substantial progress in simulating and predicting climate can be achieved if mesoscale and subsynoptic atmospheric phenomena are more realistically represented in climate models. In this study the experiments carried out with the ECMWF model (atmosphere only) are described in detail. Here, the focus is on the tropics and the Northern Hemisphere extratropics during boreal winter. The resolutions considered in Project Athena for the ECMWF model are T159 (126 km), T511 (39 km), T1279 (16 km), and T2047 (10 km). It was found that increasing horizontal resolution improves the tropical precipitation, the tropical atmospheric circulation, the frequency of occurrence of Euro-Atlantic blocking, and the representation of extratropical cyclones in large parts of the Northern Hemisphere extratropics. All of these improvements come from the increase in resolution from T159 to T511 with relatively small changes for further resolution increases to T1279 and T2047, although it should be noted that results from this very highest resolution are from a previously untested model version. Problems in simulating the Madden–Julian oscillation remain unchanged for all resolutions tested. There is some evidence that increasing horizontal resolution to T1279 leads to moderate increases in seasonal forecast skill during boreal winter in the tropics and Northern Hemisphere extratropics. Sensitivity experiments are discussed, which helps to foster a better understanding of some of the resolution dependence found for the ECMWF model in Project Athena.

Northern Hemisphere tropical cyclone (TC) activity is investigated in multiyear global climate simulations with the ECMWF Integrated Forecast System (IFS) at 10-km resolution forced by the observed records of sea surface temperature and sea ice. The results are compared to analogous simulations with the 16-, 39-, and 125-km versions of the model as well as observations. In the North Atlantic, mean TC frequency in the 10-km model is comparable to the observed frequency, whereas it is too low in the other versions. While spatial distributions of the genesis and track densities improve systematically with increasing resolution, the 10-km model displays qualitatively more realistic simulation of the track density in the western subtropical North Atlantic. In the North Pacific, the TC count tends to be too high in the west and too low in the east for all resolutions. These model errors appear to be associated with the errors in the large-scale environmental conditions that are fairly similar in this region for all model versions. The largest benefits of the 10-km simulation are the dramatically more accurate representation of the TC intensity distribution and the structure of the most intense storms. The model can generate a supertyphoon with a maximum surface wind speed of 68.4 m s−1. The life cycle of an intense TC comprises intensity fluctuations that occur in apparent connection with the variations of the eyewall/rainband structure. These findings suggest that a hydrostatic model with cumulus parameterization and of high enough resolution could be efficiently used to simulate the TC intensity response (and the associated structural changes) to future climate change.

The interactions and teleconnections between the tropical and midlatitude regions on intraseasonal time scales are an important modulator of tropical and extratropical circulation anomalies and their associated weather patterns. These interactions arise due to the impact of the tropics on the extratropics, the impact of the midlatitudes on the tropics, and two‐way interactions between the regions. Observational evidence, as well as theoretical studies with models of complexity ranging from the linear barotropic framework to intricate Earth system models, suggest the involvement of a myriad of processes and mechanisms in generating and maintaining these interconnections. At this stage, our understanding of these teleconnections is primarily a collection of concepts; a comprehensive theoretical framework has yet to be established. These intraseasonal teleconnections are increasingly recognized as an untapped source of potential subseasonal predictability. However, the complexity and diversity of mechanisms associated with these teleconnections, along with the lack of a conceptual framework to relate them, prevent this potential predictability from being translated into realized forecast skill. This review synthesizes our progress in understanding the observed characteristics of intraseasonal tropical‐extratropical interactions and their associated mechanisms, identifies the significant gaps in this understanding, and recommends new research endeavors to address the remaining challenges.

The effects of horizontal resolution and the treatment of convection on simulation of the diurnal cycle of precipitation during boreal summer are analyzed in several innovative weather and climate model integrations. The simulations include: season-long integrations of the Nonhydrostatic Icosahedral Atmospheric Model (NICAM) with explicit clouds and convection; year-long integrations of the operational Integrated Forecast System (IFS) from the European Centre for Medium-range Weather Forecasts at three resolutions (125, 39 and 16 km); seasonal simulations of the same model at 10 km resolution; and seasonal simulations of the National Center for Atmospheric Research (NCAR) low-resolution climate model with and without an embedded two-dimensional cloud-resolving model in each grid box. NICAM with explicit convection simulates best the phase of the diurnal cycle, as well as many regional features such as rainfall triggered by advancing sea breezes or high topography. However, NICAM greatly overestimates mean rainfall and the magnitude of the diurnal cycle. Introduction of an embedded cloud model within the NCAR model significantly improves global statistics of the seasonal mean and diurnal cycle of rainfall, as well as many regional features. However, errors often remain larger than for the other higher-resolution models. Increasing resolution alone has little impact on the timing of daily rainfall in IFS with parameterized convection, yet the amplitude of the diurnal cycle does improve along with the representation of mean rainfall. Variations during the day in atmospheric prognostic fields appear quite similar among models, suggesting that the distinctive treatments of model physics account for the differences in representing the diurnal cycle of precipitation.

Weather and climate variations on subseasonal to decadal time scales can have enormous social, economic, and environmental impacts, making skillful predictions on these time scales a valuable tool for decision-makers. As such, there is a growing interest in the scientific, operational, and applications communities in developing forecasts to improve our foreknowledge of extreme events. On subseasonal to seasonal (S2S) time scales, these include high-impact meteorological events such as tropical cyclones, extratropical storms, floods, droughts, and heat and cold waves. On seasonal to decadal (S2D) time scales, while the focus broadly remains similar (e.g., on precipitation, surface and upper-ocean temperatures, and their effects on the probabilities of high-impact meteorological events), understanding the roles of internal variability and externally forced variability such as anthropogenic warming in forecasts also becomes important. The S2S and S2D communities share common scientific and technical challenges. These include forecast initialization and ensemble generation; initialization shock and drift; understanding the onset of model systematic errors; bias correction, calibration, and forecast quality assessment; model resolution; atmosphere–ocean coupling; sources and expectations for predictability; and linking research, operational forecasting, and end-user needs. In September 2018 a coordinated pair of international conferences, framed by the above challenges, was organized jointly by the World Climate Research Programme (WCRP) and the World Weather Research Programme (WWRP). These conferences surveyed the state of S2S and S2D prediction, ongoing research, and future needs, providing an ideal basis for synthesizing current and emerging developments in these areas that promise to enhance future operational services. This article provides such a synthesis.

The interaction of ocean coupling and model physics in the simulation of the intraseasonal oscillation (ISO) is explored with three general circulation models: the Community Atmospheric Model, versions 3 and 4 (CAM3 and CAM4), and the superparameterized CAM3 (SPCAM3). Each is integrated coupled to an ocean model, and as an atmosphere-only model using sea surface temperatures (SSTs) from the coupled SPCAM3, which simulates a realistic ISO. For each model, the ISO is best simulated with coupling. For each SST boundary condition, the ISO is best simulated in SPCAM3.Near-surface vertical gradients of specific humidity, Dq (temperature, DT), explain ;20% (50%) of tropical Indian Ocean latent (sensible) heat flux variance, and somewhat less of west Pacific variance. In turn, local SST anomalies explain ;5% (25%) of Dq (DT) variance in coupled simulations, and less in uncoupled simulations. Ergo, latent and sensible heat fluxes are strongly controlled by wind speed fluctuations, which are largest in the coupled simulations, and represent a remote response to coupling. The moisture budget reveals that wind variability in coupled simulations increases east-of-convection midtropospheric moistening via horizontal moisture advection, which influences the direction and duration of ISO propagation.These results motivate a new conceptual model for the role of ocean feedbacks on the ISO. Indian Ocean surface fluxes help developing convection attain a magnitude capable of inducing the circulation anomalies necessary for downstream moistening and propagation. The ''processing'' of surface fluxes by model physics strongly influences the moistening details, leading to model-dependent responses to coupling.