The Max Planck Institute Grand Ensemble (MPI‐GE) is the largest ensemble of a single comprehensive climate model currently available, with 100 members for the historical simulations (1850–2005) and four forcing scenarios. It is currently the only large ensemble available that includes scenario representative concentration pathway (RCP) 2.6 and a 1% CO₂ scenario. These advantages make MPI‐GE a powerful tool. We present an overview of MPI‐GE, its components, and detail the experiments completed. We demonstrate how to separate the forced response from internal variability in a large ensemble. This separation allows the quantification of both the forced signal under climate change and the internal variability to unprecedented precision. We then demonstrate multiple ways to evaluate MPI‐GE and put observations in the context of a large ensemble, including a novel approach for comparing model internal variability with estimated observed variability. Finally, we present four novel analyses, which can only be completed using a large ensemble. First, we address whether temperature and precipitation have a pathway dependence using the forcing scenarios. Second, the forced signal of the highly noisy atmospheric circulation is computed, and different drivers are identified to be important for the North Pacific and North Atlantic regions. Third, we use the ensemble dimension to investigate the time dependency of Atlantic Meridional Overturning Circulation variability changes under global warming. Last, sea level pressure is used as an example to demonstrate how MPI‐GE can be utilized to estimate the ensemble size needed for a given scientific problem and provide insights for future ensemble projects.

Single model initial-condition large ensembles (SMILEs) are valuable tools that can be used to investigate the climate system. SMILEs allow scientists to quantify and separate the internal variability of the climate system and its response to external forcing, with different types of SMILEs appropriate to answer different scientific questions. In this editorial we first provide an introduction to SMILEs and an overview of the studies in the special issue “Large Ensemble Climate Model Simulations: Exploring Natural Variability, Change Signals and Impacts”. These studies analyse a range of different types of SMILEs including global climate models (GCMs), regionally downscaled climate models (RCMs), a hydrological model with input from a RCM SMILE, a SMILE with prescribed sea surface temperature (SST) built for event attribution, a SMILE that assimilates observed data, and an initialised regional model. These studies provide novel methods, that can be used with SMILEs. The methods published in this issue include a snapshot empirical orthogonal function analysis used to investigate El Niño–Southern Oscillation teleconnections; the partitioning of future uncertainty into model differences, internal variability, and scenario choices; a weighting scheme for multi-model ensembles that can incorporate SMILEs; and a method to identify the required ensemble size for any given problem. Studies in this special issue also focus on RCM SMILEs, with projections of the North Atlantic Oscillation and its regional impacts assessed over Europe, and an RCM SMILE intercomparison. Finally a subset of studies investigate projected impacts of global warming, with increased water flows projected for future hydrometeorological events in southern Ontario; precipitation projections over central Europe are investigated and found to be inconsistent across models in the Alps, with a continuation of past tendencies in Mid-Europe; and equatorial Asia is found to have an increase in the probability of large fire and drought events under higher levels of warming. These studies demonstrate the utility of different types of SMILEs. In the second part of this editorial we provide a perspective on how three types of SMILEs could be combined to exploit the advantages of each. To do so we use a GCM SMILE and an RCM SMILE with all forcings, as well as a naturally forced GCM SMILE (nat-GCM) over the European domain. We utilise one of the key advantages of SMILEs, precisely separating the forced response and internal variability within an individual model to investigate a variety of simple questions. Broadly we show that the GCM can be used to investigate broad-scale patterns and can be directly compared to the nat-GCM to attribute forced changes to either anthropogenic emissions or volcanoes. The RCM provides high-resolution spatial information of both the forced change and the internal variability around this change at different warming levels. By combining all three ensembles we can gain information that would not be available using a single type of SMILE alone, providing a perspective on future research that could be undertaken using these tools.

Initial-condition large ensembles with ensemble sizes ranging from 30 to 100 members have become a commonly used tool for quantifying the forced response and internal variability in various components of the climate system. However, there is no consensus on the ideal or even sufficient ensemble size for a large ensemble. Here, we introduce an objective method to estimate the required ensemble size that can be applied to any given application and demonstrate its use on the examples of global mean near-surface air temperature, local temperature and precipitation, and variability in the El Niño–Southern Oscillation (ENSO) region and central United States for the Max Planck Institute Grand Ensemble (MPI-GE). Estimating the required ensemble size is relevant not only for designing or choosing a large ensemble but also for designing targeted sensitivity experiments with a model. Where possible, we base our estimate of the required ensemble size on the pre-industrial control simulation, which is available for every model. We show that more ensemble members are needed to quantify variability than the forced response, with the largest ensemble sizes needed to detect changes in internal variability itself. Finally, we highlight that the required ensemble size depends on both the acceptable error to the user and the studied quantity.

Two large ensembles are used to quantify the extent to which internal variability can contribute to long‐term changes in El Niño‐Southern Oscillation (ENSO) characteristics. We diagnose changes that are externally forced and distinguish between multi‐model simulation results that differ by chance and those that differ due to different model physics. The range of simulated ENSO amplitude changes in the large ensemble historical simulations encompasses 90% of the Coupled Model Intercomparison Project 5 historical simulations and 80% of moderate (RCP4.5) and strong (RCP8.5) warming scenarios. When considering projected ENSO pattern changes, model differences are also important. We find that ENSO has high internal variability and that single realizations of a model can produce very different results to the ensemble mean response. Due to this variability, 30–40 ensemble members of a single model are needed to robustly compute absolute ENSO variance to a 10% error when 30‐year analysis periods are used.

Single‐model initial‐condition large ensembles are powerful tools to quantify the forced response, internal climate variability, and their evolution under global warming. Here, we present the CMIP6 version of the Max Planck Institute Grand Ensemble (MPI‐GE CMIP6) with currently 30 realizations for the historical period and five emission scenarios. The power of MPI‐GE CMIP6 goes beyond its predecessor ensemble MPI‐GE by providing high‐frequency output, the full range of emission scenarios including the highly policy‐relevant low emission scenarios SSP1‐1.9 and SSP1‐2.6, and the opportunity to compare the ensemble to complementary high‐resolution simulations. First, we describe MPI‐GE CMIP6, evaluate it with observations and reanalyzes and compare it to MPI‐GE. Then, we demonstrate with six application examples how to use the power of the ensemble to better quantify and understand present and future climate extremes, to inform about uncertainty in approaching Paris Agreement global warming limits, and to combine large ensembles and artificial intelligence. For instance, MPI‐GE CMIP6 allows us to show that the recently observed Siberian and Pacific North American heatwaves would only avoid reaching 1–2 years return periods in 2071–2100 with low emission scenarios, that recently observed European precipitation extremes are captured only by complementary high‐resolution simulations, and that 3‐hourly output projects a decreasing activity of storms in mid‐latitude oceans. Further, the ensemble is ideal for estimates of probabilities of crossing global warming limits and the irreducible uncertainty introduced by internal variability, and is sufficiently large to be used for infilling surface temperature observations with artificial intelligence.

Future changes in the El Niño–Southern Oscillation (ENSO) are uncertain, both because future projections differ between climate models and because the large internal variability of ENSO clouds the diagnosis of forced changes in observations and individual climate model simulations. By leveraging 14 single model initial-condition large ensembles (SMILEs), we robustly isolate the time-evolving response of ENSO sea surface temperature (SST) variability to anthropogenic forcing from internal variability in each SMILE. We find nonlinear changes in time in many models and considerable inter-model differences in projected changes in ENSO and the mean-state tropical Pacific zonal SST gradient. We demonstrate a linear relationship between the change in ENSO SST variability and the tropical Pacific zonal SST gradient, although forced changes in the tropical Pacific SST gradient often occur later in the 21st century than changes in ENSO SST variability, which can lead to departures from the linear relationship. Single-forcing SMILEs show a potential contribution of anthropogenic forcing (aerosols and greenhouse gases) to historical changes in ENSO SST variability, while the observed historical strengthening of the tropical Pacific SST gradient sits on the edge of the model spread for those models for which single-forcing SMILEs are available. Our results highlight the value of SMILEs for investigating time-dependent forced responses and inter-model differences in ENSO projections. The nonlinear changes in ENSO SST variability found in many models demonstrate the importance of characterizing this time-dependent behavior, as it implies that ENSO impacts may vary dramatically throughout the 21st century.