Uncertainties in historical changes and future projections of drought. Part II: model-simulated historical and future drought changes

The double Intertropical Convergence Zone (ITCZ) bias is an outstanding bias in many climate models. This work assesses the annual‐mean double‐ITCZ problem in the models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) based on several quantitative indices. Within the 46 CMIP6 models, 9 models from mainland China are evaluated as a group to verify the effort of model development from one perspective. The double‐ITCZ bias and its large intermodel spread still exist in CMIP6 models. The overall performance of the models from Chinese mainland is similar with all CMIP6 models. It is found that the top‐five models with relatively low double‐ITCZ biases can effectively restrain the frequency of deep convection and related sea surface temperature (SST) bias in the southeastern Pacific dry subsidence region, which highlights the necessity of improving convective physics in climate models. Impacts of model resolution on the double‐ITCZ problem are examined by comparing the high‐ and low‐resolution groups in CMIP6 and High Resolution Model Intercomparison Project (HighResMIP) historical experiments, respectively. Increased resolution in atmospheric models is found to be able to reduce the positive precipitation bias over the tropical southern Atlantic and improve the simulation of deep convection frequency and convective precipitation ratio there. However, the double‐ITCZ bias over the Pacific is not improved significantly by increased resolution.

The double Intertropical Convergence Zone (ITCZ) bias is an outstanding bias in many climate models. This work assesses the annual-mean double-ITCZ problem in the models participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) based on several quantitative indices. Within the forty-six CMIP6 models, nine models from mainland China are evaluated as a group to verify the effort of model development from one perspective. The double-ITCZ bias and its large inter-model spread still exist in CMIP6 models. The overall performance of the models from Chinese mainland is similar with all CMIP6 models. It is found that the top-five models with relatively low double-ITCZ biases can effectively restrain the frequency of deep convection and related sea surface temperature (SST) bias in the southeastern Pacific dry subsidence region, which highlights the necessity of improving convective physics in climate models. Impacts of model resolution on the double-ITCZ problem are examined by comparing the high- and low-resolution groups in CMIP6 and High Resolution Model Intercomparison Project (HighResMIP) historical experiments, respectively. Increased resolution in atmospheric models is found to be able to reduce the positive precipitation bias over the tropical southern Atlantic, and improve the simulation of deep convection frequency and convective precipitation ratio there. However, the double-ITCZ bias over the Pacific is not improved significantly by increased resolution.

California is one of the major uncertainty hotspots for climate change, as climate models have historically been split between projecting wetter and drier future conditions over the region. We analysed the future (mid‐century and end‐century) projections of California's winter precipitation changes from the latest Coupled Model Intercomparison Project Phase 6 (CMIP6), and studied its respective model agreement in comparison to the previous CMIP5 projections. Over northern California more than two thirds of the models in each ensemble agree on wetter future conditions. However, over southern California both ensembles show highly uncertain precipitation changes, with model projections almost equally divided between wetter or drier conditions. Projected end‐century precipitation changes range from −30% to +70% in CMIP5 and −20% to +80% in CMIP6. The CMIP6 ensemble mean changes are generally wetter and show larger model disagreement compared to CMIP5. Distribution of year‐to‐year precipitation indicates more extremely wet or dry years over southern California in CMIP6 compared to CMIP5, with some models suggesting that the five wettest years account for as much as ~55% of the 20‐year rainfall, and the five driest for as little as ~5%. Dynamically, both ensembles project weakened subsidence over Baja California that is stronger in CMIP6 than in CMIP5, in line with the wetter mean conditions in CMIP6. In the western tropical Pacific we find strengthening of the Hadley circulation in CMIP6 that is not seen in CMIP5, and more El Niño than La Niña conditions in the equatorial Pacific. More CMIP6 models also project an increase in ENSO events compared to CMIP5, and a stronger impact of ENSO on California's precipitation is found in CMIP6 than in CMIP5. These factors also contribute to larger model disagreement and more extremely wet or dry years over southern California in CMIP6.

Abstract. Anthropogenic climate change is projected to lead to ocean warming, acidification, deoxygenation, reductions in near-surface nutrients, and changes to primary production, all of which are expected to affect marine ecosystems. Here we assess projections of these drivers of environmental change over the twenty-first century from Earth system models (ESMs) participating in the Coupled Model Intercomparison Project Phase 6 (CMIP6) that were forced under the CMIP6 Shared Socioeconomic Pathways (SSPs). Projections are compared to those from the previous generation (CMIP5) forced under the Representative Concentration Pathways (RCPs). A total of 10 CMIP5 and 13 CMIP6 models are used in the two multi-model ensembles. Under the high-emission scenario SSP5-8.5, the multi-model global mean change (2080–2099 mean values relative to 1870–1899) ± the inter-model SD in sea surface temperature, surface pH, subsurface (100–600 m) oxygen concentration, euphotic (0–100 m) nitrate concentration, and depth-integrated primary production is +3.47±0.78 ∘C, -0.44±0.005, -13.27±5.28, -1.06±0.45 mmol m−3 and -2.99±9.11 %, respectively. Under the low-emission, high-mitigation scenario SSP1-2.6, the corresponding global changes are +1.42±0.32 ∘C, -0.16±0.002, -6.36±2.92, -0.52±0.23 mmol m−3, and -0.56±4.12 %. Projected exposure of the marine ecosystem to these drivers of ocean change depends largely on the extent of future emissions, consistent with previous studies. The ESMs in CMIP6 generally project greater warming, acidification, deoxygenation, and nitrate reductions but lesser primary production declines than those from CMIP5 under comparable radiative forcing. The increased projected ocean warming results from a general increase in the climate sensitivity of CMIP6 models relative to those of CMIP5. This enhanced warming increases upper-ocean stratification in CMIP6 projections, which contributes to greater reductions in upper-ocean nitrate and subsurface oxygen ventilation. The greater surface acidification in CMIP6 is primarily a consequence of the SSPs having higher associated atmospheric CO2 concentrations than their RCP analogues for the same radiative forcing. We find no consistent reduction in inter-model uncertainties, and even an increase in net primary production inter-model uncertainties in CMIP6, as compared to CMIP5.

Выполнена оценка возможных изменений среднегодовой приповерхностной температуры воздуха (ПТВ) в Дальневосточном регионе, включающем территорию и окраинные моря России, а также северо-западную часть Тихого океана, до 2099 г., для чего используются осредненные по ансамблю данные 33 моделей проекта CMIP6 (Coupled Model Intercomparison Project Phase 6), полученные в рамках четырех сценариев, отвечающих разным уровням антропогенного радиационного форсинга (слабого, умеренного и значительного). Анализируются различия между осредненными за 30-летние периоды аномалиями ПТВ. Для верификации модельных результатов проанализировано потепление, произошедшее в регионе с 1940–1969 до 1994–2023 гг., для чего использованы данные реанализа ERA5 и эксперимента Historical CMIP6. По обоим видам данных средняя ПТВ в регионе выросла на 1.1 °С: с 1940–1969 к 1994–2023 гг.; это сходство обосновывает оценки будущих изменений ПТВ по моделям CMIP6. Все сценарии SSP (Shared Socio-economic Pathways) будущего радиационного форсинга показывают приблизительно одинаковое повышение ПТВ с 1994–2023 по 2024–2053 гг., оно составляет в среднем по региону 1.2–1.5 °С. К 2070–2099 гг. средняя ПТВ в рассматриваемом регионе возрастет соответственно темпу эмиссии парниковых газов – на 1.7, 2.7, 3.8 и 4.8 °С. Как показывают данные реанализа ERA5, от 1940–1969 к 1994–2023 гг. увеличение ПТВ над морскими акваториями региона происходило весьма неравномерно: наибольшие темпы наблюдались в северной части Охотского моря (до 2 °С и более) и в прибрежных районах северо-западной части Берингова моря (до 1.0–1.2 °С). Увеличение ПТВ ослабевало в направлении с северо-запада на юго-восток, т.е. с удалением от суши, и составило 0.2–0.6 °С в северо-западной части Тихого океана. Картина потепления над морскими акваториями по данным CMIP6 выражена сильнее, чем по данным реанализа ERA5, но при этом качественно им соответствует. An assessment of possible changes in the annual mean surface air temperature (SAT) in the Far East Region (35°–65° N, 130°–180° E) is made from the present to 2099, using ensemble-averaged data from 33 CMIP6 (Coupled Model Intercomparison Project Phase 6) models obtained within the framework of four scenarios corresponding to the weak, moderate, or significant anthropogenic radiative forcing resulting from СО2 emissions. To elucidate long-term climate change, SAT averaged for 30-year periods, namely, 1994–2023, 2024–2053 and 2070–2099 are analyzed. To verify the model results, the warming that occurred in the region from the mid-20th century (1940–1969) to the early 21st century (1994–2023) is analyzed, using ERA5 data with the fine spatial resolution of 0.25°, and CMIP6 data with the coarser resolution, mostly 1.0°–2.0°. According to both data types, the regional SAT increased, on average by 1.1 °C from 1940–1969 to 1994–2023, justifying the use of forecast estimates based on the CMIP6 models in this work. All scenarios of possible radiative forcing show the similar SAT increase from the 1994–2023 to 2024–2053, on average 1.2–1.5 °C. On the contrary, by the 2070–2099, the regional SAT will increase in accordance with the emission rates on average by 1.7, 2.7, 3.8 and 4.8 °C, respectively. As for the Russian Far East land area, ERA5 and CMIP6 show similar spatial warming patterns, with the warming, on average, of 1.2 °C from 1940–1969 to 1994–2023, i.e. higher than that for the entire considered region including marine areas. From 1940–1969 to 1994–2023 negative annual mean SAT changed to positive one in some areas of the Primorsky, Khabarovsky and Kamchatksky provinces, implying the permafrost melting. According to the CMIP6 models, the land warming of 2.0–2.1 °C, 3.0–3.5 °C, 4.7–5.3 °C, and 6.1–6.6 °C is expected by the end of the 21st century for the scenarios with the different levels of radiative forcing. As shown by the ERA5 data, the SAT increase from 1940–1969 to 1994–2023 was very uneven for the marine areas: the highest rates were observed in the northern Okhotsk Sea (up to 2 °C and more) and in the coastal northwestern Bering Sea (up to 1.0–1.2 °C), which can be explained by the ice cover decrease. The SAT increase weakened in the direction from the northwest to southeast, i.e. with the distance from the land, and amounted to only 0.2–0.6 °C in the northwestern Pacific, which can be attributed to the effect of Pacific Decadal Oscillation (PDO). The coastal Okhotsk Sea off the Sakhalin Island is the only area where SAT decreased by 0.2–0.6 °C from 1940–1969 to 1994–2023, which probably can be attributed to the changes in the East Sakhalin Current transporting Amur River water southward along the coast but this suggestion should be verified. The warming pattern over the marine areas according to CMIP6 data qualitatively corresponds to that one based on ERA5 data, keeping in mind the lower resolution of the modeled data. The warming in the Northwest Pacific from the modeled data exceeds that one from ERA5, which can be explained by elimination of the PDO effects when averaging CMIP6 multi-model data.

Abstract. Poor air quality is currently responsible for large impacts on human health across the world. In addition, the air pollutants ozone (O3) and particulate matter less than 2.5 µm in diameter (PM2.5) are also radiatively active in the atmosphere and can influence Earth's climate. It is important to understand the effect of air quality and climate mitigation measures over the historical period and in different future scenarios to ascertain any impacts from air pollutants on both climate and human health. The Coupled Model Intercomparison Project Phase 6 (CMIP6) presents an opportunity to analyse the change in air pollutants simulated by the current generation of climate and Earth system models that include a representation of chemistry and aerosols (particulate matter). The shared socio-economic pathways (SSPs) used within CMIP6 encompass a wide range of trajectories in precursor emissions and climate change, allowing for an improved analysis of future changes to air pollutants. Firstly, we conduct an evaluation of the available CMIP6 models against surface observations of O3 and PM2.5. CMIP6 models consistently overestimate observed surface O3 concentrations across most regions and in most seasons by up to 16 ppb, with a large diversity in simulated values over Northern Hemisphere continental regions. Conversely, observed surface PM2.5 concentrations are consistently underestimated in CMIP6 models by up to 10 µg m−3, particularly for the Northern Hemisphere winter months, with the largest model diversity near natural emission source regions. The biases in CMIP6 models when compared to observations of O3 and PM2.5 are similar to those found in previous studies. Over the historical period (1850–2014) large increases in both surface O3 and PM2.5 are simulated by the CMIP6 models across all regions, particularly over the mid to late 20th century, when anthropogenic emissions increase markedly. Large regional historical changes are simulated for both pollutants across East and South Asia with an annual mean increase of up to 40 ppb for O3 and 12 µg m−3 for PM2.5. In future scenarios containing strong air quality and climate mitigation measures (ssp126), annual mean concentrations of air pollutants are substantially reduced across all regions by up to 15 ppb for O3 and 12 µg m−3 for PM2.5. However, for scenarios that encompass weak action on mitigating climate and reducing air pollutant emissions (ssp370), annual mean increases in both surface O3 (up 10 ppb) and PM2.5 (up to 8 µg m−3) are simulated across most regions, although, for regions like North America and Europe small reductions in PM2.5 are simulated due to the regional reduction in precursor emissions in this scenario. A comparison of simulated regional changes in both surface O3 and PM2.5 from individual CMIP6 models highlights important regional differences due to the simulated interaction of aerosols, chemistry, climate and natural emission sources within models. The projection of regional air pollutant concentrations from the latest climate and Earth system models used within CMIP6 shows that the particular future trajectory of climate and air quality mitigation measures could have important consequences for regional air quality, human health and near-term climate. Differences between individual models emphasise the importance of understanding how future Earth system feedbacks influence natural emission sources, e.g. response of biogenic emissions under climate change.

An assessment is made of the modes of interannual variability in the seasonal mean summer and winter Southern Hemisphere (SH) 500-hPa geopotential height in the twentieth century in models from the Coupled Model Intercomparison Project (CMIP) phase 5 (CMIP5) dataset. Modes of variability of both the slow (signal) and intraseasonal (noise) components in the CMIP5 models are evaluated against those estimated from reanalysis data. There is general improvement in the leading modes of the slow (signal) component in CMIP5 models compared with the CMIP phase 3 (CMIP3) dataset. The largest improvement is in the spatial structures of the modes related to El Niño–Southern Oscillation variability in SH summer. An overall score metric is significantly higher for CMIP5 over CMIP3 in both seasons. The leading modes in the intraseasonal noise component are generally well reproduced in CMIP5 models, and there are few differences from CMIP3. A new total overall score metric is used to rank the CMIP5 models over both seasons. Weighting the seasons by the relative spread of overall scores is shown to be suitable for generating multimodel ensembles for further analysis of interannual variability. In multimodel ensembles, it is found that an ensemble of size 5 or 6 is sufficient in SH summer to reproduce well the dominant modes. In contrast, about 13 models are typically are required in SH winter. It is shown that it is necessary that the selected models individually reproduce well the leading modes of the slow component.

The frequency and intensity of extreme thermal stress conditions during summer are expected to increase due to climate change. This study examines sixteen models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) that have been bias-adjusted using the quantile delta mapping method. These models provide Universal Thermal Climate Index (UTCI) for summer seasons between 1979 and 2010, which are regridded to a similar spatial grid as ERA5-HEAT (available at 0.25° × 0.25° spatial resolution) using bilinear interpolation. The evaluation compares the summertime climatology and trends of the CMIP6 multi-model ensemble (MME) mean UTCI with ERA5 data, focusing on a regional hotspot in northwest India (NWI). The Pattern Correlation Coefficient (between CMIP6 models and ERA5) values exceeding 0.9 were employed to derive the MME mean of UTCI, which was subsequently used to analyze the climatology and trends of UTCI in the CMIP6 models.The spatial climatological mean of CMIP6 MME UTCI demonstrates significant thermal stress over the NWI region, similar to ERA5. Both ERA5 and CMIP6 MME UTCI show a rising trend in thermal stress conditions over NWI. The temporal variation analysis reveals that NWI experiences higher thermal stress during the summer compared to the rest of India. The number of thermal stress days is also increasing in NWI and major Indian cities according to ERA5 and CMIP6 MME. Future climate projections under different scenarios (SSP1-2.6, SSP2-4.5, and SSP5-8.5) indicate an increasing trend in thermal discomfort conditions throughout the twenty-first century. The projected rates of increase are approximately 0.09 °C per decade, 0.26 °C per decade, and 0.56 °C per decade, respectively. Assessing the near (2022–2059) and far (2060–2100) future, all three scenarios suggest a rise in intense heat stress days (UTCI > 38 °C) in NWI. Notably, the CMIP6 models predict that NWI could reach deadly levels of heat stress under the high-emission (SSP5-8.5) scenario. The findings underscore the urgency of addressing climate change and its potential impacts on human well-being and socio-economic sectors.

As the major renewable energy, wind can greatly reduce carbon emissions. Following the “carbon neutral” strategy, wind power could help to achieve the realization of energy transformation and green development. Based on ERA5 reanalysis data and the multi-ensemble historical and scenario simulations of the Coupled Model Intercomparison Project Phase 6 (CMIP6), a variety of statistical analyses are used to evaluate the performance of CMIP6 simulating the wind speed in China. The conclusions are as follows: spatial patterns of the nine CMIP6 models are similar with ERA5, but BCC-CSM2-MR and MRI-ESM2-0 highly overestimate the wind speed in northwest China. CESM2-WACCM, NorESM2-MM, and HadGEM3-GC31-MM behave better than the other six CMIP6 models in four specific regions are chosen for detailed study. CESM2-WACCM, NorESM2-MM, and HadGEM3-GC31-MM tend to simulate a larger wind speed than ERA5 except the yearly averaged wind speed in region II and region IV. CESM2-WACCM and NorESM2-MM simulate a large monthly mean wind speed, but the value is relatively close with ERA5 in the summer. HadGEM3-GC31-MM overestimates wind speed in region I and region II from April to October, but gets closer with ERA during winter. CESM2-WACCM, NorESM2-MM, and HadGEM3-GC31-MM simulate an increasing trend in Tibetan Plateau and Xinjiang in the next 100 years, while NorESM2-MM projects rising wind speed in the eastern part of Inner Mongolia, and HadGEM3-GC31-MM simulates increasing wind speed in the northeast and central China. The future wind speed in three models is projected to decline in region I, and the value of HadGEM3-GC31-MM is much larger. In region II, wind speed simulated by three models is projected to decrease, but the wind speed from HadGEM3-GC31-MM in region III and modeled wind speed in region IV from NorESM2-MM would climb with the slope equal to 0.0001 and 0.0012, respectively. This study indicates that the CMIP6 models have certain limitations to perform realistic wind changes, but CMIP6 could provide available reference for the projection of wind in specific areas.

This is the first study to show the global Cut-off Low (COL) activity in 46 models participating in the Coupled Model Intercomparison Project Phase 5 (CMIP5) and Phase 6 (CMIP6). The COL historical simulations for the period 1979–2005 obtained from the CMIP5 and CMIP6 models and their ensembles are compared with the ERA5 reanalysis using an objective feature-tracking algorithm. The results show that the CMIP6 models simulate the spatial distribution of COLs more realistically than the CMIP5 models. Some improvements include reduced equatorward bias and underestimation over regions of high COL density. Reduced biases in CMIP6 are mainly attributed to the improved representation of the zonal wind due to the poleward shift of the subtropical jet streams. The CMIP5 models systematically underestimate the COL intensity as measured by the T42 vorticity at 250 hPa. In CMIP6, the intensity is still underestimated in summer, but overestimated in winter in part due to increased westerlies. The overestimation is enhanced by the finer spatial resolution models that identify more of the strong systems compared to coarser resolution models. Other aspects of COLs such as their temporal and lifetime distributions are modestly improved in CMIP6 compared to CMIP5. Finally, the predictive skill of climate models is evaluated using five variables and the Taylor diagram. We find that 15 out of the 20 (75%) best coupled models belong to CMIP6, and this highlights the overall improvement compared to its predecessor CMIP5. Despite this, the use of the multi-model ensemble average seems to be better in simulating COLs than individual models.

Ensembles of climate change projections created by general circulation models (GCMs) with high resolution are increasingly needed to develop adaptation strategies for regional climate change. The Meteorological Research Institute atmospheric GCM version 3.2 (MRI-AGCM3.2), which is listed in the Coupled Model Intercomparison Project phase 5 (CMIP5), has been typically run with resolutions of 60 km and 20 km. Ensembles of MRI-AGCM3.2 consist of members with multiple cumulus convection schemes and different patterns of future sea surface temperature, and are utilized together with their downscaled data; however, the limited size of the high-resolution ensemble may lead to undesirable biases and uncertainty in future climate projections that will limit its appropriateness and effectiveness for studies on climate change and impact assessments. In this study, to develop a comprehensive understanding of the regional precipitation simulated with MRI-AGCM3.2, we investigate how well MRI-AGCM3.2 simulates the present-day regional precipitation around the globe and compare the uncertainty in future precipitation changes and the change projection itself between MRI-AGCM3.2 and the CMIP5 multiple atmosphere–ocean coupled GCM (AOGCM) ensemble. MRI-AGCM3.2 reduces the bias of the regional mean precipitation obtained with the high-performing CMIP5 models, with a reduction of approximately 20% in the bias over the Tibetan Plateau through East Asia and Australia. When 26 global land regions are considered, MRI-AGCM3.2 simulates the spatial pattern and the regional mean realistically in more regions than the individual CMIP5 models. As for the future projections, in 20 of the 26 regions, the sign of annual precipitation change is identical between the 50th percentiles of the MRI-AGCM3.2 ensemble and the CMIP5 multi-model ensemble. In the other six regions around the tropical South Pacific, the differences in modeling with and without atmosphere–ocean coupling may affect the projections. The uncertainty in future changes in annual precipitation from MRI-AGCM3.2 partially overlaps the maximum–minimum uncertainty range from the full ensemble of the CMIP5 models in all regions. Moreover, on average over individual regions, the projections from MRI-AGCM3.2 spread over roughly 0.8 of the uncertainty range from the high-performing CMIP5 models compared to 0.4 of the range of the full ensemble.

Abstract. By regulating the global transport of heat, freshwater, and carbon, the Atlantic meridional overturning circulation (AMOC) serves as an important component of the climate system. During the late 20th and early 21st centuries, indirect observations and models suggest a weakening of the AMOC. Direct AMOC observations also suggest a weakening during the early 21st century but with substantial interannual variability. Long-term weakening of the AMOC has been associated with increasing greenhouse gases (GHGs), but some modeling studies suggest the build up of anthropogenic aerosols (AAs) may have offset part of the GHG-induced weakening. Here, we quantify 1900–2020 AMOC variations and assess the driving mechanisms in state-of-the-art climate models from the Coupled Model Intercomparison Project phase 6 (CMIP6). The CMIP6 forcing (GHGs, anthropogenic and volcanic aerosols, solar variability, and land use and land change) multi-model mean shows negligible AMOC changes up to ∼ 1950, followed by robust AMOC strengthening during the second half of the 20th century (∼ 1950–1990) and weakening afterwards (1990–2020). These multi-decadal AMOC variations are related to changes in North Atlantic atmospheric circulation, including an altered sea level pressure gradient, storm track activity, surface winds, and heat fluxes, which drive changes in the subpolar North Atlantic surface density flux. To further investigate these AMOC relationships, we perform a regression analysis and decompose these North Atlantic climate responses into an anthropogenic aerosol-forced component and a subsequent AMOC-related feedback. Similar to previous studies, CMIP6 GHG simulations yield robust AMOC weakening, particularly during the second half of the 20th century. Changes in natural forcings, including solar variability and volcanic aerosols, yield negligible AMOC changes. In contrast, CMIP6 AA simulations yield robust AMOC strengthening (weakening) in response to increasing (decreasing) anthropogenic aerosols. Moreover, the CMIP6 all-forcing AMOC variations and atmospheric circulation responses also occur in the CMIP6 AA simulations, which suggests these are largely driven by changes in anthropogenic aerosol emissions. More specifically, our results suggest that AMOC multi-decadal variability is initiated by North Atlantic aerosol optical thickness perturbations to net surface shortwave radiation and sea surface temperature (and hence sea surface density), which in turn affect sea level pressure gradient and surface wind and – via latent and sensible heat fluxes – sea surface density flux through its thermal component. AMOC-related feedbacks act to reinforce this aerosol-forced AMOC response, largely due to changes in sea surface salinity (and hence sea surface density), with temperature-related (and cloud-related) feedbacks acting to mute the initial response. Although aspects of the CMIP6 all-forcing multi-model mean response resembles observations, notable differences exist. This includes CMIP6 AMOC strengthening from ∼ 1950 to 1990, when the indirect estimates suggest AMOC weakening. The CMIP6 multi-model mean also underestimates the observed increase in North Atlantic ocean heat content, and although the CMIP6 North Atlantic atmospheric circulation responses – particularly the overall patterns – are similar to observations, the simulated responses are weaker than those observed, implying they are only partially externally forced. The possible causes of these differences include internal climate variability, observational uncertainties, and model shortcomings, including excessive aerosol forcing. A handful of CMIP6 realizations yield AMOC evolution since 1900 similar to the indirect observations, implying the inferred AMOC weakening from 1950 to 1990 (and even from 1930 to 1990) may have a significant contribution from internal (i.e., unforced) climate variability. Nonetheless, CMIP6 models yield robust, externally forced AMOC changes, the bulk of which are due to anthropogenic aerosols.

Based on four reanalyses or gridded data sets (ERA5, 20CR, APHRODITE and REGEN), we provide an overview of 23 Historical and 7 HighResMIP experiments’ performance from the Coupled Model Intercomparison Project Phase 6 (CMIP6) (for short, 6-Hist, HighRes) in simulating seven extreme precipitation indices over Asia defined by the Expert Team on Climate Change Detection and Indices (ETCCDI). We compare them with 28 Historical experiments in CMIP5 (5-Hist). CMIP5 and CMIP6 models are generally able to reproduce extreme precipitation’s spatial distribution and their trend patterns in comparison to the benchmark data set (APHRODITE). The overall performance of individual model is summarized by a “portrait” diagram based on four statistics for each index. We divide all 58 models into three groups (A, the top 20%; B, the median 60% and C group, the last 20%) according to MR rankings (the comprehensive ranking measure). Based on the “portrait” diagram and MR rankings, models that perform relatively well for all seven extreme precipitation indices include HadCM3, HadGEM2-AO, HadGEM2-CC and HadGEM2-ES from 5-Hist, EC-Earth3, EC-Earth3-Veg from 6-Hist and ECMWF-IFS-HR, ECMWF-IFS-LR, ECMWF-IFS-MR from HighRes. The simulated performance of CMIP6 is polarized, for the top four and the last five ranking models are both from CMIP6. Compared with the counterpart models in CMIP6 and CMIP5, the improvement of PCC (pattern correlation coefficient) is more obvious. Furthermore, the dry biases of CMIP6 (both 6-Hist and HighRes) in Southern China and India and the wet biases of CMIP6 in Tibet are reduced compared to CMIP5. This may benefit from the improvement that CMIP6 models can capture the characteristics of meridional moisture flux convergence, and improve the overestimation or underestimation of meridional and zonal specific humidity eddies compared to CMIP5.

Cold winters in Eurasia considerably affect transportation, agriculture, energy, and public health. This study utilizes 31 global climate models from the Coupled Model Intercomparison Project Phase 6 (CMIP6) and 33 CMIP5 models to evaluate the historical surface air temperature, sea level pressure, 500‐hPa geopotential height, 150‐hPa meridional and zonal wind, and polar vortex indices during cold winters. Our research quantifies the advancements of CMIP6 over CMIP5. Additionally, future changes in these variables under three different Shared Socioeconomic Pathways (SSPs), that is, SSP 1–2.6, SSP 2–4.5, and SSP 5–8.5, are projected based on 20 out of the 31 CMIP6 models. The results indicate that the multimodel ensemble means from both CMIP5 and CMIP6 effectively capture the main features of the observed Eurasian cold winters and their associated factors with good simulation agreement. The CMIP6 ensemble mean outperforms its CMIP5 counterpart, and both ensemble means (CMIP5 and CMIP6) perform better than individual CMIP6 models. Among CMIP6 models, 500‐hPa geopotential height achieves the highest simulation skill, whereas sea level pressure shows the lowest. Compared with same‐institute models from CMIP5, CMIP6 models show overall improvements with sea level pressure simulation being notably advanced. Under the three SSPs, the occurrence probability of cold winters is projected to decrease as the area and intensity indices of the polar vortex decline. Moreover, surface temperature anomalies are projected to exhibit a “warm Arctic and cold Eurasia” pattern, and the anticyclonic anomalies at 500 hPa and 150 hPa are projected to be centered at high latitudes.

Precipitation variability has great economic, social, and environmental impacts across the globe, and in particular in China. This paper evaluates the historical precipitation variability based on 20 general circulation models from the Coupled Model Intercomparison Project Phase 5 (CMIP5) archive over the 20th century relative to two observational data sets and quantifies CMIP5 improvements over CMIP3. Multimodel ensemble means and individual models are assessed. Three future emission scenarios are used (representative concentration pathways (RCP) 8.5, RCP 4.5, and RCP 2.6), and 21st century CMIP5 estimates are put into context based on the 20th century biases. We find that CMIP5 models can reproduce the spatial pattern of precipitation over China during the 20th century, which represents an improvement over CMIP3. However, the models overestimate the magnitude of seasonal and annual precipitation in most regions of China, especially along the eastern edge of the Tibetan Plateau, and underestimate summer precipitation over southeastern China. For China as a whole, CMIP5's overestimation of annual precipitation is greater than CMIP3, which can be traced back to a greater underestimation of summer precipitation in CMIP3. There is a large spread among individual models, with the greatest uncertainties in simulating summer precipitation. Trends and correlations also suggest a better agreement of CMIP5 with observations than CMIP3. Throughout the 20th century, both the observations and models show an increasing trend in precipitation over parts of northwestern China and a decreasing trend over the Tibetan Plateau. There is poor agreement in precipitation trends over the southeast and northeast regions. In general, multimodel means cannot capture the amplitude of observed multidecadal precipitation variability. In the 21st century, precipitation is generally projected to increase across all of China under all three scenarios. RCP 8.5 exhibits the largest significant trend at a rate of +1.5 mm/yr, corresponding to 16% precipitation increase by the end of the century. The RCP 2.6 scenario shows the smallest increases, at +0.5 mm/yr (6%) by 2100. The greatest increases are projected to occur over the Tibetan Plateau and eastern China in summer, suggesting an altered monsoonal circulation in the future. However, due to the uncertainties in CMIP5, future precipitation projections should be interpreted with caution.