The climate modeling community has been challenged to develop a method for improving the simulation of the Pacific-North America (PNA) teleconnection pattern in climate models. The accuracy of PNA teleconnection simulation is significantly improved by considering mesoscale convection contributions to sea surface fluxes. The variation in the PNA over the past 22 years was simulated by the Grid Atmospheric Model of IAP LASG version 1.0 (GAMIL1.0), which was guided by observational SST from January 1979 to December 2000. Results show that heating in the tropical central-eastern Pacific is simulated more realistically, and sea surface latent heat flux and precipitation anomalies are more similar to the reanalysis data when mesoscale enhance ment is considered during the parameterization scheme of sea surface turbulent fluxes in GAMIL1.0. Realistic heat ing in the tropical central-eastern Pacific in turn signifi cantly improves the simulation of interannual variation and spatial patterns of PNA.