Large ensemble coupled earth system model simulations are deemed essential to robustly quantifying spatial and temporal structures of climate variabilities and distinguishing forced responses from natural variabilities. E3SM Version 3 coupled model at nominal resolution has shown substantially improved skills in simulating atmospheric variabilities from diurnal to intra-seasonal scales and surface air temperature record over the historical era. These improvements are attributable to improved treatments of physics in individual component models as well as some targeted grid resolution increase, combinedly an effective pathway known to be able to lead to improved representation of model climatology as well as major modes of variabilities in the coupled system. The E3SMv3 based historical large ensemble simulations are being carried out using the so-called ‘macro’ initialization approach. In this work, we will present an overview of the large ensemble in characterizing the simulated climate variabilities with a focus on the major modes of variabilities and the associated teleconnections and the quantification of the forced responses during the historical period.