The zonal gradient in the equatorial Pacific sea surface temperature (SST), high in the west and low in the east, is known to be a pacemaker of global warming 1. The zonal SST gradient has strengthened since the mid-20 th century 2 , but the cause is controversial because it is not reproduced by a majority of Coupled Model Intercomparison Project Phase 5 (CMIP5) climate models 3-6 , which instead suggest a weakening of the zonal SST gradient from the past to the future 7. This discrepancy between observations and models has to be reconciled not only for attributing the past climate change but also for assessing Earth's climate sensitivity 8,9. Here we show using large ensemble (LE) simulations by four different climate models, in addition to CMIP5 models, that the intensified trend in the Pacific zonal SST gradient for 1951-2010 can be captured by some realisations, suggesting that it could arise from internal climate variability. Models and members that simulate the past strengthening of the SST gradient in CMIP5 and LE commonly exhibit reversed trends in future projections, when the rate of global-mean temperature rise is likely to amplify by 9-30% with the larger values occurring in low-emission scenarios. Atmosphere-ocean state in the tropical Pacific is characterised by the zonal contrast between high SST in the warm pool and low SST in the cold tongue, tied with the Walker circulation accompanying surface easterlies. This mean state is maintained by the well-known Bjerknes feedback 10 , which acts against radiative heating that homogenises SST zonally. Specifically, equatorial easterlies pile warm water up to the west and cause upwelling to cool the eastern Pacific. The strengthened SST gradient, in turn, amplifies the easterlies by