AbstractEl Niño–Southern Oscillation (ENSO) is the dominant mode of interannual climate fluctuations with wide-ranging socioeconomic and environmental impacts. Understanding the eastern Pacific (EP) and central Pacific (CP) El Niño response to a warmer climate is paramount, yet the role of internal climate variability in modulating their response is not clear. Using large ensembles, we find that internal variability generates a spread in the standard deviation and skewness of these two El Niño types that is similar to the spread of 17 models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) that realistically simulate ENSO diversity. Based on 40 Community Earth System Model Large Ensemble (CESM-LE) and 99 Max Planck Institute for Meteorology Grand Ensemble (MPI-GE) members, unforced variability can explain more than 90% of the historical EP and CP El Niño standard deviation and all of the ENSO skewness spread in the 17 CMIP5 models. Both CESM-LE and the selected CMIP5 models show increased EP and CP El Niño variability in a warmer climate, driven by a stronger mean vertical temperature gradient in the upper ocean and faster surface warming of the eastern equatorial Pacific. However, MPI-GE shows no agreement in EP or CP standard deviation change. This is due to weaker sensitivity to the warming signal, such that when the eastern equatorial Pacific surface warming is faster, the change in upper ocean vertical temperature gradient tends to be weaker. This highlights that individual models produce a different ENSO response in a warmer climate, and that considerable uncertainty within the CMIP5 ensemble may be caused by internal climate variability.