Abstract. Single model initial-condition large ensembles (LEs) are a useful approach to understand the roles of forced responses and internal variability in historical and future climate change. Here, we produce one of the largest ensembles thus far using the MIROC6 coupled atmosphere–ocean global climate model (MIROC6-LE). The total experimental period of MIROC6-LE is longer than 76 000 years. MIROC6-LE consists of a long preindustrial control run, 50-member historical simulations, 8 single forcing historical experiments with 10 or 50 members, 5 future scenario experiments with 50 members and 3 single forcing future experiments with 50 members. Here, we describe the experimental design. The output data of most of the experiments are freely available to the public. This dataset would be useful to a wide range of research communities. We also demonstrate some examples of initial analyses. Specifically, we confirm that the linear additivity of the forcing-response relationship holds for the 1850–2020 trends of the annual mean values and extreme indices of surface air temperature and precipitation by analyzing historical fully forced runs and the sum of single forced historical runs. To isolate historical anthropogenic signals of annual mean and extreme temperature for 2000–2020 relative to 1850–1900, ensemble sizes of 4 and 15, respectively, are sufficient in most of the world. Historical anthropogenic signals of annual mean and extreme precipitation are significant with the 50-member ensembles in 76 % and 69 % of the world, respectively. Fourteen members are sufficient to examine differences in changes in annual mean values and extreme indices of temperature and precipitation between the shared socioeconomic pathways (ssp), ssp585 and ssp126, in most of the world. Ensembles larger than 50 members are desirable for investigations of differences in annual mean and extreme precipitation changes between ssp126 and ssp119. Historical and future changes in internal variability, represented by departures from the ensemble mean, are analyzed with a focus on the El Niño/Southern Oscillation (ENSO) and global annual mean temperature and precipitation. An ensemble size of 31 is large enough to detect ENSO intensification from preindustrial conditions to 1951–2000, from 1951–2000 to 2051–2100 in all future experiments, and from low- to high-emission future scenario experiments. The single forcing historical experiments with 27 members can isolate ENSO intensification due to anthropogenic greenhouse gas and aerosol forcings. Future changes in the global mean temperature variability are discernible with 23 members under all future experiments, while 50 members are not sufficient for detecting changes in the global mean precipitation variability in ssp119 and ssp126. We also confirm that these temperature and precipitation variabilities are not precisely analyzed when detrended anomalies from the long-term averages are used due to interannual climate responses to the historical natural forcing, which highlights the importance of large ensembles for assessing internal variability.