Abstract:The Hadley circulation (HC) is the strongest and largest meridional circulation, playing a crucial role in influencing the global climate. A regime shift is observed in the boreal winter HC variability after 1980, presenting as an equatorially symmetrical pattern. This study examines the performances of 26 CMIP6 models in simulating the variability of boreal winter HC after 1980, analyzes the possible physical processes of the simulation uncertainties, and predicts its future variations. The results reveal that all models can reasonably capture the spatial structure of the climatological boreal winter HC. However, significant inter–model discrepancies are existed in the simulations of the spatial distribution of its linear trend, the dominant mode of HC (EOF1), as well as El Niño–Southern Oscillation (ENSO). By examining the spatial structure of HC EOF1 and ENSO between the simulations and observations, all the models are categorized into three groups, i.e., both the HC EOF1 and ENSO are well–simulated group (WWG), the simulation of HC EOF1 is unsatisfactory and ENSO is well–simulated group (UWG), and neither the simulations of HC EOF1 and ENSO is satisfactory (UUG). The group analysis suggests that the WWG models can well simulate the robust relationship between the HC EOF1 and ENSO, while the UWG and UUG models fail in reproducing the connection between the EOF1 and ENSO. In the WWG models, a greater variability of ENSO would associate with a stronger linkage between the HC EOF1 and ENSO, implying that the amplitude of ENSO variability determines the linkage between HC and ENSO. Conversely, the opposite situation is observed in the UWG and UUG models. The failure of UWG and UUG models in reproducing the relationship of HC EOF–ENSO is mainly due to its unsatisfactory simulation in the ENSO variability. This result shows that a model's ability in simulating the ENSO variability affects its capacity in simulating the relationship between the HC EOF1 and ENSO, which in turn impacting its performance in charactering the variability of HC. Furthermore, the HC EOF1 is projected to maintain an equatorially symmetrical structure towards the end of this century. Compared to the period 1980–2013, the variability of ENSO by century–end has intensified, leading to a stronger association between ENSO and HC EOF1 and a strengthened HC EOF1. These findings further support of the conclusion. Our work focuses on the discrepancies of the simulation of the variability in boreal winter HC and its future prediction in CMIP6, providing valuable insights for the advancement of climate models and understanding the relationship between ENSO and HC.

Eurasian snow cover variability and links to winter climate in the CMIP5 models

Snow Cover on the Tibetan Plateau and Lake Baikal Intensifies the Winter North Atlantic Oscillation

Changes in the Relationship Between the Interannual Variation of Eurasian Snow Cover and Spring SAT over Eastern Eurasia

Modeled Climate Responses to Realistic Extremes of Northern Hemisphere Spring and Summer Snow Anomalies

Modeled Northern Hemisphere Autumn and Winter Climate Responses to Realistic Tibetan Plateau and Mongolia Snow Anomalies

Observed Influences of Autumn-Early Winter Eurasian Snow Cover Anomalies on the Hemispheric Pna-Like Variability in Winter

The Warm Arctic—Cold Eurasia Pattern and Its Key Region in Winter in CMIP6 Model Simulations

Diversity of the wintertime Arctic Oscillation pattern among CMIP5 models: Role of the stratospheric polar vortex

An NAO-dominated mode of atmospheric circulation drives large decadal changes in wintertime surface climate and snow mass over Eurasia

Evaluation of snow extent and its variability in the Atmospheric Model Intercomparison Project

The cooperative effects of November Arctic sea ice and Eurasian snow cover on the Eurasian surface air temperature in January–February

Impacts of Autumn‐Winter Tibetan Plateau Snow Anomalies on North Atlantic‐Europe and Arctic Climate

What Determines the Leading Mode of Wintertime Northern Hemisphere Snow-Cover Interannual Variability?

Model Sensitivity of Tibetan Plateau Surface Potential Vorticity and the Asian Summer Monsoon System to Asian Orographic Perturbation in FGOALS-f2

Autumn Snow Cover Variability over Northern Eurasia and Roles of Atmospheric Circulation

Atmospheric response to cold wintertime Tibetan Plateau conditions over eastern Asia in climate models

Upper‐level wind mode over the North Atlantic‐Arctic Ocean leading to coherent variations of Eurasian and North American winter temperatures and weakness in the CMIP6 models

Spatio-temporal variability of winter snow cover over the Tibetan Plateau and its relation to Arctic Oscillation

The Summer North Atlantic Oscillation, Arctic sea ice, and Arctic jet Rossby wave forcing

Role of Eurasian Snow Cover in Linking Winter‐Spring Eurasian Coldness to the Autumn Arctic Sea Ice Retreat

Can CMIP6 models reproduce the influence of ENSO on the variability of boreal winter Hadley Circulation?