The sixth version of the Model for Interdisciplinary Research on Climate (MIROC), called MIROC6, was cooperatively developed by a Japanese modeling community. In the present paper, simulated mean climate, internal climate variability, and climate sensitivity in MIROC6 are evaluated and briefly summarized in comparison with the previous version of our climate model (MIROC5) and observations. The results show that the overall reproducibility of mean climate and internal climate variability in MIROC6 is better than that in MIROC5. The tropical climate systems (e.g., summertime precipitation in the western Pacific and the eastward-propagating Madden–Julian oscillation) and the midlatitude atmospheric circulation (e.g., the westerlies, the polar night jet, and troposphere–stratosphere interactions) are significantly improved in MIROC6. These improvements can be attributed to the newly implemented parameterization for shallow convective processes and to the inclusion of the stratosphere. While there are significant differences in climates and variabilities between the two models, the effective climate sensitivity of 2.6 K remains the same because the differences in radiative forcing and climate feedback tend to offset each other. With an aim towards contributing to the sixth phase of the Coupled Model Intercomparison Project, designated simulations tackling a wide range of climate science issues, as well as seasonal to decadal climate predictions and future climate projections, are currently ongoing using MIROC6.

Decadal-scale climate variations over the Pacific Ocean and its surroundings are strongly related to the so-called Pacific decadal oscillation (PDO) which is coherent with wintertime climate over North America and Asian monsoon, and have important impacts on marine ecosystems and fisheries. In a near-term climate prediction covering the period up to 2030, we require knowledge of the future state of internal variations in the climate system such as the PDO as well as the global warming signal. We perform sets of ensemble hindcast and forecast experiments using a coupled atmosphere-ocean climate model to examine the predictability of internal variations on decadal timescales, in addition to the response to external forcing due to changes in concentrations of greenhouse gases and aerosols, volcanic activity, and solar cycle variations. Our results highlight that an initialization of the upper-ocean state using historical observations is effective for successful hindcasts of the PDO and has a great impact on future predictions. Ensemble hindcasts for the 20th century demonstrate a predictive skill in the upper-ocean temperature over almost a decade, particularly around the Kuroshio-Oyashio extension (KOE) and subtropical oceanic frontal regions where the PDO signals are observed strongest. A negative tendency of the predicted PDO phase in the coming decade will enhance the rising trend in surface air-temperature (SAT) over east Asia and over the KOE region, and suppress it along the west coasts of North and South America and over the equatorial Pacific. This suppression will contribute to a slowing down of the global-mean SAT rise.

The present paper presents results of seasonal‐to‐decadal climate predictions based on a coupled climate model called the Model for Interdisciplinary Research on Climate version 6 (MIROC6) contributing to the Coupled Model Intercomparison Project Phase 6 (CMIP6). MIROC6 is initialized every year for 1960–2018 by assimilating observed ocean temperature and salinity anomalies and full fields of sea ice concentration and by prescribing atmospheric initial states from reanalysis data. The impacts of updating the system on prediction skill are then evaluated by comparing hindcast experiments between the MIROC6 prediction system and a previous system based on MIROC version 5 (MIROC5). Skill of seasonal prediction is overall improved in association with representation and initialization of El Niño/Southern Oscillation (ENSO), the Quasi‐Biennial Oscillation (QBO), and the Northern Hemisphere sea ice concentration in MIROC6. In particular, the QBO is skillfully predicted up to 3 years ahead with a maximum anomaly correlation exceeding r = 0.8. The prediction skill for the North Atlantic Oscillation in winter is also enhanced, but the prediction still suffers from model's inherent errors. On decadal timescales, MIROC6 has a larger fraction of areas of the globe with better surface temperature skill at all lead times than MIROC5, and it has predictive skill in the annual‐mean sea surface temperature (SST) in the North Atlantic and the Pacific. In particular, MIROC6 hindcasts at 2–5 years lead time are able to capture the spatial structure of SST changes in the North Pacific and the eastern tropical Pacific associated with the 1970s regime shift better than MIROC5 hindcasts.

The El Niño–Southern Oscillation (ENSO) exerts a strong influence on tropical Atlantic variability, but it is also affected by Atlantic forcing. Previous research has proposed three Atlantic precursors for ENSO: the North tropical Atlantic, the equatorial Atlantic, and the entire tropical Atlantic. However, the relative importance of these Atlantic precursors for ENSO remains unclear. Here, we present evidence from a set of multimodel partial ocean assimilation experiments that equatorial Atlantic cooling is the main contributor for weakening equatorial zonal winds in the Indo‐Pacific sector and subsequent ocean warming in the tropical Pacific. Opposite tendencies occur for a warmer equatorial Atlantic. The equatorial Atlantic affects the interbasin climate seesaw between the Atlantic and Pacific through an atmospheric zonal Wavenumber 1 pattern. However, model mean state biases and systematic errors prevent a precise assessment of the response times for the equatorial Pacific trade winds to Atlantic forcing.

The Paris Agreement calls for efforts to limit anthropogenic global warming to less than 1.5^oC above pre‐industrial levels. However, natural internal variability may exacerbate anthropogenic warming to produce temporary excursions above 1.5^oC. Such excursions would not necessarily exceed the Paris Agreement, but would provide a warning that the threshold is being approached. Here we develop a new capability to predict the probability that global temperature will exceed 1.5^oC above pre‐industrial levels in the coming five‐years. For the period 2017 to 2021 we predict a 38% and 10% chance respectively of monthly or yearly temperatures exceeding 1.5^oC, with virtually no chance of the five‐year mean being above the threshold. Our forecasts will be updated annually to provide policy makers with advanced warning of the evolving probability and duration of future warming events.

We demonstrate the significant impact of the Indian Ocean on the Pacific climate on decadal timescales by comparing two sets of data assimilation experiments (pacemaker experiments) conducted over recent decades. For the Indian Ocean of an atmosphere‐ocean coupled global climate model, we assimilate ocean temperature and salinity anomalies defined as deviations from climatology or as anomalies with the area‐averaged changes for the Indian Ocean subtracted. When decadal sea surface temperature (SST) trends are observed to be strong over the Indian Ocean, the equatorial thermocline uniformly deepens, and the model simulates the eastward tendencies of surface wind aloft. Surface winds strongly converge around the maritime continent, and the associated strengthening of the Walker circulation suppresses an increasing trend in the equatorial Pacific SST through ocean thermocline shoaling, similar to common changes associated with seasonal Indian Ocean warming.