<p>Following the protocol of the fourth phase of the Paleoclimate Modeling Intercomparison Project (PMIP4), we performed numerical experiments targeting distinctive past time periods using the Model for Interdisciplinary Research on Climate Earth System Model (MIROC-ES2L), which is an Earth System Model. Setup and basic performance of the experiments are presented.</p><p> The Last Glacial Maximum was one of the most extreme climate states during the Quaternary and conducting numerical modeling experiments of this period has long been a challenge for the paleoclimate community. We conducted a Last Glacial Maximum experiment with a long spin-up of nearly 9,000 years. Globally, there is reasonable agreement between the anomalies relative to present day derived from model climatology and those derived from proxy data archives while some regional discrepancies remain.</p><p> By changing orbital and greenhouse gas forcings, we conducted experiments for two interglacial periods: 6,000 and 127,000 years before present. Model anomalies relative to present day are qualitatively consistent with variations in solar forcing. However, anomalies in the model are smaller than those derived from proxy data archives, suggesting that processes that play a role in past interglacial climates are still missing in this state-of-the-art model. We conducted transient simulations from 850 CE to 1850 CE and from 1850 CE to 2014 CE. Cooling in the model indicates clear responses to huge volcanic eruptions, which are consistent with paleo-proxy data. The contrast between cooling during the Little Ice Age and warming during the 20th to 21st centuries is well represented at the multi-decadal time scale.</p>